45 RegisterCreatorFunction(
"APEUpwind", UpwindSolver::create,
46 "Upwind solver for the APE equation");
48 UpwindSolver::UpwindSolver() :
102 characteristic[0] = pL/2 + uL*c*rho0/2;
103 characteristic[1] = pL/2 - uL*c*rho0/2;
105 characteristic[2] = pR/2 + uR*c*rho0/2;
106 characteristic[3] = pR/2 - uR*c*rho0/2;
109 for (
int j = 0; j < 2; j++)
113 W[j] = characteristic[j];
117 W[j] = characteristic[j+2];
126 pF = rho0*u + u0*p/(c*c);
127 uF = p/rho0 + u0*u + v0*vL + w0*wL;
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)
Upwind Riemann solver.
RiemannSolverFactory & GetRiemannSolverFactory()
std::map< std::string, RSParamFuncType > m_params
Map of parameter function types.
SOLVER_UTILS_EXPORT bool CheckParams(std::string name)
Determine whether a parameter has been defined in m_params.
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...