Support SciMLBase v3

Project.toml

@@ -1,6 +1,6 @@
 name = "OrdinaryDiffEqOperatorSplitting"
 uuid = "760fc936-9fa0-4281-9c1a-468957eedc89"
-version = "0.2.3"
+version = "0.2.4"
 authors = ["Dennis Ogiermann <termi-official@users.noreply.github.com> and contributors"]
 
 [deps]
@@ -27,7 +27,7 @@ OrdinaryDiffEqTsit5 = "1.1.0"
 PrecompileTools = "1.0"
 RecursiveArrayTools = "3.39.0"
 SafeTestsets = "0.1.0"
-SciMLBase = "2.77.0"
+SciMLBase = "2.77.0, 3.1"
 TimerOutputs = "0.5.28"
 Unrolled = "0.1.5"
 julia = "1.10"

src/integrator.jl

@@ -1115,3 +1115,11 @@ end
 
 DiffEqBase.u_modified!(i::OperatorSplittingIntegrator, bool) = i.u_modified = bool
 DiffEqBase.u_modified!(i::SplitSubIntegrator, bool) = i.u_modified = bool
+
+# SciMLBase v3 renamed `u_modified!` → `derivative_discontinuity!`. On v3+,
+# also provide the overloads under the new name so callers using the new
+# API dispatch correctly instead of hitting the generic `error(...)` fallback.
+@static if isdefined(SciMLBase, :derivative_discontinuity!)
+    SciMLBase.derivative_discontinuity!(i::OperatorSplittingIntegrator, bool) = i.u_modified = bool
+    SciMLBase.derivative_discontinuity!(i::SplitSubIntegrator, bool) = i.u_modified = bool
+end

src/utils.jl

@@ -74,7 +74,13 @@ function forward_sync_internal!(u_source, child::DEIntegrator, solution_indices)
     @views usrc = u_source[solution_indices]
     sync_vectors!(child.u, usrc)
     sync_vectors!(child.uprev, child.u)
-    SciMLBase.u_modified!(child, true)
+    # SciMLBase v3 renamed this to `derivative_discontinuity!`; call the
+    # appropriate name based on which SciMLBase is loaded.
+    @static if isdefined(SciMLBase, :derivative_discontinuity!)
+        SciMLBase.derivative_discontinuity!(child, true)
+    else
+        SciMLBase.u_modified!(child, true)
+    end
     return nothing
 end