From a54026e7bd22b0d53d61a355d231558c6d04b3dd Mon Sep 17 00:00:00 2001 From: Akshay Sridhar Date: Thu, 19 Dec 2024 11:03:36 -0500 Subject: [PATCH] Reduce cfl and show eps convergence for mono4, mono5 --- examples/column/vanleer_advection.jl | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/examples/column/vanleer_advection.jl b/examples/column/vanleer_advection.jl index 570ef91673..159e4dbf46 100644 --- a/examples/column/vanleer_advection.jl +++ b/examples/column/vanleer_advection.jl @@ -116,13 +116,16 @@ for (i, stretch_fn) in enumerate(stretch_fns) if limiter_method != Operators.AlgebraicMean() && stretch_fn == Meshes.Uniform() @info "Extrema with $(limiter_method): $(extrema(q_final))" @assert maximum(q_final) <= FT(1) - @assert isapprox(maximum(q_final .- 1), FT(0), atol = FT(1e-6)) - @assert isapprox(minimum(q_final .- 0), FT(0), atol = FT(1e-6)) + @show maximum(q_final .- 1) + @show minimum(q_final .- 0) + @show abs(maximum(q_final .- 1)) + @assert abs(maximum(q_final .- 1)) <= 10*sqrt(eps(FT)) + @assert abs(minimum(q_final .- 0)) <= 10*sqrt(eps(FT)) elseif limiter_method != Operators.AlgebraicMean() @info "Extrema with $(limiter_method): $(extrema(q_final))" @assert maximum(q_final) <= FT(1) - @assert isapprox(maximum(q_final .- 1), FT(0), atol = FT(0.05)) - @assert isapprox(minimum(q_final .- 0), FT(0), atol = FT(0.05)) + @assert abs(maximum(q_final .- 1)) <= FT(0.05) + @assert abs(minimum(q_final .- 0)) <= FT(0.05) end q_analytic = pulse.(z, t₁, z₀, zₕ, z₁)