Commit d6268cf7 authored by Stéphane Adjemian's avatar Stéphane Adjemian
Browse files

Fixed bugs.

(1) Wrong calling sequence for resol routine in check.m.
(2) Replaced steady_state_0 by oo.steady_state in resol.m.
parent bcf4fbfe
......@@ -66,7 +66,7 @@ if isempty(options.qz_criterium)
options.qz_criterium = 1+1e-6;
end
[dr, info] = resol(oo.steady_state,1);
[dr,info,M,options,oo] = resol(1,M,options,oo);
oo.dr = dr;
......
......@@ -112,14 +112,14 @@ end
params0 = M.params;
% check if steady_state_0 is steady state
% check if steady_state_0 (-> oo.steady_state) is steady state
tempex = oo.exo_simul;
oo.exo_simul = repmat(oo.exo_steady_state',M.maximum_lag+M.maximum_lead+1,1);
if M.exo_det_nbr > 0
tempexdet = oo.exo_det_simul;
oo.exo_det_simul = repmat(oo.exo_det_steady_state',M.maximum_lag+M.maximum_lead+1,1);
end
steady_state = steady_state_0;
steady_state = oo.steady_state;
check1 = 0;
% testing for steadystate file
if (~options.bytecode)
......@@ -144,7 +144,7 @@ if options.steadystate_flag
end
else
% testing if steady_state_0 isn't a steady state or if we aren't computing Ramsey policy
% testing if steady_state_0 (-> oo.steady_state) isn't a steady state or if we aren't computing Ramsey policy
if options.ramsey_policy == 0
if options.linear == 0
% nonlinear models
......@@ -193,7 +193,7 @@ if check1
resid = check1 ;
else
info(1)= 20;
resid = feval(fh,steady_state_0,oo.exo_steady_state, M.params);
resid = feval(fh,oo.steady_state,oo.exo_steady_state, M.params);
end
info(2) = resid'*resid ;
return
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment