Revise bounds check for MCL

examples/tree_2d_dgsem/elixir_euler_sedov_blast_wave_MCL.jl

@@ -94,7 +94,8 @@ callbacks = CallbackSet(summary_callback,
 ###############################################################################
 # run the simulation
 
-stage_callbacks = (BoundsCheckCallback(save_errors=false),)
+output_directory = "out"
+stage_callbacks = (BoundsCheckCallback(save_errors=true, interval=100, output_directory = output_directory),)
 
 sol = Trixi.solve(ode, Trixi.SimpleSSPRK33(stage_callbacks=stage_callbacks);
                   dt=1.0, # solve needs some value here but it will be overwritten by the stepsize_callback

src/callbacks_stage/subcell_bounds_check.jl

@@ -118,8 +118,6 @@ function init_callback(callback::BoundsCheckCallback, semi, limiter::SubcellLimi
         return nothing
     end
 
-    # TODO: Revise Bounds Check for MCL
-
     @unpack output_directory = callback
     mkpath(output_directory)
     open("$output_directory/deviations.txt", "a") do f
@@ -128,7 +126,7 @@ function init_callback(callback::BoundsCheckCallback, semi, limiter::SubcellLimi
         if limiter.PressurePositivityLimiterKuzmin
             print(f, ", pressure_min")
         end
-        # No check for entropy limiting rn
+        # TODO: Bounds check for entropy limiting
         println(f)
     end
 
@@ -201,12 +199,12 @@ end
     println("─"^100)
     variables = varnames(cons2cons, semi.equations)
     for v in eachvariable(semi.equations)
-        println(variables[v], ":\n- lower bound: ", idp_bounds_delta[1, v],
-                "\n- upper bound: ", idp_bounds_delta[2, v])
+        println(variables[v], ":\n- lower bound: ", idp_bounds_delta[2, 1, v],
+                "\n- upper bound: ", idp_bounds_delta[2, 2, v])
     end
     if limiter.PressurePositivityLimiterKuzmin
-        println("pressure:\n- lower bound: ",
-                idp_bounds_delta[1, nvariables(semi.equations) + 1])
+        println("pressure:\n- positivity: ",
+                idp_bounds_delta[2, 1, nvariables(semi.equations) + 1])
     end
     println("─"^100 * "\n")