Account for flow-direction signs when accumulating fluxes.

This commit is contained in:
Bård Skaflestad 2011-10-19 22:10:22 +02:00
parent fb1316802c
commit 98abee6a79

View File

@ -418,7 +418,7 @@ init_cell_contrib(grid_t *G ,
cflx += 1 * np; cflx += 1 * np;
dcflx += 2 * np; dcflx += 2 * np;
pimpl->ratio->coeff[ conn++ ] = dt; pimpl->ratio->coeff[ conn++ ] = dt * (2*(c1 == c) - 1.0);
} }
} }
@ -516,7 +516,7 @@ assemble_sources(struct compr_src *src, struct cfs_tpfa_res_data *h)
assert (src->cell[i] >= 0 ); assert (src->cell[i] >= 0 );
assert (((size_t) src->cell[i]) < h->J->m); assert (((size_t) src->cell[i]) < h->J->m);
h->F[ src->cell[ i ] ] += src->flux[ i ]; h->F[ src->cell[ i ] ] -= src->flux[ i ];
} }
} }