From 8ad661598d973e8a1ab2486d656ee40c8c101924 Mon Sep 17 00:00:00 2001 From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr> Date: Wed, 25 Apr 2018 10:39:24 +0200 Subject: [PATCH] essais sur dt et kj --- src/scheme/FiniteVolumesDiffusion.hpp | 6 +++--- src/scheme/FiniteVolumesEulerUnknowns.hpp | 2 +- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp index 694b7f43f..c47761df2 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 d1248d2e7..70ac38072 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 ? }); } -- GitLab