diff --git a/test/interface/algebraic_interpolation.jl b/test/interface/algebraic_interpolation.jl index 362a5c4071..78bf090997 100644 --- a/test/interface/algebraic_interpolation.jl +++ b/test/interface/algebraic_interpolation.jl @@ -1,4 +1,5 @@ using OrdinaryDiffEq +using DiffEqBase using Test function rober_ip(du, u, p, t) y₁, y₂, y₃ = u @@ -32,6 +33,13 @@ ref_op = solve(prob_op, Rodas5P(), reltol = 1e-8, abstol = 1e-8) sol_ip = solve(prob_ip, FBDF(), reltol = 1e-8, abstol = 1e-8) sol_op = solve(prob_op, FBDF(), reltol = 1e-8, abstol = 1e-8) +# make sure interpolation changes don't accidentally break this test suite +# the intention is that ref uses a stiffness-aware interpolation, while sol uses hermite +@test DiffEqBase.interp_summary(sol_ip) == "3rd order Hermite" +@test DiffEqBase.interp_summary(sol_op) == "3rd order Hermite" +@test occursin("stiffness-aware", DiffEqBase.interp_summary(ref_ip)) +@test occursin("stiffness-aware", DiffEqBase.interp_summary(ref_op)) + reltol = 1e-4 abstol = 1e-4 t = 1