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

Correct a nasty bug where we were not correctly defining new conn

Browse Source 
pointers in the first cell (i.e., c==0) in add_well_connections().
  With this change, computing the number of existing connections is
  split across loop iterations, but the new approach has the added
  benefit of actually being correct.

  Preliminary testing is promising and suggests this is on the right
  track, but we still do not re-enable the MEX code.
This commit is contained in:
Bård Skaflestad 2010-08-11 22:35:20 +00:00
parent fc082bcc9a
commit 9fdbc97459
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
mex_ip_simple.c
View File

@ -285,10 +285,11 @@ add_well_connections(const mxArray *W, int nc, int nf,
tmp = mxRealloc(*conn, (pconn[nc] + neconn) * sizeof *tmp);
if (tmp != NULL) {
n = pconn[nc] - pconn[nc - 1]; /* # exisiting conns */
pconn[nc] += neconn; /* ubnd, *LAST* bin */
n = pconn[nc]; /* Preserve original ubound */
pconn[nc] += neconn; /* New ubound, *LAST* bin */
for (c = nc - 1; c > 0; c--) {
for (c = nc - 1; c >= 0; c--) {
n -= pconn[c]; /* => n == # exisiting conns */
dst = pconn[c + 1] - (n + cwork[c]);
src = pconn[c + 0];
@ -299,7 +300,7 @@ add_well_connections(const mxArray *W, int nc, int nf,
/* Set cell's pointer for new connections. */
cwork[c] = dst + n;
n = pconn[c] - pconn[c - 1];
n = pconn[c]; /* Preserve original ubound */
pconn[c] = dst;
}