44 : NonlinearSWESolver(pSession)
65 static NekDouble g = m_params[
"gravity"]();
79 hstar = 0.5 * (cL + cR) + 0.25 * (uL - uR);
87 SL = uL - cL *
sqrt(0.5 * ((hstar * hstar + hstar * hL) / (hL * hL)));
98 SR = uR + cR *
sqrt(0.5 * ((hstar * hstar + hstar * hR) / (hR * hR)));
108 huf = uL * uL * hL + 0.5 * g * hL * hL;
114 huf = uR * uR * hR + 0.5 * g * hR * hR;
119 hf = (SR * hL * uL - SL * hR * uR + 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 + SL * SR * (hR * uR - hL * uL)) / (SR - SL);
125 hvf = (SR * fL - SL * fR + SL * SR * (hR * vR - hL * vL)) / (SR - SL);
static std::string solverName
HLLSolver(const LibUtilities::SessionReaderSharedPtr &pSession)
static RiemannSolverSharedPtr create(const LibUtilities::SessionReaderSharedPtr &pSession)
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
std::shared_ptr< SessionReader > SessionReaderSharedPtr
RiemannSolverFactory & GetRiemannSolverFactory()
scalarT< T > sqrt(scalarT< T > in)