OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

solve(rock_comp): Replace manual index search with csrmatrix_elm_index().

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-19 13:55:21 +01:00
parent bc04c0fa3e
commit 3d3aaa98d1
1 changed files with 1 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/core/pressure/IncompTpfa.cpp
View File

@ -186,12 +186,7 @@ namespace Opm
// The p part goes on the diagonal, the p0 on the rhs.
for (int c = 0; c < gg->number_of_cells; ++c) {
// Find diagonal
int j = h_->A->ia[c];
for (; j < h_->A->ia[c+1]; ++j) {
if (h_->A->ja[j] == c) {
break;
}
}
size_t j = csrmatrix_elm_index(c, c, h_->A);
h_->A->sa[j] += porevol[c]*rock_comp[c]/dt;
h_->b[c] += porevol[c]*rock_comp[c]*pressure[c]/dt;
}