diff --git a/src/filter/inversion.jl b/src/filter/inversion.jl
index e7c9472c..dd4420f0 100644
--- a/src/filter/inversion.jl
+++ b/src/filter/inversion.jl
@@ -3491,6 +3491,7 @@ function filter_data_with_model(𝓂::ℳ,
     data_in_deviations::KeyedArray{Float64},
     ::Val{:second_order}, # algo
     ::Val{:inversion}; # filter
+    quadratic_matrix_equation_solver::Symbol = :schur,
     warmup_iterations::Int = 0,
     filter_algorithm::Symbol = :LagrangeNewton,
     smooth::Bool = true,
@@ -3706,6 +3707,7 @@ function filter_data_with_model(𝓂::ℳ,
                                 data_in_deviations::KeyedArray{Float64},
                                 ::Val{:pruned_second_order}, # algo
                                 ::Val{:inversion}; # filter
+                                quadratic_matrix_equation_solver::Symbol = :schur,
                                 warmup_iterations::Int = 0,
                                 filter_algorithm::Symbol = :LagrangeNewton,
                                 smooth::Bool = true,
@@ -3979,6 +3981,7 @@ function filter_data_with_model(𝓂::ℳ,
                                 data_in_deviations::KeyedArray{Float64},
                                 ::Val{:third_order}, # algo
                                 ::Val{:inversion}; # filter
+                                quadratic_matrix_equation_solver::Symbol = :schur,
                                 warmup_iterations::Int = 0,
                                 filter_algorithm::Symbol = :LagrangeNewton,
                                 smooth::Bool = true,
@@ -4289,6 +4292,7 @@ function filter_data_with_model(𝓂::ℳ,
                                 data_in_deviations::KeyedArray{Float64},
                                 ::Val{:pruned_third_order}, # algo
                                 ::Val{:inversion}; # filter
+                                quadratic_matrix_equation_solver::Symbol = :schur,
                                 warmup_iterations::Int = 0,
                                 filter_algorithm::Symbol = :LagrangeNewton,
                                 smooth::Bool = true,
diff --git a/src/filter/kalman.jl b/src/filter/kalman.jl
index e6be997e..d51369dc 100644
--- a/src/filter/kalman.jl
+++ b/src/filter/kalman.jl
@@ -559,6 +559,7 @@ function filter_data_with_model(𝓂::ℳ,
     data_in_deviations::KeyedArray{Float64},
     ::Val{:first_order}, # algo
     ::Val{:kalman}; # filter
+    quadratic_matrix_equation_solver::Symbol = :schur,
     warmup_iterations::Int = 0,
     smooth::Bool = true,
     verbose::Bool = false)