44 "LaxFriedrichs Riemann solver");
65 double hL,
double huL,
double hvL,
66 double hR,
double huR,
double hvR,
67 double &hf,
double &huf,
double &hvf)
87 NekDouble uRoe = (srL * uL + srR * uR) / srLR;
88 NekDouble cRoe = sqrt(0.5 * (cL * cL + cR * cR));
91 NekDouble S = std::max(uRoe+cRoe, std::max(uR+cR, -(uL-cL)));
100 hf = 0.5 * ((huL + huR) - sign * S * (hR -hL));
103 huf = 0.5 * ((hL * uL * uL + 0.5 * g * hL * hL +
104 hR * uR * uR + 0.5 * g * hR * hR) -
105 sign * S * (huR - huL));
108 hvf = 0.5 * ((hL * uL * vL + hR * uR * vR) -
109 sign * S * (hvR - hvL));