correction boucle programme principal pour sous cycles de la diffusion

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();