Dispatch check_bounds for dimension using u

bennibolm · May 24, 2024 · 60536e1 · 60536e1
1 parent 488b794
commit 60536e1
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/src/callbacks_stage/subcell_bounds_check_2d.jl b/src/callbacks_stage/subcell_bounds_check_2d.jl
@@ -5,7 +5,8 @@
 @muladd begin
 #! format: noindent
 
-@inline function check_bounds(u, equations, solver, cache,
+@inline function check_bounds(u::AbstractArray{<:Any, 4},
+                              equations, solver, cache,
                               limiter::SubcellLimiterIDP)
     (; local_twosided, positivity, local_onesided) = solver.volume_integral.limiter
     (; variable_bounds) = limiter.cache.subcell_limiter_coefficients
@@ -145,7 +146,8 @@ end
     return nothing
 end
 
-@inline function check_bounds(u, equations, solver, cache,
+@inline function check_bounds(u::AbstractArray{<:Any, 4},
+                              equations, solver, cache,
                               limiter::SubcellLimiterMCL)
     (; var_min, var_max) = limiter.cache.subcell_limiter_coefficients
     (; bar_states1, bar_states2, lambda1, lambda2) = limiter.cache.container_bar_states
@@ -632,7 +634,7 @@ end
     n_vars = nvariables(equations)
 
     # Print errors to output file
-    open("$output_directory/deviations.txt", "a") do f
+    open(joinpath(output_directory, "deviations.txt"), "a") do f
         print(f, iter, ", ", time)
         for v in eachvariable(equations)
             print(f, ", ", mcl_bounds_delta_local[1, v], ", ",