Commit 044d19c5 authored by MichelJuillard's avatar MichelJuillard
Browse files

fixed bug in detection of near 0/0 eigenvalues

parent d91c4ac1
...@@ -99,8 +99,9 @@ if options_.ramsey_policy ...@@ -99,8 +99,9 @@ if options_.ramsey_policy
end end
ys = oo_.steady_state; ys = oo_.steady_state;
if inst_nbr == 1 if inst_nbr == 1
nl_func = @(x) dyn_ramsey_static_(x,M_,options_,oo_,it_); nl_func = @(x) dyn_ramsey_static_(x,M_,options_,oo_,it_);
inst_val = fzero(nl_func,oo_.steady_state(k_inst)); % inst_val = fzero(nl_func,oo_.steady_state(k_inst));
inst_val = csolve(nl_func,oo_.steady_state(k_inst),'',options_.solve_tolf,100);
else else
[inst_val,info1] = dynare_solve('dyn_ramsey_static_', ... [inst_val,info1] = dynare_solve('dyn_ramsey_static_', ...
oo_.steady_state(k_inst),0, ... oo_.steady_state(k_inst),0, ...
...@@ -380,7 +381,7 @@ else % use original Dynare solver ...@@ -380,7 +381,7 @@ else % use original Dynare solver
end end
for i=1:nd for i=1:nd
if abs(ss(i,i)) < 1e-6 && abs(tt(i,i)) < 1e-6 if abs(ss(i,i)) < 1e-6 && abs(tt(i,i)) < 1e-6
info(7) = 1; info(1) = 7;
end end
end end
return return
......
...@@ -49,7 +49,7 @@ if ~noprint ...@@ -49,7 +49,7 @@ if ~noprint
error(['The Jacobian matrix evaluated at the steady state contains elements ' ... error(['The Jacobian matrix evaluated at the steady state contains elements ' ...
'that are not real or are infinite']) 'that are not real or are infinite'])
case 7 case 7
error(['One of the eigenvalue is is close to 0/0 (the absolute ' ... error(['One of the eigenvalues is is close to 0/0 (the absolute ' ...
'value of numerator and denominator is smaller than 1e-6).']) 'value of numerator and denominator is smaller than 1e-6).'])
case 19 case 19
error('The steadystate file did not compute the steady state (inconsistent deep parameters).') error('The steadystate file did not compute the steady state (inconsistent deep parameters).')
......
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