How to input a tolerance in BVP

[JMmontilla]

Is there a way to pass the ftol argument to the BVProblem? As an example, say that I want the residuals of this problem to be zero up to 25 decimals:

using BoundaryValueDiffEq, OrdinaryDiffEq
const g = 9.81
L = 1.0
tspan = (0.0,pi/2)
function simplependulum!(du,u,p,t)
    θ  = u[1]
    dθ = u[2]
    du[1] = dθ
    du[2] = -(g/p[1])*sin(θ)
end
u₀_2 = BigFloat.(["-0.47", "-4.79"]) # the initial guess
tspan = BigFloat.(tspan)
function bc3!(residual, sol, p, t)
    residual[1] = sol(pi/4)[1] + pi/2 # use the interpolation here, since indexing will be wrong for adaptive methods
    residual[2] = sol(pi/2)[1] - pi/2
end
bvp3 = BVProblem(simplependulum!, bc3!, u₀_2, tspan, [L])
sol3 = solve(bvp3, Shooting(Feagin14()), abstol = 1e-40, reltol = 1e-35)

@show sol3(pi/2)[1] - pi/2

How can I do that?

[ChrisRackauckas]

Take a look at https://github.com/JuliaDiffEq/BoundaryValueDiffEq.jl/blob/master/test/orbital.jl#L51-L54 . We should update the docs if it's not clear on how to pass a nonlinear solver. Please let us know if that should be done.

[JMmontilla]

Thanks :)