Observations must be kept constant: new matrix, not alias

DynareJulia · Jun 23, 2024 · ad5586f · ad5586f
1 parent 47193d5
commit ad5586f
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/estimation/estimation.jl b/src/estimation/estimation.jl
@@ -390,7 +390,7 @@ struct SSWs{D<:AbstractFloat,I<:Integer}
         ns = length(state_ids)
         ny = length(varobs)
         nobs = size(observations, 2)
-        Y = observations
+        Y = similar(observations)
         Z = zeros(D, (ny, ns))
         H = zeros(D, (ny, ny))
         T = zeros(D, (ns, ns))
@@ -630,7 +630,6 @@ function loglikelihood(
         context.work.observed_variables,
         dim = 1
     )
-    ssws.Y
     ns = length(ssws.state_ids)
     np = model.exogenous_nbr
     ny, nobs = size(ssws.Y)