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

Solve(rock_comp): Refactor rock-compressibility diagonal term calculation.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-03-20 11:21:02 +01:00
parent f91ec88df6
commit a8d067af18
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

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