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

trunk dynare_solve.m: for solve_algo = 2: 

* when solve1() fails on a block, immediately return with info=1 instead of trying to solve subsequent blocks
* added debug messages


git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@2078 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
sebastien 2008-09-16 17:24:11 +00:00
parent 30ffa48eec
commit e9c7351aee
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
matlab/dynare_solve.m
View File

@ -102,10 +102,17 @@ function [x,info] = dynare_solve(func,x,jacobian_flag,varargin)
[j1,j2,r,s] = dmperm(fjac);
if options_.debug
disp(['DYNARE_SOLVE (solve_algo=2): number of blocks = ' num2str(length(r))]);
end
for i=length(r)-1:-1:1
if options_.debug
disp(['DYNARE_SOLVE (solve_algo=2): solving block ' num2str(i) ', of size ' num2str(r(i+1)-r(i)) ]);
end
[x,info]=solve1(func,x,j1(r(i):r(i+1)-1),j2(r(i):r(i+1)-1),jacobian_flag,varargin{:});
if info & options_.debug
error(sprintf('Solve block = %d check = %d\n',i,info));
if info
return
end
end
fvec = feval(func,x,varargin{:});