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

v3+v4 check.m: corrected bug in rank condition message 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@821 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
michel 2006-07-12 12:31:42 +00:00
parent 59c35bc1fb
commit 60cd6de62a
1 changed files with 7 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
matlab/check.m
View File

@ -27,6 +27,7 @@ global it_
eigenvalues_ = dr.eigval; eigenvalues_ = dr.eigval;
nyf = nnz(dr.kstate(:,2)>M_.maximum_lag+1); nyf = nnz(dr.kstate(:,2)>M_.maximum_lag+1);
[m_lambda,i]=sort(abs(eigenvalues_)); [m_lambda,i]=sort(abs(eigenvalues_));
n_explod = nnz(abs(eigenvalues_) > options_.qz_criterium);
if options_.noprint == 0 if options_.noprint == 0
disp(' ') disp(' ')
@ -36,17 +37,15 @@ global it_
disp(sprintf('%16.4g %16.4g %16.4g\n',z)) disp(sprintf('%16.4g %16.4g %16.4g\n',z))
options_ = set_default_option(options_,'qz_criterium',1.000001); options_ = set_default_option(options_,'qz_criterium',1.000001);
disp(sprintf('\nThere are %d eigenvalue(s) larger than 1 in modulus ', ... disp(sprintf('\nThere are %d eigenvalue(s) larger than 1 in modulus ', ...
nnz(abs(eigenvalues_) > options_.qz_criterium))); n_explod));
disp(sprintf('for %d forward-looking variable(s)',nyf)); disp(sprintf('for %d forward-looking variable(s)',nyf));
disp(' ') disp(' ')
if info(1) == 0 if dr.rank == nyf & nyf == n_explod
if dr.rank == nyf disp('The rank condition is verified.')
disp('The rank condition is verified.') else
else disp('The rank conditions ISN''T verified!')
disp('The rank conditions ISN''T verified!')
end
disp(' ')
end end
disp(' ')
end end
% keep lambda_ for backward compatibility % keep lambda_ for backward compatibility