54 int expDim = Fwd.num_elements()-1;
59 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
62 Fwd [0][i], Fwd [1][i], 0.0,
63 Bwd [0][i], Bwd [1][i], 0.0,
64 flux[0][i], flux[1][i], hvf);
69 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
72 Fwd [0][i], Fwd [1][i], Fwd [2][i],
73 Bwd [0][i], Bwd [1][i], Bwd [2][i],
74 flux[0][i], flux[1][i], flux[2][i]);
79 ASSERTL0(
false,
"No 3D Shallow water supported.");
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)
#define ASSERTL0(condition, msg)
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_PointSolve(NekDouble hL, NekDouble huL, NekDouble hvL, NekDouble hR, NekDouble huR, NekDouble hvR, NekDouble &hf, NekDouble &huf, NekDouble &hvf)
bool m_requiresRotation
Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields...
The RiemannSolver class provides an abstract interface under which solvers for various Riemann proble...