47 APESolver::APESolver() :
70 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
73 Fwd[0][i], 0.0, Fwd[1][i], 0.0, 0.0,
74 Bwd[0][i], 0.0, Bwd[1][i], 0.0, 0.0,
75 bf[0][i], bf[1][i], bf[2][i], 0.0, 0.0,
76 flux[0][i], rhoF, flux[1][i], vF, wF);
81 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
84 Fwd[0][i], 0.0, Fwd[1][i], Fwd[2][i], 0.0,
85 Bwd[0][i], 0.0, Bwd[1][i], Bwd[2][i], 0.0,
86 bf[0][i], bf[1][i], bf[2][i], bf[3][i], 0.0,
87 flux[0][i], rhoF, flux[1][i], flux[2][i], wF);
92 for (
int i = 0; i < Fwd[0].num_elements(); ++i)
95 Fwd[0][i], 0.0, Fwd[1][i], Fwd[2][i], Fwd[3][i],
96 Bwd[0][i], 0.0, Bwd[1][i], Bwd[2][i], Bwd[3][i],
97 bf[0][i], bf[1][i], bf[2][i], bf[3][i], bf[4][i],
98 flux[0][i], rhoF, flux[1][i], flux[2][i], flux[3][i]);
115 int nTracePts = normals[0].num_elements();
116 int nDim = normals.num_elements();
119 for (
int i = 0; i < nDim + 2; i++)
125 for (
int i = 0; i < nDim; ++i)
127 baseVecLocs[0][i] = i + 2;
Array< OneD, Array< OneD, NekDouble > > GetRotBasefield()
SOLVER_UTILS_EXPORT bool CheckVectors(std::string name)
Determine whether a vector has been defined in m_vectors.
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...
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 p0, NekDouble rho0, NekDouble u0, NekDouble v0, NekDouble w0, NekDouble &pF, NekDouble &rhoF, NekDouble &uF, NekDouble &vF, NekDouble &wF)
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
std::map< std::string, RSVecFuncType > m_vectors
Map of vector function types.
SOLVER_UTILS_EXPORT void rotateToNormal(const Array< OneD, const Array< OneD, NekDouble > > &inarray, const Array< OneD, const Array< OneD, NekDouble > > &normals, const Array< OneD, const Array< OneD, NekDouble > > &vecLocs, Array< OneD, Array< OneD, NekDouble > > &outarray)
Rotate a vector field to trace normal.