swap order of element and indices for consistency

trixi-framework · May 15, 2024 · 0873271 · 0873271
1 parent 746c822
commit 0873271
Show file tree

Hide file tree

Showing 4 changed files with 4 additions and 4 deletions.
diff --git a/src/callbacks_stage/positivity_zhang_shu_dg2d.jl b/src/callbacks_stage/positivity_zhang_shu_dg2d.jl
@@ -26,7 +26,7 @@ function limiter_zhang_shu!(u, threshold::Real, variable,
         total_volume = zero(eltype(u))
         for j in eachnode(dg), i in eachnode(dg)
             volume_jacobian = abs(inv(get_inverse_jacobian(inverse_jacobian, mesh,
-                                                           element, i, j)))
+                                                           i, j, element)))
             u_node = get_node_vars(u, equations, dg, i, j, element)
             u_mean += u_node * weights[i] * weights[j] * volume_jacobian
             total_volume += weights[i] * weights[j] * volume_jacobian

diff --git a/src/callbacks_stage/positivity_zhang_shu_dg3d.jl b/src/callbacks_stage/positivity_zhang_shu_dg3d.jl
@@ -26,7 +26,7 @@ function limiter_zhang_shu!(u, threshold::Real, variable,
         total_volume = zero(eltype(u))
         for k in eachnode(dg), j in eachnode(dg), i in eachnode(dg)
             volume_jacobian = abs(inv(get_inverse_jacobian(inverse_jacobian, mesh,
-                                                           element, i, j, k)))
+                                                           i, j, k, element)))
             u_node = get_node_vars(u, equations, dg, i, j, k, element)
             u_mean += u_node * weights[i] * weights[j] * weights[k] * volume_jacobian
             total_volume += weights[i] * weights[j] * weights[k] * volume_jacobian

diff --git a/src/solvers/dgsem_structured/dg.jl b/src/solvers/dgsem_structured/dg.jl
@@ -80,7 +80,7 @@ end
                                       mesh::Union{StructuredMesh, StructuredMeshView,
                                                   UnstructuredMesh2D, P4estMesh,
                                                   T8codeMesh},
-                                      element, indices...)
+                                      indices..., element)
     return inverse_jacobian[indices..., element]
 end
 

diff --git a/src/solvers/dgsem_tree/dg.jl b/src/solvers/dgsem_tree/dg.jl
@@ -43,7 +43,7 @@ function volume_jacobian(element, mesh::TreeMesh, cache)
 end
 
 @inline function get_inverse_jacobian(inverse_jacobian, mesh::TreeMesh,
-                                      element, indices...)
+                                      indices..., element)
     return inverse_jacobian[element]
 end