correction sous cycle diffusion pb temps

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;
 	}
       }