reset n_boundaries_per_direction (#1870)

trixi-framework · Mar 15, 2024 · aa9ea20 · aa9ea20
1 parent a528083
commit aa9ea20
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