Add two-point BVP BigFloat test case

Author: ErikQQY

test/misc/bigfloat_test.jl

@@ -12,15 +12,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)
 
-    sol2 = solve(prob, RadauIIa5(), dt = 0.05)
-    @test SciMLBase.successful_retcode(sol2.retcode)
+    @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
 
-    sol3 = solve(prob, LobattoIIIa4(nested_nlsolve = true), dt = 0.05)
-    @test SciMLBase.successful_retcode(sol3.retcode)
+    function f!(du, u, p, t)
+        du[1] = u[2]
+        du[2] = 1 / p * 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)
+
+    @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