65 "Auxiliary functions for the shallow water system.");
82 "No UPWINDTYPE defined in session.");
151 ASSERTL0(
false,
"This function is not implemented for this equation.");
156 ASSERTL0(
false,
"This function is not implemented for this equation.");
175 for (
int bcRegion = 0; bcRegion <
m_fields[0]->GetBndConditions().size();
179 ->GetBndConditions()[bcRegion]
193 ->GetBndCondExpansions()[bcRegion]
196 id2 =
m_fields[0]->GetTrace()->GetPhys_Offset(
197 m_fields[0]->GetTraceMap()->GetBndCondIDToGlobalTraceID(cnt +
203 cnt +=
m_fields[0]->GetBndCondExpansions()[bcRegion]->GetExpSize();
#define ASSERTL0(condition, msg)
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
NekDouble m_g
Acceleration of gravity.
void EvaluateWaterDepth(void)
void CopyBoundaryTrace(const Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
Array< OneD, Array< OneD, NekDouble > > m_bottomSlope
virtual void v_ConservativeToPrimitive()
void EvaluateCoriolis(void)
virtual void v_PrimitiveToConservative()
bool m_constantDepth
Indicates if constant depth case.
bool m_primitive
Indicates if variables are primitive or conservative.
void v_GenerateSummary(SolverUtils::SummaryList &s) override
Print a summary of time stepping parameters.
ShallowWaterSystem(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Initialises UnsteadySystem class members.
Array< OneD, NekDouble > m_coriolis
Coriolis force.
Array< OneD, Array< OneD, NekDouble > > m_vecLocs
Array< OneD, NekDouble > m_depth
Still water depth.
static SolverUtils::EquationSystemSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
Creates an instance of this class.
void v_InitObject(bool DeclareFields=true) override
Init object for UnsteadySystem class.
static std::string className
Name of class.
~ShallowWaterSystem() override
Destructor.
int m_spacedim
Spatial dimension (>= expansion dim).
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array holding all dependent variables.
SOLVER_UTILS_EXPORT int GetExpSize()
LibUtilities::SessionReaderSharedPtr m_session
The session reader.
enum MultiRegions::ProjectionType m_projectionType
Type of projection; e.g continuous or discontinuous.
SOLVER_UTILS_EXPORT SessionFunctionSharedPtr GetFunction(std::string name, const MultiRegions::ExpListSharedPtr &field=MultiRegions::NullExpListSharedPtr, bool cache=false)
Get a SessionFunction by name.
SOLVER_UTILS_EXPORT int GetTotPoints()
Base class for unsteady solvers.
SOLVER_UTILS_EXPORT void v_GenerateSummary(SummaryList &s) override
Print a summary of time stepping parameters.
SOLVER_UTILS_EXPORT void v_InitObject(bool DeclareField=true) override
Init object for UnsteadySystem class.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
MultiRegions::Direction const DirCartesianMap[]
std::vector< std::pair< std::string, std::string > > SummaryList
EquationSystemFactory & GetEquationSystemFactory()
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
std::shared_ptr< MeshGraph > MeshGraphSharedPtr
void Neg(int n, T *x, const int incx)
Negate x = -x.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)