Speed up coverage tests

.github/workflows/ci.yml

@@ -111,7 +111,17 @@ jobs:
       - uses: julia-actions/julia-buildpkg@v1
         env:
           PYTHON: ""
-      - uses: julia-actions/julia-runtest@v1
+      - name: Run tests without coverage
+        uses: julia-actions/julia-runtest@v1
+        with:
+          coverage: false
+        env:
+          PYTHON: ""
+          TRIXI_TEST: ${{ matrix.trixi_test }}
+      - name: Run tests with coverage
+        uses: julia-actions/julia-runtest@v1
+        with:
+          coverage: true
         env:
           PYTHON: ""
           TRIXI_TEST: ${{ matrix.trixi_test }}

src/callbacks_step/glm_speed.jl

@@ -84,6 +84,9 @@ end
   # c_h is proportional to its own time step divided by the complete MHD time step
   equations.c_h = glm_scale * c_h_deltat / dt
 
+  # avoid re-evaluating possible FSAL stages
+  u_modified!(integrator, false)
+
   return nothing
 end
 

test/test_mpi.jl

@@ -42,28 +42,32 @@ const EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr.jl"),
       # Expected errors are exactly the same as in the serial test!
       l2   = [4.913300828257469e-5],
-      linf = [0.00045263895394385967])
+      linf = [0.00045263895394385967],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_amr_nonperiodic.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_nonperiodic.jl"),
       # Expected errors are exactly the same as in the serial test!
       l2   = [3.2207388565869075e-5],
-      linf = [0.0007508059772436404])
+      linf = [0.0007508059772436404],
+      coverage_override = (maxiters=6,))
   end
 
   # Linear scalar advection with AMR
   # These example files are only for testing purposes and have no practical use
   @trixi_testset "elixir_advection_amr_refine_twice.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_refine_twice.jl"),
       l2   = [0.00020547512522578292],
-      linf = [0.007831753383083506])
+      linf = [0.007831753383083506],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_amr_coarsen_twice.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_coarsen_twice.jl"),
       l2   = [0.0014321062757891826],
-      linf = [0.0253454486893413])
+      linf = [0.0253454486893413],
+      coverage_override = (maxiters=6,))
   end
 
   # Hyperbolic diffusion
@@ -106,7 +110,10 @@ const EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "
   @trixi_testset "elixir_euler_source_terms_amr_refine_coarsen.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_source_terms_amr_refine_coarsen.jl"),
       l2   = [4.8226610349853444e-5, 4.117706709270575e-5, 4.1177067092959676e-5, 0.00012205252427437389],
-      linf = [0.0003543874851490436, 0.0002973166773747593, 0.0002973166773760916, 0.001154106793870291])
+      linf = [0.0003543874851490436, 0.0002973166773747593, 0.0002973166773760916, 0.001154106793870291],
+      # Let this test run until the end to cover the time-dependent lines
+      # of the indicator and the MPI-specific AMR code.
+      coverage_override = (maxiters=10^5,))
   end
 
   @trixi_testset "elixir_euler_density_wave.jl" begin
@@ -147,7 +154,8 @@ const EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_vortex_amr.jl"),
       # Expected errors are exactly the same as in the serial test!
       l2   = [2.120552206480055e-6, 0.003281541473561042, 0.003280625257336616, 0.004645872821313438],
-      linf = [4.500266027052113e-5, 0.031765399304366726, 0.03179340562764421, 0.04563622772500864])
+      linf = [4.500266027052113e-5, 0.031765399304366726, 0.03179340562764421, 0.04563622772500864],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_vortex_shockcapturing.jl" begin

test/test_p4est_2d.jl

@@ -35,13 +35,15 @@ isdir(outdir) && rm(outdir, recursive=true)
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_solution_independent.jl"),
       # Expected errors are exactly the same as with StructuredMesh!
       l2   = [4.949660644033807e-5],
-      linf = [0.0004867846262313763])
+      linf = [0.0004867846262313763],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_amr_unstructured_flag.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_unstructured_flag.jl"),
       l2   = [0.0012766060609964525],
-      linf = [0.01750280631586159])
+      linf = [0.01750280631586159],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_restart.jl" begin
@@ -82,7 +84,8 @@ isdir(outdir) && rm(outdir, recursive=true)
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_blast_wave_amr.jl"),
       l2   = [6.32183914e-01, 3.86914231e-01, 3.86869171e-01, 1.06575688e+00],
       linf = [2.76020890e+00, 2.32659890e+00, 2.32580837e+00, 2.15778188e+00],
-      tspan = (0.0, 0.3))
+      tspan = (0.0, 0.3),
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_wall_bc_amr.jl" begin

test/test_p4est_3d.jl

@@ -35,14 +35,16 @@ isdir(outdir) && rm(outdir, recursive=true)
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr.jl"),
       # Expected errors are exactly the same as with TreeMesh!
       l2   = [9.773852895157622e-6],
-      linf = [0.0005853874124926162])
+      linf = [0.0005853874124926162],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_amr_unstructured_curved.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_unstructured_curved.jl"),
       l2   = [1.6236411810065552e-5],
       linf = [0.0010554006923731395],
-      tspan = (0.0, 1.0))
+      tspan = (0.0, 1.0),
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_cubed_sphere.jl" begin
@@ -89,14 +91,16 @@ isdir(outdir) && rm(outdir, recursive=true)
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_ec.jl"),
       l2   = [0.010380390326164493, 0.006192950051354618, 0.005970674274073704, 0.005965831290564327, 0.02628875593094754],
       linf = [0.3326911600075694, 0.2824952141320467, 0.41401037398065543, 0.45574161423218573, 0.8099577682187109],
-      tspan = (0.0, 0.2))
+      tspan = (0.0, 0.2),
+      coverage_override = (polydeg=3,)) # Prevent long compile time in CI
   end
 
   @trixi_testset "elixir_euler_sedov.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_sedov.jl"),
       l2   = [7.82070951e-02, 4.33260474e-02, 4.33260474e-02, 4.33260474e-02, 3.75260911e-01],
       linf = [7.45329845e-01, 3.21754792e-01, 3.21754792e-01, 3.21754792e-01, 4.76151527e+00],
-      tspan = (0.0, 0.3))
+      tspan = (0.0, 0.3),
+      coverage_override = (polydeg=3,)) # Prevent long compile time in CI
   end
 
   @trixi_testset "elixir_mhd_alfven_wave_nonconforming.jl" begin
@@ -118,7 +122,8 @@ isdir(outdir) && rm(outdir, recursive=true)
       linf = [0.20904054009050665, 0.18622917151105936, 0.2347957890323218, 0.19432508025509926,
               0.6858860133405615, 0.15172116633332622, 0.22432820727833747, 0.16805989780225183,
               0.000535219040687628],
-      tspan = (0.0, 0.04))
+      tspan = (0.0, 0.04),
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_source_terms_nonperiodic_hohqmesh.jl" begin

test/test_paper_self_gravitating_gas_dynamics.jl

@@ -116,7 +116,8 @@ const EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_eulergravity_sedov_blast_wave.jl"),
       l2   = [0.046315994852653024, 0.0650818006233669, 0.06508180062336677, 0.4896707211656037],
       linf = [2.3874843337593776, 4.07876384374792, 4.07876384374792, 16.23914384809855],
-      tspan = (0.0, 0.05))
+      tspan = (0.0, 0.05),
+      coverage_override = (maxiters=2,))
   end
 
   @trixi_testset "elixir_eulergravity_sedov_blast_wave.jl with ref-level=8 and no AMR" begin

test/test_special_elixirs.jl

@@ -15,39 +15,50 @@ isdir(outdir) && rm(outdir, recursive=true)
 # pathof(Trixi) returns /path/to/Trixi/src/Trixi.jl, dirname gives the parent directory
 const EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples")
 
+cmd = string(Base.julia_cmd())
+coverage = occursin("--code-coverage", cmd) && !occursin("--code-coverage=none", cmd)
 
 @testset "Special elixirs" begin
   @testset "Convergence test" begin
-    @timed_testset "tree_2d_dgsem" begin
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "tree_2d_dgsem", "elixir_advection_extended.jl"), 3, initial_refinement_level=2)
-      @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
-    end
+    if !coverage
+      @timed_testset "tree_2d_dgsem" begin
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "tree_2d_dgsem", "elixir_advection_extended.jl"), 3, initial_refinement_level=2)
+        @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
+      end
 
-    @timed_testset "structured_2d_dgsem" begin
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_2d_dgsem", "elixir_advection_extended.jl"), 3, cells_per_dimension=(5, 9))
-      @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
-    end
+      @timed_testset "structured_2d_dgsem" begin
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_2d_dgsem", "elixir_advection_extended.jl"), 3, cells_per_dimension=(5, 9))
+        @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
+      end
 
-    @timed_testset "p4est_2d_dgsem" begin
-      # Run convergence test on unrefined mesh
-      no_refine = @cfunction((p4est, which_tree, quadrant) -> Cint(0), Cint, (Ptr{Trixi.p4est_t}, Ptr{Trixi.p4est_topidx_t}, Ptr{Trixi.p4est_quadrant_t}))
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "p4est_2d_dgsem", "elixir_euler_source_terms_nonconforming_unstructured_flag.jl"), 2, refine_fn_c=no_refine)
-      @test isapprox(mean_convergence[:linf], [3.2, 3.2, 4.0, 3.7], rtol=0.05)
-    end
+      @timed_testset "p4est_2d_dgsem" begin
+        # Run convergence test on unrefined mesh
+        no_refine = @cfunction((p4est, which_tree, quadrant) -> Cint(0), Cint, (Ptr{Trixi.p4est_t}, Ptr{Trixi.p4est_topidx_t}, Ptr{Trixi.p4est_quadrant_t}))
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "p4est_2d_dgsem", "elixir_euler_source_terms_nonconforming_unstructured_flag.jl"), 2, refine_fn_c=no_refine)
+        @test isapprox(mean_convergence[:linf], [3.2, 3.2, 4.0, 3.7], rtol=0.05)
+      end
 
-    @timed_testset "structured_3d_dgsem" begin
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_3d_dgsem", "elixir_advection_basic.jl"), 2, cells_per_dimension=(7, 4, 5))
-      @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
-    end
+      @timed_testset "structured_3d_dgsem" begin
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_3d_dgsem", "elixir_advection_basic.jl"), 2, cells_per_dimension=(7, 4, 5))
+        @test isapprox(mean_convergence[:l2], [4.0], rtol=0.05)
+      end
 
-    @timed_testset "p4est_3d_dgsem" begin
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "p4est_3d_dgsem", "elixir_advection_unstructured_curved.jl"), 2, initial_refinement_level=0)
-      @test isapprox(mean_convergence[:l2], [2.7], rtol=0.05)
-    end
+      @timed_testset "p4est_3d_dgsem" begin
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "p4est_3d_dgsem", "elixir_advection_unstructured_curved.jl"), 2, initial_refinement_level=0)
+        @test isapprox(mean_convergence[:l2], [2.7], rtol=0.05)
+      end
 
-    @timed_testset "paper_self_gravitating_gas_dynamics" begin
-      mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "paper_self_gravitating_gas_dynamics", "elixir_eulergravity_convergence.jl"), 2, tspan=(0.0, 0.25), initial_refinement_level=1)
-      @test isapprox(mean_convergence[:l2], 4 * ones(4), atol=0.4)
+      @timed_testset "paper_self_gravitating_gas_dynamics" begin
+        mean_convergence = convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "paper_self_gravitating_gas_dynamics", "elixir_eulergravity_convergence.jl"), 2, tspan=(0.0, 0.25), initial_refinement_level=1)
+        @test isapprox(mean_convergence[:l2], 4 * ones(4), atol=0.4)
+      end
+    else
+      # Without coverage, just run simple convergence tests to cover
+      # the convergence test logic
+      @test_nowarn convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "tree_2d_dgsem", "elixir_advection_basic.jl"), 2, tspan=(0.0, 0.01))
+      @test_nowarn convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "tree_2d_dgsem", "elixir_advection_extended.jl"), 2, initial_refinement_level=0, tspan=(0.0, 0.1))
+      @test_nowarn convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_2d_dgsem", "elixir_advection_basic.jl"), 2, tspan=(0.0, 0.01))
+      @test_nowarn convergence_test(@__MODULE__, joinpath(EXAMPLES_DIR, "structured_2d_dgsem", "elixir_advection_extended.jl"), 2, cells_per_dimension=(1, 1), tspan=(0.0, 0.1))
     end
   end
 

test/test_structured_3d.jl

@@ -98,7 +98,9 @@ isdir(outdir) && rm(outdir, recursive=true)
               0.0013647609788548722],
       linf = [0.027510637768610846, 0.02797062834945721, 0.01274249949295704, 0.038940694415543736,
               0.02200825678588325, 0.03167600959583505, 0.021420957993862344, 0.03386589835999665,
-              0.01888303191983353])
+              0.01888303191983353],
+      # Use same polydeg as everything else to prevent long compile times in CI
+      coverage_override = (polydeg=3,))
   end
 
   @trixi_testset "elixir_mhd_alfven_wave.jl with flux_lax_friedrichs" begin
@@ -109,7 +111,9 @@ isdir(outdir) && rm(outdir, recursive=true)
       linf = [0.027719103797310463, 0.027570111789910784, 0.012561901006903103, 0.03903568568480584,
               0.021311996934554767, 0.03154849824135775, 0.020996033645485412, 0.03403185137382961,
               0.019488952445771597],
-      surface_flux = (flux_lax_friedrichs, flux_nonconservative_powell))
+      surface_flux = (flux_lax_friedrichs, flux_nonconservative_powell),
+      # Use same polydeg as everything else to prevent long compile times in CI
+      coverage_override = (polydeg=3,))
   end
 
   @trixi_testset "elixir_mhd_ec_shockcapturing.jl" begin
@@ -120,7 +124,9 @@ isdir(outdir) && rm(outdir, recursive=true)
       linf = [0.30749098250807516, 0.2679008863509767, 0.271243087484388, 0.26545396569129537,
               0.9620950892188596, 0.18163281157498123, 0.15995708312378454, 0.17918221526906408,
               0.015138346608166353],
-      tspan = (0.0, 0.25))
+      tspan = (0.0, 0.25),
+      # Use same polydeg as everything else to prevent long compile times in CI
+      coverage_override = (polydeg=3,))
   end
 end
 

test/test_tree_1d_advection.jl

@@ -18,13 +18,15 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_1
   @trixi_testset "elixir_advection_amr.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr.jl"),
       l2   = [0.3540206249507417],
-      linf = [0.9999896603382347])
+      linf = [0.9999896603382347],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_advection_amr_nonperiodic.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_amr_nonperiodic.jl"),
       l2   = [4.283508859843524e-6],
-      linf = [3.235356127918171e-5])
+      linf = [3.235356127918171e-5],
+      coverage_override = (maxiters=6,))
   end
 end
 

test/test_tree_1d_euler.jl

@@ -91,35 +91,41 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_1
   @trixi_testset "elixir_euler_sedov_blast_wave.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_sedov_blast_wave.jl"),
       l2   = [1.250005061244617, 0.06878411345533507, 0.9264328311018613],
-      linf = [2.9766770877037168, 0.16838100902295852, 2.6655773445485798])
+      linf = [2.9766770877037168, 0.16838100902295852, 2.6655773445485798],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_sedov_blast_wave_pure_fv.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_sedov_blast_wave_pure_fv.jl"),
       l2   = [1.075075094036344, 0.06766902169711514, 0.9221426570128292],
-      linf = [3.3941512671408542, 0.16862631133303882, 2.6572394126490315])
+      linf = [3.3941512671408542, 0.16862631133303882, 2.6572394126490315],
+      # Let this test run longer to cover some lines in flux_hllc
+      coverage_override = (maxiters=10^5, tspan=(0.0, 0.1)))
   end
 
   @trixi_testset "elixir_euler_sedov_blast_wave.jl with pressure" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_sedov_blast_wave.jl"),
       l2   = [1.297525985166995, 0.07964929522694145, 0.9269991156246368],
       linf = [3.1773015255764427, 0.21331831536493773, 2.6650170188241047],
       shock_indicator_variable = pressure,
-      cfl = 0.2)
+      cfl = 0.2,
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_sedov_blast_wave.jl with density" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_sedov_blast_wave.jl"),
       l2   = [1.2798798835860528, 0.07103461242058921, 0.9273792517187003],
       linf = [3.1087017048015824, 0.17734706962928956, 2.666689753470263],
       shock_indicator_variable = density,
-      cfl = 0.2)
+      cfl = 0.2,
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_euler_positivity.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_positivity.jl"),
       l2   = [1.6493820253458906, 0.19793887460986834, 0.9783506076125921],
-      linf = [4.71751203912051, 0.5272411022735763, 2.7426163947635844])
+      linf = [4.71751203912051, 0.5272411022735763, 2.7426163947635844],
+      coverage_override = (maxiters=3,))
   end
 
   @trixi_testset "elixir_euler_blast_wave.jl" begin

test/test_tree_1d_mhd.jl

@@ -40,7 +40,8 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_1
   @trixi_testset "elixir_mhd_briowu_shock_tube.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_mhd_briowu_shock_tube.jl"),
       l2   = [0.17764301067932906, 0.19693621875378622, 0.3635136528288642, 0.0, 0.3757321708837591, 8.593007507325741e-16, 0.36473438378159656, 0.0],
-      linf = [0.5601530250396535, 0.43867368105486537, 1.0960903616351099, 0.0, 1.0551794137886303, 4.107825191113079e-15, 1.5374410890043144, 0.0])
+      linf = [0.5601530250396535, 0.43867368105486537, 1.0960903616351099, 0.0, 1.0551794137886303, 4.107825191113079e-15, 1.5374410890043144, 0.0],
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_mhd_torrilhon_shock_tube.jl" begin
@@ -60,7 +61,8 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_1
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_mhd_shu_osher_shock_tube.jl"),
       l2   = [1.01126210e+00, 8.27157902e+00, 1.30882545e+00, 0.00000000e+00, 5.21930435e+01, 6.56538824e-16, 1.01022340e+00, 0.00000000e+00],
       linf = [2.87172004e+00, 2.26438057e+01, 4.16672442e+00, 0.00000000e+00, 1.35152372e+02, 3.44169138e-15, 2.83556069e+00, 0.00000000e+00],
-      tspan = (0.0, 0.2))
+      tspan = (0.0, 0.2),
+      coverage_override = (maxiters=6,))
   end
 
   @trixi_testset "elixir_mhd_shu_osher_shock_tube.jl with flipped shock direction" begin
@@ -71,7 +73,8 @@ EXAMPLES_DIR = joinpath(pathof(Trixi) |> dirname |> dirname, "examples", "tree_1
       linf = [2.92876280e+00, 2.28341581e+01, 4.11643561e+00, 0.00000000e+00, 1.36966213e+02, 1.55431223e-15, 2.80548864e+00, 0.00000000e+00],
       initial_condition = initial_condition_shu_osher_shock_tube_flipped,
       boundary_conditions=BoundaryConditionDirichlet(initial_condition_shu_osher_shock_tube_flipped),
-      tspan = (0.0, 0.2))
+      tspan = (0.0, 0.2),
+      coverage_override = (maxiters=6,))
   end
 end