diff --git a/src/solvers/dgsem_p4est/subcell_limiters_2d.jl b/src/solvers/dgsem_p4est/subcell_limiters_2d.jl index 1acdb05aa7a..032047f2c4d 100644 --- a/src/solvers/dgsem_p4est/subcell_limiters_2d.jl +++ b/src/solvers/dgsem_p4est/subcell_limiters_2d.jl @@ -12,18 +12,15 @@ function calc_bounds_twosided_interface!(var_min, var_max, variable, u, t, semi, (; neighbor_ids, node_indices) = cache.interfaces index_range = eachnode(dg) - index_end = last(index_range) for interface in eachinterface(dg, cache) # Get element and side index information on the primary element primary_element = neighbor_ids[1, interface] primary_indices = node_indices[1, interface] - primary_direction = indices2direction(primary_indices) # Get element and side index information on the secondary element secondary_element = neighbor_ids[2, interface] secondary_indices = node_indices[2, interface] - secondary_direction = indices2direction(secondary_indices) # Create the local i,j indexing i_primary_start, i_primary_step = index_to_start_step_2d(primary_indices[1], diff --git a/src/solvers/dgsem_tree/dg_2d_subcell_limiters.jl b/src/solvers/dgsem_tree/dg_2d_subcell_limiters.jl index 971181fa596..a4765a25fb7 100644 --- a/src/solvers/dgsem_tree/dg_2d_subcell_limiters.jl +++ b/src/solvers/dgsem_tree/dg_2d_subcell_limiters.jl @@ -1843,7 +1843,7 @@ end normal_direction::AbstractVector, direction, equations, dg, indices...) factor = (normal_direction[1] * u_inner[2] + normal_direction[2] * u_inner[3]) - u_normal = (factor / sum(normal_direction.^2)) * normal_direction + u_normal = (factor / sum(normal_direction .^ 2)) * normal_direction return SVector(u_inner[1], u_inner[2] - 2.0 * u_normal[1], @@ -1888,7 +1888,7 @@ end u_outer = boundary_condition.boundary_value_function(boundary_condition.outer_boundary_value_function, u_inner, normal_direction / - sum(normal_direction.^2), + sum(normal_direction .^ 2), direction, x, t, equations) return u_outer