Skip to content
Snippets Groups Projects
Verified Commit c47d91c2 authored by Sébastien Villemot's avatar Sébastien Villemot
Browse files

Block decomposition: another fix for sparse indices of static Jacobian when mfs > 0

parent 7db02b23
No related branches found
No related tags found
No related merge requests found
...@@ -690,7 +690,7 @@ StaticModel::computeChainRuleJacobian() ...@@ -690,7 +690,7 @@ StaticModel::computeChainRuleJacobian()
if (eq >= nb_recursives && var >= nb_recursives) 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_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);
} }
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment