62 string ShallowWaterSystem::className =
65 ShallowWaterSystem::create,
66 "Auxiliary functions for the shallow water system.");
69 ShallowWaterSystem::ShallowWaterSystem(
83 "No UPWINDTYPE defined in session.");
91 m_vecLocs[0][i] = 1 + i;
159 ASSERTL0(
false,
"This function is not implemented for this equation.");
164 ASSERTL0(
false,
"This function is not implemented for this equation.");
183 for(
int bcRegion = 0; bcRegion <
m_fields[0]->GetBndConditions().num_elements(); ++bcRegion)
189 for (e = 0; e <
m_fields[0]->GetBndCondExpansions()[bcRegion]
192 npts =
m_fields[0]->GetBndCondExpansions()[bcRegion]->
193 GetExp(e)->GetTotPoints();
194 id2 =
m_fields[0]->GetTrace()->GetPhys_Offset(
196 GetBndCondCoeffsToGlobalCoeffsMap(cnt+e));
201 cnt +=
m_fields[0]->GetBndCondExpansions()[bcRegion]->GetExpSize();
Array< OneD, NekDouble > m_coriolis
Coriolis force.
#define ASSERTL0(condition, msg)
virtual void v_ConservativeToPrimitive()
bool m_primitive
Indicates if variables are primitive or conservative.
Array< OneD, NekDouble > m_depth
Still water depth.
std::vector< std::pair< std::string, std::string > > SummaryList
enum MultiRegions::ProjectionType m_projectionType
Type of projection; e.g continuous or discontinuous.
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
SOLVER_UTILS_EXPORT int GetTotPoints()
virtual void v_PrimitiveToConservative()
virtual SOLVER_UTILS_EXPORT void v_GenerateSummary(SummaryList &s)
Print a summary of time stepping parameters.
Array< OneD, Array< OneD, NekDouble > > m_bottomSlope
virtual void v_InitObject()
Init object for UnsteadySystem class.
Base class for unsteady solvers.
void AddSummaryItem(SummaryList &l, const std::string &name, const std::string &value)
Adds a summary item to the summary info list.
int m_spacedim
Spatial dimension (>= expansion dim).
void CopyBoundaryTrace(const Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
virtual void v_GenerateSummary(SolverUtils::SummaryList &s)
Print a summary of time stepping parameters.
void Neg(int n, T *x, const int incx)
Negate x = -x.
virtual SOLVER_UTILS_EXPORT void v_InitObject()
Init object for UnsteadySystem class.
SOLVER_UTILS_EXPORT void EvaluateFunction(Array< OneD, Array< OneD, NekDouble > > &pArray, std::string pFunctionName, const NekDouble pTime=0.0, const int domain=0)
Evaluates a function as specified in the session file.
EquationSystemFactory & GetEquationSystemFactory()
MultiRegions::Direction const DirCartesianMap[]
void EvaluateCoriolis(void)
Array< OneD, MultiRegions::ExpListSharedPtr > m_fields
Array holding all dependent variables.
LibUtilities::SessionReaderSharedPtr m_session
The session reader.
virtual ~ShallowWaterSystem()
Destructor.
bool m_constantDepth
Indicates if constant depth case.
int m_infosteps
Number of time steps between outputting status information.
void EvaluateWaterDepth(void)
Array< OneD, Array< OneD, NekDouble > > m_vecLocs
NekDouble m_g
Acceleration of gravity.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.