Hi Johannes,
this model is non-stationary but options_.qz_criterium = 0.999999 ....
I'm on the train and can't get the original model to understand why.
Best
Michel
Johannes Pfeifer writes:
Hi,
I have a problem with the mod-file at http://www.dynare.org/phpBB3/viewtopic.php?f=1 http://www.dynare.org/phpBB3/viewtopic.php?f=1&t=10760 &t=10760
I narrowed the problem down to a problem with lyapunov_symm.m
When you try the code
load workspace_dsge_simulated_theoretical_correlation
tmp = th_autocovariances(dr,ivar,M_,options_,nodecomposition);
tmp{5,1}
with the attached mat-file, you will see that the covariance matrix at lag 5 is complex, causing the subsequent crash the user reports.
The reason is that
[vx, u] =lyapunov_symm(A,B*M_.Sigma_e*B',options_.lyapunov_fixed_point_tol,options_.q z_criterium,options_.lyapunov_complex_threshold,[],[],options_.debug);
returns an vx matrix with a big negative eigenvalue, which should never happen. However, that matrix actually solves the Lyapunov equation. It simply is not a valid covariance matrix. How should we deal with an issue like this?
Best,
Johannes
Dev mailing list Dev@dynare.org https://www.dynare.org/cgi-bin/mailman/listinfo/dev