From a2caa8b4a8e3bbc1be00c53e22f29e6878cdc4aa Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?S=C3=A9bastien=20Villemot?= Date: Tue, 24 Mar 2020 15:27:31 +0100 Subject: [PATCH] Block decomposition: in symbolic normalization, no longer insert a null pointer in dynamic Jacobian Rather use DataTree::Zero, since we use numerical zero for static and contemporanous Jacobians. --- src/ModelTree.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ModelTree.cc b/src/ModelTree.cc index c8832c7c..ae87ac98 100644 --- a/src/ModelTree.cc +++ b/src/ModelTree.cc @@ -363,7 +363,7 @@ ModelTree::computeNonSingularNormalization(jacob_map_t &contemporaneous_jacobian if (static_jacobian.find(eq_and_endo) == static_jacobian.end()) static_jacobian[eq_and_endo] = 0; if (dynamic_jacobian.find({ 0, eq_and_endo.first, eq_and_endo.second }) == dynamic_jacobian.end()) - dynamic_jacobian[{ 0, eq_and_endo.first, eq_and_endo.second }] = nullptr; + dynamic_jacobian[{ 0, eq_and_endo.first, eq_and_endo.second }] = Zero; if (contemporaneous_jacobian.find(eq_and_endo) == contemporaneous_jacobian.end()) contemporaneous_jacobian[eq_and_endo] = 0; try