reconfigure CI testing to actually test the submodules in isolation

Author: ChrisRackauckas

.github/workflows/CI.yml

@@ -24,16 +24,36 @@ jobs:
           - Downstream
           - ODEInterfaceRegression
           - Multithreading
-          - FIRK
-          - Symplectic
-          - Extrapolation
-          - Feagin
-          - StabilizedRK
-          - StabilizedIRK
-          - SSPRK
-          - LowStorageRK
-          - QPRK
-          - Linear
+
+          - OrdinaryDiffEqAdamsBashforthMoulton
+          - OrdinaryDiffEqBDF
+          - OrdinaryDiffEqDefault
+          - OrdinaryDiffEqDifferentiation
+          - OrdinaryDiffEqExplicitRK
+          - OrdinaryDiffEqExponentialRK
+          - OrdinaryDiffEqExtrapolation
+          - OrdinaryDiffEqFIRK
+          - OrdinaryDiffEqFeagin
+          - OrdinaryDiffEqFunctionMap
+          - OrdinaryDiffEqHighOrderRK
+          - OrdinaryDiffEqIMEXMultistep
+          - OrdinaryDiffEqLinear
+          - OrdinaryDiffEqLowOrderRK
+          - OrdinaryDiffEqLowStorageRK
+          - OrdinaryDiffEqNordsieck
+          - OrdinaryDiffEqPDIRK
+          - OrdinaryDiffEqPRK
+          - OrdinaryDiffEqQPRK
+          - OrdinaryDiffEqRKN
+          - OrdinaryDiffEqRosenbrock
+          - OrdinaryDiffEqSDIRK
+          - OrdinaryDiffEqSSPRK
+          - OrdinaryDiffEqStabilizedIRK
+          - OrdinaryDiffEqStabilizedRK
+          - OrdinaryDiffEqSymplecticRK
+          - OrdinaryDiffEqTsit5
+          - OrdinaryDiffEqVerner
+
         version:
           - '1'
     steps:

lib/OrdinaryDiffEqAdamsBashforthMoulton/test/runtests.jl

@@ -1,3 +1,3 @@
 using SafeTestsets
 
-@time @safetestset "Adams Variable Coefficients Tests" include("adams_tests.jl")
+@time @safetestset "Adams Variable Coefficients Tests" include("adams_tests.jl")

lib/OrdinaryDiffEqBDF/test/dae_ad_tests.jl

@@ -1,30 +1,5 @@
 using OrdinaryDiffEq, LinearAlgebra, ForwardDiff, Test
 
-function rober(du, u, p, t)
-    y₁, y₂, y₃ = u
-    k₁, k₂, k₃ = p
-    du[1] = -k₁ * y₁ + k₃ * y₂ * y₃
-    du[2] = k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2
-    du[3] = y₁ + y₂ + y₃ - 1
-    nothing
-end
-function rober(u, p, t)
-    y₁, y₂, y₃ = u
-    k₁, k₂, k₃ = p
-    [-k₁ * y₁ + k₃ * y₂ * y₃,
-        k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2,
-        y₁ + y₂ + y₃ - 1]
-end
-M = [1.0 0 0
-     0 1.0 0
-     0 0 0]
-roberf = ODEFunction(rober, mass_matrix = M)
-roberf_oop = ODEFunction{false}(rober, mass_matrix = M)
-prob_mm = ODEProblem(roberf, [1.0, 0.0, 0.2], (0.0, 1e5), (0.04, 3e7, 1e4))
-prob_mm_oop = ODEProblem(roberf_oop, [1.0, 0.0, 0.2], (0.0, 1e5), (0.04, 3e7, 1e4))
-sol = solve(prob_mm, Rodas5P(), reltol = 1e-8, abstol = 1e-8)
-sol = solve(prob_mm_oop, Rodas5P(), reltol = 1e-8, abstol = 1e-8)
-
 function f(out, du, u, p, t)
     out[1] = -p[1] * u[1] + p[3] * u[2] * u[3] - du[1]
     out[2] = +p[1] * u[1] - p[2] * u[2]^2 - p[3] * u[2] * u[3] - du[2]
@@ -47,10 +22,10 @@ sol1 = solve(prob, DFBDF(), dt = 1e-5, abstol = 1e-8, reltol = 1e-8)
 # These tests flex differentiation of the solver and through the initialization
 # To only test the solver part and isolate potential issues, set the initialization to consistent
 @testset "Inplace: $(isinplace(_prob)), DAEProblem: $(_prob isa DAEProblem), BrownBasic: $(initalg isa BrownFullBasicInit), Autodiff: $autodiff" for _prob in [
-        prob, prob_mm, prob_oop, prob_mm_oop],
+        prob, prob_oop],
     initalg in [BrownFullBasicInit(), ShampineCollocationInit()], autodiff in [true, false]
 
-    alg = _prob isa DAEProblem ? DFBDF(; autodiff) : Rodas5P(; autodiff)
+    alg = DFBDF(; autodiff)
     function f(p)
         sol = solve(remake(_prob, p = p), alg, abstol = 1e-14,
             reltol = 1e-14, initializealg = initalg)

lib/OrdinaryDiffEqBDF/test/dae_initialization_tests.jl

@@ -1,102 +1,5 @@
 using OrdinaryDiffEq, StaticArrays, LinearAlgebra, Test
 
-## Mass Matrix
-
-function rober_oop(u, p, t)
-    y₁, y₂, y₃ = u
-    k₁, k₂, k₃ = p
-    du1 = -k₁ * y₁ + k₃ * y₂ * y₃
-    du2 = k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2
-    du3 = y₁ + y₂ + y₃ - 1
-    [du1, du2, du3]
-end
-M = [1.0 0 0
-     0 1.0 0
-     0 0 0]
-f_oop = ODEFunction(rober_oop, mass_matrix = M)
-prob_mm = ODEProblem(f_oop, [1.0, 0.0, 0.0], (0.0, 1e5), (0.04, 3e7, 1e4))
-sol = solve(prob_mm, Rosenbrock23(autodiff = false), reltol = 1e-8, abstol = 1e-8)
-@test sol[1] == [1.0, 0.0, 0.0] # Ensure initialization is unchanged if it works at the start!
-sol = solve(prob_mm, Rosenbrock23(), reltol = 1e-8, abstol = 1e-8,
-    initializealg = ShampineCollocationInit())
-@test sol[1] == [1.0, 0.0, 0.0] # Ensure initialization is unchanged if it works at the start!
-
-prob_mm = ODEProblem(f_oop, [1.0, 0.0, 0.2], (0.0, 1e5), (0.04, 3e7, 1e4))
-sol = solve(prob_mm, Rosenbrock23(), reltol = 1e-8, abstol = 1e-8)
-@test sum(sol[1]) ≈ 1
-@test sol[1] ≈ [1.0, 0.0, 0.0]
-for alg in [Rosenbrock23(autodiff = false), Trapezoid()]
-    local sol
-    sol = solve(prob_mm, alg, reltol = 1e-8, abstol = 1e-8,
-        initializealg = ShampineCollocationInit())
-    @test sum(sol[1]) ≈ 1
-end
-
-function rober(du, u, p, t)
-    y₁, y₂, y₃ = u
-    k₁, k₂, k₃ = p
-    du[1] = -k₁ * y₁ + k₃ * y₂ * y₃
-    du[2] = k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2
-    du[3] = y₁ + y₂ + y₃ - 1
-    nothing
-end
-M = [1.0 0 0
-     0 1.0 0
-     0 0 0]
-f = ODEFunction(rober, mass_matrix = M)
-prob_mm = ODEProblem(f, [1.0, 0.0, 0.0], (0.0, 1e5), (0.04, 3e7, 1e4))
-sol = solve(prob_mm, Rodas5(autodiff = false), reltol = 1e-8, abstol = 1e-8)
-@test sol[1] == [1.0, 0.0, 0.0] # Ensure initialization is unchanged if it works at the start!
-sol = solve(prob_mm, Rodas5(), reltol = 1e-8, abstol = 1e-8,
-    initializealg = ShampineCollocationInit())
-@test sol[1] == [1.0, 0.0, 0.0] # Ensure initialization is unchanged if it works at the start!
-
-prob_mm = ODEProblem(f, [1.0, 0.0, 1.0], (0.0, 1e5), (0.04, 3e7, 1e4))
-sol = solve(prob_mm, Rodas5(), reltol = 1e-8, abstol = 1e-8)
-@test sum(sol[1]) ≈ 1
-@test sol[1] ≈ [1.0, 0.0, 0.0]
-
-for alg in [Rodas5(autodiff = false), Trapezoid()]
-    local sol
-    sol = solve(prob_mm, alg, reltol = 1e-8, abstol = 1e-8,
-        initializealg = ShampineCollocationInit())
-    @test sum(sol[1]) ≈ 1
-end
-
-function rober_no_p(du, u, p, t)
-    y₁, y₂, y₃ = u
-    (k₁, k₂, k₃) = (0.04, 3e7, 1e4)
-    du[1] = -k₁ * y₁ + k₃ * y₂ * y₃
-    du[2] = k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2
-    du[3] = y₁ + y₂ + y₃ - 1
-    nothing
-end
-
-function rober_oop_no_p(du, u, p, t)
-    y₁, y₂, y₃ = u
-    (k₁, k₂, k₃) = (0.04, 3e7, 1e4)
-    du1 = -k₁ * y₁ + k₃ * y₂ * y₃
-    du2 = k₁ * y₁ - k₃ * y₂ * y₃ - k₂ * y₂^2
-    du3 = y₁ + y₂ + y₃ - 1
-    [du1, du2, du3]
-end
-
-# test oop and iip ODE initialization with parameters without eltype/length
-struct UnusedParam
-end
-for f in (
-    ODEFunction(rober_no_p, mass_matrix = M), ODEFunction(rober_oop_no_p, mass_matrix = M))
-    local prob, probp
-    prob = ODEProblem(f, [1.0, 0.0, 1.0], (0.0, 1e5))
-    probp = ODEProblem(f, [1.0, 0.0, 1.0], (0.0, 1e5), UnusedParam)
-    for initializealg in (ShampineCollocationInit(), BrownFullBasicInit())
-        isapprox(init(prob, Rodas5(), abstol = 1e-10; initializealg).u,
-            init(prob, Rodas5(), abstol = 1e-10; initializealg).u)
-    end
-end
-
-## DAEProblem
-
 f = function (du, u, p, t)
     out1 = -0.04u[1] + 1e4 * u[2] * u[3] - du[1]
     out2 = +0.04u[1] - 3e7 * u[2]^2 - 1e4 * u[2] * u[3] - du[2]
@@ -246,23 +149,3 @@ for initializealg in (ShampineCollocationInit(), BrownFullBasicInit())
     @test isapprox(
         init(probp, DABDF2(); initializealg).u, init(prob, DABDF2(); initializealg).u)
 end
-
-# to test that we get the right NL solve we need a broken solver.
-struct BrokenNLSolve <: SciMLBase.AbstractNonlinearAlgorithm
-    BrokenNLSolve(; kwargs...) = new()
-end
-function SciMLBase.__solve(prob::NonlinearProblem,
-        alg::BrokenNLSolve, args...;
-        kwargs...)
-    u = fill(reinterpret(Float64, 0xDEADBEEFDEADBEEF), 3)
-    SciMLBase.build_solution(prob, alg, u, copy(u);
-        retcode = ReturnCode.Success)
-end
-function f2(u, p, t)
-    u
-end
-f = ODEFunction(f2, mass_matrix = Diagonal([1.0, 1.0, 0.0]))
-prob = ODEProblem(f, ones(3), (0.0, 1.0))
-integrator = init(prob, Rodas5P(),
-    initializealg = ShampineCollocationInit(1.0, BrokenNLSolve()))
-@test all(isequal(reinterpret(Float64, 0xDEADBEEFDEADBEEF)), integrator.u)

lib/OrdinaryDiffEqBDF/test/runtest.jl

@@ -0,0 +1,6 @@
+using SafeTestsets
+
+@time @safetestset "DAE Convergence Tests" include("dae_convergence_tests.jl")
+@time @safetestset "DAE AD Tests" include("dae_ad_tests.jl")
+@time @safetestset "DAE Event Tests" include("dae_event.jl")
+@time @safetestset "DAE Initialization Tests" include("dae_initialization_tests.jl")

lib/OrdinaryDiffEqBDF/test/runtests.jl

@@ -1,3 +1,4 @@
 using SafeTestsets
 
-@time @safetestset "Linear-Nonlinear Krylov Methods Tests" include("linear_nonlinear_krylov_tests.jl")
+@time @safetestset "Linear-Nonlinear Krylov Methods Tests" include("linear_nonlinear_krylov_tests.jl")
+@time @safetestset "Linear-Nonlinear Convergence Tests" include("linear_nonlinear_convergence_tests.jl")