58 int nTracePts = Fwd[0].size();
62 for (
int i = 0; i < nDim + 3; i++)
75 for (
int i = 0; i < nTracePts; ++i)
77 v_PointSolve(Fwd[0][i], Fwd[1][i], Fwd[2][i], 0.0, 0.0, Bwd[0][i],
78 Bwd[1][i], Bwd[2][i], 0.0, 0.0, bfFwd[0][i],
79 bfFwd[1][i], bfFwd[2][i], 0.0, 0.0, bfBwd[0][i],
80 bfBwd[1][i], bfBwd[2][i], 0.0, 0.0, flux[0][i],
81 flux[1][i], flux[2][i], vF, wF);
86 for (
int i = 0; i < nTracePts; ++i)
88 v_PointSolve(Fwd[0][i], Fwd[1][i], Fwd[2][i], Fwd[3][i], 0.0,
89 Bwd[0][i], Bwd[1][i], Bwd[2][i], Bwd[3][i], 0.0,
90 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i],
91 0.0, bfBwd[0][i], bfBwd[1][i], bfBwd[2][i],
92 bfBwd[3][i], 0.0, flux[0][i], flux[1][i], flux[2][i],
98 for (
int i = 0; i < nTracePts; ++i)
100 v_PointSolve(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 bfFwd[0][i], bfFwd[1][i], bfFwd[2][i], bfFwd[3][i],
103 bfFwd[4][i], bfBwd[0][i], bfBwd[1][i], bfBwd[2][i],
104 bfBwd[3][i], bfBwd[4][i], flux[0][i], flux[1][i],
105 flux[2][i], flux[3][i], flux[4][i]);
void GetRotBasefield(Array< OneD, Array< OneD, NekDouble > > &bfFwd, Array< OneD, Array< OneD, NekDouble > > &bfBwd)
virtual void v_PointSolve(NekDouble pL, NekDouble rhoL, NekDouble uL, NekDouble vL, NekDouble wL, NekDouble pR, NekDouble rhoR, NekDouble uR, NekDouble vR, NekDouble wR, NekDouble c0sqL, NekDouble rho0L, NekDouble u0L, NekDouble v0L, NekDouble w0L, NekDouble c0sqR, NekDouble rho0R, NekDouble u0R, NekDouble v0R, NekDouble w0R, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF)=0
LEESolver(const LibUtilities::SessionReaderSharedPtr &pSession)
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) override
bool m_requiresRotation
Indicates whether the Riemann solver requires a rotation to be applied to the velocity fields.
std::shared_ptr< SessionReader > SessionReaderSharedPtr