diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp index 957c2e77f615eef56523063555d515951f0dd5dd..b83e94cd0e83bbaefd4cade9b4e2c8fb14aa0b85 100644 --- a/src/scheme/FiniteVolumesDiffusion.hpp +++ b/src/scheme/FiniteVolumesDiffusion.hpp @@ -76,7 +76,7 @@ private: Kokkos::View<Rd**> // Fonction qui calcule F_jr computeFjr(const Kokkos::View<const Rd**>& Cjr, const Kokkos::View<const Rd*>& uj, - const Kokkos::View<const Rd*>& xr, + const Kokkos::View<const Rd*>& xj, const Kokkos::View<const double*>& kj) { const Kokkos::View<const unsigned int**>& cell_nodes = m_connectivity.cellNodes(); const Kokkos::View<const unsigned short*> cell_nb_nodes @@ -141,6 +141,7 @@ private: } Kokkos::View<Rd**> m_Fjr; + Kokkos::View<Rd**> m_Gjr; Kokkos::View<double*> m_CFL; public: @@ -149,7 +150,8 @@ public: : m_mesh_data(mesh_data), m_mesh(mesh_data.mesh()), m_connectivity(m_mesh.connectivity()), - m_Fjr("Fjr", m_mesh.numberOfCells(), m_connectivity.maxNbNodePerCell()), + m_Fjr("Fjr", m_mesh.numberOfCells(), m_connectivity.maxNbNodePerCell()), + m_Gjr("Gjr", m_mesh.numberOfCells(), m_connectivity.maxNbNodePerCell()), m_CFL("CFL", m_mesh.numberOfCells()) { ; @@ -160,7 +162,7 @@ public: KOKKOS_INLINE_FUNCTION double diffusion_dt(const Kokkos::View<const double*>& Vj, const Kokkos::View<const double*>& rhoj, - const Kokkos::View<const Rd*>& xr, + const Kokkos::View<const Rd*>& xj, const Kokkos::View<const Rd*>& kj) const { Kokkos::View<double*> dt_j("dt_j", m_mesh.numberOfCells()); @@ -174,7 +176,7 @@ public: return dt; } - // Avance la valeur des inconnues pendant un pas de temps dt // A MODIFIER + // Avance la valeur des inconnues pendant un pas de temps dt void computeNextStep(const double& t, const double& dt, UnknownsType& unknowns) { @@ -194,6 +196,7 @@ public: computeExplicitFluxes(xr, xj, rhoj, uj, Cjr); const Kokkos::View<const Rd**> Fjr = m_Fjr; + const Kokkos::View<const Rd**> Gjr = m_Gjr; const Kokkos::View<const unsigned short*> cell_nb_nodes = m_connectivity.cellNbNodes(); const Kokkos::View<const unsigned int**>& cell_nodes