45 RegisterCreatorFunction(
"LaxFriedrichs", LaxFriedrichsSolver::create,
"Lax-Friedrichs Solver");
51 LaxFriedrichsSolver::LaxFriedrichsSolver() :
92 NekDouble a_1_max = std::max(std::abs(u0 - c), std::abs(u0 + c));
95 NekDouble uFL = pL/rho + u0*uL + v0*vL + w0*wL;
100 NekDouble uFR = pR/rho + u0*uR + v0*vR + w0*wR;
105 pF = 0.5*(pFL + pFR - a_1_max*(pR - pL));
106 uF = 0.5*(uFL + uFR - a_1_max*(uR - uL));
107 vF = 0.5*(vFL + vFR - a_1_max*(vR - vL));
108 wF = 0.5*(wFL + wFR - a_1_max*(wR - wL));