From e60d438e86f24d69f1961e67f3737b5d7b2d46c1 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Charybdis=29?= Date: Fri, 28 Jun 2013 16:59:47 +0200 Subject: [PATCH] Fixed bug. --- matlab/DsgeSmoother.m | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/matlab/DsgeSmoother.m b/matlab/DsgeSmoother.m index ecee9f3ed..4b58ca34e 100644 --- a/matlab/DsgeSmoother.m +++ b/matlab/DsgeSmoother.m @@ -151,14 +151,14 @@ elseif options_.lik_init == 5 % Old diffuse Kalman filter only for th end R_tmp = R(stable, :); T_tmp = T(stable,stable); - if DynareOptions.lyapunov_fp == 1 - Pstar_tmp = lyapunov_symm(T_tmp,Q,DynareOptions.lyapunov_fixed_point_tol,DynareOptions.lyapunov_complex_threshold, 3, R_tmp); - elseif DynareOptions.lyapunov_db == 1 - Pstar_tmp = disclyap_fast(T_tmp,R_tmp*Q*R_tmp',DynareOptions.lyapunov_doubling_tol); - elseif DynareOptions.lyapunov_srs == 1 - Pstar_tmp = lyapunov_symm(T_tmp,Q,DynareOptions.lyapunov_fixed_point_tol,DynareOptions.lyapunov_complex_threshold, 4, R_tmp); + if options_.lyapunov_fp == 1 + Pstar_tmp = lyapunov_symm(T_tmp,Q,options_.lyapunov_fixed_point_tol,options_.lyapunov_complex_threshold, 3, R_tmp); + elseif options_.lyapunov_db == 1 + Pstar_tmp = disclyap_fast(T_tmp,R_tmp*Q*R_tmp',options_.lyapunov_doubling_tol); + elseif options_.lyapunov_srs == 1 + Pstar_tmp = lyapunov_symm(T_tmp,Q,options_.lyapunov_fixed_point_tol,options_.lyapunov_complex_threshold, 4, R_tmp); else - Pstar_tmp = lyapunov_symm(T_tmp,R_tmp*Q*R_tmp',DynareOptions.qz_criterium,DynareOptions.lyapunov_complex_threshold); + Pstar_tmp = lyapunov_symm(T_tmp,R_tmp*Q*R_tmp',options_.qz_criterium,options_.lyapunov_complex_threshold); end Pstar(stable, stable) = Pstar_tmp; Pinf = [];