44 "HLL Riemann solver");
65 double hL,
double huL,
double hvL,
66 double hR,
double huR,
double hvR,
67 double &hf,
double &huf,
double &hvf)
84 hstar = 0.5 * (cL + cR) + 0.25 * (uL - uR);
92 SL = uL - cL * sqrt(0.5*((hstar*hstar + hstar*hL)/(hL*hL)));
99 SR = uR + cR * sqrt(0.5*((hstar*hstar + hstar*hR)/(hR*hR)));
106 huf = uL * uL * hL + 0.5 * g * hL * hL;
112 huf = uR * uR * hR + 0.5 * g * hR * hR;
117 hf = (SR * hL * uL - SL * hR * uR +
118 SL * SR * (hR - hL)) / (SR - SL);
119 fL = uL * uL * hL + 0.5 * g * hL * hL;
120 fR = uR * uR * hR + 0.5 * g * hR * hR;
121 huf =(SR * fL - SL * fR +
122 SL * SR * (hR * uR - hL * uL)) / (SR - SL);
125 hvf =(SR * fL - SL * fR +
126 SL * SR * (hR * vR - hL * vL)) / (SR - SL);