Merge pull request #310 from SciML/qqy/twp_bigfloat

Author: ErikQQY

Project.toml

@@ -58,6 +58,7 @@ Random = "1.10"
 ReTestItems = "1.29"
 Reexport = "1.2"
 SciMLBase = "2.82"
+Sparspak = "0.3.11"
 StaticArrays = "1.9.8"
 Test = "1.10"
 julia = "1.10"
@@ -76,8 +77,9 @@ Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 ReTestItems = "817f1d60-ba6b-4fd5-9520-3cf149f6a823"
 RecursiveArrayTools = "731186ca-8d62-57ce-b412-fbd966d074cd"
+Sparspak = "e56a9233-b9d6-4f03-8d0f-1825330902ac"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [targets]
-test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "NonlinearSolveFirstOrder", "ODEInterface", "OrdinaryDiffEq", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "StaticArrays", "Test"]
+test = ["Aqua", "DiffEqDevTools", "Hwloc", "InteractiveUtils", "JET", "LinearSolve", "NonlinearSolveFirstOrder", "ODEInterface", "OrdinaryDiffEq", "Pkg", "Random", "ReTestItems", "RecursiveArrayTools", "Sparspak", "StaticArrays", "Test"]

test/misc/bigfloat_test.jl

@@ -1,5 +1,8 @@
 @testitem "BigFloat compatibility" begin
     using BoundaryValueDiffEq
+    # Need Sparspak for BigFloat
+    using Sparspak
+
     tspan = (0.0, pi / 2)
     function simplependulum!(du, u, p, t)
         θ = u[1]
@@ -12,15 +15,36 @@
         residual[2] = u(pi / 2)[1] - big(pi / 2)
     end
     u0 = BigFloat.([pi / 2, pi / 2])
-    prob = BVProblem(simplependulum!, bc!, u0, tspan)
-    sol1 = solve(prob, MIRK4(), dt = 0.05)
-    @test SciMLBase.successful_retcode(sol1.retcode)
+    multi_point_bvp = BVProblem(simplependulum!, bc!, u0, tspan)
+
+    @testset "BigFloat compatiability with Multi-point BVP" begin
+        for solver in [MIRK4(), RadauIIa5(), LobattoIIIa4(nested_nlsolve = true)]
+            sol = solve(multi_point_bvp, solver, dt = 0.05)
+            @test SciMLBase.successful_retcode(sol.retcode)
+        end
+    end
 
-    sol2 = solve(prob, RadauIIa5(), dt = 0.05)
-    @test SciMLBase.successful_retcode(sol2.retcode)
+    function f!(du, u, p, t)
+        du[1] = u[2]
+        du[2] = u[1]
+    end
+    function bca!(resid_a, u_a, p)
+        resid_a[1] = u_a[1] - 1
+    end
+    function bcb!(resid_b, u_b, p)
+        resid_b[1] = u_b[1]
+    end
+    bvp_function = BVPFunction(
+        f!, (bca!, bcb!), bcresid_prototype = (zeros(1), zeros(1)), twopoint = Val(true))
+    tspan = (0.0, 1.0)
+    two_point_bvp = BVProblem(bvp_function, BigFloat.([1.0, 0.0]), tspan)
 
-    sol3 = solve(prob, LobattoIIIa4(nested_nlsolve = true), dt = 0.05)
-    @test SciMLBase.successful_retcode(sol3.retcode)
+    @testset "BigFloat compatiability with Two-point BVP" begin
+        for solver in [MIRK4(), RadauIIa5(), LobattoIIIa4(nested_nlsolve = true)]
+            sol = solve(two_point_bvp, solver, dt = 0.05)
+            @test SciMLBase.successful_retcode(sol.retcode)
+        end
+    end
 
     function second_f!(ddu, du, u, p, t)
         ϵ = 0.1