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

factorized display of critical values.

time-shift
Marco Ratto 2017-03-23 00:01:11 +01:00
parent 76097a799f
commit 234aaacf92
1 changed files with 33 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

56
matlab/perfect-foresight-models/sim1.m
View File

@ -165,6 +165,11 @@ for iter = 1:options.simul.maxit
dy = -lin_solve( A, res, verbose );
end
end
if (any(~isreal(dy)) || any(isnan(dy)) || any(isinf(dy)) ,
if verbose
display_critical_variables(reshape(dy,[ny periods])', M);
end
end
Y(i_upd) = Y(i_upd) + dy;
end
@ -184,33 +189,12 @@ if stop
if verbose
skipline()
disp(sprintf('Total time of simulation: %s.', num2str(etime(clock,h1))))
dyy=reshape(dy,[size(endogenousvariables,1) periods])';
if ~isreal(res) || ~isreal(Y)
disp('Simulation terminated with imaginary parts in the residuals or endogenous variables.')
if any(~isreal(dy))
indx = find(any(~isreal(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Newton algorithm provided complex number for variables:')
fprintf('%s, ', endo_names{:});
skipline()
end
else
disp('Simulation terminated with NaN or Inf in the residuals or endogenous variables.')
if any(isnan(dy))
indx = find(any(isnan(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Newton algorithm provided NaN for variables:')
fprintf('%s, ',endo_names{:});
skipline()
end
if any(isinf(dy))
indx = find(any(isinf(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Newton algorithm diverged (Inf) for variables:')
fprintf('%s, ',endo_names{:});
skipline()
end
end
display_critical_variables(reshape(dy,[ny periods])', M);
disp('There is most likely something wrong with your model. Try model_diagnostics or another simulation method.')
printline(105)
end
@ -322,4 +306,30 @@ function [ x, flag, relres ] = lin_solve_robust( A, b , verbose)
flag = relres > 1e-6;
if flag ~= 0 && verbose
fprintf( 'WARNING : Failed to find a solution to the linear system\n' );
end
end
function display_critical_variables(dyy, M)
if any(isnan(dyy))
indx = find(any(isnan(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Last iteration provided NaN for the following variables:')
fprintf('%s, ',endo_names{:}),
fprintf('\n'),
end
if any(isinf(dyy))
indx = find(any(isinf(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Last iteration diverged (Inf) for the following variables:')
fprintf('%s, ',endo_names{:}),
fprintf('\n'),
end
if any(~isreal(dyy))
indx = find(any(~isreal(dyy)));
endo_names=cellstr(M.endo_names(indx,:));
disp('Last iteration provided complex number for the following variables:')
fprintf('%s, ',endo_names{:}),
fprintf('\n'),
end