Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Protected Member Functions | Protected Attributes | List of all members
Nektar::LinearSWESolver Class Reference

#include <LinearSWESolver.h>

Inheritance diagram for Nektar::LinearSWESolver:
Inheritance graph
[legend]
Collaboration diagram for Nektar::LinearSWESolver:
Collaboration graph
[legend]

Protected Member Functions

 LinearSWESolver ()
 
virtual void v_Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
 
virtual void v_ArraySolve (const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
 
virtual void v_PointSolve (NekDouble etaL, NekDouble uL, NekDouble vL, NekDouble dL, NekDouble etaR, NekDouble uR, NekDouble vR, NekDouble dR, NekDouble &etaf, NekDouble &uf, NekDouble &vf)
 
- Protected Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT RiemannSolver ()
 
void GenerateRotationMatrices (const Array< OneD, const Array< OneD, NekDouble > > &normals)
 Generate rotation matrices for 3D expansions. More...
 
void FromToRotation (Array< OneD, const NekDouble > &from, Array< OneD, const NekDouble > &to, NekDouble *mat)
 A function for creating a rotation matrix that rotates a vector from into another vector to. More...
 
SOLVER_UTILS_EXPORT void rotateToNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
 Rotate a vector field to trace normal. More...
 
SOLVER_UTILS_EXPORT void rotateFromNormal (const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
 Rotate a vector field from trace normal. More...
 
SOLVER_UTILS_EXPORT bool CheckScalars (std::string name)
 Determine whether a scalar has been defined in m_scalars. More...
 
SOLVER_UTILS_EXPORT bool CheckVectors (std::string name)
 Determine whether a vector has been defined in m_vectors. More...
 
SOLVER_UTILS_EXPORT bool CheckParams (std::string name)
 Determine whether a parameter has been defined in m_params. More...
 
SOLVER_UTILS_EXPORT bool CheckAuxScal (std::string name)
 Determine whether a scalar has been defined in m_auxScal. More...
 
SOLVER_UTILS_EXPORT bool CheckAuxVec (std::string name)
 Determine whether a vector has been defined in m_auxVec. More...
 

Protected Attributes

bool m_pointSolve
 
- Protected Attributes inherited from Nektar::SolverUtils::RiemannSolver
bool m_requiresRotation
 Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields. More...
 
std::map< std::string,
RSScalarFuncType
m_scalars
 Map of scalar function types. More...
 
std::map< std::string,
RSVecFuncType
m_vectors
 Map of vector function types. More...
 
std::map< std::string,
RSParamFuncType
m_params
 Map of parameter function types. More...
 
std::map< std::string,
RSScalarFuncType
m_auxScal
 Map of auxiliary scalar function types. More...
 
std::map< std::string,
RSVecFuncType
m_auxVec
 Map of auxiliary vector function types. More...
 
Array< OneD, Array< OneD,
NekDouble > > 
m_rotMat
 Rotation matrices for each trace quadrature point. More...
 
Array< OneD, Array< OneD,
Array< OneD, NekDouble > > > 
m_rotStorage
 Rotation storage. More...
 

Additional Inherited Members

- Public Member Functions inherited from Nektar::SolverUtils::RiemannSolver
SOLVER_UTILS_EXPORT void Solve (const int nDim, const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
 Perform the Riemann solve given the forwards and backwards spaces. More...
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetScalar (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetScalar (std::string name, RSScalarFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetVector (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetVector (std::string name, RSVecFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetParam (std::string name, FuncPointerT func, ObjectPointerT obj)
 
void SetParam (std::string name, RSParamFuncType fp)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetAuxScal (std::string name, FuncPointerT func, ObjectPointerT obj)
 
template<typename FuncPointerT , typename ObjectPointerT >
void SetAuxVec (std::string name, FuncPointerT func, ObjectPointerT obj)
 
std::map< std::string,
RSScalarFuncType > & 
GetScalars ()
 
std::map< std::string,
RSVecFuncType > & 
GetVectors ()
 
std::map< std::string,
RSParamFuncType > & 
GetParams ()
 
- Public Attributes inherited from Nektar::SolverUtils::RiemannSolver
int m_spacedim
 

Detailed Description

Definition at line 45 of file LinearSWESolver.h.

Constructor & Destructor Documentation

Nektar::LinearSWESolver::LinearSWESolver ( )
protected

Definition at line 40 of file LinearSWESolver.cpp.

References Nektar::SolverUtils::RiemannSolver::m_requiresRotation.

40  : RiemannSolver(),
41  m_pointSolve(true)
42  {
43  m_requiresRotation = true;
44 
45  }
SOLVER_UTILS_EXPORT RiemannSolver()
bool m_requiresRotation
Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields...

Member Function Documentation

virtual void Nektar::LinearSWESolver::v_ArraySolve ( const Array< OneD, const Array< OneD, NekDouble > > &  Fwd,
const Array< OneD, const Array< OneD, NekDouble > > &  Bwd,
Array< OneD, Array< OneD, NekDouble > > &  flux 
)
inlineprotectedvirtual

Definition at line 58 of file LinearSWESolver.h.

References ASSERTL0.

Referenced by v_Solve().

62  {
63  ASSERTL0(false, "This function should be defined by subclasses.");
64  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
virtual void Nektar::LinearSWESolver::v_PointSolve ( NekDouble  etaL,
NekDouble  uL,
NekDouble  vL,
NekDouble  dL,
NekDouble  etaR,
NekDouble  uR,
NekDouble  vR,
NekDouble  dR,
NekDouble etaf,
NekDouble uf,
NekDouble vf 
)
inlineprotectedvirtual

Reimplemented in Nektar::LinearAverageSolver, and Nektar::LinearHLLSolver.

Definition at line 66 of file LinearSWESolver.h.

References ASSERTL0.

Referenced by v_Solve().

70  {
71  ASSERTL0(false, "This function should be defined by subclasses.");
72  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
void Nektar::LinearSWESolver::v_Solve ( const int  nDim,
const Array< OneD, const Array< OneD, NekDouble > > &  Fwd,
const Array< OneD, const Array< OneD, NekDouble > > &  Bwd,
Array< OneD, Array< OneD, NekDouble > > &  flux 
)
protectedvirtual

Implements Nektar::SolverUtils::RiemannSolver.

Definition at line 47 of file LinearSWESolver.cpp.

References ASSERTL0, m_pointSolve, Nektar::SolverUtils::RiemannSolver::m_scalars, v_ArraySolve(), and v_PointSolve().

52  {
53 
54  // extract the forward and backward trace of the depth
55  const Array<OneD, NekDouble> &dFwd = m_scalars["depthFwd"]();
56  const Array<OneD, NekDouble> &dBwd = m_scalars["depthBwd"]();
57 
58 
59  if (m_pointSolve)
60  {
61 
62  int expDim = Fwd.num_elements()-1;
63  NekDouble vf;
64 
65  if (expDim == 1)
66  {
67  for (int i = 0; i < Fwd[0].num_elements(); ++i)
68  {
70  Fwd [0][i], Fwd [1][i], 0.0, dFwd[i],
71  Bwd [0][i], Bwd [1][i], 0.0, dBwd[i],
72  flux[0][i], flux[1][i], vf);
73  }
74  }
75  else if (expDim == 2)
76  {
77  for (int i = 0; i < Fwd[0].num_elements(); ++i)
78  {
80  Fwd [0][i], Fwd [1][i], Fwd [2][i], dFwd[i],
81  Bwd [0][i], Bwd [1][i], Bwd [2][i], dBwd[i],
82  flux[0][i], flux[1][i], flux[2][i]);
83  }
84  }
85  else if (expDim == 3)
86  {
87  ASSERTL0(false, "No 3D Shallow water supported.");
88  }
89  }
90  else
91  {
92  v_ArraySolve(Fwd, Bwd, flux);
93  }
94  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
virtual void v_ArraySolve(const Array< OneD, const Array< OneD, NekDouble > > &Fwd, const Array< OneD, const Array< OneD, NekDouble > > &Bwd, Array< OneD, Array< OneD, NekDouble > > &flux)
std::map< std::string, RSScalarFuncType > m_scalars
Map of scalar function types.
virtual void v_PointSolve(NekDouble etaL, NekDouble uL, NekDouble vL, NekDouble dL, NekDouble etaR, NekDouble uR, NekDouble vR, NekDouble dR, NekDouble &etaf, NekDouble &uf, NekDouble &vf)
double NekDouble

Member Data Documentation

bool Nektar::LinearSWESolver::m_pointSolve
protected

Definition at line 48 of file LinearSWESolver.h.

Referenced by v_Solve().