OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

In collapseJacs(), avoid explicit zeros.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2015-08-26 13:29:46 +02:00 committed by babrodtk
parent dd52d444e7
commit 51dde7cec5
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/autodiff/AutoDiffHelpers.hpp
View File

@ -464,9 +464,11 @@ collapseJacs(const AutoDiffBlock<double>& x, Matrix& jacobian)
jac1.toSparse(jac); jac1.toSparse(jac);
for (Eigen::SparseMatrix<double>::Index k = 0; k < jac.outerSize(); ++k) { for (Eigen::SparseMatrix<double>::Index k = 0; k < jac.outerSize(); ++k) {
for (Eigen::SparseMatrix<double>::InnerIterator i(jac, k); i ; ++i) { for (Eigen::SparseMatrix<double>::InnerIterator i(jac, k); i ; ++i) {
t.push_back(Tri(i.row(), if (i.value() != 0.0) {
i.col() + block_col_start, t.push_back(Tri(i.row(),
i.value())); i.col() + block_col_start,
i.value()));
}
} }
} }
block_col_start += jac.cols(); block_col_start += jac.cols();