Skip to content

Commit

Permalink
Reduce cfl and show eps convergence for mono4, mono5
Browse files Browse the repository at this point in the history
  • Loading branch information
Akshay Sridhar authored and charleskawczynski committed Dec 19, 2024
1 parent 021490b commit e0c0184
Showing 1 changed file with 7 additions and 4 deletions.
11 changes: 7 additions & 4 deletions examples/column/vanleer_advection.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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₁)
Expand Down

0 comments on commit e0c0184

Please sign in to comment.