fix tests

Author: tristanmontoya

src/solvers/dgmulti/dg_manifold_covariant.jl

@@ -97,14 +97,15 @@ function Trixi.calc_interface_flux!(cache, surface_integral::SurfaceIntegralWeak
         uP = u_face_values[idP]
         auxM = aux_face_values[idM]
         auxP = aux_face_values[idP]
+        
         # Transform uP to the same coordinate system as uM
         uP_global = contravariant2global(uP, auxP, equations)
         uP_transformed_to_M = global2contravariant(uP_global, auxM, equations)
 
         # Compute ref_index from face_node_index in order to access covariant normal vector
         # in reference element coordinates.
         ref_index = mod(face_node_index - 1, Nfq) + 1
-        normal = SVector(ntuple(k -> nrstJ[k][ref_index], NDIMS))
+        normal = SVector{NDIMS}(getindex.(nrstJ, ref_index))
 
         flux_face_values[idM] = surface_flux(uM, uP_transformed_to_M, auxM, auxM, normal,
                                              equations)

test/Project.toml

@@ -6,6 +6,7 @@ OrdinaryDiffEqLowStorageRK = "b0944070-b475-4768-8dec-fb6eb410534d"
 OrdinaryDiffEqSSPRK = "669c94d9-1f4b-4b64-b377-1aa079aa2388"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Trixi = "a7f1ee26-1774-49b1-8366-f1abc58fbfcb"
+TrixiAtmo = "c9ed1054-d170-44a9-8ee2-d5566f5d1389"
 TrixiTest = "0a316866-cbd0-4425-8bcb-08103b2c1f26"
 
 [compat]

test/test_2d_shallow_water_covariant.jl

@@ -100,17 +100,4 @@ end
     @test_allocations(TrixiAtmo.Trixi.rhs!, semi, sol, 100)
 end
 
-@trixi_testset "elixir_tri_barotropic_instability" begin
-    @test_trixi_include(joinpath(EXAMPLES_DIR,
-    "elixir_barotropic_instability.jl"),
-    l2=[41.05018196765347, 0.04598801953369521, 0.03324228006147076],
-    linf=[202.17195189961058, 0.2046503536574818, 0.14813768215260187]),
-    polydeg = 3,
-    initial_refinement_level = 1,
-    tspan = (0.0, 1.0 * SECONDS_PER_DAY)
-    # Ensure that we do not have excessive memory allocations
-    # (e.g., from type instabilities)
-    @test_allocations(TrixiAtmo.Trixi.rhs!, semi, sol, 100)
-end
-
 end # module

test/test_threaded.jl

@@ -29,7 +29,20 @@ include("test_trixiatmo.jl")
                         tspan=(0.0, 0.1))
     # Ensure that we do not have excessive memory allocations
     # (e.g., from type instabilities)
-    @test_allocations(TrixiAtmo.Trixi.rhs!, semi, sol, 2000)
+    @test_allocations(TrixiAtmo.Trixi.rhs!, semi, sol, 5000)
+end
+
+@trixi_testset "elixir_tri_barotropic_instability threaded" begin
+    @test_trixi_include(abspath(joinpath(EXAMPLES_DIR, "shallow_water/covariant",
+                                         "elixir_tri_barotropic_instability.jl")),
+                        l2=[41.05018196765347, 0.04598801953369521, 0.03324228006147076],
+                        linf=[202.17195189961058, 0.2046503536574818, 0.14813768215260187],
+                        polydeg=3,
+                        initial_refinement_level=1,
+                        tspan=(0.0, 1.0 * SECONDS_PER_DAY))
+    # Ensure that we do not have excessive memory allocations
+    # (e.g., from type instabilities)
+    @test_allocations(TrixiAtmo.Trixi.rhs!, semi, sol, 5000)
 end
 end
 end # module