Fix solve_level == solve_level_max

1 files changed, 3 insertions, 2 deletions

diff --git a/src/qms.c b/src/qms.c
index 047e4b6..e3a0990 100644
--- a/src/qms.c
+++ b/src/qms.c
@@ -1660,12 +1660,13 @@ void qms_mg_inithist(CCTK_ARGUMENTS)
     DECLARE_CCTK_PARAMETERS;
 
     const int reflevel = ctz(cctk_levfac[0]);
+    const int reflevel_solve = MIN(reflevel, solve_level);
 
     size_t grid_size = cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2];
     double dt;
 
-    dt = cctk_delta_time / (1 << MIN(reflevel, solve_level));
-    if (reflevel > solve_level_max)
+    dt = cctk_delta_time / (1 << reflevel_solve);
+    if (reflevel_solve > solve_level_max)
         dt *= 2;
 
     fprintf(stderr, "%d inithist\n", reflevel);