From 21f4dc8f2353b4a02bbf359a16ec92afec2f2184 Mon Sep 17 00:00:00 2001
From: Hendrik Ranocha <ranocha@users.noreply.github.com>
Date: Fri, 9 Sep 2022 08:04:00 +0200
Subject: [PATCH] fix for wrapped f (#148)

* fix for wrapped f

* .f

* do not fix the bug here
---
 Project.toml                     | 4 ++--
 src/SummationByPartsOperators.jl | 3 ++-
 src/diffeq_callbacks.jl          | 3 ++-
 3 files changed, 6 insertions(+), 4 deletions(-)

diff --git a/Project.toml b/Project.toml
index 1eb16e9f..187e8f93 100644
--- 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
index ff42934c..fed57068 100644
--- 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
index 14c231c9..270063ee 100644
--- 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