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);
87 ustar = 0.5*(uL + uR) + cL - cR;
93 SL = uL - cL * sqrt(0.5*((hstar*hstar + hstar*hL)/(hL*hL)));
100 SR = uR + cR * sqrt(0.5*((hstar*hstar + hstar*hR)/(hR*hR)));
107 huf = uL * uL * hL + 0.5 * g * hL * hL;
113 huf = uR * uR * hR + 0.5 * g * hR * hR;
118 hf = (SR * hL * uL - SL * hR * uR +
119 SL * SR * (hR - hL)) / (SR - SL);
120 fL = uL * uL * hL + 0.5 * g * hL * hL;
121 fR = uR * uR * hR + 0.5 * g * hR * hR;
122 huf =(SR * fL - SL * fR +
123 SL * SR * (hR * uR - hL * uL)) / (SR - SL);
126 hvf =(SR * fL - SL * fR +
127 SL * SR * (hR * vR - hL * vL)) / (SR - SL);