Make make_transition_matrices! allocate less again

nathanaelbosch · Oct 11, 2023 · c42f9f3 · c42f9f3
1 parent 51b6b02
commit c42f9f3
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/priors/common.jl b/src/priors/common.jl
@@ -73,8 +73,8 @@ function make_transition_matrices!(cache, prior::AbstractODEFilterPrior, dt)
     _Ah, _Qh = discretize(cache.prior, dt)
     copy!(Ah, _Ah)
     copy!(Qh, _Qh)
-    A .= P * Ah * PI
-    @warn "`make_transition_matrices!` currently allocates more than it should!"
+    # A = P * Ah * PI
+    _matmul!(A, P, _matmul!(A, Ah, PI))
     fast_X_A_Xt!(Q, Qh, P)
 end