874   } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}, star;  in PhysicsRiemann_IsoGas_Exact()  local
891         star.rho = rho;  in PhysicsRiemann_IsoGas_Exact()
892         star.u   = L.u - c * fl;  in PhysicsRiemann_IsoGas_Exact()
905   if (L.u - c < 0 && 0 < star.u - c) { /* 1-wave is sonic rarefaction */  in PhysicsRiemann_IsoGas_Exact()
910   } else if (star.u + c < 0 && 0 < R.u + c) { /* 2-wave is sonic rarefaction */  in PhysicsRiemann_IsoGas_Exact()
915 …} else if ((L.rho >= star.rho && L.u - c >= 0) || (L.rho < star.rho && (star.rho * star.u - L.rho …  in PhysicsRiemann_IsoGas_Exact()
918 …} else if ((star.rho <= R.rho && R.u + c <= 0) || (star.rho > R.rho && (R.rho * R.u - star.rho * s…  in PhysicsRiemann_IsoGas_Exact()
922     ustar[0] = star.rho;  in PhysicsRiemann_IsoGas_Exact()
923     ustar[1] = star.rho * star.u;  in PhysicsRiemann_IsoGas_Exact()
926   *maxspeed = MaxAbs(MaxAbs(star.u - c, star.u + c), MaxAbs(L.u - c, R.u + c));  in PhysicsRiemann_IsoGas_Exact()
1018   } L = {uL[0], uL[1] / uL[0]}, R = {uR[0], uR[1] / uR[0]}, star;  in PhysicsRiemann_Shallow_Exact()  local
1040         star.h = h;  in PhysicsRiemann_Shallow_Exact()
1041         star.u = L.u - fl;  in PhysicsRiemann_Shallow_Exact()
1060   cstar = PetscSqrtScalar(g * star.h);  in PhysicsRiemann_Shallow_Exact()
1061   if (L.u - cL < 0 && 0 < star.u - cstar) { /* 1-wave is sonic rarefaction */  in PhysicsRiemann_Shallow_Exact()
1066   } else if (star.u + cstar < 0 && 0 < R.u + cR) { /* 2-wave is sonic rarefaction */  in PhysicsRiemann_Shallow_Exact()
1071 …} else if ((L.h >= star.h && L.u - c >= 0) || (L.h < star.h && (star.h * star.u - L.h * L.u) / (st…  in PhysicsRiemann_Shallow_Exact()
1074 …} else if ((star.h <= R.h && R.u + c <= 0) || (star.h > R.h && (R.h * R.u - star.h * star.h) / (R.…  in PhysicsRiemann_Shallow_Exact()
1078     ustar[0] = star.h;  in PhysicsRiemann_Shallow_Exact()
1079     ustar[1] = star.h * star.u;  in PhysicsRiemann_Shallow_Exact()
1082   *maxspeed = MaxAbs(MaxAbs(star.u - cstar, star.u + cstar), MaxAbs(L.u - cL, R.u + cR));  in PhysicsRiemann_Shallow_Exact()