diff --git a/src/Utils/matrix_exponential.jl b/src/Utils/matrix_exponential.jl
index 5d7c22bb0c..3faaed6eb1 100644
--- a/src/Utils/matrix_exponential.jl
+++ b/src/Utils/matrix_exponential.jl
@@ -30,8 +30,8 @@ end
 
 end end  # quote / load_exponentialutilities
 
-function _expmv(::Val{:ExponentialUtilities}, t, A, b)
-    return ExponentialUtilities.expv(t, A, b)
+function _expmv(::Val{:ExponentialUtilities}, t, A, b; kwargs...)
+    return ExponentialUtilities.expv(t, A, b; kwargs...)
 end
 
 # =======
@@ -47,11 +47,11 @@ end
 
 end end  # quote / load_expokit
 
-function _expmv(::Val{:Expokit}, t, A, b::Vector)
-    return Expokit.expmv(t, A, b)
+function _expmv(::Val{:Expokit}, t, A, b::Vector; kwargs...)
+    return Expokit.expmv(t, A, b; kwargs...)
 end
 
-function _expmv(backend::Val{:Expokit}, t, A, b)
+function _expmv(backend::Val{:Expokit}, t, A, b; kwargs...)
     b = Vector(b)  # Expokit requires a dense vector
-    return _expmv(backend, t, A, b)
+    return _expmv(backend, t, A, b; kwargs...)
 end