up

SciML · Dec 4, 2024 · a9f2106 · a9f2106
1 parent a9fdfd6
commit a9f2106
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 6 deletions.
diff --git a/src/systems/diffeqs/abstractodesystem.jl b/src/systems/diffeqs/abstractodesystem.jl
@@ -534,15 +534,15 @@ function SciMLBase.BVProblem{iip, specialize}(sys::AbstractODESystem, u0map = []
 
     # Define the boundary conditions.
     bc = if iip 
-        (residual, u, p, t) -> (residual .= u[1] - _u0)
+        (residual, u, p, t) -> (residual .= u[1] .- _u0)
     else
         (u, p, t) -> (u[1] - _u0)
     end
 
     return BVProblem{iip}(f, bc, u0, tspan, p; kwargs1..., kwargs...)
 end
 
-get_callback(prob::BVProblem) = prob.kwargs[:callback]
+get_callback(prob::BVProblem) = error("BVP solvers do not support callbacks.")
 
 """
 ```julia

diff --git a/test/bvproblem.jl b/test/bvproblem.jl
@@ -23,8 +23,8 @@ sol2 = solve(bvp, MIRK4(), dt = 0.01);
 op = ODEProblem(lotkavolterra, u0map, tspan, parammap)
 osol = solve(op, Vern9())
 
-@test isapprox(sol.u[end],osol.u[end]; atol = 0.001)
-@test isapprox(sol2.u[end],osol.u[end]; atol = 0.001)
+@test isapprox(sol.u[end],osol.u[end]; atol = 0.01)
+@test isapprox(sol2.u[end],osol.u[end]; atol = 0.01)
 @test sol.u[1] == [1., 2.]
 @test sol2.u[1] == [1., 2.]
 
@@ -50,7 +50,7 @@ sol2 = solve(bvp2, MIRK4(), dt = 0.01);
 op = ODEProblem(pend, u0map, tspan, parammap)
 osol = solve(op, Vern9())
 
-@test sol.u[end] ≈ osol.u[end]
+@test isapprox(sol.u[end], osol.u[end]; atol = 0.01)
 @test sol.u[1] == [π/2, π/2]
-@test sol2.u[end] ≈ osol.u[end]
+@test isapprox(sol2.u[end], osol.u[end]; atol = 0.01)
 @test sol2.u[1] == [π/2, π/2]
diff --git a/test/runtests.jl b/test/runtests.jl
@@ -80,6 +80,7 @@ end
             @safetestset "NonlinearSystem Test" include("nonlinearsystem.jl")
             @safetestset "PDE Construction Test" include("pde.jl")
             @safetestset "JumpSystem Test" include("jumpsystem.jl")
+            @safetestset "BVProblem Test" include("bvproblem.jl")
             @safetestset "print_tree" include("print_tree.jl")
             @safetestset "Constraints Test" include("constraints.jl")
         end