fix for wrapped f (#148)

* fix for wrapped f * .f * do not fix the bug here
ranocha · Sep 9, 2022 · 21f4dc8 · 21f4dc8 · ranocha · Sep 9, 2022
1 parent 1c5858a
commit 21f4dc8
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 4 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "SummationByPartsOperators"
 uuid = "9f78cca6-572e-554e-b819-917d2f1cf240"
 author = ["Hendrik Ranocha"]
-version = "0.5.21"
+version = "0.5.22"
 
 [deps]
 ArgCheck = "dce04be8-c92d-5529-be00-80e4d2c0e197"
@@ -29,7 +29,7 @@ PolynomialBases = "0.4.5"
 RecursiveArrayTools = "2.11"
 Reexport = "0.2, 1.0"
 Requires = "0.5.2, 1"
-SciMLBase = "1.21"
+SciMLBase = "1.51"
 StaticArrays = "1.0"
 UnPack = "1.0"
 Unrolled = "0.1.3"

diff --git a/src/SummationByPartsOperators.jl b/src/SummationByPartsOperators.jl
@@ -57,7 +57,8 @@ using StaticArrays: SVector, StaticVector, FieldVector
 using UnPack: @unpack
 using Unrolled: @unroll
 
-@reexport using SciMLBase: DiscreteCallback, ODEProblem, SecondOrderODEProblem
+@reexport using SciMLBase: SciMLBase, DiscreteCallback, ODEProblem,
+                           SecondOrderODEProblem
 
 import LinearAlgebra: mul!
 @reexport using PolynomialBases

diff --git a/src/diffeq_callbacks.jl b/src/diffeq_callbacks.jl
@@ -4,7 +4,8 @@ function SavingCallback(semi::AbstractSemidiscretization; kwargs...)
   T = eltype(semi.derivative)
 
   save_func = (u,t,integrator) -> integrate(u->ScalarIntegralQuantities(u,u^2),
-                                              u, integrator.f.f)
+                                            u,
+                                            SciMLBase.unwrapped_f(integrator.f.f))
   saved_values = SavedValues(T, ScalarIntegralQuantities{T})
   SavingCallback(save_func, saved_values; kwargs...)
 end