From 7ded49e6a0f4c50b379da08767a24103ffef3fff Mon Sep 17 00:00:00 2001
From: Fanny CHOPOT <fanny.chopot.ocre@cea.fr>
Date: Fri, 1 Jun 2018 10:49:12 +0200
Subject: [PATCH] correction sous cycle diffusion pb temps

---
 src/main.cpp | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/src/main.cpp b/src/main.cpp
index 728770738..dff800896 100644
--- a/src/main.cpp
+++ b/src/main.cpp
@@ -189,18 +189,22 @@ int main(int argc, char *argv[])
       // ETAPE 2 DU SPLITTING - DIFFUSION
       
       double dt_diff = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj, cj);
+      double t_diff = t-dt_euler;
  
       if (dt_euler <= dt_diff) {
 	dt_diff = dt_euler;
-	finite_volumes_diffusion.computeNextStep(t, dt_diff, unknowns);
+	finite_volumes_diffusion.computeNextStep(t_diff, dt_diff, unknowns);
       } else {
-	double t_diff = t-dt_euler;
 	while (t > t_diff) {
-	  finite_volumes_diffusion.computeNextStep(t_diff, dt_diff, unknowns);
 	  dt_diff = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj, cj);
 	  if (t_diff+dt_diff > t) {
 	    dt_diff = t-t_diff;
 	  }
+	  finite_volumes_diffusion.computeNextStep(t_diff, dt_diff, unknowns);
+	  //dt_diff = 0.4*finite_volumes_diffusion.diffusion_dt(rhoj, kj, cj);
+	  //if (t_diff+dt_diff > t) {
+	  // dt_diff = t-t_diff;
+	  //}
 	  t_diff += dt_diff;
 	}
       }
-- 
GitLab