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

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.
issue#70
Sébastien Villemot 2020-03-24 15:27:31 +01:00
parent 5d6e6338f4
commit a2caa8b4a8
No known key found for this signature in database
GPG Key ID: 2CECE9350ECEBE4A
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/ModelTree.cc
View File

@ -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