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

Exit if initial guess is a solution... 

Without entering in the nonlinear solvers.
fix-nonlinear-solvers
Stéphane Adjemian (Ryûk) 2022-04-04 15:35:32 +02:00
parent b236027365
commit 56863d6e2f
Signed by: stepan
GPG Key ID: 295C1FE89E17EB3C
1 changed files with 9 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
matlab/dynare_solve.m
View File

@ -92,13 +92,13 @@ end
if jacobian_flag
[fvec, fjac] = feval(f, x, arguments{:});
wrong_initial_guess_flag = false;
if ~all(isfinite(fvec)) || any(isinf(fjac(:))) || any(isnan((fjac(:)))) ...
|| any(~isreal(fvec)) || any(~isreal(fjac(:)))
if max(abs(fvec)) < tolf %return if initial value solves problem
if ~all(isfinite(fvec)) || any(isinf(fjac(:))) || any(isnan((fjac(:)))) || any(~isreal(fvec)) || any(~isreal(fjac(:)))
if max(abs(fvec)) < tolf
% return if initial value solves problem
exitflag = -1;
return;
end
disp_verbose('Randomize initial guess...',options.verbosity)
disp_verbose('Randomize initial guess...', options.verbosity)
% Let's try random numbers for the variables initialized with the default value.
wrong_initial_guess_flag = true;
% First try with positive numbers.
@ -129,6 +129,11 @@ if jacobian_flag
else
fvec = feval(f, x, arguments{:});
fjac = zeros(nn, nn);
if max(abs(fvec)) < tolf
% return if initial value solves problem
exitflag = -1;
return;
end
wrong_initial_guess_flag = false;
if ~all(isfinite(fvec))
% Let's try random numbers for the variables initialized with the default value.