From 85bdf32d305f3ccc23585ac0711c3eb10c76c06a Mon Sep 17 00:00:00 2001
From: HOCH PHILIPPE <philippe.hoch@gmail.com>
Date: Sun, 18 Aug 2024 12:42:03 +0200
Subject: [PATCH] nettoyage Flux Rusanov, v1 et v2 et Roe v2

---
 ...RoeViscousFormEulerianCompositeSolver_v2.cpp | 17 +----------------
 src/scheme/RusanovEulerianCompositeSolver.cpp   | 17 +----------------
 .../RusanovEulerianCompositeSolver_v2.cpp       | 17 +----------------
 3 files changed, 3 insertions(+), 48 deletions(-)

diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp
index 2aab892d8..cfe2edb3f 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 d8db36aa2..dd0f65402 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 d490202f5..bd7a407b0 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);
             }
           });
-- 
GitLab