diff --git a/src/main.cpp b/src/main.cpp index 397f552a84dafbf4b23d2b04ce1e8c73d2d02f6c..a0acdd5f753a05c9c65ab6bbfcfeefb8d894709f 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 447311e00b6936a0527a042b07c2222e1e6eb43d..39539dbae4fa7e4a6f6ea1d45da0d5ccfb0c3ea9 100644 --- a/src/mesh/Mesh.hpp +++ b/src/mesh/Mesh.hpp @@ -85,7 +85,7 @@ public: } }); } - */ + */ ~Mesh() {