stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

made evaluate_steady_state file consistent with evaluation of auxiliary 

variables in <fname>_static.m file
time-shift
Michel Juillard 2016-11-20 15:40:08 +01:00
parent e14a215556
commit 73e0f9d591
1 changed files with 3 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
matlab/evaluate_steady_state.m
View File

@ -224,13 +224,13 @@ function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadysta
[ys,check] = dynare_solve(@static_problem,...
ys_init,...
options, exo_ss, params,...
M.orig_endo_nbr,...
M.endo_nbr,...
static_model);
if check && options.debug
[ys,check,fvec,fjac] = dynare_solve(@static_problem,...
ys_init,...
options, exo_ss, params,...
M.orig_endo_nbr,...
M.endo_nbr,...
static_model);
[infrow,infcol]=find(isinf(fjac) | isnan(fjac));
if ~isempty(infrow)
@ -277,12 +277,7 @@ function [ys,params,info] = evaluate_steady_state(ys_init,M,options,oo,steadysta
[infrow,infcol]=find(isinf(jacob) | isnan(jacob));
fprintf('\nSTEADY: The Jacobian contains Inf or NaN. The problem arises from: \n\n')
for ii=1:length(infrow)
if infcol(ii)<=M.orig_endo_nbr
fprintf('STEADY: Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n',infrow(ii),deblank(M.endo_names(infcol(ii),:)),deblank(M.endo_names(infcol(ii),:)),ys_init(infcol(ii)))
else %auxiliary vars
orig_var_index=M.aux_vars(1,infcol(ii)-M.orig_endo_nbr).orig_index;
fprintf('STEADY: Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n',infrow(ii),deblank(M.endo_names(orig_var_index,:)),deblank(M.endo_names(orig_var_index,:)),ys_init(infcol(ii)))
end
fprintf('STEADY: Derivative of Equation %d with respect to Variable %s (initial value of %s: %g) \n',infrow(ii),deblank(M.endo_names(infcol(ii),:)),deblank(M.endo_names(infcol(ii),:)),ys_init(infcol(ii)))
end
fprintf('Check whether your model is truly linear. Put "resid(1);" before "steady;" to see the problematic equations.\n')
end