From bbdcdc89a28ecedc2f8215113bb82649fc6f3bc4 Mon Sep 17 00:00:00 2001 From: Johannes Pfeifer Date: Mon, 28 Nov 2022 16:19:05 +0100 Subject: [PATCH] sim1.m: use operation that preserves sparsity --- matlab/perfect-foresight-models/sim1.m | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/matlab/perfect-foresight-models/sim1.m b/matlab/perfect-foresight-models/sim1.m index 7f9fc84c6..8f0f65a90 100644 --- a/matlab/perfect-foresight-models/sim1.m +++ b/matlab/perfect-foresight-models/sim1.m @@ -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)