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

Implement prototype TPFA discretisation for local problems.

Browse Source 
Contingent upon preprocessor symbol "USE_MIM_IP_TPFA" and only a
mimetic two-point discretisation (not cell-centred systems).  Hide
existing discretisation behind preprocessor symbol "USE_MIM_IP_SIMPLE".
This commit is contained in:
Bård Skaflestad 2011-06-21 14:02:39 +02:00
parent e15bdb5d51
commit aa4f77e0db
1 changed files with 44 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

44
src/coarse_sys.c
View File

@ -470,7 +470,11 @@ coarse_sys_meta_construct(grid_t *g, const int *p,
return m;
}
#define USE_MIM_IP_SIMPLE 0
#define USE_MIM_IP_TPFA 1
#define USE_MIM_IP_QFAMILY 0
#if USE_MIM_IP_SIMPLE
/* ---------------------------------------------------------------------- */
static double *
compute_fs_ip(grid_t *g, const double *perm,
@ -491,6 +495,46 @@ compute_fs_ip(grid_t *g, const double *perm,
return Binv;
}
#elif USE_MIM_IP_TPFA
#include "trans_tpfa.h"
/* ---------------------------------------------------------------------- */
static double *
compute_fs_ip(grid_t *g, const double *perm,
const struct coarse_sys_meta *m)
/* ---------------------------------------------------------------------- */
{
size_t c, nc, nconn, p1, p2, i, j;
double *Binv, *htrans;
nc = g->number_of_cells;
Binv = malloc(m->sum_ngconn2 * sizeof *Binv);
htrans = malloc(g->cell_facepos[ nc ] * sizeof *htrans);
if ((Binv != NULL) && (htrans != NULL)) {
tpfa_htrans_compute(g, perm, htrans);
for (c = p2 = 0; c < nc; c++) {
p1 = g->cell_facepos[c + 0] ;
nconn = g->cell_facepos[c + 1] - p1;
for (i = 0; i < nconn; i++) {
Binv[p2 + i*(nconn + 1)] = htrans[p1 + i];
for (j = (i + 1) % nconn; j != i; j = (j + 1) % nconn) {
Binv[p2 + i*nconn + j] = 0.0;
}
}
p2 += nconn * nconn;
}
}
free(htrans);
return Binv;
}
#endif
/* Create basis function weighting source term (unsigned) based on