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

Fix indexing bug in well handling code for incremental (Newton) pressure assembly.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-10 15:31:51 +02:00
parent f93cbf3c4a
commit 7eb2ed764e
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/tpfa/ifs_tpfa.c
View File

@ -760,7 +760,7 @@ ifs_tpfa_assemble_comprock_increment(struct UnstructuredGrid *G ,
struct ifs_tpfa_data *h )
/* ---------------------------------------------------------------------- */
{
int c, w, system_singular, ok;
int c, w, wdof, system_singular, ok;
size_t j;
double *v;
@ -782,7 +782,8 @@ ifs_tpfa_assemble_comprock_increment(struct UnstructuredGrid *G ,
}
if (F->W != NULL) {
for (w = 0; w < F->W->number_of_wells; ++w) {
h->b[c] += -v[c];
wdof = G->number_of_cells + w;
h->b[wdof] += -v[wdof];
}
}
}