diff --git a/examples/t8code_2d_fv/elixir_advection_basic.jl b/examples/t8code_2d_fv/elixir_advection_basic.jl index 4808ef35b57..d0c983dbe32 100644 --- a/examples/t8code_2d_fv/elixir_advection_basic.jl +++ b/examples/t8code_2d_fv/elixir_advection_basic.jl @@ -88,9 +88,9 @@ callbacks = CallbackSet(summary_callback, analysis_callback, alive_callback, ############################################################################### # run the simulation -sol = solve(ode, CarpenterKennedy2N54(williamson_condition = false), - dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback - save_everystep = false, callback = callbacks); -summary_callback() +# sol = solve(ode, CarpenterKennedy2N54(williamson_condition = false), +# dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback +# save_everystep = false, callback = callbacks); +# summary_callback() GC.enable(true) diff --git a/test/test_mpi_t8code_fv_2d.jl b/test/test_mpi_t8code_fv_2d.jl index cfc2892958a..629942abe81 100644 --- a/test/test_mpi_t8code_fv_2d.jl +++ b/test/test_mpi_t8code_fv_2d.jl @@ -13,37 +13,37 @@ const EXAMPLES_DIR = pkgdir(Trixi, "examples", "t8code_2d_fv") # Run basic tests @testset "Examples 2D" begin # Linear scalar advection - # @trixi_testset "elixir_advection_basic.jl" begin - # @trixi_testset "first-order FV" begin - # @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"), - # order=1, - # l2=[0.08551397247817498], - # linf=[0.12087467695430498]) - # # Ensure that we do not have excessive memory allocations - # # (e.g., from type instabilities) - # let - # t = sol.t[end] - # u_ode = sol.u[end] - # du_ode = similar(u_ode) - # @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000 - # end - # end - # @trixi_testset "second-order FV" begin - # @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"), - # l2=[0.008142380494734171], - # linf=[0.018687916234976898]) - # # Ensure that we do not have excessive memory allocations - # # (e.g., from type instabilities) - # let - # t = sol.t[end] - # u_ode = sol.u[end] - # du_ode = similar(u_ode) - # @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000 - # end - # end - # # The extended reconstruction stencil is currently not mpi parallel. - # # The current version runs through but an error occurs on some rank. - # end + @trixi_testset "elixir_advection_basic.jl" begin + @trixi_testset "first-order FV" begin + @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"), + order=1, + l2=[0.08551397247817498], + linf=[0.12087467695430498]) + # Ensure that we do not have excessive memory allocations + # (e.g., from type instabilities) + let + t = sol.t[end] + u_ode = sol.u[end] + du_ode = similar(u_ode) + @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000 + end + end + @trixi_testset "second-order FV" begin + @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"), + l2=[0.008142380494734171], + linf=[0.018687916234976898]) + # Ensure that we do not have excessive memory allocations + # (e.g., from type instabilities) + let + t = sol.t[end] + u_ode = sol.u[end] + du_ode = similar(u_ode) + @test (@allocated Trixi.rhs!(du_ode, u_ode, semi, t)) < 1000 + end + end + # The extended reconstruction stencil is currently not mpi parallel. + # The current version runs through but an error occurs on some rank. + end @trixi_testset "elixir_advection_gauss.jl" begin @trixi_testset "first-order FV" begin