diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp index 694b7f43fc7afe7788b45b67c1025e41eb56d770..c47761df296c5562ae78a3870ae339fa47a959a3 100644 --- a/src/scheme/FiniteVolumesDiffusion.hpp +++ b/src/scheme/FiniteVolumesDiffusion.hpp @@ -166,12 +166,12 @@ public: const Kokkos::View<const double*>& Vj = m_mesh_data.Vj(); Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){ - //dt_j[j]= rhoj(j)*Vj(j)*(2./(kj(j+1) + 2*kj(j) + kj(j-1))) - // * std::min(Dj(j),Dj(j-1)); + //dt_j[j]= 0.5*rhoj(j)*Vj(j)*(2./(kj(j+1) + 2*kj(j) + kj(j-1))) + // * std::min(Dj(j),Dj(j-1)); // ATTENTION : kj != 0 ! // * std::min(xj(j+1)-xj(j), xj(j)-xj(j-1)); // * std::min((Dxj(j+1),Cjr(j,1)) + (xj(j),Cjr(j,0)), //(xj(j),Cjr(j,1)) + (xj(j-1),Cjr(j,0)) ); - dt_j[j] = 0.0001; // le signal vient d ici ! + dt_j[j] = 0.0001; }); double dt = std::numeric_limits<double>::max(); diff --git a/src/scheme/FiniteVolumesEulerUnknowns.hpp b/src/scheme/FiniteVolumesEulerUnknowns.hpp index d1248d2e7535c4ed4910e160281e4be6483cc38b..70ac380723b33a1d7ac6cef47455fba30756d196 100644 --- a/src/scheme/FiniteVolumesEulerUnknowns.hpp +++ b/src/scheme/FiniteVolumesEulerUnknowns.hpp @@ -173,7 +173,7 @@ public: }); Kokkos::parallel_for(m_mesh.numberOfCells(), KOKKOS_LAMBDA(const int& j){ - m_kj[j] = 0.; // Par quoi initialiser k ? + m_kj[j] = 0.1; // Par quoi initialiser k ? }); }