From e85fa7e1cae760fea9058f63eb3eca7395d59dbe Mon Sep 17 00:00:00 2001 From: Tristan Montoya Date: Mon, 19 Aug 2024 20:29:58 +0200 Subject: [PATCH] modify calc_node_coordinates! to handle P4estMesh{2} with 3D coordinates --- src/solvers/dgsem_p4est/containers_2d.jl | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/solvers/dgsem_p4est/containers_2d.jl b/src/solvers/dgsem_p4est/containers_2d.jl index 236d7d24c06..a302266153f 100644 --- a/src/solvers/dgsem_p4est/containers_2d.jl +++ b/src/solvers/dgsem_p4est/containers_2d.jl @@ -43,7 +43,8 @@ function calc_node_coordinates!(node_coordinates, # places and the additional information passed to the compiler makes them faster # than native `Array`s. tmp1 = StrideArray(undef, real(mesh), - StaticInt(2), static_length(nodes), static_length(mesh.nodes)) + StaticInt(size(mesh.tree_node_coordinates, 1)), + static_length(nodes), static_length(mesh.nodes)) matrix1 = StrideArray(undef, real(mesh), static_length(nodes), static_length(mesh.nodes)) matrix2 = similar(matrix1)