Merge branch 'master' into qqy/more_AD

Author: ErikQQY

Project.toml

@@ -42,21 +42,22 @@ BoundaryValueDiffEqFIRK = "1"
 BoundaryValueDiffEqMIRK = "1"
 BoundaryValueDiffEqMIRKN = "1"
 BoundaryValueDiffEqShooting = "1"
-DiffEqBase = "6.161.0"
+DiffEqBase = "6.167"
 DiffEqDevTools = "2.45"
 FastClosures = "0.3.2"
 ForwardDiff = "0.10.38, 1"
 Hwloc = "3.3"
 InteractiveUtils = "<0.0.1, 1"
 JET = "0.9"
 LinearAlgebra = "1.10"
+NonlinearSolveFirstOrder = "1.3"
 ODEInterface = "0.5"
-OrdinaryDiffEq = "6.90.1"
+OrdinaryDiffEq = "6.93"
 Pkg = "1.10.0"
 Random = "1.10"
 ReTestItems = "1.29"
 Reexport = "1.2"
-SciMLBase = "2.71"
+SciMLBase = "2.82"
 StaticArrays = "1.9.8"
 Test = "1.10"
 julia = "1.10"

docs/Project.toml

@@ -19,6 +19,16 @@ Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 SimpleBoundaryValueDiffEq = "be0294bd-f90f-4760-ac4e-3421ce2b2da0"
 
+[sources]
+BoundaryValueDiffEq = {path = "../"}
+BoundaryValueDiffEqAscher = {path = "../lib/BoundaryValueDiffEqAscher"}
+BoundaryValueDiffEqCore = {path = "../lib/BoundaryValueDiffEqCore"}
+BoundaryValueDiffEqFIRK = {path = "../lib/BoundaryValueDiffEqFIRK"}
+BoundaryValueDiffEqMIRK = {path = "../lib/BoundaryValueDiffEqMIRK"}
+BoundaryValueDiffEqMIRKN = {path = "../lib/BoundaryValueDiffEqMIRKN"}
+BoundaryValueDiffEqShooting = {path = "../lib/BoundaryValueDiffEqShooting"}
+
+
 [compat]
 ADTypes = "1.9"
 BoundaryValueDiffEq = "5.13.0"

docs/src/assets/Project.toml

@@ -19,6 +19,16 @@ Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 SimpleBoundaryValueDiffEq = "be0294bd-f90f-4760-ac4e-3421ce2b2da0"
 
+[sources]
+BoundaryValueDiffEq = {path = "../"}
+BoundaryValueDiffEqAscher = {path = "../lib/BoundaryValueDiffEqAscher"}
+BoundaryValueDiffEqCore = {path = "../lib/BoundaryValueDiffEqCore"}
+BoundaryValueDiffEqFIRK = {path = "../lib/BoundaryValueDiffEqFIRK"}
+BoundaryValueDiffEqMIRK = {path = "../lib/BoundaryValueDiffEqMIRK"}
+BoundaryValueDiffEqMIRKN = {path = "../lib/BoundaryValueDiffEqMIRKN"}
+BoundaryValueDiffEqShooting = {path = "../lib/BoundaryValueDiffEqShooting"}
+
+
 [compat]
 ADTypes = "1.9"
 BoundaryValueDiffEq = "5.13.0"

docs/src/basics/error_control.md

@@ -25,6 +25,7 @@ DefectControl
 GlobalErrorControl
 SequentialErrorControl
 HybridErrorControl
+NoErrorControl
 ```
 
 While we can achieve global error control in different ways, we can use different methods to estimate the global error:

docs/src/basics/solve.md

@@ -2,8 +2,7 @@
 
 ## Iteration Controls
 
-  - `abstol::Number`: The absolute tolerance. Defaults to `1e-3`.
-  - `reltol::Number`: The relative tolerance. Defaults to `1e-3`.
+  - `abstol::Number`: The absolute tolerance. Defaults to `1e-6`.
   - `defect_threshold`: Monitor of the size of defect norm. Defaults to `0.1`.
   - `odesolve_kwargs`: OrdinaryDiffEq.jl solvers kwargs for passing to ODE solving in shooting methods.
   - `nlsolve_kwargs`: NonlinearSolve.jl solvers kwargs for passing to nonlinear solving in collocation methods and shootingn methods.

docs/src/devdocs/internal_interfaces.md

@@ -3,5 +3,5 @@
 ## Solvers
 
 ```@docs
-BoundaryValueDiffEqCore.BoundaryValueDiffEqAlgorithm
+BoundaryValueDiffEqCore.AbstractBoundaryValueDiffEqAlgorithm
 ```

docs/src/tutorials/solve_nlls_bvp.md

@@ -35,7 +35,7 @@ end
 tspan = (0.0, 100.0)
 u0 = [0.0, 1.0]
 prob = BVProblem(BVPFunction(f!, bc!; bcresid_prototype = zeros(3)), u0, tspan)
-sol = solve(prob, MIRK4(), dt = 0.01)
+sol = solve(prob, MIRK4(), dt = 0.01, abstol = 1e-3)
 plot(sol)
 ```
 

ext/BoundaryValueDiffEqODEInterfaceExt.jl

@@ -1,8 +1,7 @@
 module BoundaryValueDiffEqODEInterfaceExt
 
 using BoundaryValueDiffEq: BVPM2, BVPSOL, COLNEW
-using BoundaryValueDiffEqCore: BoundaryValueDiffEqAlgorithm, __extract_u0,
-                               __initial_guess_length, __extract_mesh,
+using BoundaryValueDiffEqCore: __extract_u0, __initial_guess_length, __extract_mesh,
                                __flatten_initial_guess, __get_bcresid_prototype,
                                __has_initial_guess, __initial_guess
 using SciMLBase: SciMLBase, BVProblem, TwoPointBVProblem, ReturnCode

lib/BoundaryValueDiffEqAscher/src/BoundaryValueDiffEqAscher.jl

@@ -2,12 +2,15 @@ module BoundaryValueDiffEqAscher
 
 using ADTypes: ADTypes, AutoSparse, AutoForwardDiff
 using AlmostBlockDiagonals: AlmostBlockDiagonals, IntermediateAlmostBlockDiagonal
-using BoundaryValueDiffEqCore: BVPJacobianAlgorithm, __extract_problem_details,
-                               concrete_jacobian_algorithm, __Fix3,
-                               __concrete_nonlinearsolve_algorithm, __default_nonsparse_ad,
-                               __internal_nlsolve_problem, BoundaryValueDiffEqAlgorithm,
-                               __vec, __vec_f, __vec_f!, __vec_bc, __vec_bc!,
-                               __extract_mesh, get_dense_ad, __get_bcresid_prototype
+
+using BoundaryValueDiffEqCore: AbstractBoundaryValueDiffEqAlgorithm,
+                               AbstractBoundaryValueDiffEqCache, BVPJacobianAlgorithm,
+                               __extract_problem_details, concrete_jacobian_algorithm,
+                               __Fix3, __concrete_nonlinearsolve_algorithm,
+                               __internal_nlsolve_problem, __vec, __vec_f, __vec_f!,
+                               __vec_bc, __vec_bc!, __extract_mesh, get_dense_ad,
+                               __get_bcresid_prototype, __split_kwargs, __default_nonsparse_ad
+
 using ConcreteStructs: @concrete
 using DiffEqBase: DiffEqBase
 using DifferentiationInterface: DifferentiationInterface, Constant, prepare_jacobian

lib/BoundaryValueDiffEqAscher/src/algorithms.jl

@@ -1,4 +1,4 @@
-abstract type AbstractAscher <: BoundaryValueDiffEqAlgorithm end
+abstract type AbstractAscher <: AbstractBoundaryValueDiffEqAlgorithm end
 
 for stage in (1, 2, 3, 4, 5, 6, 7)
     alg = Symbol("Ascher$(stage)")