From a38f0fe285a43f0e22b542a0abfd14cdc3bb3631 Mon Sep 17 00:00:00 2001 From: HOCH PHILIPPE <philippe.hoch@gmail.com> Date: Mon, 31 Mar 2025 16:39:39 +0200 Subject: [PATCH] Correction of some Boundary Condition: Edge vers Face --- src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp | 6 +++--- src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp | 6 +++--- src/scheme/RusanovEulerianCompositeSolver.cpp | 6 +++--- src/scheme/RusanovEulerianCompositeSolver_o2.cpp | 6 +++--- src/scheme/RusanovEulerianCompositeSolver_v2.cpp | 6 +++--- src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp | 6 +++--- 6 files changed, 18 insertions(+), 18 deletions(-) diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp index ce30fc35..d56dc2a5 100644 --- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp +++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2.cpp @@ -352,7 +352,7 @@ class RoeViscousFormEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -459,7 +459,7 @@ class RoeViscousFormEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -579,7 +579,7 @@ class RoeViscousFormEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; diff --git a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp index d22e940a..e707e3f8 100644 --- a/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp +++ b/src/scheme/RoeViscousFormEulerianCompositeSolver_v2_o2.cpp @@ -360,7 +360,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -467,7 +467,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -587,7 +587,7 @@ class RoeViscousFormEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; diff --git a/src/scheme/RusanovEulerianCompositeSolver.cpp b/src/scheme/RusanovEulerianCompositeSolver.cpp index a13b34fb..cda7e625 100644 --- a/src/scheme/RusanovEulerianCompositeSolver.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver.cpp @@ -245,7 +245,7 @@ class RusanovEulerianCompositeSolver Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -481,7 +481,7 @@ class RusanovEulerianCompositeSolver Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -588,7 +588,7 @@ class RusanovEulerianCompositeSolver Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; diff --git a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp index 867e6c4e..7d53c2f8 100644 --- a/src/scheme/RusanovEulerianCompositeSolver_o2.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver_o2.cpp @@ -253,7 +253,7 @@ class RusanovEulerianCompositeSolver_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -489,7 +489,7 @@ class RusanovEulerianCompositeSolver_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -596,7 +596,7 @@ class RusanovEulerianCompositeSolver_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp index cde234af..f367cb16 100644 --- a/src/scheme/RusanovEulerianCompositeSolver_v2.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver_v2.cpp @@ -246,7 +246,7 @@ class RusanovEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -483,7 +483,7 @@ class RusanovEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -590,7 +590,7 @@ class RusanovEulerianCompositeSolver_v2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; diff --git a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp index 7386f4e8..7a22ea9c 100644 --- a/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp +++ b/src/scheme/RusanovEulerianCompositeSolver_v2_o2.cpp @@ -253,7 +253,7 @@ class RusanovEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; @@ -488,7 +488,7 @@ class RusanovEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; Rdxd MatriceProj(identity); MatriceProj -= tensorProduct(normal, normal); @@ -595,7 +595,7 @@ class RusanovEulerianCompositeSolver_v2_o2 Rd vectorSym(zero); for (size_t dim = 0; dim < Dimension; ++dim) - vectorSym[dim] = stateEdge[face_cell_id][face_local_number_in_cell][1 + dim]; + vectorSym[dim] = stateFace[face_cell_id][face_local_number_in_cell][1 + dim]; vectorSym -= dot(vectorSym, normal) * normal; -- GitLab