From 7400c2c1431c1c4e24df638f3811605cdd1a715f Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= Date: Tue, 8 Jul 2014 16:10:09 +0200 Subject: [PATCH] evaluate_steady_state.m: fix bugs when static and dynamic models differ. Closes #661 --- matlab/evaluate_steady_state.m | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/matlab/evaluate_steady_state.m b/matlab/evaluate_steady_state.m index 38aee0a32..dcaa56133 100644 --- a/matlab/evaluate_steady_state.m +++ b/matlab/evaluate_steady_state.m @@ -128,12 +128,12 @@ function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadysta % Evaluate residual of *dynamic* model using the steady state % computed on the *static* one z = repmat(ys,1,M.maximum_lead + M.maximum_lag + 1); - zx = repmat([oo.exo_simul oo.exo_det_simul],M.maximum_lead + M.maximum_lag + 1, 1); + zx = repmat([exo_ss'], M.maximum_lead + M.maximum_lag + 1, 1); if options.bytecode [chck, r, junk]= bytecode('dynamic','evaluate', z, zx, M.params, ys, 1); mexErrCheck('bytecode', chck); elseif options.block - [r, data] = feval([M.fname '_dynamic'], z', zx, M.params, ys, M.maximum_lag+1, data); + [r, oo.dr] = feval([M.fname '_dynamic'], z', zx, M.params, ys, M.maximum_lag+1, oo.dr); else iyv = M.lead_lag_incidence'; iyr0 = find(iyv(:));