From 9d047dff05c7976bede70337e8f74a76855916b3 Mon Sep 17 00:00:00 2001 From: huiyuxie Date: Tue, 3 Sep 2024 15:55:57 -1000 Subject: [PATCH] Change for 1D --- src/solvers/dg_1d.jl | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/src/solvers/dg_1d.jl b/src/solvers/dg_1d.jl index 003452c..233d8cb 100644 --- a/src/solvers/dg_1d.jl +++ b/src/solvers/dg_1d.jl @@ -245,9 +245,6 @@ function prolong_boundaries_kernel!(boundaries_u, u, neighbor_ids, neighbor_side side = neighbor_sides[k] @inbounds begin - # boundaries_u[1, j, k] = u[j, size(u, 2), element] * isequal(side, 1) - # boundaries_u[2, j, k] = u[j, 1, element] * (1 - isequal(side, 1)) - boundaries_u[1, j, k] = u[j, size(u, 2), element] * (2 - side) # set to 0 instead of NaN boundaries_u[2, j, k] = u[j, 1, element] * (side - 1) # set to 0 instead of NaN end @@ -272,7 +269,7 @@ function boundary_flux_kernel!(surface_flux_values, boundaries_u, node_coordinat orientation = orientations[boundary] u_ll, u_rr = get_surface_node_vars(boundaries_u, equations, boundary) - u_inner = isequal(side, 1) * u_ll + (1 - isequal(side, 1)) * u_rr + u_inner = (2 - side) * u_ll + (side - 1) * u_rr x = get_node_coords(node_coordinates, equations, boundary) # TODO: Improve this part