OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-12 00:45:36 -06:00
Code Issues Packages Projects Releases Wiki Activity

Fix remaining bugs for incompressible flows driven by sources.

Browse Source 
In particular:

  - Zero sufficient portion of the ->mat_row to hold all connections
    of a cell in addition to the accumulation term.

  - Don't write the residual into the accumulation term of the
    Jacobian matrix row.
This commit is contained in:
Bård Skaflestad 2011-10-20 13:59:15 +02:00
parent 9392b3700f
commit afdb88678f
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/cfs_tpfa_residual.c
View File

@ -466,20 +466,20 @@ compute_cell_contrib(grid_t *G ,
/* Jacobian row */
vector_zero(1 + nconn, pimpl->ratio->mat_row);
vector_zero(1 + (G->cell_facepos[c + 1] - G->cell_facepos[c]),
pimpl->ratio->mat_row);
/* t2 <- A \ ((dA/dp) * t1) */
matvec(np, np, dAc, pimpl->ratio->t1, pimpl->ratio->t2);
solve_linear_systems(np, 1, pimpl->ratio, pimpl->ratio->t2);
dF1 = dF2 = 0.0;
dF2 = 0.0;
for (p = 0; p < np; p++) {
dF1 += pimpl->ratio->t1[ p ];
dF2 += pimpl->ratio->t2[ p ];
}
pimpl->is_incomp = pimpl->is_incomp && (! (fabs(dF2) > 0));
pimpl->ratio->mat_row[ 0 ] = dF1 - dF2;
pimpl->ratio->mat_row[ 0 ] = - dF2;
/* Accumulate inter-cell Jacobian contributions */
dv = pimpl->ratio->linsolve_buffer + (1 + nconn)*np;