From c47d91c20b8553e1d8b318a9d328d0763e6dc43e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= <sebastien@dynare.org>
Date: Thu, 19 Jan 2023 10:03:37 +0100
Subject: [PATCH] Block decomposition: another fix for sparse indices of static
Jacobian when mfs > 0
---
src/StaticModel.cc | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/StaticModel.cc b/src/StaticModel.cc
index 5c30f2ad..f3cfc2da 100644
--- a/src/StaticModel.cc
+++ b/src/StaticModel.cc
@@ -690,7 +690,7 @@ StaticModel::computeChainRuleJacobian()
if (eq >= nb_recursives && var >= nb_recursives)
blocks_jacobian_sparse_column_major_order[blk].emplace(pair{eq-nb_recursives, var-nb_recursives}, d1);
}
- blocks_jacobian_sparse_colptr[blk] = computeCSCColPtr(blocks_jacobian_sparse_column_major_order[blk], size);
+ blocks_jacobian_sparse_colptr[blk] = computeCSCColPtr(blocks_jacobian_sparse_column_major_order[blk], blocks[blk].mfs_size);
}
}
}
--
GitLab