Allow to dispatch default_options() on integrator method and add custom options for SPARK methods.

Author: michakraus

src/integrators/integrator.jl

@@ -48,7 +48,7 @@ function GeometricIntegrator(
         integratormethod::GeometricMethod,
         solvermethod::SolverMethod,
         iguess::Extrapolation;
-        options = default_options(),
+        options = default_options(integratormethod),
         method = initmethod(integratormethod, problem),
         caches = CacheDict(problem, method),
         solver = initsolver(solvermethod, options, method, caches)

src/integrators/solvers.jl

@@ -1,7 +1,7 @@
 
 struct NoSolver <: SolverMethod end
 
-default_options() = Options(
+default_options(method = nothing) = Options(
     min_iterations = 1,
     x_abstol = 8eps(),
     f_abstol = 8eps(),

src/spark/abstract.jl

@@ -5,22 +5,28 @@ abstract type LSPARKMethod <: LDAEMethod end
 abstract type PSPARKMethod <: PDAEMethod end
 
 const AbstractSPARKMethod = Union{HSPARKMethod, ISPARKMethod, LSPARKMethod, PSPARKMethod}
-const AbstractSPARKProblem{DT<:Number, TT<:Real} = 
-    Union{HDAEProblem{DT,TT}, 
+const AbstractSPARKProblem{DT<:Number, TT<:Real} =
+    Union{HDAEProblem{DT,TT},
           IDAEProblem{DT,TT},
           LDAEProblem{DT,TT},
           PDAEProblem{DT,TT}}
 
-Integrators.default_solver(::AbstractSPARKMethod) = Newton()
 Integrators.default_iguess(::AbstractSPARKMethod) = HermiteExtrapolation()
+Integrators.default_solver(::AbstractSPARKMethod) = Newton()
+Integrators.default_options(::AbstractSPARKMethod) = Options(
+    x_reltol = 8eps(),
+    x_suctol = 2eps(),
+    f_abstol = 8eps(),
+    f_reltol = 8eps(),
+    f_suctol = 2eps(),
+)
 
 
 nstages(method::AbstractSPARKMethod) = nstages(tableau(method))
 pstages(method::AbstractSPARKMethod) = pstages(tableau(method))
 eachstage(method::AbstractSPARKMethod) = eachstage(tableau(method))
 hasnullvector(method::AbstractSPARKMethod) = hasnullvector(tableau(method))
 
-
 # create nonlinear solver
 function initsolver(::Newton, config::Options, method::AbstractSPARKMethod, caches::CacheDict)
     NewtonSolver(zero(nlsolution(caches)), zero(nlsolution(caches)); linesearch = Backtracking(), config = config)