OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-18 21:43:27 -06:00
Code Issues Packages Projects Releases Wiki Activity

comprock_increment(): Don't proceed to mult_csr() unless OK.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-05-10 22:29:41 +02:00
parent a4b32311b2
commit a865f01478
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/core/pressure/tpfa/ifs_tpfa.c
View File

@ -766,15 +766,15 @@ ifs_tpfa_assemble_comprock_increment(struct UnstructuredGrid *G ,
assemble_incompressible(G, F, trans, gpress, h, &system_singular, &ok); assemble_incompressible(G, F, trans, gpress, h, &system_singular, &ok);
v = h->pimpl->work;
mult_csr_matrix(h->A, prev_pressure, v);
/* We want to solve a Newton step for the residual /* We want to solve a Newton step for the residual
* (porevol(pressure)-porevol(initial_pressure))/dt + residual_for_imcompressible * (porevol(pressure)-porevol(initial_pressure))/dt + residual_for_imcompressible
* *
*/ */
if (ok) { if (ok) {
v = h->pimpl->work;
mult_csr_matrix(h->A, prev_pressure, v);
for (c = 0; c < G->number_of_cells; ++c) { for (c = 0; c < G->number_of_cells; ++c) {
j = csrmatrix_elm_index(c, c, h->A); j = csrmatrix_elm_index(c, c, h->A);
h->A->sa[j] += porevol[c] * rock_comp[c] / dt; h->A->sa[j] += porevol[c] * rock_comp[c] / dt;