diff options
| -rw-r--r-- | test/h5test.h | 15 | ||||
| -rw-r--r-- | test/hyperslab.c | 6 |
2 files changed, 20 insertions, 1 deletions
diff --git a/test/h5test.h b/test/h5test.h index 640a882..8187bcf 100644 --- a/test/h5test.h +++ b/test/h5test.h @@ -117,6 +117,21 @@ extern MPI_Info h5_io_info_g; /* MPI INFO object for IO */ /* set alarms to N seconds if N > 0, else use default alarm_seconds. */ #define ALARM_SET(N) HDalarm((N)>0 ? N : alarm_seconds) +/* + * The methods to compare the equality of floating-point values: + * 1. XXX_ABS_EQUAL - check if the difference is smaller than the + * Epsilon value. The Epsilon values, FLT_EPSILON, DBL_EPSILON, + * and LDBL_EPSILON, are defined by compiler in float.h. + * 2. To be defined later. + */ +#define FLT_ABS_EQUAL(X,Y) (fabsf(X-Y)<FLT_EPSILON) +#define DBL_ABS_EQUAL(X,Y) (fabs(X-Y)<DBL_EPSILON) +#define LDBL_ABS_EQUAL(X,Y) (fabsl(X-Y)<LDBL_EPSILON) + +/*#define FP_ABS_UNEQUAL(X,Y,T) ((T==1 && fabsf(X-Y)>FLT_EPSILON) || \ + (T==2 && fabs(X-Y)>DBL_EPSILON) || \ + (T==3 && fabsl(X-Y)>LDBL_EPSILON))*/ + #ifdef __cplusplus extern "C" { #endif diff --git a/test/hyperslab.c b/test/hyperslab.c index 0a46faf..16bd1bb 100644 --- a/test/hyperslab.c +++ b/test/hyperslab.c @@ -663,7 +663,11 @@ test_multifill(size_t nx) for (i = 0; i < nx; i++) { if (dst[i].left != 3333333) { sprintf(s, "bad dst[%lu].left", (unsigned long)i); - } else if (dst[i].mid != fill.mid) { + } else if (!DBL_ABS_EQUAL(dst[i].mid, fill.mid)) { + /* Check if two DOUBLE values are equal. If their difference + * is smaller than the EPSILON value for double, they are + * considered equal. See the definition in h5test.h. + */ sprintf(s, "bad dst[%lu].mid", (unsigned long)i); } else if (dst[i].right != 4444444) { sprintf(s, "bad dst[%lu].right", (unsigned long)i); |