Bug fix (solve_algo={12,14}).

Removed threshold for detecting the non zero elements in the rows of
the Jacobian matrix. Using tolf as a threshold parameter for identifying
the non zero elements leaded (not systematically) the algorithm to not
reevaluate the residuals of the dynamic model while necessary.
time-shift
Stéphane Adjemian (Charybdis) 2020-03-11 12:38:59 +01:00
parent 492b0fffea
commit a4519f6b86
Signed by: stepan
GPG Key ID: 295C1FE89E17EB3C
1 changed files with 1 additions and 1 deletions

View File

@ -243,7 +243,7 @@ elseif ismember(options.solve_algo, [2, 12, 4, 14])
end
end
[j1,j2,r,s] = dmperm(fjac);
JAC = abs(fjac(j1,j2))>tolf;
JAC = abs(fjac(j1,j2))>0;
if options.debug
disp(['DYNARE_SOLVE (solve_algo=2|4|12|14): number of blocks = ' num2str(length(r))]);
end