summaryrefslogtreecommitdiffstats
path: root/test
diff options
context:
space:
mode:
authorRobb Matzke <matzke@llnl.gov>1998-07-06 23:02:32 (GMT)
committerRobb Matzke <matzke@llnl.gov>1998-07-06 23:02:32 (GMT)
commitcf419d05ab65b423bc4ce395972e84522ca3eec5 (patch)
tree6f689a9e56f4e251f51aacec02a067b9f8d232a1 /test
parent4d3351f1d777bc3cb891443b6165a87251d4e643 (diff)
downloadhdf5-cf419d05ab65b423bc4ce395972e84522ca3eec5.zip
hdf5-cf419d05ab65b423bc4ce395972e84522ca3eec5.tar.gz
hdf5-cf419d05ab65b423bc4ce395972e84522ca3eec5.tar.bz2
[svn-r454] Changes since 19980702
---------------------- ./src/H5T.c Prints statistics for no-op conversions also. For now, hardware floating point conversions may or may not raise SIGFPE depending on the hardware, operating system, etc. Software conversions never raise SIGFPE and use +Inf or -Inf for overflow. ./test/dtypes.c Catches SIGFPE and causes the test to be skipped. Better test for NaN. ./config/irix5.3 Removed the -U_POSIX_SOURCE because it was removed from the main makefiles. ./bin/trace ./src/H5S.c Fixed a typo that prevented tracing info from being added to new API functions.
Diffstat (limited to 'test')
-rw-r--r--test/.distdep4
-rw-r--r--test/dtypes.c211
2 files changed, 161 insertions, 54 deletions
diff --git a/test/.distdep b/test/.distdep
index a2fcc62..5706161 100644
--- a/test/.distdep
+++ b/test/.distdep
@@ -182,7 +182,9 @@ dtypes.o: \
../src/H5Tpkg.h \
../src/H5HGprivate.h \
../src/H5Fprivate.h \
- ../src/H5private.h
+ ../src/H5private.h \
+ ../src/H5Tprivate.h \
+ ../src/H5Gprivate.h
hyperslab.o: \
hyperslab.c \
../src/H5private.h \
diff --git a/test/dtypes.c b/test/dtypes.c
index 2cab029..e9e8f06 100644
--- a/test/dtypes.c
+++ b/test/dtypes.c
@@ -8,12 +8,13 @@
* Purpose: Tests the data type interface (H5T)
*/
#include <assert.h>
-#include <float.h>
#include <hdf5.h>
#include <math.h>
#include <signal.h>
#include <stdio.h>
#include <stdlib.h>
+#include <sys/types.h>
+#include <sys/wait.h>
#include <unistd.h>
#define H5T_PACKAGE
@@ -52,6 +53,39 @@ typedef enum flt_t {
FLT_FLOAT, FLT_DOUBLE, FLT_LDOUBLE, FLT_OTHER
} flt_t;
+/*
+ * Some machines generate SIGFPE on floating point overflows. According to
+ * the Posix standard, we cannot assume that we can continue from such a
+ * signal. Therefore, if the following constant is defined then tests that
+ * might raise SIGFPE are executed in a child process.
+ */
+#define HANDLE_SIGFPE
+
+
+/*-------------------------------------------------------------------------
+ * Function: fpe_handler
+ *
+ * Purpose: Exit with 255
+ *
+ * Return: void
+ *
+ * Programmer: Robb Matzke
+ * Monday, July 6, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+#ifdef HANDLE_SIGFPE
+static void
+fpe_handler(int __unused__ signo)
+{
+ puts(" -SKIP-");
+ puts(" Test skipped due to SIGFPE from probable overflow.");
+ exit(255);
+}
+#endif
+
/*-------------------------------------------------------------------------
* Function: cleanup
@@ -675,6 +709,74 @@ test_conv_int (void)
/*-------------------------------------------------------------------------
+ * Function: my_isnan
+ *
+ * Purpose: Determines whether VAL points to NaN.
+ *
+ * Return: TRUE or FALSE
+ *
+ * Programmer: Robb Matzke
+ * Monday, July 6, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+static int
+my_isnan(flt_t type, void *val)
+{
+ int retval;
+ char s[256];
+
+ switch (type) {
+ case FLT_FLOAT:
+ retval = (*((float*)val)!=*((float*)val));
+ break;
+
+ case FLT_DOUBLE:
+ retval = (*((double*)val)!=*((double*)val));
+ break;
+
+ case FLT_LDOUBLE:
+ retval = (*((long double*)val)!=*((long double*)val));
+ break;
+
+ default:
+ return 0;
+ }
+
+ /*
+ * Sometimes NaN==NaN (e.g., DEC Alpha) so we try to print it and see if
+ * the result contains a NaN string.
+ */
+ if (!retval) {
+ switch (type) {
+ case FLT_FLOAT:
+ sprintf(s, "%g", *((float*)val));
+ break;
+
+ case FLT_DOUBLE:
+ sprintf(s, "%g", *((double*)val));
+ break;
+
+ case FLT_LDOUBLE:
+ sprintf(s, "%Lg", *((long double*)val));
+ break;
+
+ default:
+ return 0;
+ }
+
+ if (!strstr(s, "NaN") || !strstr(s, "NAN") || !strstr(s, "nan")) {
+ retval = 1;
+ }
+ }
+
+ return retval;
+}
+
+
+/*-------------------------------------------------------------------------
* Function: test_conv_flt_1
*
* Purpose: Test conversion of random floating point values from SRC to
@@ -714,6 +816,39 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
size_t i, j, k; /*counters */
int endian; /*machine endianess */
+#ifdef HANDLE_SIGFPE
+ pid_t child_pid; /*process ID of child */
+ int status; /*child exit status */
+
+ /*
+ * Some systems generage SIGFPE during floating point overflow and we
+ * cannot assume that we can continue from such a signal. Therefore, we
+ * fork here and let the child run the test and return the number of
+ * failures with the exit status.
+ */
+ fflush(stdout);
+ fflush(stderr);
+ if ((child_pid=fork())<0) {
+ perror("fork");
+ return 1;
+ } else if (child_pid>0) {
+ while (child_pid!=waitpid(child_pid, &status, 0)) /*void*/;
+ if (WIFEXITED(status) && 255==WEXITSTATUS(status)) {
+ return 0; /*child exit after catching SIGFPE*/
+ } else if (WIFEXITED(status)) {
+ return WEXITSTATUS(status);
+ } else {
+ puts(" Child didn't exit normally.");
+ return 1;
+ }
+ }
+
+ /*
+ * The remainder of this function is executed only by the child.
+ */
+ signal(SIGFPE,fpe_handler);
+#endif
+
/* What are the names of the source and destination types */
if (H5Tequal(src, H5T_NATIVE_FLOAT)) {
src_type_name = "float";
@@ -785,6 +920,9 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
/* Check the software results against the hardware */
for (j=0; j<nelmts; j++) {
+ hw_f = 911.0;
+ hw_d = 911.0;
+ hw_ld = 911.0;
/* The hardware conversion */
if (FLT_FLOAT==src_type) {
@@ -800,14 +938,7 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
}
} else if (FLT_DOUBLE==src_type) {
if (FLT_FLOAT==dst_type) {
- /* Watch out for that FPE on overflow! */
- if (((double*)saved)[j] > FLT_MAX) {
- hw_f = FLT_MAX;
- } else if (((double*)saved)[j] < -FLT_MAX) {
- hw_f = -FLT_MAX;
- } else {
- hw_f = ((double*)saved)[j];
- }
+ hw_f = ((double*)saved)[j];
hw = (unsigned char*)&hw_f;
} else if (FLT_DOUBLE==dst_type) {
hw_d = ((double*)saved)[j];
@@ -818,24 +949,10 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
}
} else {
if (FLT_FLOAT==dst_type) {
- /* Watch out for that FPE on overflow! */
- if (((long double*)saved)[j] > FLT_MAX) {
- hw_f = FLT_MAX;
- } else if (((long double*)saved)[j] < -FLT_MAX) {
- hw_f = -FLT_MAX;
- } else {
- hw_f = ((long double*)saved)[j];
- }
+ hw_f = ((long double*)saved)[j];
hw = (unsigned char*)&hw_f;
} else if (FLT_DOUBLE==dst_type) {
- /* Watch out for that FPE! */
- if (((long double*)saved)[j] > DBL_MAX) {
- hw_d = DBL_MAX;
- } else if (((long double*)saved)[j] < -DBL_MAX) {
- hw_d = -DBL_MAX;
- } else {
- hw_d = ((long double*)saved)[j];
- }
+ hw_d = ((long double*)saved)[j];
hw = (unsigned char*)&hw_d;
} else {
hw_ld = ((long double*)saved)[j];
@@ -857,16 +974,16 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
* pattern for NaN by setting the significand to all ones.
*/
if (FLT_FLOAT==dst_type &&
- ((float*)buf)[j]!=((float*)buf)[j] &&
- ((float*)hw)[0]!=((float*)hw)[0]) {
+ my_isnan(dst_type, (float*)buf+j) &&
+ my_isnan(dst_type, hw)) {
continue;
} else if (FLT_DOUBLE==dst_type &&
- ((double*)buf)[j]!=((double*)buf)[j] &&
- ((double*)hw)[0]!=((double*)hw)[0]) {
+ my_isnan(dst_type, (double*)buf+j) &&
+ my_isnan(dst_type, hw)) {
continue;
} else if (FLT_LDOUBLE==dst_type &&
- ((long double*)buf)[j]!=((long double*)buf)[j] &&
- ((long double*)hw)[0]!=((long double*)hw)[0]) {
+ my_isnan(dst_type, (long double*)buf+j) &&
+ my_isnan(dst_type, hw)) {
continue;
}
#endif
@@ -877,16 +994,7 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
* hardware conversions on some machines return NaN instead of
* overflowing to +Inf or -Inf or underflowing to +0 or -0.
*/
- if (FLT_FLOAT==dst_type &&
- *((float*)hw)!=*((float*)hw)) {
- continue;
- } else if (FLT_DOUBLE==dst_type &&
- *((double*)hw)!=*((double*)hw)) {
- continue;
- } else if (FLT_LDOUBLE==dst_type &&
- *((long double*)hw)!=*((long double*)hw)) {
- continue;
- }
+ if (my_isnan(dst_type, hw)) continue;
#endif
#if 1
@@ -970,12 +1078,20 @@ test_conv_flt_1 (const char *name, hid_t src, hid_t dst)
done:
if (buf) free (buf);
if (saved) free (saved);
+#ifdef HANDLE_SIGFPE
+ exit(MIN((int)fails_all_tests, 254));
+#else
return (int)fails_all_tests;
+#endif
error:
if (buf) free (buf);
if (saved) free (saved);
- return (int)MIN(1, fails_all_tests);
+#ifdef HANDLE_SIGFPE
+ exit(MIN(MAX((int)fails_all_tests, 1), 254));
+#else
+ return MAX((int)fails_all_tests, 1);
+#endif
}
@@ -1000,10 +1116,6 @@ main(void)
{
unsigned long nerrors = 0;
-#if 0
- signal(SIGFPE,SIG_IGN);
-#endif
-
/* Set the error handler */
H5Eset_auto (display_error_cb, NULL);
@@ -1015,12 +1127,6 @@ main(void)
nerrors += test_named ()<0 ? 1 : 0;
nerrors += test_conv_int ()<0 ? 1 : 0;
-#ifndef LATER
- /*
- * NOT READY FOR TESTING YET BECAUSE SOME SYSTEMS GENERATE A SIGFPE WHEN
- * AN OVERFLOW OCCURS CASTING A DOUBLE TO A FLOAT.
- */
-#else
/* Test degenerate cases */
nerrors += test_conv_flt_1("noop", H5T_NATIVE_FLOAT, H5T_NATIVE_FLOAT);
nerrors += test_conv_flt_1("noop", H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE);
@@ -1038,7 +1144,6 @@ main(void)
nerrors += test_conv_flt_1("sw", H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE);
nerrors += test_conv_flt_1("sw", H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT);
nerrors += test_conv_flt_1("sw", H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE);
-#endif
if (nerrors) {
printf("***** %lu FAILURE%s! *****\n",