From 6c60c3ec4a4737fcabddc876e0043fa3e8c82b98 Mon Sep 17 00:00:00 2001 From: Simon Candelaresi <10759273+SimonCan@users.noreply.github.com> Date: Mon, 11 Dec 2023 18:07:35 +0000 Subject: [PATCH 1/3] Update src/semidiscretization/semidiscretization_coupled.jl Co-authored-by: Erik Faulhaber <44124897+efaulhaber@users.noreply.github.com> --- src/semidiscretization/semidiscretization_coupled.jl | 1 + 1 file changed, 1 insertion(+) diff --git a/src/semidiscretization/semidiscretization_coupled.jl b/src/semidiscretization/semidiscretization_coupled.jl index 38c1f9f34cb..2681e8ceb36 100644 --- a/src/semidiscretization/semidiscretization_coupled.jl +++ b/src/semidiscretization/semidiscretization_coupled.jl @@ -499,6 +499,7 @@ function mesh_equations_solver_cache(other_semi_index, i, semi_, semi_tuple...) if i == other_semi_index return mesh_equations_solver_cache(semi_) else + # Walk through semidiscretizations until we find `i` mesh_equations_solver_cache(other_semi_index, i + 1, semi_tuple...) end end From 388a30990ff50d7872e6440ce95b7438e66c997c Mon Sep 17 00:00:00 2001 From: Simon Candelaresi <10759273+SimonCan@users.noreply.github.com> Date: Mon, 11 Dec 2023 18:09:05 +0000 Subject: [PATCH 2/3] Update src/semidiscretization/semidiscretization_coupled.jl Co-authored-by: Erik Faulhaber <44124897+efaulhaber@users.noreply.github.com> --- src/semidiscretization/semidiscretization_coupled.jl | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/src/semidiscretization/semidiscretization_coupled.jl b/src/semidiscretization/semidiscretization_coupled.jl index 2681e8ceb36..8c97f155169 100644 --- a/src/semidiscretization/semidiscretization_coupled.jl +++ b/src/semidiscretization/semidiscretization_coupled.jl @@ -545,9 +545,10 @@ function copy_to_coupled_boundary!(boundary_condition::BoundaryConditionCoupled{ element_id = linear_indices[i_cell, j_cell] for element_id in eachnode(solver) - x = @view cache.elements.node_coordinates[:, i_node, j_node, - linear_indices[i_cell, j_cell]] - u_node = u[:, i_node, j_node, linear_indices[i_cell, j_cell]] + x = get_node_coords(node_coordinates, equations, solver, i_node, j_node, + linear_indices[i_cell, j_cell]) + u_node = get_node_vars(u, equations, solver, i_node, j_node, + linear_indices[i_cell, j_cell]] converted_u = coupling_converter(x, u_node) u_boundary[:, element_id, cell] = converted_u i_node += i_node_step From 0e7fd1960e227ec8c634908f2c3fd541ecaa6544 Mon Sep 17 00:00:00 2001 From: Simon Candelaresi <10759273+SimonCan@users.noreply.github.com> Date: Mon, 11 Dec 2023 18:09:30 +0000 Subject: [PATCH 3/3] Update src/semidiscretization/semidiscretization_coupled.jl Co-authored-by: Erik Faulhaber <44124897+efaulhaber@users.noreply.github.com> --- src/semidiscretization/semidiscretization_coupled.jl | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/semidiscretization/semidiscretization_coupled.jl b/src/semidiscretization/semidiscretization_coupled.jl index 8c97f155169..f03fae7b574 100644 --- a/src/semidiscretization/semidiscretization_coupled.jl +++ b/src/semidiscretization/semidiscretization_coupled.jl @@ -550,7 +550,9 @@ function copy_to_coupled_boundary!(boundary_condition::BoundaryConditionCoupled{ u_node = get_node_vars(u, equations, solver, i_node, j_node, linear_indices[i_cell, j_cell]] converted_u = coupling_converter(x, u_node) - u_boundary[:, element_id, cell] = converted_u + @inbounds for i in eachindex(converted_u) + u_boundary[i, element_id, cell] = converted_u[i] + end i_node += i_node_step j_node += j_node_step end