use linear interpolations (#62)

* use linear interpolation for MPE * format * fix undefined reference * test also MVector and out-of-place arrays * in-place MPE * update MPRK22 - needs to be checked * format * fix convergence tests * bump version
SKopecz · Apr 9, 2024 · 4cc77e3 · 4cc77e3 · ranocha · Apr 9, 2024
1 parent 290fd8d
commit 4cc77e3
Show file tree

Hide file tree

Showing 5 changed files with 207 additions and 65 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "PositiveIntegrators"
 uuid = "d1b20bf0-b083-4985-a874-dc5121669aa5"
 authors = ["Stefan Kopecz, Hendrik Ranocha, and contributors"]
-version = "0.1.6-pre"
+version = "0.1.6"
 
 [deps]
 FastBroadcast = "7034ab61-46d4-4ed7-9d0f-46aef9175898"

diff --git a/src/PositiveIntegrators.jl b/src/PositiveIntegrators.jl
@@ -24,20 +24,24 @@ using SymbolicIndexingInterface
 
 using LinearSolve: LinearSolve, LinearProblem, LUFactorization
 
-import SciMLBase: __has_mass_matrix, __has_analytic, __has_tgrad,
+using SciMLBase: DEFAULT_OBSERVED
+import SciMLBase: interp_summary,
+                  __has_mass_matrix, __has_analytic, __has_tgrad,
                   __has_jac, __has_jvp, __has_vjp, __has_jac_prototype,
                   __has_sparsity, __has_Wfact, __has_Wfact_t,
                   __has_paramjac, __has_syms, __has_indepsym, __has_paramsyms,
-                  __has_observed, DEFAULT_OBSERVED, __has_colorvec,
-                  __has_sys
-
-import OrdinaryDiffEq: @cache,
-                       OrdinaryDiffEqAdaptiveAlgorithm, calculate_residuals,
-                       calculate_residuals!, False,
-                       OrdinaryDiffEqMutableCache, OrdinaryDiffEqConstantCache,
+                  __has_observed, __has_colorvec, __has_sys
+
+using OrdinaryDiffEq: @cache,
+                      DEFAULT_PRECS,
+                      OrdinaryDiffEqAdaptiveAlgorithm,
+                      OrdinaryDiffEqConstantCache, OrdinaryDiffEqMutableCache,
+                      False,
+                      recursivefill!, _vec, wrapprecs, dolinsolve
+import OrdinaryDiffEq: alg_order, isfsal,
+                       calculate_residuals, calculate_residuals!,
                        alg_cache, initialize!, perform_step!,
-                       recursivefill!, _vec, DEFAULT_PRECS, wrapprecs,
-                       dolinsolve, alg_order
+                       _ode_interpolant, _ode_interpolant!
 
 # 2. Export functionality defining the public API
 export PDSFunction, PDSProblem