Merge branch 'master' into qqy/quality

Author: ErikQQY

Project.toml

@@ -67,15 +67,17 @@ DiffEqDevTools = "f3b72e0c-5b89-59e1-b016-84e28bfd966d"
 Hwloc = "0e44f5e4-bd66-52a0-8798-143a42290a1d"
 InteractiveUtils = "b77e0a4c-d291-57a0-90e8-8db25a27a240"
 JET = "c3a54625-cd67-489e-a8e7-0a5a0ff4e31b"
+LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
 NonlinearSolveFirstOrder = "5959db7a-ea39-4486-b5fe-2dd0bf03d60d"
 ODEInterface = "54ca160b-1b9f-5127-a996-1867f4bc2a2c"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 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", "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"]

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
 
 ```julia
-BoundaryValueDiffEqCore.BoundaryValueDiffEqAlgorithm
+BoundaryValueDiffEqCore.AbstractBoundaryValueDiffEqAlgorithm
 ```

docs/src/solvers/firk.md

@@ -12,9 +12,11 @@ solve(prob::BVProblem, alg, dt; kwargs...)
 solve(prob::TwoPointBVProblem, alg, dt; kwargs...)
 ```
 
-!!! note "Nested nonlinear solving in FIRK methods"
-    
-    When encountered with large BVP system, setting `nested_nlsolve` to `true` enables FIRK methods to use nested nonlinear solving for the implicit FIRK step instead of solving as a part of the global residual(when default as `nested_nlsolve=false`),
+## Nested nonlinear solving in FIRK methods
+
+When working with large boundary value problems, especially those involving stiff systems, computational efficiency and solver robustness become critical concerns. To improve the efficiency of FIRK methods on large BVPs, we can use nested nonlinear solving to obtain the implicit FIRK step instead of solving them as part of the global residual. In BoundaryValueDiffEq.jl, we can set `nested_nlsolve` as `true` to enable FIRK methods to compute the implicit FIRK steps using nested nonlinear solving(default option in FIRK methods is `nested_nlsolve=false`).
+
+Moreover, the nested nonlinear problem solver can be finely tuned to meet specific accuracy requirements by providing detailed keyword arguments through the `nested_nlsolve_kwargs` option in any FIRK solver, for example, `RadauIIa5(; nested_nlsolve = true, nested_nlsolve_kwargs = (; abstol = 1e-6, reltol = 1e-6))`, where `nested_nlsolve_kwargs` can be any common keyword arguments in NonlinearSolve.jl, see [Common Solver Options in NonlinearSolve.jl](https://docs.sciml.ai/NonlinearSolve/stable/basics/solve/).
 
 ## Full List of Methods
 

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,14 @@ module BoundaryValueDiffEqAscher
 
 using ADTypes: ADTypes, AutoSparse, AutoForwardDiff
 using AlmostBlockDiagonals: AlmostBlockDiagonals, IntermediateAlmostBlockDiagonal
-using BoundaryValueDiffEqCore: BVPJacobianAlgorithm, __extract_problem_details,
-                               concrete_jacobian_algorithm, __Fix3,
-                               __concrete_nonlinearsolve_algorithm,
-                               __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
+
 using ConcreteStructs: @concrete
 using DiffEqBase: DiffEqBase
 using DifferentiationInterface: DifferentiationInterface, Constant, prepare_jacobian