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

model_diagnostics: remove workaround for Octave 

The rank() and null() functions accept a tolerance argument since at least
Octave 3.8.
silicon
Sébastien Villemot 2023-03-16 16:02:10 +01:00
parent 441173ec5a
commit 060c694066
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
matlab/model_diagnostics.m
View File

@ -183,7 +183,7 @@ for b=1:nb
display_problematic_vars_Jacobian(imagrow,imagcol,M,dr.ys,'static','MODEL_DIAGNOSTICS: ')
end
try
if isoctave || matlab_ver_less_than('9.12') || isempty(options.jacobian_tolerance)
if (~isoctave && matlab_ver_less_than('9.12')) || isempty(options.jacobian_tolerance)
rank_jacob = rank(jacob); %can sometimes fail
else
rank_jacob = rank(jacob,options.jacobian_tolerance); %can sometimes fail
@ -198,7 +198,7 @@ for b=1:nb
'singular'])
disp(['MODEL_DIAGNOSTICS: there is ' num2str(n_vars_jacob-rank_jacob) ...
' collinear relationships between the variables and the equations'])
if isoctave || matlab_ver_less_than('9.12') || isempty(options.jacobian_tolerance)
if (~isoctave && matlab_ver_less_than('9.12')) || isempty(options.jacobian_tolerance)
ncol = null(jacob);
else
ncol = null(jacob,options.jacobian_tolerance); %can sometimes fail
@ -221,7 +221,7 @@ for b=1:nb
fprintf('%s\n',endo_names{k})
end
end
if isoctave || matlab_ver_less_than('9.12') || isempty(options.jacobian_tolerance)
if (~isoctave && matlab_ver_less_than('9.12')) || isempty(options.jacobian_tolerance)
neq = null(jacob'); %can sometimes fail
else
neq = null(jacob',options.jacobian_tolerance); %can sometimes fail