From 5f78d383ef1717527a232e86ca5426ed3d797965 Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Thu, 17 May 2018 15:19:50 +0200
Subject: [PATCH] correction boucle programme principal pour sous cycles de la
 diffusion

---
 src/main.cpp      | 10 ++++++----
 src/mesh/Mesh.hpp |  2 +-
 2 files changed, 7 insertions(+), 5 deletions(-)

diff --git a/src/main.cpp b/src/main.cpp
index 397f552a8..a0acdd5f7 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -173,12 +173,14 @@ int main(int argc, char *argv[])
 	finite_volumes_diffusion.computeNextStep(t, dt_diff, unknowns);
 	t += dt_euler;
       } else {
-	while (dt_euler > dt_diff) {
-	  finite_volumes_diffusion.computeNextStep(t, dt_diff, unknowns);
+	double t_diff = t + dt_diff;
+	while (t + dt_euler > t_diff) {
+	  finite_volumes_diffusion.computeNextStep(t_diff, dt_diff, unknowns);
 	  dt_diff = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj);
-	  std::cout << dt_diff << '\n';
+	  t_diff += dt_diff;
+	  std::cout << t_diff << '\n';
 	}
-	t += dt_diff;
+	t = t_diff;
       }
       
       block_eos.updatePandCFromRhoE();    
diff --git a/src/mesh/Mesh.hpp b/src/mesh/Mesh.hpp
index 447311e00..39539dbae 100644
--- a/src/mesh/Mesh.hpp
+++ b/src/mesh/Mesh.hpp
@@ -85,7 +85,7 @@ public:
   	}
       });
   }
-  */  
+  */
 
   ~Mesh()
   {
-- 
GitLab