Skip to content

Commit

Permalink
Remove initial_condition; Disable advection_basic tests
Browse files Browse the repository at this point in the history
  • Loading branch information
bennibolm committed Oct 15, 2024
1 parent de7b527 commit 92e1d3b
Show file tree
Hide file tree
Showing 3 changed files with 75 additions and 90 deletions.
2 changes: 1 addition & 1 deletion src/solvers/fv_t8code/fv.jl
Original file line number Diff line number Diff line change
Expand Up @@ -160,7 +160,7 @@ function create_cache(mesh::T8codeMesh, equations::AbstractEquations, solver::FV
end

function rhs!(du, u, t, mesh::T8codeMesh, equations,
initial_condition, boundary_conditions, source_terms::Source,
boundary_conditions, source_terms::Source,
solver::FV, cache) where {Source}
# Reset du
@trixi_timeit timer() "reset ∂u/∂t" du.=zero(eltype(du))
Expand Down
77 changes: 31 additions & 46 deletions test/test_mpi_t8code_fv.jl
Original file line number Diff line number Diff line change
Expand Up @@ -13,52 +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.
# @trixi_testset "second-order FV, extended reconstruction stencil" begin
# @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"),
# order=2,
# extended_reconstruction_stencil=true,
# l2=[0.020331012873518642],
# linf=[0.05571209803860677])
# # 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
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
Expand Down
86 changes: 43 additions & 43 deletions test/test_t8code_fv.jl
Original file line number Diff line number Diff line change
Expand Up @@ -44,49 +44,49 @@ mkdir(outdir)
# end
# 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
@trixi_testset "second-order FV, extended reconstruction stencil" begin
@test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"),
extended_reconstruction_stencil=true,
l2=[0.028436326251639936],
linf=[0.08696815845435057])
# 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
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
# @trixi_testset "second-order FV, extended reconstruction stencil" begin
# @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_advection_basic.jl"),
# extended_reconstruction_stencil=true,
# l2=[0.028436326251639936],
# linf=[0.08696815845435057])
# # 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
# end

@trixi_testset "elixir_advection_gauss.jl" begin
@trixi_testset "first-order FV" begin
Expand Down

0 comments on commit 92e1d3b

Please sign in to comment.