62 int expDim = Fwd.num_elements()-1;
67 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
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);
77 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
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]);
87 ASSERTL0(
false,
"No 3D Shallow water supported.");
#define ASSERTL0(condition, msg)
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)
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)
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...