55 int nvariables = Fwd.num_elements();
62 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
65 Fwd [0][i], Fwd [1][i], 0.0, 0.0, Fwd [2][i],
66 Bwd [0][i], Bwd [1][i], 0.0, 0.0, Bwd [2][i],
67 flux[0][i], flux[1][i], rhouf, rhovf, flux[2][i]);
72 if (nvariables == expDim+2)
74 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
77 Fwd [0][i], Fwd [1][i], Fwd [2][i], 0.0, Fwd [3][i],
78 Bwd [0][i], Bwd [1][i], Bwd [2][i], 0.0, Bwd [3][i],
79 flux[0][i], flux[1][i], flux[2][i], rhovf, flux[3][i]);
83 if (nvariables > expDim+2)
85 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
88 Fwd [0][i], Fwd [1][i], Fwd [2][i], 0.0, Fwd [3][i], Fwd [4][i],
89 Bwd [0][i], Bwd [1][i], Bwd [2][i], 0.0, Bwd [3][i], Bwd [4][i],
90 flux[0][i], flux[1][i], flux[2][i], rhovf, flux[3][i], flux[4][i]);
97 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
100 Fwd [0][i], Fwd [1][i], Fwd [2][i], Fwd [3][i], Fwd [4][i],
101 Bwd [0][i], Bwd [1][i], Bwd [2][i], Bwd [3][i], Bwd [4][i],
102 flux[0][i], flux[1][i], flux[2][i], flux[3][i], flux[4][i]);
104 if (nvariables > expDim+2)
106 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
109 Fwd [0][i], Fwd [1][i], Fwd [2][i], Fwd [3][i], Fwd [4][i], Fwd [5][i],
110 Bwd [0][i], Bwd [1][i], Bwd [2][i], Bwd [3][i], Bwd [4][i], Bwd [5][i],
111 flux[0][i], flux[1][i], flux[2][i], flux[3][i], flux[4][i], flux[5][i]);
virtual void v_PointSolve(NekDouble rhoL, NekDouble rhouL, NekDouble rhovL, NekDouble rhowL, NekDouble EL, NekDouble rhoR, NekDouble rhouR, NekDouble rhovR, NekDouble rhowR, NekDouble ER, NekDouble &rhof, NekDouble &rhouf, NekDouble &rhovf, NekDouble &rhowf, NekDouble &Ef)
virtual void v_PointSolveVisc(NekDouble rhoL, NekDouble rhouL, NekDouble rhovL, NekDouble rhowL, NekDouble EL, NekDouble EpsL, NekDouble rhoR, NekDouble rhouR, NekDouble rhovR, NekDouble rhowR, NekDouble ER, NekDouble EpsR, NekDouble &rhof, NekDouble &rhouf, NekDouble &rhovf, NekDouble &rhowf, NekDouble &Ef, NekDouble &Epsf)
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_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...