829 …ubar = (uL[1] / PetscSqrtScalar(uL[0]) + uR[1] / PetscSqrtScalar(uR[0])) / (PetscSqrtScalar(uL[0])…  in PhysicsRiemann_IsoGas_Roe()
884       fl  = (L.rho < rho) ? (rho - L.rho) / PetscSqrtScalar(L.rho * rho) /* shock */  in PhysicsRiemann_IsoGas_Exact()
886       fr  = (R.rho < rho) ? (rho - R.rho) / PetscSqrtScalar(R.rho * rho) /* shock */  in PhysicsRiemann_IsoGas_Exact()
895 …dfl = (L.rho < rho) ? 1 / PetscSqrtScalar(L.rho * rho) * (1 - 0.5 * (rho - L.rho) / rho) : 1 / rho;  in PhysicsRiemann_IsoGas_Exact()
896 …dfr = (R.rho < rho) ? 1 / PetscSqrtScalar(R.rho * rho) * (1 - 0.5 * (rho - R.rho) / rho) : 1 / rho;  in PhysicsRiemann_IsoGas_Exact()
1023   cL = PetscSqrtScalar(g * L.h);  in PhysicsRiemann_Shallow_Exact()
1024   cR = PetscSqrtScalar(g * R.h);  in PhysicsRiemann_Shallow_Exact()
1033 …   fl  = (L.h < h) ? PetscSqrtScalar(0.5 * g * (h * h - L.h * L.h) * (1 / L.h - 1 / h)) /* shock */  in PhysicsRiemann_Shallow_Exact()
1034 …             : 2 * PetscSqrtScalar(g * h) - 2 * PetscSqrtScalar(g * L.h);         /* rarefaction */  in PhysicsRiemann_Shallow_Exact()
1035 …   fr  = (R.h < h) ? PetscSqrtScalar(0.5 * g * (h * h - R.h * R.h) * (1 / R.h - 1 / h)) /* shock */  in PhysicsRiemann_Shallow_Exact()
1036 …             : 2 * PetscSqrtScalar(g * h) - 2 * PetscSqrtScalar(g * R.h);         /* rarefaction */  in PhysicsRiemann_Shallow_Exact()
1050 …(L.h < h) ? 0.5 / fl * 0.5 * g * (-L.h * L.h / (h * h) - 1 + 2 * h / L.h) : PetscSqrtScalar(g / h);  in PhysicsRiemann_Shallow_Exact()
1051 …(R.h < h) ? 0.5 / fr * 0.5 * g * (-R.h * R.h / (h * h) - 1 + 2 * h / R.h) : PetscSqrtScalar(g / h);  in PhysicsRiemann_Shallow_Exact()
1060   cstar = PetscSqrtScalar(g * star.h);  in PhysicsRiemann_Shallow_Exact()
1064     ufan[1] = PetscSqrtScalar(g * ufan[0]) * ufan[0];  in PhysicsRiemann_Shallow_Exact()
1069     ufan[1] = -PetscSqrtScalar(g * ufan[0]) * ufan[0];  in PhysicsRiemann_Shallow_Exact()
1098 …s         = PetscMax(PetscAbs(L.u) + PetscSqrtScalar(g * L.h), PetscAbs(R.u) + PetscSqrtScalar(g *…  in PhysicsRiemann_Shallow_Rusanov()
1111   c            = PetscSqrtScalar(u[0] * phys->gravity);  in PhysicsCharacteristic_Shallow()