From fbac03415e1bb8bf16de148e874fc030c0b4c287 Mon Sep 17 00:00:00 2001 From: Marco Ratto <marco.ratto@jrc.ec.europa.eu> Date: Wed, 1 Feb 2012 15:38:59 +0100 Subject: [PATCH] bug fix when STEADY_STATE instances are used in model definition. --- matlab/getH.m | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/matlab/getH.m b/matlab/getH.m index 2b505c86b..864d6319b 100644 --- a/matlab/getH.m +++ b/matlab/getH.m @@ -31,8 +31,9 @@ yy0=oo_.dr.ys(I); % yy0 = [ yy0; oo_.dr.ys(find(M_.lead_lag_incidence(j,:)))]; % end dyssdtheta=zeros(length(oo_.dr.ys),M_.param_nbr); +d2yssdtheta=zeros(length(oo_.dr.ys),M_.param_nbr,M_.param_nbr); df = feval([M_.fname,'_params_derivs'],yy0, oo_.exo_steady_state', ... - M_.params, oo_.dr.ys, 1, dyssdtheta); + M_.params, oo_.dr.ys, 1, dyssdtheta,d2yssdtheta); [residual, gg1] = feval([M_.fname,'_static'],oo_.dr.ys, oo_.exo_steady_state', M_.params); dyssdtheta = -gg1\df; if any(any(isnan(dyssdtheta))), @@ -47,7 +48,7 @@ if any(any(isnan(dyssdtheta))), dyssdtheta = -U(:,k+1:end)*(T\U(:,k+1:end)')*df; end [df, gp] = feval([M_.fname,'_params_derivs'],yy0, oo_.exo_steady_state', ... - M_.params, oo_.dr.ys, 1, dyssdtheta); + M_.params, oo_.dr.ys, 1, dyssdtheta,d2yssdtheta); [residual, g1, g2 ] = feval([M_.fname,'_dynamic'],yy0, oo_.exo_steady_state', ... M_.params, oo_.dr.ys, 1); -- GitLab