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

sim1.m: use operation that preserves sparsity

unit-tests
Johannes Pfeifer 2022-11-28 16:19:05 +01:00
parent ba9b0e9829
commit bbdcdc89a2
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
matlab/perfect-foresight-models/sim1.m
View File

@ -68,8 +68,9 @@ for iter = 1:options.simul.maxit
% A is the stacked Jacobian with period x equations alongs the rows and
% periods times variables (in declaration order) along the columns
if options.debug && iter==1
row=find(all(A==0,2));
column=find(all(A==0,1));
[row,col]=find(A);
row=setdiff(1:periods*ny,row);
column=setdiff(1:periods*ny,col);
if ~isempty(row) || ~isempty(column)
fprintf('The stacked Jacobian is singular. The problem derives from:\n')
if ~isempty(row)