diff --git a/src/rules/continuous_transition/W.jl b/src/rules/continuous_transition/W.jl index 99b5fc8c9..5159a6a13 100644 --- a/src/rules/continuous_transition/W.jl +++ b/src/rules/continuous_transition/W.jl @@ -9,12 +9,10 @@ function compute_delta(my, Vy, mx, Vx, Vyx, mA, Va, ma, Fs) G₆ = zeros(eltype(ma), dy, dy) mamat = ma * ma' - Y = similar(Ex_xx) Z = similar(Ex_xx) @inbounds for (i, j) in Iterators.product(1:dy, 1:dy) - mul!(Y, Ex_xx, Fs[j]) mul!(Z, Fs[i]', Y) diff --git a/src/rules/continuous_transition/a.jl b/src/rules/continuous_transition/a.jl index ca20ddd18..1ce867a61 100644 --- a/src/rules/continuous_transition/a.jl +++ b/src/rules/continuous_transition/a.jl @@ -23,7 +23,6 @@ for i in 1:dy xi += Fs[i]' * Vxymxy * mW[:, i] for j in 1:dy - mul!(Y, Vxmx, Fs[j]) mul!(W, Fs[i]', Y, mW[j, i], 1) diff --git a/src/rules/continuous_transition/marginals.jl b/src/rules/continuous_transition/marginals.jl index d55f854d2..3d1b9ac9b 100644 --- a/src/rules/continuous_transition/marginals.jl +++ b/src/rules/continuous_transition/marginals.jl @@ -26,10 +26,8 @@ function continuous_tranition_marginal(m_y::MultivariateNormalDistributionsFamil Ξ = Wx - Y = similar(Va) for (i, j) in Iterators.product(1:dy, 1:dy) - mul!(Y, Va, Fs[i]') mul!(Ξ, Fs[j], Y, mW[j, i], 1) diff --git a/src/rules/continuous_transition/x.jl b/src/rules/continuous_transition/x.jl index 4ffdfc24a..79a32a6cc 100644 --- a/src/rules/continuous_transition/x.jl +++ b/src/rules/continuous_transition/x.jl @@ -18,7 +18,6 @@ Y = similar(Va) for (i, j) in Iterators.product(1:dy, 1:dy) - mul!(Y, Va, Fs[i]') mul!(Ξ, Fs[j], Y, mW[j, i], 1)