44 std::string LaxFriedrichsSolver::solverName =
46 RegisterCreatorFunction(
"LaxFriedrichs", LaxFriedrichsSolver::create,
47 "Lax-Friedrichs Solver");
53 LaxFriedrichsSolver::LaxFriedrichsSolver() :
96 NekDouble a_1_max = std::max(std::abs(u0 - c), std::abs(u0 + c));
99 NekDouble uFL = pL/rho0 + u0*uL + v0*vL + w0*wL;
104 NekDouble uFR = pR/rho0 + u0*uR + v0*vR + w0*wR;
109 pF = 0.5*(pFL + pFR - a_1_max*(pR/(c*c) - pL/(c*c)));
110 uF = 0.5*(uFL + uFR - a_1_max*(uR - uL));
111 vF = 0.5*(vFL + vFR - a_1_max*(vR - vL));
112 wF = 0.5*(wFL + wFR - a_1_max*(wR - 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)
Lax-Friedrichs 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...