summaryrefslogtreecommitdiffstats
path: root/test
diff options
context:
space:
mode:
authorRaymond Lu <songyulu@hdfgroup.org>2005-04-05 20:25:21 (GMT)
committerRaymond Lu <songyulu@hdfgroup.org>2005-04-05 20:25:21 (GMT)
commit138abc32798ac465e600965cf63aba760d968f88 (patch)
tree76071c4ef3e331df158b83ce05bcc65c92ff6124 /test
parentfdbeca1316b313c319cfe3abe4ca254ad2095f1a (diff)
downloadhdf5-138abc32798ac465e600965cf63aba760d968f88.zip
hdf5-138abc32798ac465e600965cf63aba760d968f88.tar.gz
hdf5-138abc32798ac465e600965cf63aba760d968f88.tar.bz2
[svn-r10550] Purpose: Splitting dtypes.c
Description: Yesterday, dtypes.c was split into two programs, dtypes.c and dt_atomic.c. After dt_atomic.c was added, dtypes.c remained the same as before waiting for the daily test to pass. While dt_atomic.c turned out fine today, the test in dt_atomic.c was taken out from dtypes.c. There's also some minor change in dt_atomic.c. Platforms tested: h5committest
Diffstat (limited to 'test')
-rw-r--r--test/dt_atomic.c104
-rw-r--r--test/dtypes.c4667
2 files changed, 129 insertions, 4642 deletions
diff --git a/test/dt_atomic.c b/test/dt_atomic.c
index e1386f3..f60be7c 100644
--- a/test/dt_atomic.c
+++ b/test/dt_atomic.c
@@ -69,6 +69,20 @@ typedef enum dtype_t {
FLT_LDOUBLE, OTHER
} dtype_t;
+/* Skip overflow tests if non-zero */
+static int skip_overflow_tests_g = 0;
+
+/*
+ * Although we check whether a floating point overflow generates a SIGFPE and
+ * turn off overflow tests in that case, it might still be possible for an
+ * overflow condition to occur. Once a SIGFPE is raised the program cannot
+ * be allowed to continue (cf. Posix signals) so in order to recover from a
+ * SIGFPE we run tests that might generate one in a child process.
+ */
+#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID)
+#define HANDLE_SIGFPE
+#endif
+
/* Don't use hardware conversions if set */
static int without_hardware_g = 0;
@@ -314,6 +328,7 @@ static int without_hardware_g = 0;
free(value); \
}
+void some_dummy_func(float x);
static hbool_t overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits);
static int my_isnan(dtype_t type, void *val);
@@ -433,6 +448,92 @@ except_func(int except_type, hid_t UNUSED src_id, hid_t UNUSED dst_id, void UNUS
/*-------------------------------------------------------------------------
+ * Function: some_dummy_func
+ *
+ * Purpose: A dummy function to help check for overflow.
+ *
+ * Note: DO NOT DECLARE THIS FUNCTION STATIC OR THE COMPILER MIGHT
+ * PROMOTE ARGUMENT `x' TO DOUBLE AND DEFEAT THE OVERFLOW
+ * CHECKING.
+ *
+ * Return: void
+ *
+ * Programmer: Robb Matzke
+ * Tuesday, July 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+void
+some_dummy_func(float x)
+{
+ char s[128];
+ sprintf(s, "%g", x);
+}
+
+
+/*-------------------------------------------------------------------------
+ * Function: generates_sigfpe
+ *
+ * Purpose: Determines if SIGFPE is generated from overflows. We must be
+ * able to fork() and waitpid() in order for this test to work
+ * properly. Sets skip_overflow_tests_g to non-zero if they
+ * would generate SIGBUS, zero otherwise.
+ *
+ * Programmer: Robb Matzke
+ * Tuesday, July 21, 1998
+ *
+ * Modifications:
+ *
+ *-------------------------------------------------------------------------
+ */
+static void
+generates_sigfpe(void)
+{
+#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID)
+ pid_t pid;
+ int status;
+ size_t i, j;
+ double d;
+ unsigned char *dp = (unsigned char*)&d;
+ float f;
+
+ HDfflush(stdout);
+ HDfflush(stderr);
+ if ((pid=fork())<0) {
+ HDperror("fork");
+ HDexit(1);
+ } else if (0==pid) {
+ for (i=0; i<2000; i++) {
+ for (j=0; j<sizeof(double); j++) dp[j] = HDrand();
+ f = (float)d;
+ some_dummy_func((float)f);
+ }
+ HDexit(0);
+ }
+
+ while (pid!=waitpid(pid, &status, 0))
+ /*void*/;
+ if (WIFEXITED(status) && 0==WEXITSTATUS(status)) {
+ HDputs("Floating-point overflow cases will be tested.");
+ skip_overflow_tests_g = FALSE;
+ } else if (WIFSIGNALED(status) && SIGFPE==WTERMSIG(status)) {
+ HDputs("Floating-point overflow cases cannot be safely tested.");
+ skip_overflow_tests_g = TRUE;
+ /* delete the core dump file that SIGFPE may have created */
+ HDunlink("core");
+ }
+#else
+ HDputs("Cannot determine if floating-point overflows generate a SIGFPE;");
+ HDputs("assuming yes.");
+ HDputs("Overflow cases will not be tested.");
+ skip_overflow_tests_g = TRUE;
+#endif
+}
+
+
+/*-------------------------------------------------------------------------
* Function: test_derived_flt
*
* Purpose: Tests user-define and query functions of floating-point types.
@@ -4383,6 +4484,9 @@ main(void)
* for user-defined integer types */
nerrors += test_derived_integer();
+ /* Does floating point overflow generate a SIGFPE? */
+ generates_sigfpe();
+
/* Test degenerate cases */
nerrors += run_fp_tests("noop");
diff --git a/test/dtypes.c b/test/dtypes.c
index f7d9c8f..ba18b91 100644
--- a/test/dtypes.c
+++ b/test/dtypes.c
@@ -32,9 +32,6 @@
/* For test_compound_10 */
#define ARRAY_DIM 4
-/* Define if you want to see a count of overflows */
-#undef SHOW_OVERFLOWS
-
/* Epsilon for floating-point comparisons */
#define FP_EPSILON 0.000001
@@ -89,276 +86,9 @@ typedef enum dtype_t {
FLT_LDOUBLE, OTHER
} dtype_t;
-/* Count the number of overflows */
-#ifdef SHOW_OVERFLOWS
-static int noverflows_g = 0;
-#endif
-
-/* Skip overflow tests if non-zero */
-static int skip_overflow_tests_g = 0;
-
-/* Don't use hardware conversions if set */
-static int without_hardware_g = 0;
-
/* Count opaque conversions */
static int num_opaque_conversions_g = 0;
-/*
- * Although we check whether a floating point overflow generates a SIGFPE and
- * turn off overflow tests in that case, it might still be possible for an
- * overflow condition to occur. Once a SIGFPE is raised the program cannot
- * be allowed to continue (cf. Posix signals) so in order to recover from a
- * SIGFPE we run tests that might generate one in a child process.
- */
-#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID)
-/*#define HANDLE_SIGFPE*/
-#endif
-
-/* Allocates memory aligned on a certain boundary. */
-#define aligned_malloc(Z) ((void*)((char*)HDmalloc(ALIGNMENT+Z)+ALIGNMENT))
-#define aligned_free(M) HDfree((char*)(M)-ALIGNMENT)
-
-/* Initialize source buffer of integer for integer->integer and integer->floating-point conversion test.
- * This algorithm is mainly to avoid any casting and comparison between source and destination types
- * for compiler, because we're testing conversions. */
-#define INIT_INTEGER(TYPE, SRC_MAX, SRC_MIN, SRC_SIZE, DST_SIZE, SRC_PREC, BUF, SAVED, NELMTS) \
-{ \
- unsigned char *buf_p, *saved_p; \
- int n; \
- TYPE value1 = 1; \
- TYPE value2 = 0; \
- \
- /* Allocate buffers */ \
- NELMTS=(SRC_PREC-1)*3+1; \
- BUF = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- SAVED = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- \
- buf_p = BUF; \
- saved_p = SAVED; \
- \
- /*positive values, ascending order. VALUE1 starts from 00000001, to 00000010, until 10000000*/ \
- /*VALUE2 ascends from 00000000, to 00000011, 00000111,..., until 11111111.*/ \
- for(n=0; n<SRC_PREC-1; n++) { \
- if(value1<=SRC_MAX && value1>=SRC_MIN) { \
- memcpy(buf_p, &value1, SRC_SIZE); \
- memcpy(saved_p, &value1, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- if(value2<=SRC_MAX && value2>=SRC_MIN) { \
- memcpy(buf_p, &value2, SRC_SIZE); \
- memcpy(saved_p, &value2, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- \
- value1 <<= 1; \
- value2 = (value1 - 1) | value1; \
- } \
- \
- /* negative values for signed; descending positive values for unsigned */ \
- /* VALUE2 descends from 11111111 to 11111110, 11111100, ..., until 10000000. */ \
- for(n=0; n<SRC_PREC; n++) { \
- if(value2<=SRC_MAX && value2>=SRC_MIN) { \
- memcpy(buf_p, &value2, SRC_SIZE); \
- memcpy(saved_p, &value2, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- value2 <<= 1; \
- } \
-}
-
-/* Change a buffer's byte order from big endian to little endian. It's mainly for library's
- * bit operations which handle only little endian order.
- */
-#define CHANGE_ORDER(EBUF, EORDER, ESIZE) \
-{ \
- if (H5T_ORDER_BE==EORDER) { \
- int m; \
- unsigned char mediator; \
- size_t half_size = ESIZE/2; \
- for (m=0; m<half_size; m++) { \
- mediator = EBUF[ESIZE-(m+1)]; \
- EBUF[ESIZE-(m+1)] = EBUF[m]; \
- EBUF[m] = mediator; \
- } \
- } \
-}
-
-/* Allocate buffer and initialize it with floating-point normalized and denormalized values.
- * It's for conversion test of floating-point as the source.
- */
-#define INIT_FP(TYPE, SRC_MAX, SRC_MIN, SRC_MAX_10_EXP, SRC_MIN_10_EXP, SRC_MANT_DIG, SRC_SIZE, \
- SRC_PREC, SRC_ORDR, DST_SIZE, BUF, SAVED, NELMTS) \
-{ \
- unsigned char *buf_p, *saved_p; \
- unsigned char *tmp1, *tmp2; \
- size_t num_norm, factor, n; \
- TYPE value1, value2; \
- TYPE multiply; \
- \
- /*Determine the number of normalized values and increment pace. The values start from \
- *minimal normalized value and are multiplied by MULTIPLY each step until reach to maximal \
- *normalized value.*/ \
- if(SRC_MAX_10_EXP<100) { /*for float*/ \
- factor = 0; \
- multiply = 10; \
- } else if(SRC_MAX_10_EXP>=100 && SRC_MAX_10_EXP<400) { /*for double*/ \
- factor = 2; \
- multiply = 10000; \
- } else { /*for long double*/ \
- factor = 3; \
- multiply = 100000000; \
- } \
- \
- /*The number of values if multiplied by 10 for each step.*/ \
- num_norm = (SRC_MAX_10_EXP - SRC_MIN_10_EXP); \
- /*Reduce the number of values by 2^factor. MULTIPLY=10^(2^factor). Using this algorithm \
- *instead of arithmatic operation to avoid any conversion*/ \
- num_norm >>= factor; \
- \
- /*Total number of values*/ \
- NELMTS = 2 * /*both positive and negative*/ \
- (num_norm + /*number of normalized values*/ \
- 1 + /*maximal normalized value*/ \
- SRC_MANT_DIG - 1); /*number of denormalized values*/ \
- \
- /* Allocate buffers */ \
- BUF = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- SAVED = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- tmp1 = (unsigned char*)malloc(SRC_SIZE); \
- tmp2 = (unsigned char*)malloc(SRC_SIZE); \
- \
- buf_p = BUF; \
- saved_p = SAVED; \
- \
- /*Normalized values*/ \
- value1 = SRC_MIN; \
- value2 = -SRC_MIN; \
- for(n=0; n<num_norm; n++) { \
- if(value1<SRC_MAX) { /*positive*/ \
- memcpy(buf_p, &value1, SRC_SIZE); \
- memcpy(saved_p, &value1, SRC_SIZE); \
- value1 *= multiply; \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- if(value2>-SRC_MAX) { /*negative*/ \
- memcpy(buf_p, &value2, SRC_SIZE); \
- memcpy(saved_p, &value2, SRC_SIZE); \
- value2 *= multiply; \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- } \
- \
- value1 = SRC_MAX; /*maximal value*/ \
- memcpy(buf_p, &value1, SRC_SIZE); \
- memcpy(saved_p, &value1, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- \
- value2 = -SRC_MAX; /*negative value*/ \
- memcpy(buf_p, &value2, SRC_SIZE); \
- memcpy(saved_p, &value2, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- \
- /*Denormalized values. Exponent is 0. Let mantissa starts from 00000001, 00000011, \
- *00000111,..., until 11111111.*/ \
- memset(tmp1, 0, SRC_SIZE); \
- memset(tmp2, 0, SRC_SIZE); \
- H5T_bit_set (tmp2, SRC_PREC-1, 1, TRUE); /*the negative value*/ \
- for(n=0; n<SRC_MANT_DIG-1; n++) { \
- H5T_bit_set (tmp1, n, 1, TRUE); /*turn on 1 bit each time*/ \
- CHANGE_ORDER(tmp1, SRC_ORDR, SRC_SIZE); /*change order for big endian*/ \
- memcpy(buf_p, tmp1, SRC_SIZE); \
- memcpy(saved_p, tmp1, SRC_SIZE); \
- CHANGE_ORDER(tmp1, SRC_ORDR, SRC_SIZE); /*change back the order for bit operation*/ \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- \
- /*negative values*/ \
- H5T_bit_set (tmp2, n, 1, TRUE); \
- CHANGE_ORDER(tmp2, SRC_ORDR, SRC_SIZE); \
- memcpy(buf_p, tmp2, SRC_SIZE); \
- memcpy(saved_p, tmp2, SRC_SIZE); \
- CHANGE_ORDER(tmp2, SRC_ORDR, SRC_SIZE); \
- buf_p += SRC_SIZE; \
- saved_p += SRC_SIZE; \
- } \
- free(tmp1); \
- free(tmp2); \
-}
-
-/* Allocate buffer and initialize it with floating-point special values, +/-0, +/-infinity,
- * +/-QNaN, +/-SNaN. It's for conversion test of floating-point as the source.
- */
-#define INIT_FP_SPECIAL(SRC_SIZE, SRC_PREC, SRC_ORDR, SRC_MANT_DIG, DST_SIZE, \
- BUF, SAVED, NELMTS) \
-{ \
- unsigned char *buf_p, *saved_p; \
- unsigned char *value; \
- int n; \
- \
- /*Total number of values*/ \
- NELMTS = 2 * /*both positive and negative*/ \
- 4; /*infinity, SNaN, QNaN */ \
- \
- /* Allocate buffers */ \
- BUF = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- SAVED = (unsigned char*)aligned_malloc( NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- value = (unsigned char*)calloc(SRC_SIZE, sizeof(unsigned char)); \
- \
- buf_p = BUF; \
- saved_p = SAVED; \
- \
- /* +0 */ \
- H5T_bit_set(value, 0, SRC_PREC, FALSE); \
- memcpy(buf_p, value, SRC_SIZE*sizeof(unsigned char)); \
- buf_p += SRC_SIZE; \
- \
- for(n=0; n<2; n++) { \
- if(n==1) { \
- memset(value, 0, SRC_SIZE*sizeof(unsigned char)); \
- /* -0 */ \
- H5T_bit_set(value, SRC_PREC-1, 1, TRUE); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE);/*change order for big endian*/ \
- memcpy(buf_p, value, SRC_SIZE*sizeof(unsigned char)); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE);/*change back the order for bit operation*/ \
- buf_p += SRC_SIZE; \
- } \
- \
- /* +/-infinity */ \
- H5T_bit_set(value, SRC_MANT_DIG-1, SRC_PREC-SRC_MANT_DIG, TRUE); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change order for big endian*/ \
- memcpy(buf_p, value, SRC_SIZE*sizeof(unsigned char)); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change back the order for bit operation*/ \
- buf_p += SRC_SIZE; \
- \
- /* +/-SNaN */ \
- H5T_bit_set(value, 0, 1, TRUE); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change order for big endian*/ \
- memcpy(buf_p, value, SRC_SIZE*sizeof(unsigned char)); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change back the order for bit operation*/ \
- buf_p += SRC_SIZE; \
- \
- /* +/-QNaN */ \
- H5T_bit_set(value, SRC_MANT_DIG-2, 1, TRUE); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change order for big endian*/ \
- memcpy(buf_p, value, SRC_SIZE*sizeof(unsigned char)); \
- CHANGE_ORDER(value, SRC_ORDR, SRC_SIZE); /*change back the order for bit operation*/ \
- buf_p += SRC_SIZE; \
- } \
- \
- memcpy(SAVED, BUF, NELMTS*MAX(SRC_SIZE, DST_SIZE)); \
- free(value); \
-}
-
-void some_dummy_func(float x);
-static hbool_t overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits);
-static int my_isnan(dtype_t type, void *val);
static int opaque_check(int tag_it);
static herr_t convert_opaque(hid_t UNUSED st, hid_t UNUSED dt,
H5T_cdata_t *cdata,
@@ -368,194 +98,6 @@ static herr_t convert_opaque(hid_t UNUSED st, hid_t UNUSED dt,
/*-------------------------------------------------------------------------
- * Function: fpe_handler
- *
- * Purpose: Exit with 255
- *
- * Return: void
- *
- * Programmer: Robb Matzke
- * Monday, July 6, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static void
-fpe_handler(int UNUSED signo)
-{
- SKIPPED();
- HDputs(" Test skipped due to SIGFPE.");
-#ifndef HANDLE_SIGFPE
- HDputs(" Remaining tests could not be run.");
- HDputs(" Please turn off SIGFPE on overflows and try again.");
-#endif
- HDexit(255);
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: overflow_handler
- *
- * Purpose: Gets called for all data type conversion overflows.
- *
- * Return: Success: 0
- *
- * Failure: -1
- *
- * Programmer: Robb Matzke
- * Tuesday, July 7, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-#ifdef SHOW_OVERFLOWS
-static herr_t
-overflow_handler(hid_t UNUSED src_id, hid_t UNUSED dst_id,
- void UNUSED *src_buf, void UNUSED *dst_buf)
-{
- noverflows_g++;
- return -1;
-}
-#endif
-
-
-/*-------------------------------------------------------------------------
- * Function: except_func
- *
- * Purpose: Gets called for all data type conversion exceptions.
- *
- * Return: H5T_CONV_ABORT: -1
- *
- * H5T_CONV_UNHANDLED 0
- *
- * H5T_CONV_HANDLED 1
- *
- * Programmer: Raymond Lu
- * April 19, 2004
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static H5T_conv_ret_t
-except_func(int except_type, hid_t UNUSED src_id, hid_t UNUSED dst_id, void UNUSED *src_buf,
- void *dst_buf, void *user_data)
-{
- H5T_conv_ret_t ret = H5T_CONV_HANDLED;
-
- if(except_type == H5T_CONV_EXCEPT_RANGE_HI)
- /*only test integer case*/
- *(int*)dst_buf = *(int*)user_data;
- else if(except_type == H5T_CONV_EXCEPT_RANGE_LOW)
- /*only test integer case*/
- *(int*)dst_buf = *(int*)user_data;
- else if(except_type == H5T_CONV_EXCEPT_TRUNCATE)
- ret = H5T_CONV_UNHANDLED;
- else if(except_type == H5T_CONV_EXCEPT_PRECISION)
- ret = H5T_CONV_UNHANDLED;
- else if(except_type == H5T_CONV_EXCEPT_PINF)
- /*only test integer case*/
- *(int*)dst_buf = *(int*)user_data;
- else if(except_type == H5T_CONV_EXCEPT_NINF)
- /*only test integer case*/
- *(int*)dst_buf = *(int*)user_data;
- else if(except_type == H5T_CONV_EXCEPT_NAN)
- /*only test integer case*/
- *(int*)dst_buf = *(int*)user_data;
-
- return ret;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: some_dummy_func
- *
- * Purpose: A dummy function to help check for overflow.
- *
- * Note: DO NOT DECLARE THIS FUNCTION STATIC OR THE COMPILER MIGHT
- * PROMOTE ARGUMENT `x' TO DOUBLE AND DEFEAT THE OVERFLOW
- * CHECKING.
- *
- * Return: void
- *
- * Programmer: Robb Matzke
- * Tuesday, July 21, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-void
-some_dummy_func(float x)
-{
- char s[128];
- sprintf(s, "%g", x);
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: generates_sigfpe
- *
- * Purpose: Determines if SIGFPE is generated from overflows. We must be
- * able to fork() and waitpid() in order for this test to work
- * properly. Sets skip_overflow_tests_g to non-zero if they
- * would generate SIGBUS, zero otherwise.
- *
- * Programmer: Robb Matzke
- * Tuesday, July 21, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static void
-generates_sigfpe(void)
-{
-#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID)
- pid_t pid;
- int status;
- size_t i, j;
- double d;
- unsigned char *dp = (unsigned char*)&d;
- float f;
-
- HDfflush(stdout);
- HDfflush(stderr);
- if ((pid=fork())<0) {
- HDperror("fork");
- HDexit(1);
- } else if (0==pid) {
- for (i=0; i<2000; i++) {
- for (j=0; j<sizeof(double); j++) dp[j] = HDrand();
- f = (float)d;
- some_dummy_func((float)f);
- }
- HDexit(0);
- }
-
- while (pid!=waitpid(pid, &status, 0))
- /*void*/;
- if (WIFEXITED(status) && 0==WEXITSTATUS(status)) {
- HDputs("Floating-point overflow cases will be tested.");
- skip_overflow_tests_g = FALSE;
- } else if (WIFSIGNALED(status) && SIGFPE==WTERMSIG(status)) {
- HDputs("Floating-point overflow cases cannot be safely tested.");
- skip_overflow_tests_g = TRUE;
- /* delete the core dump file that SIGFPE may have created */
- HDunlink("core");
- }
-#else
- HDputs("Cannot determine if floating-point overflows generate a SIGFPE;");
- HDputs("assuming yes.");
- HDputs("Overflow cases will not be tested.");
- skip_overflow_tests_g = TRUE;
-#endif
-}
-
-
-/*-------------------------------------------------------------------------
* Function: reset_hdf5
*
* Purpose: Reset the hdf5 library. This causes statistics to be printed
@@ -574,10 +116,6 @@ static void
reset_hdf5(void)
{
h5_reset();
-#ifdef SHOW_OVERFLOWS
- H5Tset_overflow(overflow_handler);
-#endif
- if (without_hardware_g) h5_no_hwconv();
#ifdef TEST_ALIGNMENT
SET_ALIGNMENT(SCHAR, H5_SIZEOF_CHAR);
SET_ALIGNMENT(UCHAR, H5_SIZEOF_CHAR);
@@ -626,7 +164,7 @@ test_classes(void)
hid_t memb_id; /* Compound member datatype */
H5T_class_t memb_cls;
H5T_class_t tcls;
- int nmembs, i;
+ unsigned int nmembs, i;
TESTING("H5Tget_class()");
@@ -2311,19 +1849,23 @@ test_compound_11(void)
/* Verify converted buffer is correct */
for(u=0; u<NTESTELEM; u++) {
if(((big_t *)buf_orig)[u].d1!=((little_t *)buf)[u].d1) {
- printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
+ printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].i1!=((little_t *)buf)[u].i1) {
- printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
+ printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].s1==NULL || ((little_t *)buf)[u].s1==NULL) {
- printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
else if(HDstrcmp(((big_t *)buf_orig)[u].s1,((little_t *)buf)[u].s1)) {
- printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
HDfree(((little_t *)buf)[u].s1);
@@ -2347,19 +1889,23 @@ test_compound_11(void)
/* Verify converted buffer is correct */
for(u=0; u<NTESTELEM; u++) {
if(((big_t *)buf_orig)[u].d1!=((little_t *)buf)[u].d1) {
- printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
+ printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].i1!=((little_t *)buf)[u].i1) {
- printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
+ printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].s1==NULL || ((little_t *)buf)[u].s1==NULL) {
- printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
else if(HDstrcmp(((big_t *)buf_orig)[u].s1,((little_t *)buf)[u].s1)) {
- printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
HDfree(((little_t *)buf)[u].s1);
@@ -2377,19 +1923,23 @@ test_compound_11(void)
/* Verify converted buffer is correct */
for(u=0; u<NTESTELEM; u++) {
if(((big_t *)buf_orig)[u].d1!=((little_t *)buf)[u].d1) {
- printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
+ printf("Error, line #%d: buf_orig[%u].d1=%f, buf[%u].d1=%f\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].d1,(unsigned)u,((little_t *)buf)[u].d1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].i1!=((little_t *)buf)[u].i1) {
- printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
+ printf("Error, line #%d: buf_orig[%u].i1=%d, buf[%u].i1=%d\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].i1,(unsigned)u,((little_t *)buf)[u].i1);
TEST_ERROR
} /* end if */
if(((big_t *)buf_orig)[u].s1==NULL || ((little_t *)buf)[u].s1==NULL) {
- printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%p, buf[%u].s1=%p\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
else if(HDstrcmp(((big_t *)buf_orig)[u].s1,((little_t *)buf)[u].s1)) {
- printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,(unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
+ printf("Error, line #%d: buf_orig[%u].s1=%s, buf[%u].s1=%s\n",__LINE__,
+ (unsigned)u,((big_t *)buf_orig)[u].s1,(unsigned)u,((little_t *)buf)[u].s1);
TEST_ERROR
} /* end if */
HDfree(((little_t *)buf)[u].s1);
@@ -2753,753 +2303,6 @@ test_query(void)
/*-------------------------------------------------------------------------
- * Function: test_derived_flt
- *
- * Purpose: Tests user-define and query functions of floating-point types.
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Raymond Lu
- * Thursday, Jan 6, 2005
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_derived_flt(void)
-{
- hid_t file=-1, tid1=-1, tid2=-1;
- hid_t dxpl_id=-1;
- char filename[1024];
- size_t precision, spos, epos, esize, mpos, msize, size, ebias;
- int offset;
- size_t src_size, dst_size;
- unsigned char *buf=NULL, *saved_buf=NULL;
- int *aligned=NULL;
- int endian; /*endianess */
- size_t nelmts = NTESTELEM;
- unsigned int fails_this_test = 0;
- const size_t max_fails=40; /*max number of failures*/
- char str[256]; /*message string */
- unsigned int i, j;
-
- TESTING("user-define and query functions of floating-point types");
-
- /* Create File */
- h5_fixname(FILENAME[6], H5P_DEFAULT, filename, sizeof filename);
- if((file=H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT))<0) {
- H5_FAILED();
- printf("Can't create file\n");
- goto error;
- }
-
- if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0) {
- H5_FAILED();
- printf("Can't create data transfer property list\n");
- goto error;
- }
-
- if((tid1 = H5Tcopy(H5T_IEEE_F64LE))<0) {
- H5_FAILED();
- printf("Can't copy data type\n");
- goto error;
- }
-
- if((tid2 = H5Tcopy(H5T_IEEE_F32LE))<0) {
- H5_FAILED();
- printf("Can't copy data type\n");
- goto error;
- }
-
- /*------------------------------------------------------------------------
- * 1st floating-point type
- * size=7 byte, precision=42 bits, offset=3 bits, mantissa size=31 bits,
- * mantissa position=3, exponent size=10 bits, exponent position=34,
- * exponent bias=511. It can be illustrated in little-endian order as
- *
- * 6 5 4 3 2 1 0
- * ???????? ???SEEEE EEEEEEMM MMMMMMMM MMMMMMMM MMMMMMMM MMMMM???
- *
- * To create a new floating-point type, the following properties must be
- * set in the order of
- * set fields -> set offset -> set precision -> set size.
- * All these properties must be set before the type can function. Other
- * properties can be set anytime. Derived type size cannot be expanded
- * bigger than original size but can be decreased. There should be no
- * holes among the significant bits. Exponent bias usually is set
- * 2^(n-1)-1, where n is the exponent size.
- *-----------------------------------------------------------------------*/
- if(H5Tset_fields(tid1, 44, 34, 10, 3, 31)<0) {
- H5_FAILED();
- printf("Can't set fields\n");
- goto error;
- }
- if(H5Tset_offset(tid1, 3)<0) {
- H5_FAILED();
- printf("Can't set offset\n");
- goto error;
- }
- if(H5Tset_precision(tid1, 42)<0) {
- H5_FAILED();
- printf("Can't set precision 1\n");
- goto error;
- }
- if(H5Tset_size(tid1, 7)<0) {
- H5_FAILED();
- printf("Can't set size\n");
- goto error;
- }
-
- if(H5Tset_ebias(tid1, 511)<0) {
- H5_FAILED();
- printf("Can't set exponent bias\n");
- goto error;
- }
- if(H5Tset_pad(tid1, H5T_PAD_ZERO, H5T_PAD_ZERO)<0) {
- H5_FAILED();
- printf("Can't set padding\n");
- goto error;
- }
-
- if(H5Tcommit(file, "new float type 1", tid1)<0) {
- H5_FAILED();
- printf("Can't set inpad\n");
- goto error;
- }
- if(H5Tclose(tid1)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if((tid1 = H5Topen(file, "new float type 1"))<0) {
- H5_FAILED();
- printf("Can't open datatype\n");
- goto error;
- }
- if(H5Tget_fields(tid1, &spos, &epos, &esize, &mpos, &msize)<0) {
- H5_FAILED();
- printf("Can't get fields\n");
- goto error;
- }
- if(spos!=44 || epos!=34 || esize!=10 || mpos!=3 || msize!=31) {
- H5_FAILED();
- printf("Wrong field values\n");
- goto error;
- }
-
- if((precision = H5Tget_precision(tid1))!=42) {
- H5_FAILED();
- printf("Can't get precision or wrong precision\n");
- goto error;
- }
- if((offset = H5Tget_offset(tid1))!=3) {
- H5_FAILED();
- printf("Can't get offset or wrong offset\n");
- goto error;
- }
- if((size = H5Tget_size(tid1))!=7) {
- H5_FAILED();
- printf("Can't get size or wrong size\n");
- goto error;
- }
- if((ebias = H5Tget_ebias(tid1))!=511) {
- H5_FAILED();
- printf("Can't get exponent bias or wrong bias\n");
- goto error;
- }
-
- /* Convert data from native integer to the 1st derived floating-point type.
- * Then convert data from the floating-point type back to native integer.
- * Compare the final data with the original data.
- */
- src_size = H5Tget_size(H5T_NATIVE_INT);
- endian = H5Tget_order(H5T_NATIVE_INT);
- buf = (unsigned char*)malloc(nelmts*(MAX(src_size, size)));
- saved_buf = (unsigned char*)malloc(nelmts*src_size);
- aligned = (int*)malloc(src_size);
-
- for(i=0; i<nelmts*src_size; i++)
- buf[i] = saved_buf[i] = HDrand();
-
- /* Convert data from native integer to derived floating-point type.
- * The mantissa is big enough to retain the integer's precision. */
- if(H5Tconvert(H5T_NATIVE_INT, tid1, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
- /* Convert data from the derived floating-point type back to native integer. */
- if(H5Tconvert(tid1, H5T_NATIVE_INT, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
-
- /* Are the values still the same?*/
- for(i=0; i<nelmts; i++) {
- for(j=0; j<src_size; j++)
- if(buf[i*src_size+j]!=saved_buf[i*src_size+j])
- break;
- if(j==src_size)
- continue; /*no error*/
-
- /* Print errors */
- if (0==fails_this_test++) {
- sprintf(str, "\nTesting random sw int -> derived floating-point conversions");
- printf("%-70s", str);
- HDfflush(stdout);
- H5_FAILED();
- }
- printf(" test %u elmt %u: \n", 1, (unsigned)i);
-
- printf(" src = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", saved_buf[i*src_size+ENDIAN(src_size, j)]);
-
- HDmemcpy(aligned, saved_buf+i*sizeof(int), sizeof(int));
- printf(" %29d\n", *aligned);
-
- printf(" dst = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", buf[i*src_size+ENDIAN(src_size, j)]);
-
- HDmemcpy(aligned, buf+i*sizeof(int), sizeof(int));
- printf(" %29d\n", *aligned);
-
- if (fails_this_test>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto error;
- }
- }
-
- fails_this_test = 0;
- if(buf) free(buf);
- if(saved_buf) free(saved_buf);
- if(aligned) free(aligned);
- buf = NULL;
- saved_buf = NULL;
- aligned = NULL;
-
- /*--------------------------------------------------------------------------
- * 2nd floating-point type
- * size=3 byte, precision=24 bits, offset=0 bits, mantissa size=16 bits,
- * mantissa position=0, exponent size=7 bits, exponent position=16, exponent
- * bias=63. It can be illustrated in little-endian order as
- *
- * 2 1 0
- * SEEEEEEE MMMMMMMM MMMMMMMM
- *--------------------------------------------------------------------------*/
- if(H5Tset_fields(tid2, 23, 16, 7, 0, 16)<0) {
- H5_FAILED();
- printf("Can't set fields\n");
- goto error;
- }
- if(H5Tset_offset(tid2, 0)<0) {
- H5_FAILED();
- printf("Can't set offset\n");
- goto error;
- }
- if(H5Tset_precision(tid2, 24)<0) {
- H5_FAILED();
- printf("Can't set precision 2\n");
- goto error;
- }
- if(H5Tset_size(tid2, 3)<0) {
- H5_FAILED();
- printf("Can't set size\n");
- goto error;
- }
- if(H5Tset_ebias(tid2, 63)<0) {
- H5_FAILED();
- printf("Can't set size\n");
- goto error;
- }
- if(H5Tset_pad(tid2, H5T_PAD_ZERO, H5T_PAD_ZERO)<0) {
- H5_FAILED();
- printf("Can't set padding\n");
- goto error;
- }
-
- if(H5Tcommit(file, "new float type 2", tid2)<0) {
- H5_FAILED();
- printf("Can't set inpad\n");
- goto error;
- }
- if(H5Tclose(tid2)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if((tid2 = H5Topen(file, "new float type 2"))<0) {
- H5_FAILED();
- printf("Can't open datatype\n");
- goto error;
- }
- if(H5Tget_fields(tid2, &spos, &epos, &esize, &mpos, &msize)<0) {
- H5_FAILED();
- printf("Can't get fields\n");
- goto error;
- }
- if(spos!=23 || epos!=16 || esize!=7 || mpos!=0 || msize!=16) {
- H5_FAILED();
- printf("Wrong field values\n");
- goto error;
- }
-
- if((precision = H5Tget_precision(tid2))!=24) {
- H5_FAILED();
- printf("Can't get precision or wrong precision\n");
- goto error;
- }
- if((offset = H5Tget_offset(tid2))!=0) {
- H5_FAILED();
- printf("Can't get offset or wrong offset\n");
- goto error;
- }
- if((size = H5Tget_size(tid2))!=3) {
- H5_FAILED();
- printf("Can't get size or wrong size\n");
- goto error;
- }
- if((ebias = H5Tget_ebias(tid2))!=63) {
- H5_FAILED();
- printf("Can't get exponent bias or wrong bias\n");
- goto error;
- }
-
- /* Convert data from the 2nd to the 1st derived floating-point type.
- * Then convert data from the 1st type back to the 2nd type.
- * Compare the final data with the original data.
- */
- src_size = H5Tget_size(tid2);
- dst_size = H5Tget_size(tid1);
- endian = H5Tget_order(tid2);
- buf = (unsigned char*)malloc(nelmts*(MAX(src_size, dst_size)));
- saved_buf = (unsigned char*)malloc(nelmts*src_size);
-
- for(i=0; i<nelmts*src_size; i++)
- buf[i] = saved_buf[i] = HDrand();
-
- /* Convert data from the 2nd to the 1st derived floating-point type.
- * The mantissa and exponent of the 2nd type are big enough to retain
- * the precision and exponent power. */
- if(H5Tconvert(tid2, tid1, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
- /* Convert data from the 1st back to the 2nd derived floating-point type. */
- if(H5Tconvert(tid1, tid2, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
-
- /* Are the values still the same?*/
- for(i=0; i<nelmts; i++) {
- for(j=0; j<src_size; j++)
- if(buf[i*src_size+j]!=saved_buf[i*src_size+j])
- break;
- if(j==src_size)
- continue; /*no error*/
-
- /* If original value is NaN(exponent bits are all ones, 11..11),
- * the library simply sets all mantissa bits to ones. So don't
- * compare values in this case.
- */
- if((buf[i*src_size+2]==0x7f && saved_buf[i*src_size+2]==0x7f) ||
- (buf[i*src_size+2]==0xff && saved_buf[i*src_size+2]==0xff))
- continue;
-
- /* Print errors */
- if (0==fails_this_test++) {
- sprintf(str, "\nTesting random sw derived floating-point -> derived floating-point conversions");
- printf("%-70s", str);
- HDfflush(stdout);
- H5_FAILED();
- }
- printf(" test %u elmt %u: \n", 1, (unsigned)i);
-
- printf(" src = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", saved_buf[i*src_size+ENDIAN(src_size, j)]);
- printf("\n");
-
- printf(" dst = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", buf[i*src_size+ENDIAN(src_size, j)]);
- printf("\n");
-
- if (fails_this_test>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto error;
- }
- }
-
- if (buf) free(buf);
- if (saved_buf) free(saved_buf);
-
- if(H5Tclose(tid1)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if(H5Tclose(tid2)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if(H5Pclose(dxpl_id)<0) {
- H5_FAILED();
- printf("Can't close property list\n");
- goto error;
- }
-
- if(H5Fclose(file)<0) {
- H5_FAILED();
- printf("Can't close file\n");
- goto error;
- } /* end if */
-
- PASSED();
- reset_hdf5(); /*print statistics*/
-
- return 0;
-
- error:
- if (buf) free(buf);
- if (saved_buf) free(saved_buf);
- if (aligned) free(aligned);
- HDfflush(stdout);
- H5E_BEGIN_TRY {
- H5Tclose (tid1);
- H5Tclose (tid2);
- H5Pclose (dxpl_id);
- H5Fclose (file);
- } H5E_END_TRY;
- reset_hdf5(); /*print statistics*/
- return MAX((int)fails_this_test, 1);
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: test_derived_integer
- *
- * Purpose: Tests user-define and query functions of integer types.
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Raymond Lu
- * Saturday, Jan 29, 2005
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_derived_integer(void)
-{
- hid_t file=-1, tid1=-1, tid2=-1;
- hid_t dxpl_id=-1;
- char filename[1024];
- size_t precision, size;
- int offset;
- H5T_order_t order;
- H5T_sign_t sign;
- size_t src_size, dst_size;
- unsigned char *buf=NULL, *saved_buf=NULL;
- int *aligned=NULL;
- int endian; /*endianess */
- size_t nelmts = NTESTELEM;
- unsigned int fails_this_test = 0;
- const size_t max_fails=40; /*max number of failures*/
- char str[256]; /*message string */
- unsigned int i, j;
-
- TESTING("user-define and query functions of integer types");
-
- /* Create File */
- h5_fixname(FILENAME[7], H5P_DEFAULT, filename, sizeof filename);
- if((file=H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT))<0) {
- H5_FAILED();
- printf("Can't create file\n");
- goto error;
- }
-
- if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0) {
- H5_FAILED();
- printf("Can't create data transfer property list\n");
- goto error;
- }
-
- if((tid1 = H5Tcopy(H5T_STD_I32LE))<0) {
- H5_FAILED();
- printf("Can't copy data type\n");
- goto error;
- }
-
- if((tid2 = H5Tcopy(H5T_STD_U64LE))<0) {
- H5_FAILED();
- printf("Can't copy data type\n");
- goto error;
- }
-
- /*--------------------------------------------------------------------------
- * 1st integer type
- * size=3 byte, precision=24 bits, offset=0 bits, order=big endian.
- * It can be illustrated in big-endian order as
- *
- * 0 1 2
- * SIIIIIII IIIIIIII IIIIIIII
- *
- * There's no specific order for these functions to define the attributes
- * of a new integer type, H5Tset_precision, H5Tset_offset, H5Tset_size,
- * H5Tset_order, H5Tset_pad, H5Tset_sign.
- *--------------------------------------------------------------------------*/
- if(H5Tset_offset(tid1,0)<0) {
- H5_FAILED();
- printf("Can't set offset\n");
- goto error;
- }
-
- if(H5Tset_size(tid1, 3)<0) {
- H5_FAILED();
- printf("Can't set size\n");
- goto error;
- }
-
- if(H5Tset_precision(tid1,24)<0) {
- H5_FAILED();
- printf("Can't set precision\n");
- goto error;
- }
-
- if(H5Tset_order(tid1, H5T_ORDER_BE)<0) {
- H5_FAILED();
- printf("Can't set order\n");
- goto error;
- }
-
- if(H5Tcommit(file, "new integer type 1", tid1)<0) {
- H5_FAILED();
- printf("Can't commit data type\n");
- goto error;
- }
-
- if(H5Tclose(tid1)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if((tid1 = H5Topen(file, "new integer type 1"))<0) {
- H5_FAILED();
- printf("Can't open datatype\n");
- goto error;
- }
-
- if((precision = H5Tget_precision(tid1))!=24) {
- H5_FAILED();
- printf("Can't get precision or wrong precision\n");
- goto error;
- }
- if((offset = H5Tget_offset(tid1))!=0) {
- H5_FAILED();
- printf("Can't get offset or wrong offset\n");
- goto error;
- }
- if((size = H5Tget_size(tid1))!=3) {
- H5_FAILED();
- printf("Can't get size or wrong size\n");
- goto error;
- }
- if((order = H5Tget_order(tid1))!=H5T_ORDER_BE) {
- H5_FAILED();
- printf("Can't get order or wrong order\n");
- goto error;
- }
-
- /*--------------------------------------------------------------------------
- * 2nd integer type
- * size=8 byte, precision=48 bits, offset=10 bits, order=little endian.
- * It can be illustrated in little-endian order as
- *
- * 7 6 5 4 3 2 1 0
- * ??????SI IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIII?? ????????
- *--------------------------------------------------------------------------*/
- if(H5Tset_precision(tid2,48)<0) {
- H5_FAILED();
- printf("Can't set precision\n");
- goto error;
- }
-
- if(H5Tset_offset(tid2,10)<0) {
- H5_FAILED();
- printf("Can't set offset\n");
- goto error;
- }
-
- if(H5Tset_sign(tid2,H5T_SGN_2)<0) {
- H5_FAILED();
- printf("Can't set offset\n");
- goto error;
- }
-
- if(H5Tcommit(file, "new integer type 2", tid2)<0) {
- H5_FAILED();
- printf("Can't commit data type\n");
- goto error;
- }
-
- if(H5Tclose(tid2)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if((tid2 = H5Topen(file, "new integer type 2"))<0) {
- H5_FAILED();
- printf("Can't open datatype\n");
- goto error;
- }
-
- if((precision = H5Tget_precision(tid2))!=48) {
- H5_FAILED();
- printf("Can't get precision or wrong precision\n");
- goto error;
- }
- if((offset = H5Tget_offset(tid2))!=10) {
- H5_FAILED();
- printf("Can't get offset or wrong offset\n");
- goto error;
- }
- if((size = H5Tget_size(tid2))!=8) {
- H5_FAILED();
- printf("Can't get size or wrong size\n");
- goto error;
- }
- if((sign = H5Tget_sign(tid2))!=H5T_SGN_2) {
- H5_FAILED();
- printf("Can't get sign or wrong sign\n");
- goto error;
- }
-
- /* Convert data from the 1st to the 2nd derived integer type.
- * Then convert data from the 2nd type back to the 1st type.
- * Compare the final data with the original data.
- */
- src_size = H5Tget_size(tid1);
- dst_size = H5Tget_size(tid2);
- endian = H5Tget_order(tid1);
- buf = (unsigned char*)malloc(nelmts*(MAX(src_size, dst_size)));
- saved_buf = (unsigned char*)malloc(nelmts*src_size);
-
- for(i=0; i<nelmts*src_size; i++)
- buf[i] = saved_buf[i] = HDrand();
-
- /* Convert data from the 1st to the 2nd derived integer type.
- * The precision of the 2nd type are big enough to retain
- * the 1st type's precision. */
- if(H5Tconvert(tid1, tid2, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
- /* Convert data from the 2nd back to the 1st derived integer type. */
- if(H5Tconvert(tid2, tid1, nelmts, buf, NULL, dxpl_id)<0) {
- H5_FAILED();
- printf("Can't convert data\n");
- goto error;
- }
-
- /* Are the values still the same?*/
- for(i=0; i<nelmts; i++) {
- for(j=0; j<src_size; j++)
- if(buf[i*src_size+j]!=saved_buf[i*src_size+j])
- break;
- if(j==src_size)
- continue; /*no error*/
-
- /* Print errors */
- if (0==fails_this_test++) {
- sprintf(str, "\nTesting random sw derived integer -> derived integer conversions");
- printf("%-70s", str);
- HDfflush(stdout);
- H5_FAILED();
- }
- printf(" test %u elmt %u: \n", 1, (unsigned)i);
-
- printf(" src = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", saved_buf[i*src_size+ENDIAN(src_size, j)]);
- printf("\n");
-
- printf(" dst = ");
- for (j=0; j<src_size; j++)
- printf(" %02x", buf[i*src_size+ENDIAN(src_size, j)]);
- printf("\n");
-
- if (fails_this_test>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto error;
- }
- }
-
- if(H5Tclose(tid1)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if(H5Tclose(tid2)<0) {
- H5_FAILED();
- printf("Can't close datatype\n");
- goto error;
- }
-
- if(H5Pclose(dxpl_id)<0) {
- H5_FAILED();
- printf("Can't close property list\n");
- goto error;
- }
-
- if(H5Fclose(file)<0) {
- H5_FAILED();
- printf("Can't close file\n");
- goto error;
- } /* end if */
-
- PASSED();
- reset_hdf5(); /*print statistics*/
-
- return 0;
-
- error:
- if (buf) free(buf);
- if (saved_buf) free(saved_buf);
- if (aligned) free(aligned);
- HDfflush(stdout);
- H5E_BEGIN_TRY {
- H5Tclose (tid1);
- H5Tclose (tid2);
- H5Pclose (dxpl_id);
- H5Fclose (file);
- } H5E_END_TRY;
- reset_hdf5(); /*print statistics*/
- return MAX((int)fails_this_test, 1);
-}
-
-
-/*-------------------------------------------------------------------------
* Function: test_transient
*
* Purpose: Tests transient data types.
@@ -4555,3377 +3358,6 @@ opaque_check(int tag_it)
/*-------------------------------------------------------------------------
- * Function: test_conv_int
- *
- * Purpose: Test atomic number conversions.
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Robb Matzke
- * Wednesday, June 10, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_conv_int (void)
-{
- unsigned char byte[4];
-
- /*---------------------------------------------------------------------
- * Test some specific overflow/underflow cases.
- *---------------------------------------------------------------------
- */
- TESTING("integer overflow conversions");
-
- /* (unsigned)0x80000000 -> (unsigned)0xffff */
- byte[0] = byte[1] = byte[2] = 0;
- byte[3] = 0x80;
- if (H5Tconvert (H5T_STD_U32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0xff || byte[1]!=0xff) {
- H5_FAILED();
- printf(" src: 0x80000000 unsigned\n");
- printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]);
- printf(" ans: 0xffff unsigned\n");
- goto error;
- }
-
- /* (unsigned)0xffffffff -> (signed)0x7fff */
- byte[0] = byte[1] = byte[2] = byte[3] = 0xff;
- if (H5Tconvert (H5T_STD_U32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0xff || byte[1]!=0x7f) {
- H5_FAILED();
- printf(" src: 0xffffffff unsigned\n");
- printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]);
- printf(" ans: 0x7fff signed\n");
- goto error;
- }
-
- /* (signed)0xffffffff -> (unsigned)0x0000 */
- byte[0] = byte[1] = byte[2] = byte[3] = 0xff;
- if (H5Tconvert (H5T_STD_I32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0x00 || byte[1]!=0x00) {
- H5_FAILED();
- printf(" src: 0xffffffff signed\n");
- printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]);
- printf(" ans: 0x0000 unsigned\n");
- goto error;
- }
-
- /* (signed)0x7fffffff -> (unsigned)0xffff */
- byte[0] = byte[1] = byte[2] = 0xff;
- byte[3] = 0x7f;
- if (H5Tconvert (H5T_STD_I32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0xff || byte[1]!=0xff) {
- H5_FAILED();
- printf(" src: 0x7fffffff signed\n");
- printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]);
- printf(" ans: 0xffff unsigned\n");
- goto error;
- }
-
- /* (signed)0x7fffffff -> (signed)0x7fff */
- byte[0] = byte[1] = byte[2] = 0xff;
- byte[3] = 0x7f;
- if (H5Tconvert (H5T_STD_I32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0xff || byte[1]!=0x7f) {
- H5_FAILED();
- printf(" src: 0x7fffffff signed\n");
- printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]);
- printf(" ans: 0x7fff signed\n");
- goto error;
- }
-
- /* (signed)0xbfffffff -> (signed)0x8000 */
- byte[0] = byte[1] = byte[2] = 0xff;
- byte[3] = 0xbf;
- if (H5Tconvert (H5T_STD_I32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) {
- goto error;
- }
- if (byte[0]!=0x00 || byte[1]!=0x80) {
- H5_FAILED();
- printf(" src: 0xbfffffff signed\n");
- printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]);
- printf(" ans: 0x8000 signed\n");
- goto error;
- }
-
- PASSED();
- reset_hdf5();
- return 0;
-
- error:
- reset_hdf5();
- return 1;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: test_conv_int_1
- *
- * Purpose: Test conversion of random integer values from SRC to DST.
- * These types should be any combination of:
- *
- * H5T_NATIVE_SCHAR H5T_NATIVE_UCHAR
- * H5T_NATIVE_SHORT H5T_NATIVE_USHORT
- * H5T_NATIVE_INT H5T_NATIVE_UINT
- * H5T_NATIVE_LONG H5T_NATIVE_ULONG
- * H5T_NATIVE_LLONG H5T_NATIVE_ULLONG
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Robb Matzke
- * Monday, November 16, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_conv_int_1(const char *name, hid_t src, hid_t dst)
-{
- const size_t ntests=NTESTS; /*number of tests */
- size_t nelmts=0; /*num values per test */
- const size_t max_fails=8; /*max number of failures*/
- size_t fails_all_tests=0; /*number of failures */
- size_t fails_this_test; /*fails for this test */
- char str[256]; /*hello string */
- dtype_t src_type, dst_type; /*data types */
- const char *src_type_name=NULL; /*source type name */
- const char *dst_type_name=NULL; /*destination type name */
- int endian; /*machine endianess */
- size_t src_size, dst_size; /*type sizes */
- unsigned char *buf=NULL; /*buffer for conversion */
- unsigned char *saved=NULL; /*original values */
- size_t i, j, k; /*counters */
- unsigned char *hw=NULL; /*hardware conv result */
- unsigned char src_bits[32]; /*src value in LE order */
- unsigned char dst_bits[32]; /*dest value in LE order*/
- size_t src_nbits; /*source length in bits */
- size_t dst_nbits; /*dst length in bits */
- H5T_sign_t src_sign; /*source sign type */
- H5T_sign_t dst_sign; /*dst sign type */
- void *aligned=NULL; /*aligned temp buffer */
- signed char hw_char;
- unsigned char hw_uchar;
- short hw_short;
- unsigned short hw_ushort;
- int hw_int;
- unsigned hw_uint;
- long hw_long;
- unsigned long hw_ulong;
- long_long hw_llong;
- unsigned long_long hw_ullong;
-
- /* What are the names of the source and destination types */
- if (H5Tequal(src, H5T_NATIVE_SCHAR)) {
- src_type_name = "signed char";
- src_type = INT_SCHAR;
- } else if (H5Tequal(src, H5T_NATIVE_UCHAR)) {
- src_type_name = "unsigned char";
- src_type = INT_UCHAR;
- } else if (H5Tequal(src, H5T_NATIVE_SHORT)) {
- src_type_name = "short";
- src_type = INT_SHORT;
- } else if (H5Tequal(src, H5T_NATIVE_USHORT)) {
- src_type_name = "unsigned short";
- src_type = INT_USHORT;
- } else if (H5Tequal(src, H5T_NATIVE_INT)) {
- src_type_name = "int";
- src_type = INT_INT;
- } else if (H5Tequal(src, H5T_NATIVE_UINT)) {
- src_type_name = "unsigned int";
- src_type = INT_UINT;
- } else if (H5Tequal(src, H5T_NATIVE_LONG)) {
- src_type_name = "long";
- src_type = INT_LONG;
- } else if (H5Tequal(src, H5T_NATIVE_ULONG)) {
- src_type_name = "unsigned long";
- src_type = INT_ULONG;
- } else if (H5Tequal(src, H5T_NATIVE_LLONG)) {
- src_type_name = "long long";
- src_type = INT_LLONG;
- } else if (H5Tequal(src, H5T_NATIVE_ULLONG)) {
- src_type_name = "unsigned long long";
- src_type = INT_ULLONG;
- } else {
- src_type_name = "UNKNOWN";
- src_type = OTHER;
- }
-
- if (H5Tequal(dst, H5T_NATIVE_SCHAR)) {
- dst_type_name = "signed char";
- dst_type = INT_SCHAR;
- } else if (H5Tequal(dst, H5T_NATIVE_UCHAR)) {
- dst_type_name = "unsigned char";
- dst_type = INT_UCHAR;
- } else if (H5Tequal(dst, H5T_NATIVE_SHORT)) {
- dst_type_name = "short";
- dst_type = INT_SHORT;
- } else if (H5Tequal(dst, H5T_NATIVE_USHORT)) {
- dst_type_name = "unsigned short";
- dst_type = INT_USHORT;
- } else if (H5Tequal(dst, H5T_NATIVE_INT)) {
- dst_type_name = "int";
- dst_type = INT_INT;
- } else if (H5Tequal(dst, H5T_NATIVE_UINT)) {
- dst_type_name = "unsigned int";
- dst_type = INT_UINT;
- } else if (H5Tequal(dst, H5T_NATIVE_LONG)) {
- dst_type_name = "long";
- dst_type = INT_LONG;
- } else if (H5Tequal(dst, H5T_NATIVE_ULONG)) {
- dst_type_name = "unsigned long";
- dst_type = INT_ULONG;
- } else if (H5Tequal(dst, H5T_NATIVE_LLONG)) {
- dst_type_name = "long long";
- dst_type = INT_LLONG;
- } else if (H5Tequal(dst, H5T_NATIVE_ULLONG)) {
- dst_type_name = "unsigned long long";
- dst_type = INT_ULLONG;
- } else {
- dst_type_name = "UNKNOWN";
- dst_type = OTHER;
- }
-
- /* Sanity checks */
- if (OTHER==src_type || OTHER==dst_type) {
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- printf("%-70s", str);
- H5_FAILED();
- HDputs(" Unknown data type.");
- goto error;
- }
-
-
- /* Some information about datatypes */
- endian = H5Tget_order(H5T_NATIVE_INT);
- src_size = H5Tget_size(src);
- dst_size = H5Tget_size(dst);
- src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */
- dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */
- src_sign = H5Tget_sign(src);
- dst_sign = H5Tget_sign(dst);
- aligned = HDmalloc(sizeof(long_long));
-#ifdef SHOW_OVERFLOWS
- noverflows_g = 0;
-#endif
-
- /* Allocate and initialize the source buffer through macro INIT_INTEGER. The BUF
- * will be used for the conversion while the SAVED buffer will be
- * used for the comparison later.
- */
- if(src_type == INT_SCHAR) {
- INIT_INTEGER(signed char, SCHAR_MAX, SCHAR_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_UCHAR) {
- INIT_INTEGER(unsigned char, UCHAR_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_SHORT) {
- INIT_INTEGER(short, SHRT_MAX, SHRT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_USHORT) {
- INIT_INTEGER(unsigned short, USHRT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_INT) {
- INIT_INTEGER(int, INT_MAX, INT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_UINT) {
- INIT_INTEGER(unsigned int, UINT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_LONG) {
- INIT_INTEGER(long, LONG_MAX, LONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_ULONG) {
- INIT_INTEGER(unsigned long, ULONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_LLONG) {
- INIT_INTEGER(long_long, LLONG_MAX, LLONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_ULLONG) {
- INIT_INTEGER(unsigned long_long, ULLONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else
- goto error;
-
- /* The tests */
- for (i=0; i<ntests; i++) {
- if (ntests>1) {
- sprintf(str, "Testing random %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- } else {
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- }
- printf("%-70s", str);
- HDfflush(stdout);
- fails_this_test=0;
-
- /* Perform the conversion */
- if (H5Tconvert(src, dst, nelmts, buf, NULL, H5P_DEFAULT)<0)
- goto error;
-
- /* Check the results from the library against hardware */
- for (j=0; j<nelmts; j++) {
- if (INT_SCHAR==dst_type) {
- hw = (unsigned char*)&hw_char;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_char = (char)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_char = (char)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_char = (char)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_char = (char)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_char = (char)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_char = (char)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_char = (char)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_char = (char)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_char = (char)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_char = (char)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_UCHAR==dst_type) {
- hw = (unsigned char*)&hw_uchar;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_uchar = (unsigned char)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_uchar = (unsigned char)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_uchar = (unsigned char)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_uchar = (unsigned char)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_uchar = (unsigned char)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_uchar = (unsigned char)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_uchar = (unsigned char)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_uchar = (unsigned char)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_uchar = (unsigned char)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_uchar = (unsigned char)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_SHORT==dst_type) {
- hw = (unsigned char*)&hw_short;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_short = (short)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_short = (short)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_short = (short)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_short = (short)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_short = (short)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_short = (short)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_short = (short)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_short = (short)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_short = (short)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_short = (short)(*((unsigned long_long*)aligned));
- break;
-
- default:
- break;
- }
- } else if (INT_USHORT==dst_type) {
- hw = (unsigned char*)&hw_ushort;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_ushort = (unsigned short)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_ushort = (unsigned short)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_ushort = (unsigned short)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_ushort = (unsigned short)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_ushort = (unsigned short)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_ushort = (unsigned short)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_ushort = (unsigned short)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_ushort = (unsigned short)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_ushort = (unsigned short)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_ushort = (unsigned short)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_INT==dst_type) {
- hw = (unsigned char*)&hw_int;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_int = (int)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_int = (int)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_int = (int)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_int = (int)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_int = (int)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_int = (int)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_int = (int)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_int = (int)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_int = (int)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_int = (int)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_UINT==dst_type) {
- hw = (unsigned char*)&hw_uint;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_uint = (unsigned int)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_uint = (unsigned int)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_uint = (unsigned int)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_uint = (unsigned int)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_uint = (unsigned int)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_uint = (unsigned int)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_uint = (unsigned int)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_uint = (unsigned int)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_uint = (unsigned int)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_uint = (unsigned int)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_LONG==dst_type) {
- hw = (unsigned char*)&hw_long;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_long = (long int)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_long = (long int)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_long = (long int)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_long = (long int)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_long = (long int)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_long = (long int)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_long = (long int)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_long = (long int)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_long = (long int)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_long = (long int)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_ULONG==dst_type) {
- hw = (unsigned char*)&hw_ulong;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_ulong = (unsigned long)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_ulong = (unsigned long)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_ulong = (unsigned long)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_ulong = (unsigned long)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_ulong = (unsigned long)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_ulong = (unsigned long)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_ulong = (unsigned long)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_ulong = (unsigned long)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_ulong = (unsigned long)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_ulong = (unsigned long)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_LLONG==dst_type) {
- hw = (unsigned char*)&hw_llong;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_llong = (long_long)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_llong = (long_long)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_llong = (long_long)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_llong = (long_long)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_llong = (long_long)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_llong = (long_long)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_llong = (long_long)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_llong = (long_long)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_llong = (long_long)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_llong = (long_long)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_ULLONG==dst_type) {
- hw = (unsigned char*)&hw_ullong;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- hw_ullong = (unsigned long_long)(*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_ullong = (unsigned long_long)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_ullong = (unsigned long_long)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_ullong = (unsigned long_long)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_ullong = (unsigned long_long)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_ullong = (unsigned long_long)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_ullong = (unsigned long_long)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_ullong = (unsigned long_long)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_ullong = (unsigned long_long)(*((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_ullong = (unsigned long_long)(*((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
- }
-
- /* Make certain that there isn't some weird number of destination bits */
- assert(dst_nbits%8==0);
-
- /* Are the two results the same? */
- for (k=(dst_size-(dst_nbits/8)); k<dst_size; k++)
- if (buf[j*dst_size+k]!=hw[k])
- break;
- if (k==dst_size)
- continue; /*no error*/
-
- /*
- * Convert the source and destination values to little endian
- * order so we can use the HDF5 bit vector operations to test
- * certain things. These routines have already been tested by
- * the `bittests' program.
- */
- for (k=0; k<src_size; k++)
- src_bits[src_size-(k+1)] = saved[j*src_size+ENDIAN(src_size, k)];
-
- for (k=0; k<dst_size; k++)
- dst_bits[dst_size-(k+1)] = buf[j*dst_size+ENDIAN(dst_size, k)];
-
- /*
- * Hardware usually doesn't handle overflows too gracefully. The
- * hardware conversion result during overflows is usually garbage
- * so we must handle those cases differetly when checking results.
- */
- if (H5T_SGN_2==src_sign && H5T_SGN_2==dst_sign) {
- if (src_nbits>dst_nbits) {
- if(0==H5T_bit_get_d(src_bits, src_nbits-1, 1) &&
- H5T_bit_find(src_bits, dst_nbits-1, (src_nbits-dst_nbits),
- H5T_BIT_MSB, 1)>=0) {
- /*
- * Source is positive and the magnitude is too large for
- * the destination. The destination should be set to the
- * maximum possible value: 0x7f...f
- */
- if (0==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) &&
- H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- } else if (1==H5T_bit_get_d(src_bits, src_nbits-1, 1) &&
- H5T_bit_find(src_bits, 0, src_nbits-1, H5T_BIT_MSB,
- 0)+1>=(ssize_t)dst_nbits) {
- /*
- * Source is negative but the magnitude is too large for
- * the destination. The destination should be set to the
- * smallest possible value: 0x80...0
- */
- if (1==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) &&
- H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 1)<0)
- continue; /*no error*/
- }
- } else if(src_nbits<dst_nbits) {
- /* Source is smaller than the destination */
- if(0==H5T_bit_get_d(src_bits, src_nbits-1, 1)) {
- /*
- * Source is positive, so the excess bits in the
- * destination should be set to 0's.
- */
- if (0==H5T_bit_get_d(dst_bits, src_nbits-1, 1) &&
- H5T_bit_find(dst_bits, src_nbits, dst_nbits-src_nbits, H5T_BIT_LSB, 1)<0)
- continue; /*no error*/
- } else {
- /*
- * Source is negative, so the excess bits in the
- * destination should be set to 1's.
- */
- if (1==H5T_bit_get_d(dst_bits, src_nbits-1, 1) &&
- H5T_bit_find(dst_bits, src_nbits, dst_nbits-src_nbits, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- }
- }
- } else if (H5T_SGN_2==src_sign && H5T_SGN_NONE==dst_sign) {
- if (H5T_bit_get_d(src_bits, src_nbits-1, 1)) {
- /*
- * The source is negative so the result should be zero.
- * The source is negative if the most significant bit is
- * set. The destination is zero if all bits are zero.
- */
- if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 1)<0)
- continue; /*no error*/
- } else if (src_nbits>dst_nbits &&
- H5T_bit_find(src_bits, dst_nbits-1,
- src_nbits-dst_nbits, H5T_BIT_LSB, 1)>=0) {
- /*
- * The source is a value with a magnitude too large for
- * the destination. The destination should be the
- * largest possible value: 0xff...f
- */
- if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- }
- } else if (H5T_SGN_NONE==src_sign && H5T_SGN_2==dst_sign) {
- if (src_nbits>=dst_nbits &&
- H5T_bit_find(src_bits, dst_nbits-1, (src_nbits-dst_nbits)+1,
- H5T_BIT_LSB, 1)>=0) {
- /*
- * The source value has a magnitude that is larger than
- * the destination can handle. The destination should be
- * set to the largest possible positive value: 0x7f...f
- */
- if (0==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) &&
- H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- }
- } else {
- if (src_nbits>dst_nbits &&
- H5T_bit_find(src_bits, dst_nbits, src_nbits-dst_nbits,
- H5T_BIT_LSB, 1)>=0) {
- /*
- * The unsigned source has a value which is too large for
- * the unsigned destination. The destination should be
- * set to the largest possible value: 0xff...f
- */
- if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- }
- }
-
- /* Print errors */
- if (0==fails_this_test++)
- H5_FAILED();
- printf(" test %u elmt %u\n", (unsigned)i+1, (unsigned)j);
-
- printf(" src = ");
- for (k=0; k<src_size; k++)
- printf(" %02x", saved[j*src_size+ENDIAN(src_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)dst_size-(ssize_t)src_size)), "");
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(char), sizeof(char));
- printf(" %29d\n", (int)*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- printf(" %29u\n", (unsigned)*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- printf(" %29hd\n", *((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- printf(" %29hu\n", *((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- printf(" %29d\n", *((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- printf(" %29u\n", *((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- printf(" %29ld\n", *((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- printf(" %29lu\n", *((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
-
- printf(" dst = ");
- for (k=0; k<dst_size; k++)
- printf(" %02x", buf[j*dst_size+ENDIAN(dst_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- switch (dst_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, buf+j*sizeof(char), sizeof(char));
- printf(" %29d\n", (int)*((char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, buf+j*sizeof(unsigned char), sizeof(unsigned char));
- printf(" %29u\n", (unsigned)*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, buf+j*sizeof(short), sizeof(short));
- printf(" %29hd\n", *((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, buf+j*sizeof(unsigned short), sizeof(unsigned short));
- printf(" %29hu\n", *((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, buf+j*sizeof(int), sizeof(int));
- printf(" %29d\n", *((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, buf+j*sizeof(unsigned), sizeof(unsigned));
- printf(" %29u\n", *((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, buf+j*sizeof(long), sizeof(long));
- printf(" %29ld\n", *((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, buf+j*sizeof(unsigned long), sizeof(unsigned long));
- printf(" %29lu\n", *((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, buf+j*sizeof(long_long), sizeof(long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, buf+j*sizeof(long_long), sizeof(unsigned long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)aligned));
- break;
- default:
- break;
- }
-
- printf(" ans = ");
- for (k=0; k<dst_size; k++)
- printf(" %02x", hw[ENDIAN(dst_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- switch (dst_type) {
- case INT_SCHAR:
- printf(" %29d\n", (int)*((signed char*)hw));
- break;
- case INT_UCHAR:
- printf(" %29u\n", (unsigned)*((unsigned char*)hw));
- break;
- case INT_SHORT:
- printf(" %29hd\n", *((short*)hw));
- break;
- case INT_USHORT:
- printf(" %29hu\n", *((unsigned short*)hw));
- break;
- case INT_INT:
- printf(" %29d\n", *((int*)hw));
- break;
- case INT_UINT:
- printf(" %29u\n", *((unsigned*)hw));
- break;
- case INT_LONG:
- printf(" %29ld\n", *((long*)hw));
- break;
- case INT_ULONG:
- printf(" %29lu\n", *((unsigned long*)hw));
- break;
- case INT_LLONG:
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)hw));
- break;
- case INT_ULLONG:
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)hw));
- break;
- default:
- break;
- }
-
- if (++fails_all_tests>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto done;
- }
- }
- PASSED();
- }
-#ifdef SHOW_OVERFLOWS
- if (noverflows_g>0) {
- printf(" %d overflow%s in previous test\n",
- noverflows_g, 1==noverflows_g?"":"s");
- }
-#endif
-
-done:
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
- reset_hdf5(); /*print statistics*/
- return (int)fails_all_tests;
-
-error:
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
- reset_hdf5(); /*print statistics*/
- return MAX((int)fails_all_tests, 1);
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: test_conv_int_2
- *
- * Purpose: Tests overlap calculates in H5T_conv_i_i(), which should be
- * the same as for H5T_conv_f_f() and H5T_conv_s_s().
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Robb Matzke
- * Friday, April 30, 1999
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_conv_int_2(void)
-{
- int i, j;
- hid_t src_type, dst_type;
- char buf[32*100];
-
- printf("%-70s", "Testing overlap calculations");
- HDfflush(stdout);
-
- HDmemset(buf, 0, sizeof buf);
- for (i=1; i<=32; i++) {
- for (j=1; j<=32; j++) {
-
- /* Source type */
- src_type = H5Tcopy(H5T_NATIVE_CHAR);
- H5Tset_size(src_type, (size_t)i);
-
- /* Destination type */
- dst_type = H5Tcopy(H5T_NATIVE_CHAR);
- H5Tset_size(dst_type, (size_t)j);
-
- /*
- * Conversion. If overlap calculations aren't right then an
- * assertion will fail in H5T_conv_i_i()
- */
- H5Tconvert(src_type, dst_type, 100, buf, NULL, H5P_DEFAULT);
- H5Tclose(src_type);
- H5Tclose(dst_type);
- }
- }
- PASSED();
- return 0;
-}
-
-
-/*-------------------------------------------------------------------------
- * 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(dtype_t type, void *val)
-{
- int retval;
- char s[256];
-
- if (FLT_FLOAT==type) {
- float x;
- HDmemcpy(&x, val, sizeof(float));
- retval = (x!=x);
- } else if (FLT_DOUBLE==type) {
- double x;
- HDmemcpy(&x, val, sizeof(double));
- retval = (x!=x);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (FLT_LDOUBLE==type) {
- long double x;
- HDmemcpy(&x, val, sizeof(long double));
- retval = (x!=x);
-#endif
- } else {
- 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) {
- if (FLT_FLOAT==type) {
- float x;
- HDmemcpy(&x, val, sizeof(float));
- sprintf(s, "%g", x);
- } else if (FLT_DOUBLE==type) {
- double x;
- HDmemcpy(&x, val, sizeof(double));
- sprintf(s, "%g", x);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (FLT_LDOUBLE==type) {
- long double x;
- HDmemcpy(&x, val, sizeof(long double));
- sprintf(s, "%Lg", x);
-#endif
- } else {
- return 0;
- }
- if (HDstrstr(s, "NaN") || HDstrstr(s, "NAN") || HDstrstr(s, "nan"))
- retval = 1;
- }
-
- return retval;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: test_conv_flt_1
- *
- * Purpose: Test conversion of random floating point values from SRC to
- * DST. These types should be H5T_NATIVE_FLOAT,
- * H5T_NATIVE_DOUBLE, or H5T_NATIVE_LDOUBLE.
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Robb Matzke
- * Tuesday, June 23, 1998
- *
- * Modifications:
- * Albert Cheng, Apr 16, 2004
- * Check for underflow condition. If the src number is
- * smaller than the dst MIN float number, consider it okay
- * if the converted sw and hw dst are both less than or
- * equal to the dst MIN float number.
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_conv_flt_1 (const char *name, hbool_t run_special, hid_t src, hid_t dst)
-{
- dtype_t src_type, dst_type; /*data types */
- const size_t ntests=NTESTS; /*number of tests */
- size_t nelmts=0; /*num values per test */
- const size_t max_fails=8; /*max number of failures*/
- size_t fails_all_tests=0; /*number of failures */
- size_t fails_this_test; /*fails for this test */
- const char *src_type_name = NULL; /*source type name */
- const char *dst_type_name = NULL; /*destination type name */
- size_t src_size, dst_size; /*type sizes */
- unsigned char *buf = NULL; /*buffer for conversion */
- unsigned char *saved = NULL; /*original values */
- char str[256]; /*hello string */
- float hw_f; /*hardware-converted */
- double hw_d; /*hardware-converted */
- void *aligned=NULL; /*aligned buffer */
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- long double hw_ld; /*hardware-converted */
-#endif
- unsigned char *hw=NULL; /*ptr to hardware-conv'd*/
- int underflow; /*underflow occurred */
- int uflow=0; /*underflow debug counters*/
- size_t i, j, k; /*counters */
- int endian; /*machine endianess */
- size_t dst_ebias; /* Destination type's exponent bias */
- size_t src_epos; /* Source type's exponent position */
- size_t src_esize; /* Source type's exponent size */
- size_t dst_epos; /* Destination type's exponent position */
- size_t dst_esize; /* Destination type's exponent size */
- size_t dst_msize; /* Destination type's mantissa size */
- size_t src_nbits; /* source length in bits */
- size_t dst_nbits; /* dst length in bits */
-
-#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.
- */
- HDfflush(stdout);
- HDfflush(stderr);
- if ((child_pid=fork())<0) {
- HDperror("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 {
- HDputs(" Child didn't exit normally.");
- return 1;
- }
- }
-#endif
-
- /*
- * The remainder of this function is executed only by the child if
- * HANDLE_SIGFPE is defined.
- */
- HDsignal(SIGFPE,fpe_handler);
-
- /* What are the names of the source and destination types */
- if (H5Tequal(src, H5T_NATIVE_FLOAT)) {
- src_type_name = "float";
- src_type = FLT_FLOAT;
- } else if (H5Tequal(src, H5T_NATIVE_DOUBLE)) {
- src_type_name = "double";
- src_type = FLT_DOUBLE;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (H5Tequal(src, H5T_NATIVE_LDOUBLE)) {
- src_type_name = "long double";
- src_type = FLT_LDOUBLE;
-#endif
- } else {
- src_type_name = "UNKNOWN";
- src_type = OTHER;
- }
-
- if (H5Tequal(dst, H5T_NATIVE_FLOAT)) {
- dst_type_name = "float";
- dst_type = FLT_FLOAT;
- } else if (H5Tequal(dst, H5T_NATIVE_DOUBLE)) {
- dst_type_name = "double";
- dst_type = FLT_DOUBLE;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (H5Tequal(dst, H5T_NATIVE_LDOUBLE)) {
- dst_type_name = "long double";
- dst_type = FLT_LDOUBLE;
-#endif
- } else {
- dst_type_name = "UNKNOWN";
- dst_type = OTHER;
- }
-
- /* Sanity checks */
- if(sizeof(float)==sizeof(double))
- HDputs("Sizeof(float)==sizeof(double) - some tests may not be sensible.");
- if (OTHER==src_type || OTHER==dst_type) {
- if(!strcmp(name, "noop"))
- sprintf(str, "Testing %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else if(run_special)
- sprintf(str, "Testing special %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else
- sprintf(str, "Testing regular %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
-
- printf("%-70s", str);
- H5_FAILED();
- HDputs(" Unknown data type.");
- goto error;
- }
-
- /* Get "interesting" values */
- src_size = H5Tget_size(src);
- dst_size = H5Tget_size(dst);
- src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */
- dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */
- dst_ebias=H5Tget_ebias(dst);
- H5Tget_fields(src,NULL,&src_epos,&src_esize,NULL,NULL);
- H5Tget_fields(dst,NULL,&dst_epos,&dst_esize,NULL,&dst_msize);
- endian = H5Tget_order(H5T_NATIVE_FLOAT);
-
- /* Allocate buffers */
- aligned = HDmalloc(MAX(sizeof(long double), sizeof(double)));
-#ifdef SHOW_OVERFLOWS
- noverflows_g = 0;
-#endif
-
- /* Allocate and initialize the source buffer through macro INIT_FP or INIT_FP_SPECIAL.
- * The BUF will be used for the conversion while the SAVED buffer will be used for
- * the comparison later. INIT_FP will fill in the buffer with regular values like
- * normalized and denormalized values; INIT_FP_SPECIAL will fill with special values
- * like infinity, NaN.
- */
- switch (run_special) {
- case FALSE:
- if(src_type == FLT_FLOAT) {
- INIT_FP(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, FLT_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
- } else if(src_type == FLT_DOUBLE) {
- INIT_FP(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, DBL_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if(src_type == FLT_LDOUBLE) {
- INIT_FP(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, LDBL_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
-#endif
- } else
- goto error;
-
- break;
- case TRUE:
- if(src_type == FLT_FLOAT) {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, FLT_MANT_DIG, dst_size,
- buf, saved, nelmts);
- } else if(src_type == FLT_DOUBLE) {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, DBL_MANT_DIG, dst_size,
- buf, saved, nelmts);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if(src_type == FLT_LDOUBLE) {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, LDBL_MANT_DIG, dst_size,
- buf, saved, nelmts);
-#endif
- } else
- goto error;
-
- break;
- default:
- goto error;
- }
-
- for (i=0; i<ntests; i++) {
- /*
- * If it looks like it might take a long time then print a progress
- * report between each test.
- */
- if (ntests>1) {
- if(!strcmp(name, "noop"))
- sprintf(str, "Testing %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- else if(run_special)
- sprintf(str, "Testing special %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- else
- sprintf(str, "Testing regular %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- } else {
- if(!strcmp(name, "noop"))
- sprintf(str, "Testing %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- else if(run_special)
- sprintf(str, "Testing special %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else
- sprintf(str, "Testing regular %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- }
- printf("%-70s", str);
- HDfflush(stdout);
- fails_this_test = 0;
-
- /* Perform the conversion in software */
- if (H5Tconvert(src, dst, nelmts, buf, NULL, H5P_DEFAULT)<0)
- goto error;
-
- /* Check the software results against the hardware */
- for (j=0; j<nelmts; j++) {
- underflow = 0;
- hw_f = 911.0;
- hw_d = 911.0;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- hw_ld = 911.0;
-#endif
-
- /* The hardware conversion */
- /* Check for underflow when src is a "larger" float than dst.*/
- if (FLT_FLOAT==src_type) {
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- if (FLT_FLOAT==dst_type) {
- hw_f = *((float*)aligned);
- hw = (unsigned char*)&hw_f;
- } else if (FLT_DOUBLE==dst_type) {
- hw_d = *((float*)aligned);
- hw = (unsigned char*)&hw_d;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- hw_ld = *((float*)aligned);
- hw = (unsigned char*)&hw_ld;
-#endif
- }
- } else if (FLT_DOUBLE==src_type) {
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- if (FLT_FLOAT==dst_type) {
- hw_f = (float)(*((double*)aligned));
- hw = (unsigned char*)&hw_f;
- underflow = HDfabs(*((double*)aligned)) < FLT_MIN;
- } else if (FLT_DOUBLE==dst_type) {
- hw_d = *((double*)aligned);
- hw = (unsigned char*)&hw_d;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- hw_ld = *((double*)aligned);
- hw = (unsigned char*)&hw_ld;
-#endif
- }
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- if (FLT_FLOAT==dst_type) {
- hw_f = *((long double*)aligned);
- hw = (unsigned char*)&hw_f;
- underflow = HDfabsl(*((long double*)aligned)) < FLT_MIN;
- } else if (FLT_DOUBLE==dst_type) {
- hw_d = *((long double*)aligned);
- hw = (unsigned char*)&hw_d;
- underflow = HDfabsl(*((long double*)aligned)) < DBL_MIN;
- } else {
- hw_ld = *((long double*)aligned);
- hw = (unsigned char*)&hw_ld;
- }
-#endif
- }
- if (underflow){
- uflow++;
- }
-
- /* For Intel machines, the size of "long double" is 12 bytes, precision
- * is 80 bits; for Intel IA64 and AMD processors, the size of "long double"
- * is 16 bytes, precision is 80 bits. During hardware conversion, the
- * last few unused bytes may have garbage in them. Clean them out with
- * 0s before compare the values.
- */
- if(endian==H5T_ORDER_LE && dst_type==FLT_LDOUBLE) {
- int q;
- for(q=dst_nbits/8; q<dst_size; q++) {
- buf[j*dst_size+q] = 0x00;
- hw[q] = 0x00;
- }
- }
-
- /* Are the two results the same? */
- for (k=(dst_size-(dst_nbits/8)); k<dst_size; k++)
- if (buf[j*dst_size+k]!=hw[k])
- break;
- if (k==dst_size)
- continue; /*no error*/
-
- /*
- * Assume same if both results are NaN. There are many NaN bit
- * patterns and the software doesn't attemt to emulate the
- * hardware in this regard. Instead, software uses a single bit
- * pattern for NaN by setting the significand to all ones.
- */
- if (FLT_FLOAT==dst_type &&
- my_isnan(dst_type, buf+j*sizeof(float)) &&
- my_isnan(dst_type, hw)) {
- continue;
- } else if (FLT_DOUBLE==dst_type &&
- my_isnan(dst_type, buf+j*sizeof(double)) &&
- my_isnan(dst_type, hw)) {
- continue;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (FLT_LDOUBLE==dst_type &&
- my_isnan(dst_type, buf+j*sizeof(long double)) &&
- my_isnan(dst_type, hw)) {
- continue;
-#endif
- }
-
- /*
- * Assume same if hardware result is NaN. This is because the
- * hardware conversions on some machines return NaN instead of
- * overflowing to +Inf or -Inf or underflowing to +0 or -0.
- */
- if (my_isnan(dst_type, hw))
- continue;
-
- /*
- * Instead of matching down to the bit, just make sure the
- * exponents are the same and the mantissa is the same to a
- * certain precision. This is needed on machines that don't
- * round as expected.
- * If the src number is smaller than the dst MIN float number,
- * consider it okay if the converted sw and hw dst are both
- * less than or equal to the dst MIN float number.
- */
- {
- double check_mant[2];
- int check_expo[2];
-
- if (FLT_FLOAT==dst_type) {
- float x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(float));
- if (underflow &&
- HDfabsf(x) <= FLT_MIN && HDfabsf(hw_f) <= FLT_MIN)
- continue; /* all underflowed, no error */
- check_mant[0] = HDfrexpf(x, check_expo+0);
- check_mant[1] = HDfrexpf(hw_f, check_expo+1);
- } else if (FLT_DOUBLE==dst_type) {
- double x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(double));
- if (underflow &&
- HDfabs(x) <= DBL_MIN && HDfabs(hw_d) <= DBL_MIN)
- continue; /* all underflowed, no error */
- check_mant[0] = HDfrexp(x, check_expo+0);
- check_mant[1] = HDfrexp(hw_d, check_expo+1);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- long double x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(long double));
- /* dst is largest float, no need to check underflow. */
- check_mant[0] = HDfrexpl(x, check_expo+0);
- check_mant[1] = HDfrexpl(hw_ld, check_expo+1);
-#endif
- }
-#ifdef H5_CONVERT_DENORMAL_FLOAT
- /* Special check for denormalized values */
- if(check_expo[0]<(-(int)dst_ebias) || check_expo[1]<(-(int)dst_ebias)) {
- int expo_diff=check_expo[0]-check_expo[1];
- int valid_bits=(int)((dst_ebias+dst_msize)+MIN(check_expo[0],check_expo[1]))-1;
- double epsilon=1.0;
-
- /* Re-scale the mantissas based on any exponent difference */
- if(expo_diff!=0)
- check_mant[0] = HDldexp(check_mant[0],expo_diff);
-
- /* Compute the proper epsilon */
- epsilon=HDldexp(epsilon,-valid_bits);
-
- /* Check for "close enough" fit with scaled epsilon value */
- if (HDfabs(check_mant[0]-check_mant[1])<=epsilon)
- continue;
- } /* end if */
- else {
- if (check_expo[0]==check_expo[1] &&
- HDfabs(check_mant[0]-check_mant[1])<FP_EPSILON)
- continue;
- } /* end else */
-#else /* H5_CONVERT_DENORMAL_FLOAT */
- {
- hssize_t expo; /*exponent */
- uint8_t tmp[32];
-
- assert(src_size<=sizeof(tmp));
- if(endian==H5T_ORDER_LE)
- HDmemcpy(tmp,&saved[j*src_size],src_size);
- else
- for (k=0; k<src_size; k++)
- tmp[k]=saved[j*src_size+(src_size-(k+1))];
- expo = H5T_bit_get_d(tmp, src_epos, src_esize);
- if(expo==0)
- continue; /* Denormalized floating-point value detected */
- else {
- assert(dst_size<=sizeof(tmp));
- if(endian==H5T_ORDER_LE)
- HDmemcpy(tmp,&buf[j*dst_size],dst_size);
- else
- for (k=0; k<dst_size; k++)
- tmp[k]=buf[j*dst_size+(dst_size-(k+1))];
- expo = H5T_bit_get_d(tmp, dst_epos, dst_esize);
- if(expo==0)
- continue; /* Denormalized floating-point value detected */
- else {
- if (check_expo[0]==check_expo[1] &&
- HDfabs(check_mant[0]-check_mant[1])<FP_EPSILON)
- continue;
- } /* end else */
- } /* end else */
- }
-#endif /* H5_CONVERT_DENORMAL_FLOAT */
- }
-
- if (0==fails_this_test++)
- H5_FAILED();
- printf(" test %u, elmt %u\n", (unsigned)i+1, (unsigned)j);
-
- printf(" src =");
- for (k=0; k<src_size; k++)
- printf(" %02x", saved[j*src_size+ENDIAN(src_size,k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)dst_size-(ssize_t)src_size)), "");
- if (FLT_FLOAT==src_type) {
- float x;
- HDmemcpy(&x, &saved[j*src_size], sizeof(float));
- printf(" %29.20e\n", x);
- } else if (FLT_DOUBLE==src_type) {
- double x;
- HDmemcpy(&x, &saved[j*src_size], sizeof(double));
- printf(" %29.20e\n", x);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- long double x;
- HDmemcpy(&x, &saved[j*src_size], sizeof(long double));
- HDfprintf(stdout," %29.20Le\n", x);
-#endif
- }
-
- printf(" dst =");
- for (k=0; k<dst_size; k++)
- printf(" %02x", buf[j*dst_size+ENDIAN(dst_size,k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- if (FLT_FLOAT==dst_type) {
- float x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(float));
- printf(" %29.20e\n", x);
- } else if (FLT_DOUBLE==dst_type) {
- double x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(double));
- printf(" %29.20e\n", x);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else {
- long double x;
- HDmemcpy(&x, &buf[j*dst_size], sizeof(long double));
- HDfprintf(stdout," %29.20Le\n", x);
-#endif
- }
-
- printf(" ans =");
- for (k=0; k<dst_size; k++)
- printf(" %02x", hw[ENDIAN(dst_size,k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- if (FLT_FLOAT==dst_type)
- printf(" %29.20e\n", hw_f);
- else if (FLT_DOUBLE==dst_type)
- printf(" %29.20e\n", hw_d);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- else
- HDfprintf(stdout," %29.20Le\n", hw_ld);
-#endif
-
- if (++fails_all_tests>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto done;
- }
- }
- PASSED();
- }
-#ifdef SHOW_OVERFLOWS
- if (noverflows_g>0)
- printf(" %d overflow%s in previous test\n",
- noverflows_g, 1==noverflows_g?"":"s");
-#endif
-
- done:
-#ifdef AKCDEBUG
- printf("uflow=%d, fails_all_tests=%d\n", uflow, fails_all_tests);
-#endif
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
-#ifdef HANDLE_SIGFPE
- HDexit(MIN((int)fails_all_tests, 254));
-#else
- reset_hdf5();
- return (int)fails_all_tests;
-#endif
-
- error:
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
-#ifdef HANDLE_SIGFPE
- HDexit(MIN(MAX((int)fails_all_tests, 1), 254));
-#else
- reset_hdf5();
- return MAX((int)fails_all_tests, 1);
-#endif
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: test_conv_int_float
- *
- * Purpose: Test conversion between random integer and float values
- * from SRC to DST. These types should be any combination of:
- *
- * H5T_NATIVE_SCHAR H5T_NATIVE_FLOAT
- * H5T_NATIVE_SHORT H5T_NATIVE_DOUBLE
- * H5T_NATIVE_INT H5T_NATIVE_LDOUBLE
- * H5T_NATIVE_LONG
- * H5T_NATIVE_LLONG
- *
- * Return: Success: 0
- *
- * Failure: number of errors
- *
- * Programmer: Raymond Lu
- * Thursday, November 6, 2003
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-test_conv_int_float(const char *name, hbool_t run_special, hid_t src, hid_t dst)
-{
- hid_t dxpl_id; /*dataset transfer property list*/
- int fill_value=9; /*fill value for conversion exception*/
- H5T_conv_except_func_t op; /*returned callback function for conversion exception*/
- void *user_data; /*returned pointer to user data passed in to the callback*/
- hbool_t except_set = FALSE; /*whether user's exception handling is set*/
- const size_t ntests=NTESTS; /*number of tests */
- size_t nelmts=0; /*num values per test */
- const size_t max_fails=40; /*max number of failures*/
- size_t fails_all_tests=0; /*number of failures */
- size_t fails_this_test; /*fails for this test */
- char str[256]; /*hello string */
- dtype_t src_type; /*data types */
- dtype_t dst_type; /*data types */
- const char *src_type_name=NULL; /*source type name */
- const char *dst_type_name=NULL; /*destination type name */
- int endian; /*machine endianess */
- size_t src_size, dst_size; /*type sizes */
- unsigned char *buf=NULL; /*buffer for conversion */
- unsigned char *saved=NULL; /*original values */
- size_t i, j, k; /*counters */
- unsigned char *hw=NULL; /*hardware conv result */
- unsigned char src_bits[32]; /*src value in LE order */
- unsigned char dst_bits[32]; /*dest value in LE order*/
- size_t src_nbits; /*source length in bits */
- size_t dst_nbits; /*dst length in bits */
- void *aligned=NULL; /*aligned temp buffer */
- float hw_float=0;
- double hw_double=0;
- long double hw_ldouble=0;
- signed char hw_schar=0;
- unsigned char hw_uchar=0;
- short hw_short=0;
- unsigned short hw_ushort=0;
- int hw_int=0;
- unsigned hw_uint=0;
- long hw_long=0;
- unsigned long hw_ulong=0;
- long_long hw_llong=0;
- unsigned long_long hw_ullong=0;
-
- /* What is the name of the source type */
- if (H5Tequal(src, H5T_NATIVE_SCHAR)) {
- src_type_name = "signed char";
- src_type = INT_SCHAR;
- } else if (H5Tequal(src, H5T_NATIVE_UCHAR)) {
- src_type_name = "unsigned char";
- src_type = INT_UCHAR;
- } else if (H5Tequal(src, H5T_NATIVE_SHORT)) {
- src_type_name = "short";
- src_type = INT_SHORT;
- } else if (H5Tequal(src, H5T_NATIVE_USHORT)) {
- src_type_name = "unsigned short";
- src_type = INT_USHORT;
- } else if (H5Tequal(src, H5T_NATIVE_INT)) {
- src_type_name = "int";
- src_type = INT_INT;
- } else if (H5Tequal(src, H5T_NATIVE_UINT)) {
- src_type_name = "unsigned int";
- src_type = INT_UINT;
- } else if (H5Tequal(src, H5T_NATIVE_LONG)) {
- src_type_name = "long";
- src_type = INT_LONG;
- } else if (H5Tequal(src, H5T_NATIVE_ULONG)) {
- src_type_name = "unsigned long";
- src_type = INT_ULONG;
- } else if (H5Tequal(src, H5T_NATIVE_LLONG)) {
- src_type_name = "long long";
- src_type = INT_LLONG;
- } else if (H5Tequal(src, H5T_NATIVE_ULLONG)) {
- src_type_name = "unsigned long long";
- src_type = INT_ULLONG;
- } else if (H5Tequal(src, H5T_NATIVE_FLOAT)) {
- src_type_name = "float";
- src_type = FLT_FLOAT;
- } else if (H5Tequal(src, H5T_NATIVE_DOUBLE)) {
- src_type_name = "double";
- src_type = FLT_DOUBLE;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (H5Tequal(src, H5T_NATIVE_LDOUBLE)) {
- src_type_name = "long double";
- src_type = FLT_LDOUBLE;
-#endif
- } else {
- src_type_name = "UNKNOWN";
- src_type = OTHER;
- }
-
- /* What is the name of the destination type */
- if (H5Tequal(dst, H5T_NATIVE_SCHAR)) {
- dst_type_name = "signed char";
- dst_type = INT_SCHAR;
- } else if (H5Tequal(dst, H5T_NATIVE_UCHAR)) {
- dst_type_name = "unsigned char";
- dst_type = INT_UCHAR;
- } else if (H5Tequal(dst, H5T_NATIVE_SHORT)) {
- dst_type_name = "short";
- dst_type = INT_SHORT;
- } else if (H5Tequal(dst, H5T_NATIVE_USHORT)) {
- dst_type_name = "unsigned short";
- dst_type = INT_USHORT;
- } else if (H5Tequal(dst, H5T_NATIVE_INT)) {
- dst_type_name = "int";
- dst_type = INT_INT;
- } else if (H5Tequal(dst, H5T_NATIVE_UINT)) {
- dst_type_name = "unsigned int";
- dst_type = INT_UINT;
- } else if (H5Tequal(dst, H5T_NATIVE_LONG)) {
- dst_type_name = "long";
- dst_type = INT_LONG;
- } else if (H5Tequal(dst, H5T_NATIVE_ULONG)) {
- dst_type_name = "unsigned long";
- dst_type = INT_ULONG;
- } else if (H5Tequal(dst, H5T_NATIVE_LLONG)) {
- dst_type_name = "long long";
- dst_type = INT_LLONG;
- } else if (H5Tequal(dst, H5T_NATIVE_ULLONG)) {
- dst_type_name = "unsigned long long";
- dst_type = INT_ULLONG;
- } else if (H5Tequal(dst, H5T_NATIVE_FLOAT)) {
- dst_type_name = "float";
- dst_type = FLT_FLOAT;
- } else if (H5Tequal(dst, H5T_NATIVE_DOUBLE)) {
- dst_type_name = "double";
- dst_type = FLT_DOUBLE;
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if (H5Tequal(dst, H5T_NATIVE_LDOUBLE)) {
- dst_type_name = "long double";
- dst_type = FLT_LDOUBLE;
-#endif
- } else {
- dst_type_name = "UNKNOWN";
- dst_type = OTHER;
- }
-
- /* Sanity checks */
- if (OTHER==src_type || OTHER==dst_type) {
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- printf("%-70s", str);
- H5_FAILED();
- HDputs(" Unknown data type.");
- goto error;
- }
-
- if ((INT_SCHAR==src_type || INT_UCHAR==src_type || INT_SHORT==src_type ||
- INT_USHORT==src_type || INT_INT==src_type || INT_UINT==src_type ||
- INT_LONG==src_type || INT_ULONG==src_type || INT_LLONG==src_type ||
- INT_ULLONG==src_type) &&
- (FLT_FLOAT!=dst_type && FLT_DOUBLE!=dst_type && FLT_LDOUBLE!=dst_type)) {
- if(!run_special)
- sprintf(str, "Testing regular %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else
- sprintf(str, "Testing special %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- printf("%-70s", str);
- H5_FAILED();
- HDputs(" 1. Not an integer-float conversion.");
- goto error;
- }
-
- if ((FLT_FLOAT==src_type || FLT_DOUBLE==src_type || FLT_LDOUBLE==src_type)
- && (INT_SCHAR!=dst_type && INT_UCHAR!=dst_type && INT_SHORT!=dst_type
- && INT_USHORT!=dst_type && INT_INT!=dst_type && INT_UINT!=dst_type
- && INT_LONG!=dst_type && INT_ULONG!=dst_type && INT_LLONG!=dst_type
- && INT_ULLONG!=dst_type)) {
- if(!run_special)
- sprintf(str, "Testing regular %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else
- sprintf(str, "Testing special %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- printf("%-70s", str);
- H5_FAILED();
- HDputs(" 2. Not a float-integer conversion.");
- goto error;
- }
-
- /* Some information about datatypes */
- endian = H5Tget_order(H5T_NATIVE_INT);
- src_size = H5Tget_size(src);
- dst_size = H5Tget_size(dst);
- src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */
- dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */
- aligned = HDmalloc(MAX(sizeof(long double), sizeof(long_long)));
-#ifdef SHOW_OVERFLOWS
- noverflows_g = 0;
-#endif
-
- /* This is for some Linux systems where long double has the size
- * 12 bytes but precision is 10 bytes. The 2 unused bytes may
- * have garbage causing wrong value comparison.
- */
- HDmemset(&hw_ldouble, 0, sizeof(long double));
-
- /* Create a dataset transfer property list and datatype conversion
- * exception handler function and pass in fill value. This is mainly
- * for NetCDF compatibility, which requests fill in fill value when
- * conversion exception happens. We only test (unsigned) int - float
- * and float - (unsigned) int conversions, which should cover more cases.
- */
- if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0)
- goto error;
-
- if((src_type == INT_INT && dst_type == FLT_FLOAT) ||
- (src_type == INT_UINT && dst_type == FLT_FLOAT) ||
- (src_type == FLT_FLOAT && dst_type == INT_UINT) ||
- (src_type == FLT_FLOAT && dst_type == INT_INT)) {
- if(H5Pset_type_conv_cb(dxpl_id, except_func, &fill_value)<0)
- goto error;
- else
- except_set = TRUE;
-
- if(H5Pget_type_conv_cb(dxpl_id, &op, &user_data)<0)
- goto error;
-
- if(op != except_func || *(int*)user_data != fill_value)
- goto error;
- }
-
- /* Allocate and initialize the source buffer through macro INIT_INTEGER if the source is integer,
- * INIT_FP if floating-point. The BUF will be used for the conversion while the SAVED buffer will be
- * used for the comparison later.
- */
- if(src_type == INT_SCHAR) {
- INIT_INTEGER(signed char, SCHAR_MAX, SCHAR_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_UCHAR) {
- INIT_INTEGER(unsigned char, UCHAR_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_SHORT) {
- INIT_INTEGER(short, SHRT_MAX, SHRT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_USHORT) {
- INIT_INTEGER(unsigned short, USHRT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_INT) {
- INIT_INTEGER(int, INT_MAX, INT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_UINT) {
- INIT_INTEGER(unsigned int, UINT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_LONG) {
- INIT_INTEGER(long, LONG_MAX, LONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_ULONG) {
- INIT_INTEGER(unsigned long, ULONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_LLONG) {
- INIT_INTEGER(long_long, LLONG_MAX, LLONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == INT_ULLONG) {
- INIT_INTEGER(unsigned long_long, ULLONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts);
- } else if(src_type == FLT_FLOAT) {
- if(!run_special) {
- INIT_FP(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, FLT_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
- } else {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, FLT_MANT_DIG, dst_size, buf, saved, nelmts);
- }
- } else if(src_type == FLT_DOUBLE) {
- if(!run_special) {
- INIT_FP(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, DBL_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
- } else {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, DBL_MANT_DIG, dst_size, buf, saved, nelmts);
- }
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- } else if(src_type == FLT_LDOUBLE) {
- if(!run_special) {
- INIT_FP(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, LDBL_MANT_DIG,
- src_size, src_nbits, endian, dst_size, buf, saved, nelmts);
- } else {
- INIT_FP_SPECIAL(src_size, src_nbits, endian, LDBL_MANT_DIG, dst_size, buf, saved, nelmts);
- }
-#endif
- } else
- goto error;
-
- /* The tests */
- for (i=0; i<ntests; i++) {
- if (ntests>1) {
- if(!run_special)
- sprintf(str, "Testing regular %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- else
- sprintf(str, "Testing special %s %s -> %s conversions (test %d/%d)",
- name, src_type_name, dst_type_name, (int)i+1, (int)ntests);
- } else {
- if(!run_special)
- sprintf(str, "Testing regular %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- else
- sprintf(str, "Testing special %s %s -> %s conversions",
- name, src_type_name, dst_type_name);
- }
- printf("%-70s", str);
- HDfflush(stdout);
- fails_this_test=0;
-
- /* Perform the conversion */
- if (H5Tconvert(src, dst, nelmts, buf, NULL, dxpl_id)<0)
- goto error;
-
- /* Check the results from the library against hardware */
- for (j=0; j<nelmts; j++) {
-
- if (FLT_FLOAT==dst_type) {
- hw = (unsigned char*)&hw_float;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(signed char), sizeof(signed char));
- hw_float = (float)(*((signed char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_float = (float)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_float = (float)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_float = (float)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_float = (float)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_float = (float)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_float = (float)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_float = (float)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_float = (float)(*((long_long*)aligned));
- break;
-#ifdef H5_ULLONG_TO_FP_CAST_WORKS
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_float = (float)(*((unsigned long_long*)aligned));
- break;
-#endif /* H5_ULLONG_TO_FP_CAST_WORKS */
- default:
- break;
- }
- } else if (FLT_DOUBLE==dst_type) {
- hw = (unsigned char*)&hw_double;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(signed char), sizeof(signed char));
- hw_double = (double)(*((signed char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_double = (double)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_double = (double)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_double = (double)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_double = (double)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_double = (double)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_double = (double)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_double = (double)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_double = (double)(*((long_long*)aligned));
- break;
-#ifdef H5_ULLONG_TO_FP_CAST_WORKS
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_double = (double)(*((unsigned long_long*)aligned));
- break;
-#endif /* H5_ULLONG_TO_FP_CAST_WORKS */
- default:
- break;
- }
- } else if (FLT_LDOUBLE==dst_type) {
- hw = (unsigned char*)&hw_ldouble;
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(signed char), sizeof(signed char));
- hw_ldouble = (long double)(*((signed char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- hw_ldouble = (long double)(*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- hw_ldouble = (long double)(*((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- hw_ldouble = (long double)(*((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- hw_ldouble = (long double)(*((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- hw_ldouble = (long double)(*((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- hw_ldouble = (long double)(*((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- hw_ldouble = (long double)(*((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- hw_ldouble = (long double)(*((long_long*)aligned));
- break;
-#ifdef H5_ULLONG_TO_FP_CAST_WORKS
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- hw_ldouble = (long double)(*((unsigned long_long*)aligned));
- break;
-#endif /* H5_ULLONG_TO_FP_CAST_WORKS */
- default:
- break;
- }
- } else if (INT_SCHAR==dst_type) {
- hw = (unsigned char*)&hw_schar;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_schar = (signed char)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_schar = (signed char)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_schar = (signed char)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_UCHAR==dst_type) {
- hw = (unsigned char*)&hw_uchar;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_uchar = (unsigned char)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_uchar = (unsigned char)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_uchar = (unsigned char)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_SHORT==dst_type) {
- hw = (unsigned char*)&hw_short;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_short = (short)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_short = (short)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_short = (short)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_USHORT==dst_type) {
- hw = (unsigned char*)&hw_ushort;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_ushort = (unsigned short)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_ushort = (unsigned short)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_ushort = (unsigned short)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_INT==dst_type) {
- hw = (unsigned char*)&hw_int;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_int = (int)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_int = (int)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_int = (int)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_UINT==dst_type) {
- hw = (unsigned char*)&hw_uint;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_uint = (unsigned int)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_uint = (unsigned int)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_uint = (unsigned int)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_LONG==dst_type) {
- hw = (unsigned char*)&hw_long;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_long = (long)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_long = (long)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_long = (long)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_ULONG==dst_type) {
- hw = (unsigned char*)&hw_ulong;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_ulong = (unsigned long)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_ulong = (unsigned long)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_ulong = (unsigned long)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_LLONG==dst_type) {
- hw = (unsigned char*)&hw_llong;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_llong = (long_long)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_llong = (long_long)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_llong = (long_long)(*((long double*)aligned));
- break;
- default:
- break;
- }
- } else if (INT_ULLONG==dst_type) {
- hw = (unsigned char*)&hw_ullong;
- switch (src_type) {
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- hw_ullong = (unsigned long_long)(*((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- hw_ullong = (unsigned long_long)(*((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- hw_ullong = (unsigned long_long)(*((long double*)aligned));
- break;
- default:
- break;
- }
- }
-
- /* Make certain that there isn't some weird number of destination bits */
- assert(dst_nbits%8==0);
-
- /* For Intel machines, the size of "long double" is 12 bytes, precision
- * is 80 bits; for AMD processors, the size of "long double" is 16 bytes,
- * precision is 80 bits. During hardware conversion, the last few unused
- * bytes may have garbage in them. Clean them out with 0s before compare
- * the values.
- */
- if(endian==H5T_ORDER_LE && dst_type==FLT_LDOUBLE) {
- int q;
- for(q=dst_nbits/8; q<dst_size; q++) {
- buf[j*dst_size+q] = 0x00;
- }
- }
-
- /* Are the two results the same? */
- for (k=(dst_size-(dst_nbits/8)); k<dst_size; k++)
- if (buf[j*dst_size+k]!=hw[k])
- break;
- if (k==dst_size)
- continue; /*no error*/
-
- /* Test library's handling of NaN:
- * Hardware usually doesn't handle NaN too gracefully. The
- * hardware conversion result is usually garbage so we must handle
- * those cases differetly when checking results.
- *
- * Test user's exception handler when NaN happens:
- * Try to follow the except_func callback function to check if the
- * desired value was set.
- */
- if (run_special && (FLT_FLOAT==src_type || FLT_DOUBLE==src_type || FLT_LDOUBLE==src_type)) {
- if (my_isnan(src_type, saved+j*src_size)) {
- if(!except_set) {
- if(!strcmp(name, "sw")) {
- /*Although bit operation should be on little-endian order, this special
- *case should be fine because we're searching non-zero bit only.*/
- if (H5T_bit_find(buf+j*dst_size, 0, 8*dst_size-1, H5T_BIT_LSB, 1)<0)
- continue; /*no error. Software conversion set integer to 0.*/
- } else if(!strcmp(name, "hw")) {
- continue; /*don't compare because hardware conversion may fill in garbage*/
- }
- } else {
- /* fill_value is small so we only want to know the 1st byte is set. Also
- * consider endianess here. */
- if(endian==H5T_ORDER_LE)
- if(buf[j*dst_size] == fill_value)
- continue; /*no error*/
- else
- if(buf[(j+1)*dst_size-1] == fill_value)
- continue; /*no error*/
- }
- }
- }
-
- /*
- * Convert the source and destination values to little endian
- * order so we can use the HDF5 bit vector operations to test
- * certain things. These routines have already been tested by
- * the `bittests' program.
- */
- for (k=0; k<src_size; k++)
- src_bits[src_size-(k+1)] = saved[j*src_size+ENDIAN(src_size, k)];
-
- for (k=0; k<dst_size; k++)
- dst_bits[dst_size-(k+1)] = buf[j*dst_size+ENDIAN(dst_size, k)];
-
- /* Test library's default overflow handling:
- * Hardware usually doesn't handle overflows too gracefully. The
- * hardware conversion result during overflows is usually garbage
- * so we must handle those cases differetly when checking results.
- *
- * Test user's exception handler when overflows:
- * Try to follow the except_func callback function to check if the
- * desired value was set.
- */
- if ((FLT_FLOAT==src_type || FLT_DOUBLE==src_type || FLT_LDOUBLE==src_type)
- && (INT_SCHAR==dst_type || INT_SHORT==dst_type || INT_INT==dst_type
- || INT_LONG==dst_type || INT_LLONG==dst_type)) {
- if(0==H5T_bit_get_d(src_bits, src_nbits-1, 1) &&
- overflows(src_bits, src, dst_nbits-1)) {
- /*
- * Source is positive and the magnitude is too large for
- * the destination. The destination should be set to the
- * maximum possible value: 0x7f...f
- */
- if(!except_set) {
- if (0==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) &&
- H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- } else {
- /* fill_value is small so we know only the 1st byte is set */
- if (dst_bits[0] == fill_value)
- continue; /*no error*/
- }
- } else if (1==H5T_bit_get_d(src_bits, src_nbits-1, 1) &&
- overflows(src_bits, src, dst_nbits-1)) {
- /*
- * Source is negative but the magnitude is too large for
- * the destination. The destination should be set to the
- * smallest possible value: 0x80...0
- */
- if(!except_set) {
- if (1==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) &&
- H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 1)<0)
- continue; /*no error*/
- } else {
- if (dst_bits[0] == fill_value)
- continue; /*no error*/
- }
- }
- }
-
- if ((FLT_FLOAT==src_type || FLT_DOUBLE==src_type || FLT_LDOUBLE==src_type)
- && (INT_UCHAR==dst_type || INT_USHORT==dst_type || INT_UINT==dst_type
- || INT_ULONG==dst_type || INT_ULLONG==dst_type)) {
- if (H5T_bit_get_d(src_bits, src_nbits-1, 1)) {
- /*
- * The source is negative so the result should be zero.
- * The source is negative if the most significant bit is
- * set. The destination is zero if all bits are zero.
- */
- if(!except_set) {
- if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 1)<0)
- continue; /*no error*/
- } else {
- if (dst_bits[0] == fill_value)
- continue; /*no error*/
- }
- } else if (overflows(src_bits, src, dst_nbits)) {
- /*
- * The source is a value with a magnitude too large for
- * the destination. The destination should be the
- * largest possible value: 0xff...f
- */
- if(!except_set) {
- if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 0)<0)
- continue; /*no error*/
- } else {
- if (dst_bits[0] == fill_value)
- continue; /*no error*/
- }
- }
- }
-
-/* On some machines (notably the SGI and Solaris 64-bit machines) unsigned long
- * values are not converted to float or double values correctly, they are
- * consistently off by the lowest bit being rounded oppositely to our
- * software conversion routines output. So, on those machines, we allow
- * the converted value to be +/- 1 from the machine's value. -QAK
- */
-#ifndef H5_SW_ULONG_TO_FP_BOTTOM_BIT_WORKS
- if(dst_size==sizeof(unsigned)) {
- unsigned tmp_s, tmp_h;
- HDmemcpy(&tmp_s,&buf[j*dst_size],sizeof(unsigned));
- HDmemcpy(&tmp_h,&hw[0],sizeof(unsigned));
- if((tmp_s+1)==tmp_h || (tmp_s-1)==tmp_h)
- continue; /*no error*/
- } /* end if */
- else if (dst_size==sizeof(unsigned long)) {
- unsigned long tmp_s, tmp_h;
- HDmemcpy(&tmp_s,&buf[j*dst_size],sizeof(unsigned long));
- HDmemcpy(&tmp_h,&hw[0],sizeof(unsigned long));
- if((tmp_s+1)==tmp_h || (tmp_s-1)==tmp_h)
- continue; /*no error*/
- } /* end if */
- else if (dst_size==sizeof(unsigned long_long)) {
- unsigned long_long tmp_s, tmp_h;
- HDmemcpy(&tmp_s,&buf[j*dst_size],sizeof(unsigned long_long));
- HDmemcpy(&tmp_h,&hw[0],sizeof(unsigned long_long));
- if((tmp_s+1)==tmp_h || (tmp_s-1)==tmp_h)
- continue; /*no error*/
- } /* end if */
-#endif /* end H5_ULONG_FP_BOTTOM_BIT_WORKS */
-
-/* For PGI compiler on Linux, during conversion from 'float' or 'double' to
- * 'unsigned long long', round-up happens when the fraction of float-point
- * value is greater than 0.5. So we allow the converted value to be off by 1.
- */
-#ifndef H5_FP_TO_ULLONG_BOTTOM_BIT_WORKS
- if((src_type==FLT_FLOAT || src_type==FLT_DOUBLE) && dst_type==INT_ULLONG) {
- unsigned long_long tmp_s, tmp_h;
- HDmemcpy(&tmp_s,&buf[j*dst_size],sizeof(unsigned long_long));
- HDmemcpy(&tmp_h,&hw[0],sizeof(unsigned long_long));
- if((tmp_s+1)==tmp_h)
- continue; /*no error*/
- }
-#endif /*end H5_FP_TO_ULLONG_BOTTOM_BIT_WORKS*/
-
-/* For GNU compilers on FreeBSD(sleipnir), during conversion from 'unsigned long long'
- * to 'long double', the last 2 bytes of mantissa are lost. But this loss seems
- * acceptable. We allow it to go through instead of fail it. Sometimes, there's roundup
- * to the 3rd last byte of mantissa. So we only try to compare all but the last 3 bytes.
- */
-#ifndef H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS
- if(src_type==INT_ULLONG && dst_type==FLT_LDOUBLE) {
- long double tmp_s, tmp_h;
- HDmemcpy(&tmp_s,&buf[j*dst_size],sizeof(long double));
- HDmemcpy(&tmp_h,&hw[0],sizeof(long double));
- /*Don't compare the last 3 bytes of mantissa*/
- if(!HDmemcmp(&tmp_s+4, &tmp_h+4, sizeof(long double)-4))
- continue; /*no error*/
- }
-#endif /*end H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS*/
-
- /* Print errors */
- if (0==fails_this_test++)
- H5_FAILED();
- printf(" test %u elmt %u: \n", (unsigned)i+1, (unsigned)j);
-
- printf(" src = ");
- for (k=0; k<src_size; k++)
- printf(" %02x", saved[j*src_size+ENDIAN(src_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)dst_size-(ssize_t)src_size)), "");
- switch (src_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, saved+j*sizeof(signed char), sizeof(signed char));
- printf(" %29d\n", (int)*((signed char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, saved+j*sizeof(unsigned char), sizeof(unsigned char));
- printf(" %29u\n", (unsigned)*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, saved+j*sizeof(short), sizeof(short));
- printf(" %29hd\n", *((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned short), sizeof(unsigned short));
- printf(" %29hu\n", *((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, saved+j*sizeof(int), sizeof(int));
- printf(" %29d\n", *((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, saved+j*sizeof(unsigned), sizeof(unsigned));
- printf(" %29u\n", *((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, saved+j*sizeof(long), sizeof(long));
- printf(" %29ld\n", *((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long), sizeof(unsigned long));
- printf(" %29lu\n", *((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, saved+j*sizeof(long_long), sizeof(long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, saved+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)aligned));
- break;
- case FLT_FLOAT:
- HDmemcpy(aligned, saved+j*sizeof(float), sizeof(float));
- printf(" %29f\n", *((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(double), sizeof(double));
- printf(" %29f\n", *((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, saved+j*sizeof(long double), sizeof(long double));
- printf(" %29Lf\n", *((long double*)aligned));
- break;
- case OTHER:
- break;
- }
-
- printf(" dst = ");
- for (k=0; k<dst_size; k++)
- printf(" %02x", buf[j*dst_size+ENDIAN(dst_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- switch (dst_type) {
- case INT_SCHAR:
- HDmemcpy(aligned, buf+j*sizeof(signed char), sizeof(signed char));
- printf(" %29d\n", (int)*((signed char*)aligned));
- break;
- case INT_UCHAR:
- HDmemcpy(aligned, buf+j*sizeof(unsigned char), sizeof(unsigned char));
- printf(" %29u\n", (unsigned)*((unsigned char*)aligned));
- break;
- case INT_SHORT:
- HDmemcpy(aligned, buf+j*sizeof(short), sizeof(short));
- printf(" %29hd\n", *((short*)aligned));
- break;
- case INT_USHORT:
- HDmemcpy(aligned, buf+j*sizeof(unsigned short), sizeof(unsigned short));
- printf(" %29hu\n", *((unsigned short*)aligned));
- break;
- case INT_INT:
- HDmemcpy(aligned, buf+j*sizeof(int), sizeof(int));
- printf(" %29d\n", *((int*)aligned));
- break;
- case INT_UINT:
- HDmemcpy(aligned, buf+j*sizeof(unsigned), sizeof(unsigned));
- printf(" %29u\n", *((unsigned*)aligned));
- break;
- case INT_LONG:
- HDmemcpy(aligned, buf+j*sizeof(long), sizeof(long));
- printf(" %29ld\n", *((long*)aligned));
- break;
- case INT_ULONG:
- HDmemcpy(aligned, buf+j*sizeof(unsigned long), sizeof(unsigned long));
- printf(" %29lu\n", *((unsigned long*)aligned));
- break;
- case INT_LLONG:
- HDmemcpy(aligned, buf+j*sizeof(long_long), sizeof(long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)aligned));
- break;
- case INT_ULLONG:
- HDmemcpy(aligned, buf+j*sizeof(unsigned long_long), sizeof(unsigned long_long));
- HDfprintf(stdout," %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)aligned));
- break;
- case FLT_FLOAT:
- HDmemcpy(aligned, buf+j*sizeof(float), sizeof(float));
- printf(" %29f\n", *((float*)aligned));
- break;
- case FLT_DOUBLE:
- HDmemcpy(aligned, buf+j*sizeof(double), sizeof(double));
- printf(" %29f\n", *((double*)aligned));
- break;
- case FLT_LDOUBLE:
- HDmemcpy(aligned, buf+j*sizeof(long double), sizeof(long double));
- printf(" %29Lf\n", *((long double*)aligned));
- break;
- case OTHER:
- break;
- }
-
- printf(" ans = ");
- for (k=0; k<dst_size; k++)
- printf(" %02x", hw[ENDIAN(dst_size, k)]);
- printf("%*s", (int)(3*MAX(0, (ssize_t)src_size-(ssize_t)dst_size)), "");
- switch (dst_type) {
- case INT_SCHAR:
- printf(" %29d\n", (int)*((signed char*)hw));
- break;
- case INT_UCHAR:
- printf(" %29u\n", (unsigned)*((unsigned char*)hw));
- break;
- case INT_SHORT:
- printf(" %29hd\n", *((short*)hw));
- break;
- case INT_USHORT:
- printf(" %29hu\n", *((unsigned short*)hw));
- break;
- case INT_INT:
- printf(" %29d\n", *((int*)hw));
- break;
- case INT_UINT:
- printf(" %29u\n", *((unsigned int*)hw));
- break;
- case INT_LONG:
- printf(" %29ld\n", *((long*)hw));
- break;
- case INT_ULONG:
- printf(" %29lu\n", *((unsigned long*)hw));
- break;
- case INT_LLONG:
- printf(" %29"H5_PRINTF_LL_WIDTH"d\n", *((long_long*)hw));
- break;
- case INT_ULLONG:
- printf(" %29"H5_PRINTF_LL_WIDTH"u\n", *((unsigned long_long*)hw));
- break;
- case FLT_FLOAT:
- printf(" %29f\n", *((float*)hw));
- break;
- case FLT_DOUBLE:
- printf(" %29f\n", *((double*)hw));
- break;
- case FLT_LDOUBLE:
- printf(" %29Lf\n", *((long double*)hw));
- break;
- case OTHER:
- break;
- }
-
- if (++fails_all_tests>=max_fails) {
- HDputs(" maximum failures reached, aborting test...");
- goto done;
- }
- }
- PASSED();
- }
-#ifdef SHOW_OVERFLOWS
- if (noverflows_g>0) {
- printf(" %d overflow%s in previous test\n",
- noverflows_g, 1==noverflows_g?"":"s");
- }
-#endif
-
- done:
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
- reset_hdf5(); /*print statistics*/
- return (int)fails_all_tests;
-
- error:
- if (buf) aligned_free(buf);
- if (saved) aligned_free(saved);
- if (aligned) HDfree(aligned);
- HDfflush(stdout);
- reset_hdf5(); /*print statistics*/
- return MAX((int)fails_all_tests, 1);
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: overflows
- *
- * Purpose: When convert from float or double to any integer type,
- * check if overflow occurs.
- *
- *
- * Return: TRUE: overflow happens
- *
- * FALSE: no overflow
- *
- * Programmer: Raymond Lu
- * Monday, Nov 17, 2003
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static hbool_t
-overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits)
-{
- hbool_t ret_value=FALSE;
- hsize_t expt;
- size_t mant_digits=0, expt_digits=0, bias=0;
- size_t epos, mpos;
- size_t src_prec=0; /*source type precision in bits*/
- H5T_norm_t norm;
- ssize_t indx;
- unsigned char bits[32], mant_bits[32];
-
- HDmemset(bits, 0, 32);
- HDmemset(mant_bits, 0, 32);
-
- /*
- * Sometimes, type size isn't equal to the precision like Linux's "long
- * double", where size is 96 bits and precision is 80 bits.
- */
-
- src_prec = H5Tget_precision(src_id);
- H5Tget_fields(src_id, NULL, &epos, &expt_digits, &mpos, &mant_digits);
- bias = H5Tget_ebias(src_id);
- norm = H5Tget_norm(src_id);
-
- HDmemcpy(bits, origin_bits, src_prec/8+1);
-
- /*Check for special cases: +Inf, -Inf*/
- if (H5T_bit_find (bits, mpos, mant_digits, H5T_BIT_LSB, TRUE)<0) {
- if (H5T_bit_find (bits, epos, expt_digits, H5T_BIT_LSB, FALSE)<0) {
- ret_value=TRUE;
- goto done;
- }
- } else if (H5T_NORM_NONE==norm && H5T_bit_find (bits, mpos, mant_digits-1,
- H5T_BIT_LSB, TRUE)<0 && H5T_bit_find (bits, epos, expt_digits,
- H5T_BIT_LSB, FALSE)<0) {
- /*This is a special case for the source of no implied mantissa bit.
- *If the exponent bits are all 1s and only the 1st bit of mantissa
- *is set to 1. It's infinity. The Intel-Linux "long double" is this case.*/
- ret_value=TRUE;
- goto done;
- }
-
- /* get exponent */
- expt = H5T_bit_get_d(bits, mant_digits, expt_digits) - bias;
-
- if(expt>=(dst_num_bits-1)) {
- ret_value=TRUE;
- goto done;
- }
-
- /* get significand */
- H5T_bit_copy (mant_bits, 0, bits, 0, mant_digits);
-
-
- /* restore implicit bit if normalization is implied*/
- if(norm == H5T_NORM_IMPLIED) {
- H5T_bit_inc(mant_bits, mant_digits, 1);
- mant_digits++;
- }
-
- /* shift significand */
- H5T_bit_shift (mant_bits, (expt-expt_digits), 0, 32*8);
-
- indx = H5T_bit_find(mant_bits, 0, 32*8, H5T_BIT_MSB, 1);
-
- if((size_t)indx>=dst_num_bits)
- ret_value=TRUE;
-
-done:
- return ret_value;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: run_integer_tests
- *
- * Purpose: Runs all integer tests.
- *
- * Return: Number of errors
- *
- * Programmer: Robb Matzke
- * Tuesday, November 24, 1998
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-run_integer_tests(const char *name)
-{
- int nerrors = 0;
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SCHAR, H5T_NATIVE_ULLONG);
-#endif
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UCHAR, H5T_NATIVE_ULLONG);
-#endif
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_SHORT, H5T_NATIVE_ULLONG);
-#endif
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_USHORT, H5T_NATIVE_ULLONG);
-#endif
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_INT, H5T_NATIVE_ULLONG);
-#endif
-
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_INT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_ULONG);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_UINT, H5T_NATIVE_ULLONG);
-#endif
-
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_UINT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_ULONG);
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LONG, H5T_NATIVE_ULLONG);
-#endif
-#endif
-
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_UINT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_LONG);
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULONG, H5T_NATIVE_ULLONG);
-#endif
-#endif
-
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_ULONG);
-#endif
- nerrors += test_conv_int_1(name, H5T_NATIVE_LLONG, H5T_NATIVE_ULLONG);
-#endif
-
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_INT);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_UINT);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_LONG);
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_ULONG);
-#endif
- nerrors += test_conv_int_1(name, H5T_NATIVE_ULLONG, H5T_NATIVE_LLONG);
-#endif
-
- return nerrors;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: run_fp_tests
- *
- * Purpose: Runs all floating-point tests.
- *
- * Return: Number of errors
- *
- * Programmer: Raymond Lu
- * Tuesday, March 22, 2005
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-run_fp_tests(const char *name)
-{
- int nerrors = 0;
-
- if(!strcmp(name, "noop")) {
- nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_FLOAT);
- nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE);
- nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LDOUBLE);
- goto done;
- }
-
- /*Test normalized and denormalized values. FALSE indicates non-special values.*/
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE);
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT);
- nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE);
-#endif
-
- /*Test special values, +/-0, +/-infinity, +/-QNaN, +/-SNaN.*/
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE);
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT);
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT);
- nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE);
-#endif
-
-done:
- return nerrors;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: run_int_fp_conv
- *
- * Purpose: Runs all integer-float tests.
- *
- * Return: Number of errors
- *
- * Programmer: Raymond Lu
- * Monday, November 10, 2003
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-run_int_fp_conv(const char *name)
-{
- int nerrors = 0;
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_DOUBLE);
-
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_DOUBLE);
-
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_DOUBLE);
-#endif
-
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_DOUBLE);
-
-#ifdef H5_ULLONG_TO_FP_CAST_WORKS
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_FLOAT);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_DOUBLE);
-#else /* H5_ULLONG_TO_FP_CAST_WORKS */
- {
- char str[256]; /*hello string */
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "unsigned long long", "float");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to compiler not handling conversion.");
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "unsigned long long", "double");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to compiler not handling conversion.");
- }
-#endif /* H5_ULLONG_TO_FP_CAST_WORKS */
-#endif
-
-#if H5_SW_INTEGER_TO_LDOUBLE_WORKS
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_LDOUBLE);
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_LDOUBLE);
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_LDOUBLE);
-#endif
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_LDOUBLE);
-#if H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS
- nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_LDOUBLE);
-#else /* H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS */
- {
- char str[256]; /*hello string */
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "unsigned long long", "long double");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to compiler not handling conversion.");
- }
-#endif /* H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS */
-#endif
-#endif
-#else /*H5_SW_INTEGER_TO_LDOUBLE_WORKS*/
- {
- char str[256]; /*string */
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "all integers", "long double");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
- }
-#endif /*H5_SW_INTEGER_TO_LDOUBLE_WORKS*/
-
- return nerrors;
-}
-
-
-/*-------------------------------------------------------------------------
- * Function: run_fp_int_conv
- *
- * Purpose: Runs all float-integer tests.
- *
- * Return: Number of errors
- *
- * Programmer: Raymond Lu
- * Monday, November 10, 2003
- *
- * Modifications:
- *
- *-------------------------------------------------------------------------
- */
-static int
-run_fp_int_conv(const char *name)
-{
- int nerrors = 0;
- hbool_t special;
- int i;
-
- for(i=0; i<2; i++) {
- if(i==0)
- special = FALSE;
- else
- special = TRUE;
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_SCHAR);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_UCHAR);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_SHORT);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_USHORT);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_INT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_INT);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_UINT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_UINT);
-
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LONG);
-
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_ULONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULONG);
-#endif
-
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- if(!strcmp(name, "hw")) { /* Hardware conversion */
- /* Windows .NET 2003 doesn't work for hardware conversion of this case.
- * .NET should define this macro H5_HW_FP_TO_LLONG_NOT_WORKS. */
-#ifndef H5_HW_FP_TO_LLONG_NOT_WORKS
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG);
-#endif /*H5_HW_FP_TO_LLONG_NOT_WORKS*/
- } else { /* Software conversion */
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG);
- }
-#ifdef H5_FP_TO_ULLONG_RIGHT_MAXIMUM
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_ULLONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULLONG);
-#else /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/
- {
- char str[256]; /*hello string */
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "float", "unsigned long long");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
-
- sprintf(str, "Testing random %s %s -> %s conversions",
- name, "double", "unsigned long long");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
- }
-#endif /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/
-#endif
-
-#if H5_SW_LDOUBLE_TO_INTEGER_WORKS
-#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SCHAR);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UCHAR);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SHORT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_USHORT);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_INT);
-#if H5_CV_LDOUBLE_TO_UINT_WORKS
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UINT);
-#else /*H5_CV_LDOUBLE_TO_UINT_WORKS*/
- {
- char str[256]; /*string */
-
- sprintf(str, "Testing %s %s -> %s conversions",
- name, "long double", "unsigned int");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
- }
-#endif /*H5_CV_LDOUBLE_TO_UINT_WORKS*/
-#if H5_SIZEOF_LONG!=H5_SIZEOF_INT
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LONG);
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULONG);
-#endif
-
-#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LLONG);
-#ifdef H5_FP_TO_ULLONG_RIGHT_MAXIMUM
- nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULLONG);
-#else /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/
- {
- char str[256]; /*string */
-
- sprintf(str, "Testing %s %s -> %s conversions",
- name, "long double", "unsigned long long");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
- }
-#endif /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/
-#endif
-#endif
-#else /*H5_SW_LDOUBLE_TO_INTEGER_WORKS*/
- {
- char str[256]; /*hello string */
-
- sprintf(str, "Testing %s %s -> %s conversions",
- name, "long double", "all integers");
- printf("%-70s", str);
- SKIPPED();
- HDputs(" Test skipped due to hardware conversion error.");
- }
-#endif /*H5_SW_LDOUBLE_TO_INTEGER_WORKS*/
- }
-
- return nerrors;
-}
-
-
-/*-------------------------------------------------------------------------
* Function: test_encode
*
* Purpose: Tests functions of encoding and decoding data type.
@@ -8379,60 +3811,11 @@ main(void)
nerrors += test_compound_10();
nerrors += test_compound_11();
nerrors += test_compound_12();
- nerrors += test_conv_int ();
nerrors += test_conv_enum_1();
nerrors += test_conv_enum_2();
nerrors += test_conv_bitfield();
nerrors += test_opaque();
- /* Test user-define, query functions and software conversion
- * for user-defined floating-point types */
- nerrors += test_derived_flt();
-
- /* Test user-define, query functions and software conversion
- * for user-defined integer types */
- nerrors += test_derived_integer();
-
- /* Does floating point overflow generate a SIGFPE? */
- generates_sigfpe();
-
- /* Test degenerate cases */
- nerrors += run_fp_tests("noop");
-
- /* Test hardware floating-point conversion functions */
- nerrors += run_fp_tests("hw");
-
- /* Test hardware integer conversion functions */
- nerrors += run_integer_tests("hw");
-
- /* Test hardware integer-float conversion functions */
- nerrors += run_int_fp_conv("hw");
-
- /* Test hardware float-integer conversion functions */
- nerrors += run_fp_int_conv("hw");
-
- /*----------------------------------------------------------------------
- * Software tests
- *----------------------------------------------------------------------
- */
- without_hardware_g = TRUE;
- reset_hdf5();
-
- /* Test software floating-point conversion functions */
- nerrors += run_fp_tests("sw");
-
- /* Test software integer conversion functions */
- nerrors += test_conv_int_2();
- nerrors += run_integer_tests("sw");
-
- /* Test software float-integer conversion functions */
- nerrors += run_fp_int_conv("sw");
-
- /* Test software integer-float conversion functions */
- nerrors += run_int_fp_conv("sw");
-
- reset_hdf5();
-
if (nerrors) {
printf("***** %lu FAILURE%s! *****\n",
nerrors, 1==nerrors?"":"S");