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

Made floating point tests equal to test in old code

Browse Source
This commit is contained in:
Bjørn Spjelkavik 2009-08-14 11:12:15 +00:00
parent 987d59c92c
commit 756454f0fc
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
uniquepoints.c
View File

@ -90,16 +90,14 @@ static int uniquify(int n, double *list, double tolerance)
double val = list[pos++];/* Keep first value */
for (i=1; i<n; ++i){
if (val + tolerance < list [i]){
if (list [i] - val > tolerance){
val = list[i];
list[pos++] = val;
}
}
/* Keep last value (one way or the other...) */
if (val + tolerance < list[n-1]){
list[pos-1] = list[n-1];
}
return pos;
}
@ -141,12 +139,12 @@ static int assignPointNumbers(int begin,
}
/* Find next k such that zlist[k] < z[i] < zlist[k+1] */
while ((k < end) && (zlist[k] + tolerance < z[i])){
while ((k < end) && (z[i] - zlist[k] > tolerance)){
k++;
}
/* assert (k < len && z[i] - zlist[k] <= tolerance) */
if ((k == end) || ( zlist[k] + tolerance < z[i])){
if ((k == end) || ( z[i] - zlist[k] > tolerance)){
fprintf(stderr, "Cannot associate zcorn values with given list\n");
fprintf(stderr, "of z-coordinates to given tolerance\n");
return 0;