reset n_boundaries_per_direction

trixi-framework · Mar 12, 2024 · 738de75 · 738de75
1 parent f235619
commit 738de75
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/solvers/dgsem_tree/containers_2d.jl b/src/solvers/dgsem_tree/containers_2d.jl
@@ -421,6 +421,8 @@ end
 function init_boundaries!(boundaries, elements, mesh::TreeMesh2D)
     # Exit early if there are no boundaries to initialize
     if nboundaries(boundaries) == 0
+        # In this case n_boundaries_per_direction still needs to be reset!
+        boundaries.n_boundaries_per_direction = SVector(0, 0, 0, 0)
         return nothing
     end