This commit is contained in:
Atgeirr Flø Rasmussen 2012-03-20 12:12:21 +01:00
commit 39e55a34cf

View File

@ -193,8 +193,9 @@ namespace Opm
for (int c = 0; c < gg->number_of_cells; ++c) {
// Find diagonal
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;
double d = porevol[c] * rock_comp[c] / dt;
h_->A->sa[j] += d;
h_->b[c] += d * pressure[c];
}
}