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

Cosmetic changes.

time-shift
Stéphane Adjemian (Charybdis) 2014-09-10 16:30:17 +02:00
parent 308ce25031
commit ce97a27c9d
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
matlab/dynare_solve.m
View File

@ -8,7 +8,7 @@ function [x,info] = dynare_solve(func,x,jacobian_flag,varargin)
% jacobian_flag=1: jacobian given by the 'func' function % jacobian_flag=1: jacobian given by the 'func' function
% jacobian_flag=0: jacobian obtained numerically % jacobian_flag=0: jacobian obtained numerically
% varargin: list of arguments following jacobian_flag % varargin: list of arguments following jacobian_flag
% %
% OUTPUTS % OUTPUTS
% x: solution % x: solution
% info=1: the model can not be solved % info=1: the model can not be solved
@ -47,7 +47,7 @@ if jacobian_flag
M=evalin('base','M_'); %get variable names from workspace M=evalin('base','M_'); %get variable names from workspace
fprintf('\nSTEADY: The Jacobian contains Inf or NaN. The problem arises from: \n\n') fprintf('\nSTEADY: The Jacobian contains Inf or NaN. The problem arises from: \n\n')
display_problematic_vars_Jacobian(infrow,infcol,M,x,'static','STEADY: ') display_problematic_vars_Jacobian(infrow,infcol,M,x,'static','STEADY: ')
error('An element of the Jacobian is not finite or NaN') error('An element of the Jacobian is not finite or NaN')
end end
else else
fvec = feval(func,x,varargin{:}); fvec = feval(func,x,varargin{:});
@ -62,7 +62,7 @@ if ~isempty(i)
disp(i') disp(i')
disp('Please check for example') disp('Please check for example')
disp(' i) if all parameters occurring in these equations are defined') disp(' i) if all parameters occurring in these equations are defined')
disp(' ii) that no division by an endogenous variable initialized to 0 occurs') disp(' ii) that no division by an endogenous variable initialized to 0 occurs')
info = 1; info = 1;
x = NaN(size(fvec)); x = NaN(size(fvec));
return; return;
@ -81,8 +81,8 @@ if options_.solve_algo == 0
options=optimset('fsolve'); options=optimset('fsolve');
options.MaxFunEvals = 50000; options.MaxFunEvals = 50000;
options.MaxIter = options_.steady.maxit; options.MaxIter = options_.steady.maxit;
options.TolFun=tolf; options.TolFun = tolf;
options.Display = 'iter'; options.Display = 'iter';
if jacobian_flag if jacobian_flag
options.Jacobian = 'on'; options.Jacobian = 'on';
else else
@ -102,7 +102,7 @@ if options_.solve_algo == 0
exitval = 3; exitval = 3;
end; end;
end end
if exitval > 0 if exitval > 0
info = 0; info = 0;
else else
@ -135,7 +135,7 @@ elseif options_.solve_algo == 2 || options_.solve_algo == 4
end end
[j1,j2,r,s] = dmperm(fjac); [j1,j2,r,s] = dmperm(fjac);
if options_.debug if options_.debug
disp(['DYNARE_SOLVE (solve_algo=2|4): number of blocks = ' num2str(length(r))]); disp(['DYNARE_SOLVE (solve_algo=2|4): number of blocks = ' num2str(length(r))]);
end end