diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp index 2aab892d8df4533b7d2aef934add050cfe926dce..cfe2edb3ff33fac4480fd3c199c259137d4ccaed 100644 --- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp +++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp @@ -594,21 +594,6 @@ class RoeViscousFormEulerianCompositeSolver_v2 } public: - double - EvaluateMaxEigenValueInGivenUnitDirection(const double& rho_mean, - const Rd& U_mean, - const double& E_mean, - const double& c_mean, - const Rd& normal) const - { - const double norme_normal = l2norm(normal); - Rd unit_normal = normal; - unit_normal *= 1. / norme_normal; - const double uscaln = dot(U_mean, unit_normal); - - return std::max(std::fabs(uscaln - c_mean), std::fabs(uscaln + c_mean)); - } - inline double pression(const double rho, const double epsilon, const double gam) const { @@ -1631,7 +1616,7 @@ class RoeViscousFormEulerianCompositeSolver_v2 const unsigned int R = edge_local_number_in_its_cells[k]; SumGje += Gje[K][R]; } - // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2norm(SumGje)); + // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2Norm(SumGje)); MaxErrorConservationEdge[l] = l2Norm(SumGje); } }); diff --git a/src/scheme/RusanovEulerianCompositeSolver.cpp b/src/scheme/RusanovEulerianCompositeSolver.cpp index d8db36aa23080b4c29b2ea3e60910964a8f14520..dd0f654023b7894fa9104f9a00f1371ee7bec577 100644 --- a/src/scheme/RusanovEulerianCompositeSolver.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver.cpp @@ -592,21 +592,6 @@ class RusanovEulerianCompositeSolver } public: - double - EvaluateMaxEigenValueInGivenUnitDirection(const double& rho_mean, - const Rd& U_mean, - const double& E_mean, - const double& c_mean, - const Rd& normal) const - { - const double norme_normal = l2norm(normal); - Rd unit_normal = normal; - unit_normal *= 1. / norme_normal; - const double uscaln = dot(U_mean, unit_normal); - - return std::max(std::fabs(uscaln - c_mean), std::fabs(uscaln + c_mean)); - } - inline double pression(const double rho, const double epsilon, const double gam) const { @@ -1209,7 +1194,7 @@ class RusanovEulerianCompositeSolver const unsigned int R = edge_local_number_in_its_cells[k]; SumGje += Gje[K][R]; } - // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2norm(SumGje)); + // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2Norm(SumGje)); MaxErrorConservationEdge[l] = l2Norm(SumGje); } }); diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp index d490202f5f433aaa13867f97d30134f38d2be5fd..bd7a407b009b1c0cbe6e1f9ec0fa27849138fd60 100644 --- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp @@ -592,21 +592,6 @@ class RusanovEulerianCompositeSolver_v2 } public: - double - EvaluateMaxEigenValueInGivenUnitDirection(const double& rho_mean, - const Rd& U_mean, - const double& E_mean, - const double& c_mean, - const Rd& normal) const - { - const double norme_normal = l2norm(normal); - Rd unit_normal = normal; - unit_normal *= 1. / norme_normal; - const double uscaln = dot(U_mean, unit_normal); - - return std::max(std::fabs(uscaln - c_mean), std::fabs(uscaln + c_mean)); - } - inline double pression(const double rho, const double epsilon, const double gam) const { @@ -1107,7 +1092,7 @@ class RusanovEulerianCompositeSolver_v2 const unsigned int R = edge_local_number_in_its_cells[k]; SumGje += Gje[K][R]; } - // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2norm(SumGje)); + // MaxErrorConservationEdge = std::max(MaxErrorConservationEdge, l2Norm(SumGje)); MaxErrorConservationEdge[l] = l2Norm(SumGje); } });