/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Copyright by The HDF Group. * * All rights reserved. * * * * This file is part of HDF5. The full HDF5 copyright notice, including * * terms governing use, modification, and redistribution, is contained in * * the COPYING file, which can be found at the root of the source code * * distribution tree, or in https://www.hdfgroup.org/licenses. * * If you do not have access to either file, you may request a copy from * * help@hdfgroup.org. * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ /* * Purpose: Tests the data type interface (H5T) */ #include "h5test.h" /* Number of elements in each random test */ #define NTESTELEM 10000 /* Epsilon for floating-point comparisons */ #define FP_EPSILON 0.000001F /* * Offset from aligned memory returned by malloc(). This can be used to test * that type conversions handle non-aligned buffers correctly. */ #define ALIGNMENT 1 /* * Define if you want to test alignment code on a machine that doesn't * normally require alignment. When set, all native data types must be aligned * on a byte boundary equal to the data size. */ #define TEST_ALIGNMENT /* Alignment test stuff */ #ifdef TEST_ALIGNMENT #define H5T_FRIEND /*suppress error about including H5Tpkg */ #include "H5Tpkg.h" #endif #define SET_ALIGNMENT(TYPE, VAL) H5T_NATIVE_##TYPE##_ALIGN_g = MAX(H5T_NATIVE_##TYPE##_ALIGN_g, VAL) static const char *FILENAME[] = {"dt_arith1", "dt_arith2", NULL}; /* * Count up or down depending on whether the machine is big endian or little * endian. If local variable `endian' is H5T_ORDER_BE then the result will * be I, otherwise the result will be Z-(I+1). */ #define ENDIAN(Z, I, E) (H5T_ORDER_BE == (E) ? (I) : (Z) - ((I) + 1)) typedef enum dtype_t { INT_SCHAR, INT_UCHAR, INT_SHORT, INT_USHORT, INT_INT, INT_UINT, INT_LONG, INT_ULONG, INT_LLONG, INT_ULLONG, FLT_FLOAT, FLT_DOUBLE, FLT_LDOUBLE, OTHER } dtype_t; /* * 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 /* * Decide what values of floating-point number we want to test. They are * 1 - normalized; 2 - denormalized; 3 - special. */ #define TEST_NOOP 0 #define TEST_NORMAL 1 #define TEST_DENORM 2 #define TEST_SPECIAL 3 /* Temporary buffer sizes */ #define TMP_BUF_DIM1 32 #define TMP_BUF_DIM2 100 /* Don't use hardware conversions if set */ static int without_hardware_g = 0; /* Allocates memory aligned on a certain boundary. */ #define aligned_malloc(Z) ((void *)((char *)malloc(ALIGNMENT + Z) + ALIGNMENT)) #define aligned_free(M) free((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) \ do { \ unsigned char *buf_p, *saved_p; \ unsigned int n; \ TYPE value1 = 1; \ TYPE value2 = 0; \ \ /* Allocate buffers */ \ NELMTS = SRC_PREC * 3; \ BUF = (unsigned char *)aligned_malloc(NELMTS * MAX(SRC_SIZE, DST_SIZE)); \ SAVED = (unsigned char *)aligned_malloc(NELMTS * MAX(SRC_SIZE, DST_SIZE)); \ memset(BUF, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ memset(SAVED, 0, 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; n++) { \ { \ memcpy(buf_p, &value1, SRC_SIZE); \ memcpy(saved_p, &value1, SRC_SIZE); \ buf_p += SRC_SIZE; \ saved_p += SRC_SIZE; \ } \ { \ memcpy(buf_p, &value2, SRC_SIZE); \ memcpy(saved_p, &value2, SRC_SIZE); \ buf_p += SRC_SIZE; \ saved_p += SRC_SIZE; \ } \ \ if (n < SRC_PREC - 2) { \ value1 = (TYPE)(value1 << 1); \ value2 = (TYPE)((value1 - 1) | value1); \ } \ else if (n == SRC_PREC - 2) { /*to avoid overflow of negative values for signed integer*/ \ value1 = (TYPE)(value1 << 1); \ value2 = (TYPE)((~value1) | 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 - 1; n++) { \ { \ memcpy(buf_p, &value2, SRC_SIZE); \ memcpy(saved_p, &value2, SRC_SIZE); \ buf_p += SRC_SIZE; \ saved_p += SRC_SIZE; \ } \ if (n < SRC_PREC - 1) \ value2 = (TYPE)(value2 << 1); \ } \ } while (0) /* 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) \ do { \ unsigned int m; \ if (H5T_ORDER_BE == EORDER) { \ 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; \ } \ } \ else if (H5T_ORDER_VAX == EORDER) { \ unsigned char mediator1, mediator2; \ for (m = 0; m < ESIZE; m += 4) { \ mediator1 = EBUF[m]; \ mediator2 = EBUF[m + 1]; \ \ EBUF[m] = EBUF[(ESIZE - 2) - m]; \ EBUF[m + 1] = EBUF[(ESIZE - 1) - m]; \ \ EBUF[(ESIZE - 2) - m] = mediator1; \ EBUF[(ESIZE - 1) - m] = mediator2; \ } \ } \ } while (0) /* Allocate buffer and initialize it with floating-point normalized values. * It's for conversion test of floating-point as the source. */ #define INIT_FP_NORM(TYPE, SRC_MAX, SRC_MIN, SRC_MAX_10_EXP, SRC_MIN_10_EXP, SRC_SIZE, DST_SIZE, BUF, SAVED, \ NELMTS) \ do { \ unsigned char *buf_p, *saved_p; \ 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 arithmetic 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*/ \ \ /* Allocate buffers */ \ BUF = (unsigned char *)aligned_malloc(NELMTS * MAX(SRC_SIZE, DST_SIZE)); \ SAVED = (unsigned char *)aligned_malloc(NELMTS * MAX(SRC_SIZE, DST_SIZE)); \ memset(BUF, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ memset(SAVED, 0, NELMTS *MAX(SRC_SIZE, DST_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; \ } while (0) /* Allocate buffer and initialize it with floating-point denormalized values. * It's for conversion test of floating-point as the source. */ #define INIT_FP_DENORM(TYPE, SRC_MANT_DIG, SRC_SIZE, SRC_PREC, SRC_ORDR, DST_SIZE, BUF, SAVED, NELMTS) \ do { \ unsigned char *buf_p, *saved_p; \ unsigned char *tmp1, *tmp2; \ size_t n; \ \ /*Total number of values*/ \ NELMTS = 2 * /*both positive and negative*/ \ (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)); \ memset(BUF, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ memset(SAVED, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ \ tmp1 = (unsigned char *)calloc((size_t)1, (size_t)SRC_SIZE); \ tmp2 = (unsigned char *)calloc((size_t)1, (size_t)SRC_SIZE); \ \ buf_p = BUF; \ saved_p = SAVED; \ \ /*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, (size_t)1, true); /*the negative value*/ \ for (n = 0; n < SRC_MANT_DIG - 1; n++) { \ H5T__bit_set(tmp1, n, (size_t)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, (size_t)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); \ } while (0) /* 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) \ do { \ unsigned char *buf_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)); \ memset(BUF, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ memset(SAVED, 0, NELMTS *MAX(SRC_SIZE, DST_SIZE)); \ value = (unsigned char *)calloc(SRC_SIZE, sizeof(unsigned char)); \ \ buf_p = BUF; \ \ /* +0 */ \ H5T__bit_set(value, (size_t)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, (size_t)(SRC_PREC - 1), (size_t)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, (size_t)(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, (size_t)0, (size_t)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, (size_t)(SRC_MANT_DIG - 2), (size_t)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); \ } while (0) static bool 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 my_isinf(int endian, const unsigned char *val, size_t size, size_t mpos, size_t msize, size_t epos, size_t esize); /*------------------------------------------------------------------------- * Function: fpe_handler * * Purpose: Exit with 255 * * Return: void * *------------------------------------------------------------------------- */ static void fpe_handler(int H5_ATTR_UNUSED signo) { SKIPPED(); puts(" Test skipped due to SIGFPE."); #ifndef HANDLE_SIGFPE puts(" Remaining tests could not be run."); puts(" Please turn off SIGFPE on overflows and try again."); #endif exit(255); } /*------------------------------------------------------------------------- * Function: reset_hdf5 * * Purpose: Reset the hdf5 library. This causes statistics to be printed * and counters to be reset. * * Return: void * *------------------------------------------------------------------------- */ static void reset_hdf5(void) { h5_reset(); if (without_hardware_g) h5_no_hwconv(); #ifdef TEST_ALIGNMENT SET_ALIGNMENT(SCHAR, H5_SIZEOF_CHAR); SET_ALIGNMENT(UCHAR, H5_SIZEOF_CHAR); SET_ALIGNMENT(SHORT, H5_SIZEOF_SHORT); SET_ALIGNMENT(USHORT, H5_SIZEOF_SHORT); SET_ALIGNMENT(INT, H5_SIZEOF_INT); SET_ALIGNMENT(UINT, H5_SIZEOF_INT); SET_ALIGNMENT(LONG, H5_SIZEOF_LONG); SET_ALIGNMENT(ULONG, H5_SIZEOF_LONG); SET_ALIGNMENT(LLONG, H5_SIZEOF_LONG_LONG); SET_ALIGNMENT(ULLONG, H5_SIZEOF_LONG_LONG); SET_ALIGNMENT(FLOAT, H5_SIZEOF_FLOAT); SET_ALIGNMENT(DOUBLE, H5_SIZEOF_DOUBLE); SET_ALIGNMENT(LDOUBLE, H5_SIZEOF_LONG_DOUBLE); #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 * *------------------------------------------------------------------------- */ static H5T_conv_ret_t except_func(H5T_conv_except_t except_type, hid_t H5_ATTR_UNUSED src_id, hid_t H5_ATTR_UNUSED dst_id, void H5_ATTR_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: test_hard_query * * Purpose: Tests H5Tcompiler_conv() for querying whether a conversion is * a hard one. * * Return: Success: 0 * * Failure: number of errors * *------------------------------------------------------------------------- */ static int test_hard_query(void) { TESTING("query functions of compiler conversion"); /* Verify the conversion from int to float is a hard conversion. */ if (H5Tcompiler_conv(H5T_NATIVE_INT, H5T_NATIVE_FLOAT) != true) { H5_FAILED(); printf("Can't query conversion function\n"); goto error; } /* Unregister the hard conversion from int to float. Verify the conversion * is a soft conversion. */ H5Tunregister(H5T_PERS_HARD, NULL, H5T_NATIVE_INT, H5T_NATIVE_FLOAT, (H5T_conv_t)((void (*)(void))H5T__conv_int_float)); if (H5Tcompiler_conv(H5T_NATIVE_INT, H5T_NATIVE_FLOAT) != false) { H5_FAILED(); printf("Can't query conversion function\n"); goto error; } /* Register the hard conversion from int to float. Verify the conversion * is a hard conversion. */ H5Tregister(H5T_PERS_HARD, "int_flt", H5T_NATIVE_INT, H5T_NATIVE_FLOAT, (H5T_conv_t)((void (*)(void))H5T__conv_int_float)); if (H5Tcompiler_conv(H5T_NATIVE_INT, H5T_NATIVE_FLOAT) != true) { H5_FAILED(); printf("Can't query conversion function\n"); goto error; } PASSED(); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); return 0; error: /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); return 1; } /*------------------------------------------------------------------------- * Function: expt_handle * * Purpose: Gets called from test_particular_fp_integer() for data type * conversion exceptions. * * Return: H5T_CONV_HANDLED 1 * *------------------------------------------------------------------------- */ static H5T_conv_ret_t expt_handle(H5T_conv_except_t except_type, hid_t H5_ATTR_UNUSED src_id, hid_t H5_ATTR_UNUSED dst_id, void H5_ATTR_UNUSED *src_buf, void *dst_buf, void *user_data) { signed char fill_value1 = 7; int fill_value2 = 13; if (except_type == H5T_CONV_EXCEPT_RANGE_HI || except_type == H5T_CONV_EXCEPT_RANGE_LOW || except_type == H5T_CONV_EXCEPT_TRUNCATE) { if (*(bool *)user_data) *(signed char *)dst_buf = fill_value1; else *(int *)dst_buf = fill_value2; } /* end if */ return H5T_CONV_HANDLED; } /*------------------------------------------------------------------------- * Function: test_particular_fp_integer * * Purpose: Tests hard conversions from floating numbers to integers in * a special situation when the source is "float" and assigned * the value of "INT_MAX". A compiler may do roundup making * this value "INT_MAX+1". When this float value is casted to * int, overflow happens. This test makes sure the library * returns exception in this situation. * * Also verifies the library handles conversion from double to * signed char correctly when the value of double is SCHAR_MAX. * The test makes sure the signed char doesn't overflow. * * This test is mainly for netCDF's request. * * Return: Success: 0 * * Failure: number of errors *------------------------------------------------------------------------- */ static int test_particular_fp_integer(void) { hid_t dxpl_id; bool flag; double src_d = (double)SCHAR_MAX; signed char dst_c; unsigned char *buf1 = NULL, *buf2 = NULL; unsigned char *saved_buf1 = NULL, *saved_buf2 = NULL; size_t src_size1, src_size2; size_t dst_size1, dst_size2; float src_f = (float)INT_MAX; int dst_i; int fill_value = 13; int endian; /*endianness */ unsigned int fails_this_test = 0; size_t j; TESTING("hard particular floating number -> integer conversions"); if ((dxpl_id = H5Pcreate(H5P_DATASET_XFER)) < 0) { H5_FAILED(); printf("Can't create data transfer property list\n"); goto error; } /* Test conversion from double (the value is SCHAR_MAX) to signed char. */ endian = H5Tget_order(H5T_NATIVE_DOUBLE); src_size1 = H5Tget_size(H5T_NATIVE_DOUBLE); dst_size1 = H5Tget_size(H5T_NATIVE_SCHAR); buf1 = (unsigned char *)calloc((size_t)1, (size_t)MAX(src_size1, dst_size1)); saved_buf1 = (unsigned char *)calloc((size_t)1, (size_t)MAX(src_size1, dst_size1)); memcpy(buf1, &src_d, src_size1); memcpy(saved_buf1, &src_d, src_size1); /* Register exception handling function and signal the destination is "signed char". */ flag = 1; if (H5Pset_type_conv_cb(dxpl_id, expt_handle, &flag) < 0) { H5_FAILED(); printf("Can't register conversion callback\n"); goto error; } /* Do conversion */ if (H5Tconvert(H5T_NATIVE_DOUBLE, H5T_NATIVE_SCHAR, (size_t)1, buf1, NULL, dxpl_id) < 0) { H5_FAILED(); printf("Can't convert data\n"); goto error; } memcpy(&dst_c, buf1, dst_size1); /* Print errors */ if (dst_c != SCHAR_MAX) { double x = 0.0; signed char y; if (0 == fails_this_test++) H5_FAILED(); printf(" test double to signed char:\n"); printf(" src = "); for (j = 0; j < src_size1; j++) printf(" %02x", saved_buf1[ENDIAN(src_size1, j, endian)]); memcpy(&x, saved_buf1, src_size1); printf(" %29.20e\n", x); printf(" dst = "); for (j = 0; j < dst_size1; j++) printf(" %02x", buf1[ENDIAN(dst_size1, j, endian)]); memcpy(&y, buf1, dst_size1); printf(" %29d\n", y); } /* Test conversion from float (the value is INT_MAX) to int. */ src_size2 = H5Tget_size(H5T_NATIVE_FLOAT); dst_size2 = H5Tget_size(H5T_NATIVE_INT); buf2 = (unsigned char *)calloc((size_t)1, (size_t)MAX(src_size2, dst_size2)); saved_buf2 = (unsigned char *)calloc((size_t)1, (size_t)MAX(src_size2, dst_size2)); memcpy(buf2, &src_f, src_size2); memcpy(saved_buf2, &src_f, src_size2); /* signal exception handling function that the destination is "int". */ flag = 0; /* Do conversion */ if (H5Tconvert(H5T_NATIVE_FLOAT, H5T_NATIVE_INT, (size_t)1, buf2, NULL, dxpl_id) < 0) { H5_FAILED(); printf("Can't convert data\n"); goto error; } memcpy(&dst_i, buf2, dst_size2); /* Print errors */ if (dst_i != fill_value) { float x = 0.0F; int y; if (0 == fails_this_test++) H5_FAILED(); printf(" test float to int:\n"); printf(" src = "); for (j = 0; j < src_size2; j++) printf(" %02x", saved_buf2[ENDIAN(src_size2, j, endian)]); memcpy(&x, saved_buf2, src_size2); printf(" %29.20e\n", (double)x); printf(" dst = "); for (j = 0; j < dst_size2; j++) printf(" %02x", buf2[ENDIAN(dst_size2, j, endian)]); memcpy(&y, buf2, dst_size2); printf(" %29d\n", y); } if (fails_this_test) goto error; if (H5Pclose(dxpl_id) < 0) { H5_FAILED(); printf("Can't close property list\n"); goto error; } if (buf1) free(buf1); if (buf2) free(buf2); if (saved_buf1) free(saved_buf1); if (saved_buf2) free(saved_buf2); PASSED(); return 0; error: fflush(stdout); H5E_BEGIN_TRY { H5Pclose(dxpl_id); } H5E_END_TRY if (buf1) free(buf1); if (buf2) free(buf2); if (saved_buf1) free(saved_buf1); if (saved_buf2) free(saved_buf2); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ return MAX((int)fails_this_test, 1); } /*------------------------------------------------------------------------- * Function: test_derived_flt * * Purpose: Tests user-define and query functions of floating-point types. * * Return: Success: 0 * * Failure: number of errors * *------------------------------------------------------------------------- */ static int test_derived_flt(void) { hid_t file = H5I_INVALID_HID, tid1 = H5I_INVALID_HID, tid2 = H5I_INVALID_HID; hid_t dxpl_id = H5I_INVALID_HID; char filename[1024]; size_t spos, epos, esize, mpos, msize, size; size_t src_size, dst_size; unsigned char *buf = NULL, *saved_buf = NULL; int *aligned = NULL; int endian; /*endianness */ 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[0], 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, (size_t)44, (size_t)34, (size_t)10, (size_t)3, (size_t)31) < 0) { H5_FAILED(); printf("Can't set fields\n"); goto error; } if (H5Tset_offset(tid1, (size_t)3) < 0) { H5_FAILED(); printf("Can't set offset\n"); goto error; } if (H5Tset_precision(tid1, (size_t)42) < 0) { H5_FAILED(); printf("Can't set precision 1\n"); goto error; } if (H5Tset_size(tid1, (size_t)7) < 0) { H5_FAILED(); printf("Can't set size\n"); goto error; } if (H5Tset_ebias(tid1, (size_t)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 (H5Tcommit2(file, "new float type 1", tid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 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 = H5Topen2(file, "new float type 1", H5P_DEFAULT)) < 0) FAIL_PUTS_ERROR("Can't open datatype"); 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 (H5Tget_precision(tid1) != 42) { H5_FAILED(); printf("Can't get precision or wrong precision\n"); goto error; } if (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 (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); memset(buf, 0, nelmts * MAX(src_size, size)); memset(saved_buf, 0, nelmts * src_size); aligned = (int *)calloc((size_t)1, src_size); for (i = 0; i < nelmts * src_size; i++) buf[i] = saved_buf[i] = (unsigned char)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++) { snprintf(str, sizeof(str), "\nTesting random sw derived floating-point -> derived floating-point conversions"); printf("%-70s", str); fflush(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, endian)]); memcpy(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, endian)]); memcpy(aligned, buf + i * sizeof(int), sizeof(int)); printf(" %29d\n", *aligned); if (fails_this_test >= max_fails) { puts(" maximum failures reached, aborting test..."); goto error; } } fails_this_test = 0; free(buf); free(saved_buf); 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, (size_t)23, (size_t)16, (size_t)7, (size_t)0, (size_t)16) < 0) { H5_FAILED(); printf("Can't set fields\n"); goto error; } if (H5Tset_offset(tid2, (size_t)0) < 0) { H5_FAILED(); printf("Can't set offset\n"); goto error; } if (H5Tset_precision(tid2, (size_t)24) < 0) { H5_FAILED(); printf("Can't set precision 2\n"); goto error; } if (H5Tset_size(tid2, (size_t)3) < 0) { H5_FAILED(); printf("Can't set size\n"); goto error; } if (H5Tset_ebias(tid2, (size_t)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 (H5Tcommit2(file, "new float type 2", tid2, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 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 = H5Topen2(file, "new float type 2", H5P_DEFAULT)) < 0) FAIL_PUTS_ERROR("Can't open datatype"); 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 (H5Tget_precision(tid2) != 24) { H5_FAILED(); printf("Can't get precision or wrong precision\n"); goto error; } if (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 (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); memset(buf, 0, nelmts * MAX(src_size, dst_size)); memset(saved_buf, 0, nelmts * src_size); for (i = 0; i < nelmts * src_size; i++) buf[i] = saved_buf[i] = (unsigned char)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++) { snprintf(str, sizeof(str), "\nTesting random sw derived floating-point -> derived floating-point conversions"); printf("%-70s", str); fflush(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, endian)]); printf("\n"); printf(" dst = "); for (j = 0; j < src_size; j++) printf(" %02x", buf[i * src_size + ENDIAN(src_size, j, endian)]); printf("\n"); if (fails_this_test >= max_fails) { puts(" 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(); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ return 0; error: if (buf) free(buf); if (saved_buf) free(saved_buf); if (aligned) free(aligned); fflush(stdout); H5E_BEGIN_TRY { H5Tclose(tid1); H5Tclose(tid2); H5Pclose(dxpl_id); H5Fclose(file); } H5E_END_TRY /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); 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 * *------------------------------------------------------------------------- */ static int test_derived_integer(void) { hid_t file = H5I_INVALID_HID, tid1 = H5I_INVALID_HID, tid2 = H5I_INVALID_HID; hid_t dxpl_id = H5I_INVALID_HID; char filename[1024]; size_t src_size, dst_size; unsigned char *buf = NULL, *saved_buf = NULL; int endian; /*endianness */ 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[1], 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, (size_t)0) < 0) { H5_FAILED(); printf("Can't set offset\n"); goto error; } if (H5Tset_size(tid1, (size_t)3) < 0) { H5_FAILED(); printf("Can't set size\n"); goto error; } if (H5Tset_precision(tid1, (size_t)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 (H5Tcommit2(file, "new integer type 1", tid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 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 = H5Topen2(file, "new integer type 1", H5P_DEFAULT)) < 0) FAIL_PUTS_ERROR("Can't open datatype"); if (H5Tget_precision(tid1) != 24) { H5_FAILED(); printf("Can't get precision or wrong precision\n"); goto error; } if (H5Tget_offset(tid1) != 0) { H5_FAILED(); printf("Can't get offset or wrong offset\n"); goto error; } if (H5Tget_size(tid1) != 3) { H5_FAILED(); printf("Can't get size or wrong size\n"); goto error; } if (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, (size_t)48) < 0) { H5_FAILED(); printf("Can't set precision\n"); goto error; } if (H5Tset_offset(tid2, (size_t)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 (H5Tcommit2(file, "new integer type 2", tid2, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 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 = H5Topen2(file, "new integer type 2", H5P_DEFAULT)) < 0) FAIL_PUTS_ERROR("Can't open datatype"); if (H5Tget_precision(tid2) != 48) { H5_FAILED(); printf("Can't get precision or wrong precision\n"); goto error; } if (H5Tget_offset(tid2) != 10) { H5_FAILED(); printf("Can't get offset or wrong offset\n"); goto error; } if (H5Tget_size(tid2) != 8) { H5_FAILED(); printf("Can't get size or wrong size\n"); goto error; } if (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); memset(buf, 0, nelmts * MAX(src_size, dst_size)); memset(saved_buf, 0, nelmts * src_size); for (i = 0; i < nelmts * src_size; i++) buf[i] = saved_buf[i] = (unsigned char)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++) { snprintf(str, sizeof(str), "\nTesting random sw derived integer -> derived integer conversions"); printf("%-70s", str); fflush(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, endian)]); printf("\n"); printf(" dst = "); for (j = 0; j < src_size; j++) printf(" %02x", buf[i * src_size + ENDIAN(src_size, j, endian)]); printf("\n"); if (fails_this_test >= max_fails) { puts(" 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 */ free(buf); free(saved_buf); PASSED(); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ return 0; error: if (buf) free(buf); if (saved_buf) free(saved_buf); fflush(stdout); H5E_BEGIN_TRY { H5Tclose(tid1); H5Tclose(tid2); H5Pclose(dxpl_id); H5Fclose(file); } H5E_END_TRY /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ return MAX((int)fails_this_test, 1); } /*------------------------------------------------------------------------- * Function: test_conv_int_1 * * Purpose: Test conversion of 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 * *------------------------------------------------------------------------- */ static int test_conv_int_1(const char *name, hid_t src, hid_t dst) { 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 endianness */ size_t src_size, dst_size; /*type sizes */ unsigned char *buf = NULL; /*buffer for conversion */ unsigned char *saved = NULL; /*original values */ size_t 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) { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); H5_FAILED(); puts(" Unknown data type."); goto error; } else { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); fflush(stdout); fails_this_test = 0; } /* 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 = calloc((size_t)1, sizeof(long long)); /* 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; /* 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: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_char = (signed char)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_char = (signed char)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_char = (signed char)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_char = (signed char)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_char = (signed char)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_char = (signed char)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_char = (signed char)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_char = (signed char)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_char = (signed char)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_char = (signed char)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_UCHAR == dst_type) { hw = (unsigned char *)&hw_uchar; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_uchar = (unsigned char)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_uchar = (unsigned char)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_uchar = (unsigned char)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_uchar = (unsigned char)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_uchar = (unsigned char)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_uchar = (unsigned char)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_uchar = (unsigned char)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_uchar = (unsigned char)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_uchar = (unsigned char)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_uchar = (unsigned char)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_SHORT == dst_type) { hw = (unsigned char *)&hw_short; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(char), sizeof(char)); hw_short = (short)(*((char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_short = (short)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_short = (short)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_short = (short)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_short = (short)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_short = (short)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_short = (short)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_short = (short)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_short = (short)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_short = (short)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_USHORT == dst_type) { hw = (unsigned char *)&hw_ushort; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_ushort = (unsigned short)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_ushort = (unsigned short)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_ushort = (unsigned short)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_ushort = (unsigned short)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_ushort = (unsigned short)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_ushort = (unsigned short)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_ushort = (unsigned short)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_ushort = (unsigned short)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_ushort = (unsigned short)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_ushort = (unsigned short)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_INT == dst_type) { hw = (unsigned char *)&hw_int; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_int = (int)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_int = (int)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_int = (int)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_int = (int)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_int = (int)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_int = (int)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_int = (int)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_int = (int)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_int = (int)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_int = (int)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_UINT == dst_type) { hw = (unsigned char *)&hw_uint; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_uint = (unsigned int)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_uint = (unsigned int)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_uint = (unsigned int)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_uint = (unsigned int)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_uint = (unsigned int)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_uint = (unsigned int)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_uint = (unsigned int)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_uint = (unsigned int)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_uint = (unsigned int)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_uint = (unsigned int)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_LONG == dst_type) { hw = (unsigned char *)&hw_long; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_long = (long int)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_long = (long int)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_long = (long int)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_long = (long int)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_long = (long int)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_long = (long int)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_long = (long int)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_long = (long int)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_long = (long int)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_long = (long int)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_ULONG == dst_type) { hw = (unsigned char *)&hw_ulong; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_ulong = (unsigned long)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_ulong = (unsigned long)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_ulong = (unsigned long)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_ulong = (unsigned long)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_ulong = (unsigned long)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_ulong = (unsigned long)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_ulong = (unsigned long)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_ulong = (unsigned long)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_ulong = (unsigned long)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_ulong = (unsigned long)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_LLONG == dst_type) { hw = (unsigned char *)&hw_llong; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(char), sizeof(char)); hw_llong = (long long)(*((char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_llong = (long long)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_llong = (long long)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_llong = (long long)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_llong = (long long)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_llong = (long long)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_llong = (long long)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_llong = (long long)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_llong = (long long)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_llong = (long long)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_ULLONG == dst_type) { hw = (unsigned char *)&hw_ullong; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_ullong = (unsigned long long)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_ullong = (unsigned long long)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_ullong = (unsigned long long)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_ullong = (unsigned long long)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_ullong = (unsigned long long)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_ullong = (unsigned long long)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_ullong = (unsigned long long)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_ullong = (unsigned long long)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_ullong = (unsigned long long)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_ullong = (unsigned long long)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); 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, endian)]; for (k = 0; k < dst_size; k++) dst_bits[dst_size - (k + 1)] = buf[j * dst_size + ENDIAN(dst_size, k, endian)]; /* * 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, (size_t)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, (size_t)1) && H5T__bit_find(dst_bits, (size_t)0, dst_nbits - 1, H5T_BIT_LSB, 0) < 0) continue; /*no error*/ } else if (1 == H5T__bit_get_d(src_bits, src_nbits - 1, (size_t)1) && H5T__bit_find(src_bits, (size_t)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, (size_t)1) && H5T__bit_find(dst_bits, (size_t)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, (size_t)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, (size_t)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, (size_t)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, (size_t)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, (size_t)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, (size_t)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, (size_t)1) && H5T__bit_find(dst_bits, (size_t)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, (size_t)0, dst_nbits, H5T_BIT_LSB, 0) < 0) continue; /*no error*/ } } /* Print errors */ if (0 == fails_this_test++) H5_FAILED(); printf(" elmt %u\n", (unsigned)j); printf(" src = "); for (k = 0; k < src_size; k++) printf(" %02x", saved[j * src_size + ENDIAN(src_size, k, endian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)dst_size - (ssize_t)src_size)), ""); switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); printf(" %29d\n", (int)*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); printf(" %29u\n", (unsigned)*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); printf(" %29hd\n", *((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); printf(" %29hu\n", *((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); printf(" %29d\n", *((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); printf(" %29u\n", *((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); printf(" %29ld\n", *((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); printf(" %29lu\n", *((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); fprintf(stdout, " %29lld\n", *((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); fprintf(stdout, " %29llu\n", *((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } printf(" dst = "); for (k = 0; k < dst_size; k++) printf(" %02x", buf[j * dst_size + ENDIAN(dst_size, k, endian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)src_size - (ssize_t)dst_size)), ""); switch (dst_type) { case INT_SCHAR: memcpy(aligned, buf + j * sizeof(signed char), sizeof(signed char)); printf(" %29d\n", (int)*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, buf + j * sizeof(unsigned char), sizeof(unsigned char)); printf(" %29u\n", (unsigned)*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, buf + j * sizeof(short), sizeof(short)); printf(" %29hd\n", *((short *)aligned)); break; case INT_USHORT: memcpy(aligned, buf + j * sizeof(unsigned short), sizeof(unsigned short)); printf(" %29hu\n", *((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, buf + j * sizeof(int), sizeof(int)); printf(" %29d\n", *((int *)aligned)); break; case INT_UINT: memcpy(aligned, buf + j * sizeof(unsigned), sizeof(unsigned)); printf(" %29u\n", *((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, buf + j * sizeof(long), sizeof(long)); printf(" %29ld\n", *((long *)aligned)); break; case INT_ULONG: memcpy(aligned, buf + j * sizeof(unsigned long), sizeof(unsigned long)); printf(" %29lu\n", *((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, buf + j * sizeof(long long), sizeof(long long)); fprintf(stdout, " %29lld\n", *((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, buf + j * sizeof(long long), sizeof(unsigned long long)); fprintf(stdout, " %29llu\n", *((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } printf(" ans = "); for (k = 0; k < dst_size; k++) printf(" %02x", hw[ENDIAN(dst_size, k, endian)]); 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 *)((void *)hw))); break; case INT_UCHAR: printf(" %29u\n", (unsigned)*((unsigned char *)((void *)hw))); break; case INT_SHORT: printf(" %29hd\n", *((short *)((void *)hw))); break; case INT_USHORT: printf(" %29hu\n", *((unsigned short *)((void *)hw))); break; case INT_INT: printf(" %29d\n", *((int *)((void *)hw))); break; case INT_UINT: printf(" %29u\n", *((unsigned *)((void *)hw))); break; case INT_LONG: printf(" %29ld\n", *((long *)((void *)hw))); break; case INT_ULONG: printf(" %29lu\n", *((unsigned long *)((void *)hw))); break; case INT_LLONG: fprintf(stdout, " %29lld\n", *((long long *)((void *)hw))); break; case INT_ULLONG: fprintf(stdout, " %29llu\n", *((unsigned long long *)((void *)hw))); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } if (++fails_all_tests >= max_fails) { puts(" maximum failures reached, aborting test..."); puts(" (dst is library's conversion output. ans is compiler's conversion output.)"); goto done; } } PASSED(); done: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ return (int)fails_all_tests; error: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); 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 * *------------------------------------------------------------------------- */ static int test_conv_int_2(void) { int i, j; hid_t src_type, dst_type; char *buf; printf("%-70s", "Testing overlap calculations"); fflush(stdout); buf = (char *)calloc(TMP_BUF_DIM1, TMP_BUF_DIM2); assert(buf); for (i = 1; i <= TMP_BUF_DIM1; i++) { for (j = 1; j <= TMP_BUF_DIM1; 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, (size_t)TMP_BUF_DIM2, buf, NULL, H5P_DEFAULT); H5Tclose(src_type); H5Tclose(dst_type); } } PASSED(); free(buf); return 0; } /*------------------------------------------------------------------------- * Function: my_isnan * * Purpose: Determines whether VAL points to NaN. * * Return: true or false * *------------------------------------------------------------------------- */ static int my_isnan(dtype_t type, void *val) { int retval = 0; char s[256]; if (FLT_FLOAT == type) { float x = 0.0F; memcpy(&x, val, sizeof(float)); retval = isnan(x); } else if (FLT_DOUBLE == type) { double x = 0.0; memcpy(&x, val, sizeof(double)); retval = isnan(x); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else if (FLT_LDOUBLE == type) { long double x = 0.0L; memcpy(&x, val, sizeof(long double)); retval = isnan(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 = 0.0F; memcpy(&x, val, sizeof(float)); snprintf(s, sizeof(s), "%g", (double)x); } else if (FLT_DOUBLE == type) { double x = 0.0; memcpy(&x, val, sizeof(double)); snprintf(s, sizeof(s), "%g", x); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else if (FLT_LDOUBLE == type) { long double x = 0.0L; memcpy(&x, val, sizeof(long double)); snprintf(s, sizeof(s), "%Lg", x); #endif } else { return 0; } if (strstr(s, "NaN") || strstr(s, "NAN") || strstr(s, "nan")) retval = 1; } return retval; } /*------------------------------------------------------------------------- * Function: my_isinf * * Purpose: Determines whether VAL points to +/-infinity. * * Return: true or false * *------------------------------------------------------------------------- */ static int my_isinf(int endian, const unsigned char *val, size_t size, size_t mpos, size_t msize, size_t epos, size_t esize) { unsigned char *bits; int retval = 0; size_t i; bits = (unsigned char *)calloc((size_t)1, size); for (i = 0; i < size; i++) bits[size - (i + 1)] = *(val + ENDIAN(size, i, endian)); if (H5T__bit_find(bits, mpos, msize, H5T_BIT_LSB, 1) < 0 && H5T__bit_find(bits, epos, esize, H5T_BIT_LSB, 0) < 0) retval = 1; free(bits); return retval; } /*------------------------------------------------------------------------- * Function: test_conv_flt_1 * * Purpose: Test conversion of 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 * *------------------------------------------------------------------------- */ static int test_conv_flt_1(const char *name, int run_test, hid_t src, hid_t dst) { dtype_t src_type, dst_type; /*data types */ 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 */ void *aligned = NULL; /*aligned buffer */ float hw_f; /*hardware-converted */ double hw_d; /*hardware-converted */ #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 overflow = 0; /*overflow occurred */ size_t j, k; /*counters */ int sendian; /* source type endianness */ int dendian; /* Destination type endianness */ 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_mpos; /* Destination type's mantissa position */ 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 generate SIGFPE during floating point overflow and we * cannot assume that we can continue from such a signal. Therefore, we * fork here and let the child run the test and return the number of * failures with the exit status. */ fflush(stdout); fflush(stderr); if ((child_pid = fork()) < 0) { perror("fork"); return 1; } else if (child_pid > 0) { while (child_pid != waitpid(child_pid, &status, 0)) /*void*/ ; if (WIFEXITED(status) && 255 == WEXITSTATUS(status)) { return 0; /*child exit after catching SIGFPE*/ } else if (WIFEXITED(status)) { return WEXITSTATUS(status); } else if (WIFSIGNALED(status)) { snprintf(str, sizeof(str), " Child caught signal %d.", WTERMSIG(status)); puts(str); return 1; /*child exit after catching non-SIGFPE signal */ } else { puts(" 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)) puts("Sizeof(float)==sizeof(double) - some tests may not be sensible."); if (OTHER == src_type || OTHER == dst_type) { if (!strcmp(name, "noop")) snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_SPECIAL) snprintf(str, sizeof(str), "Testing %s special %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_NORMAL) snprintf(str, sizeof(str), "Testing %s normalized %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_DENORM) snprintf(str, sizeof(str), "Testing %s denormalized %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); H5_FAILED(); puts(" Unknown data type."); goto error; } else { if (!strcmp(name, "noop")) snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_SPECIAL) snprintf(str, sizeof(str), "Testing %s special %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_NORMAL) snprintf(str, sizeof(str), "Testing %s normalized %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_DENORM) snprintf(str, sizeof(str), "Testing %s denormalized %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); fflush(stdout); fails_this_test = 0; } /* 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, &dst_mpos, &dst_msize); sendian = H5Tget_order(src); dendian = H5Tget_order(dst); /* Allocate buffers */ aligned = calloc((size_t)1, MAX(sizeof(long double), sizeof(double))); /* Allocate and initialize the source buffer through macro INIT_FP_NORM 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_NORM 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_test) { case TEST_NOOP: case TEST_NORMAL: if (src_type == FLT_FLOAT) { INIT_FP_NORM(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); } else if (src_type == FLT_DOUBLE) { INIT_FP_NORM(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else if (src_type == FLT_LDOUBLE) { INIT_FP_NORM(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); #endif } else goto error; break; case TEST_DENORM: if (src_type == FLT_FLOAT) { INIT_FP_DENORM(float, FLT_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); } else if (src_type == FLT_DOUBLE) { INIT_FP_DENORM(double, DBL_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else if (src_type == FLT_LDOUBLE) { INIT_FP_DENORM(long double, LDBL_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); #endif } else goto error; break; case TEST_SPECIAL: if (src_type == FLT_FLOAT) { INIT_FP_SPECIAL(src_size, src_nbits, sendian, FLT_MANT_DIG, dst_size, buf, saved, nelmts); } else if (src_type == FLT_DOUBLE) { INIT_FP_SPECIAL(src_size, src_nbits, sendian, 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, sendian, LDBL_MANT_DIG, dst_size, buf, saved, nelmts); #endif } else goto error; break; default: goto error; } /* 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.0F; hw_d = 911.0; #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE hw_ld = 911.0L; #endif /* The hardware conversion */ /* Check for underflow when src is a "larger" float than dst.*/ if (FLT_FLOAT == src_type) { memcpy(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 = (double)*((float *)aligned); hw = (unsigned char *)&hw_d; #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else { hw_ld = (long double)*((float *)aligned); hw = (unsigned char *)&hw_ld; #endif } } else if (FLT_DOUBLE == src_type) { memcpy(aligned, saved + j * sizeof(double), sizeof(double)); if (FLT_FLOAT == dst_type) { hw_f = (float)(*((double *)aligned)); hw = (unsigned char *)&hw_f; underflow = fabs(*((double *)aligned)) < (double)FLT_MIN; overflow = fabs(*((double *)aligned)) > (double)FLT_MAX; } 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 = (long double)*((double *)aligned); hw = (unsigned char *)&hw_ld; #endif } #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else { memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); if (FLT_FLOAT == dst_type) { hw_f = (float)*((long double *)aligned); hw = (unsigned char *)&hw_f; underflow = fabsl(*((long double *)aligned)) < (long double)FLT_MIN; overflow = fabsl(*((long double *)aligned)) > (long double)FLT_MAX; } else if (FLT_DOUBLE == dst_type) { hw_d = (double)*((long double *)aligned); hw = (unsigned char *)&hw_d; underflow = fabsl(*((long double *)aligned)) < (long double)DBL_MIN; overflow = fabsl(*((long double *)aligned)) > (long double)DBL_MAX; } else { hw_ld = *((long double *)aligned); hw = (unsigned char *)&hw_ld; } #endif } /* 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 (sendian == H5T_ORDER_LE && dst_type == FLT_LDOUBLE) { size_t 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 attempt 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. * If overflow happens when the src value is greater than * the maximum dst value, the library assign INFINITY to dst. * This might be different from what the compiler does, i.e. * the SGI compiler assigns the dst's maximal value. */ { double check_mant[2]; int check_expo[2]; if (FLT_FLOAT == dst_type) { float x = 0.0F; memcpy(&x, &buf[j * dst_size], sizeof(float)); if (underflow && fabsf(x) <= FLT_MIN && fabsf(hw_f) <= FLT_MIN) continue; /* all underflowed, no error */ if (overflow && my_isinf(dendian, buf + j * sizeof(float), dst_size, dst_mpos, dst_msize, dst_epos, dst_esize)) continue; /* all overflowed, no error */ check_mant[0] = (double)frexpf(x, check_expo + 0); check_mant[1] = (double)frexpf(hw_f, check_expo + 1); } else if (FLT_DOUBLE == dst_type) { double x = 0.0; memcpy(&x, &buf[j * dst_size], sizeof(double)); if (underflow && fabs(x) <= DBL_MIN && fabs(hw_d) <= DBL_MIN) continue; /* all underflowed, no error */ if (overflow && my_isinf(dendian, buf + j * sizeof(double), dst_size, dst_mpos, dst_msize, dst_epos, dst_esize)) continue; /* all overflowed, no error */ check_mant[0] = frexp(x, check_expo + 0); check_mant[1] = frexp(hw_d, check_expo + 1); #if (H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE) } else { long double x = 0.0L; memcpy(&x, &buf[j * dst_size], sizeof(long double)); /* dst is largest float, no need to check underflow. */ check_mant[0] = (double)frexpl(x, check_expo + 0); check_mant[1] = (double)frexpl(hw_ld, check_expo + 1); #endif } /* 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) + (size_t)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] = ldexp(check_mant[0], expo_diff); /* Compute the proper epsilon */ epsilon = ldexp(epsilon, -valid_bits); /* Check for "close enough" fit with scaled epsilon value */ if (fabs(check_mant[0] - check_mant[1]) <= epsilon) continue; } /* end if */ else { if (check_expo[0] == check_expo[1] && fabs(check_mant[0] - check_mant[1]) < (double)FP_EPSILON) continue; } /* end else */ } if (0 == fails_this_test++) { if (run_test == TEST_NOOP || run_test == TEST_NORMAL) { H5_FAILED(); } else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) { H5_WARNING(); } } printf(" elmt %u\n", (unsigned)j); printf(" src ="); for (k = 0; k < src_size; k++) printf(" %02x", saved[j * src_size + ENDIAN(src_size, k, sendian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)dst_size - (ssize_t)src_size)), ""); if (FLT_FLOAT == src_type) { float x = 0.0F; memcpy(&x, &saved[j * src_size], sizeof(float)); printf(" %29.20e\n", (double)x); } else if (FLT_DOUBLE == src_type) { double x = 0.0; memcpy(&x, &saved[j * src_size], sizeof(double)); printf(" %29.20e\n", x); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else { long double x = 0.0L; memcpy(&x, &saved[j * src_size], sizeof(long double)); fprintf(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, dendian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)src_size - (ssize_t)dst_size)), ""); if (FLT_FLOAT == dst_type) { float x = 0.0F; memcpy(&x, &buf[j * dst_size], sizeof(float)); printf(" %29.20e\n", (double)x); } else if (FLT_DOUBLE == dst_type) { double x = 0.0; memcpy(&x, &buf[j * dst_size], sizeof(double)); printf(" %29.20e\n", x); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else { long double x = 0.0L; memcpy(&x, &buf[j * dst_size], sizeof(long double)); fprintf(stdout, " %29.20Le\n", x); #endif } printf(" ans ="); for (k = 0; k < dst_size; k++) printf(" %02x", hw[ENDIAN(dst_size, k, dendian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)src_size - (ssize_t)dst_size)), ""); if (FLT_FLOAT == dst_type) printf(" %29.20e\n", (double)hw_f); else if (FLT_DOUBLE == dst_type) printf(" %29.20e\n", hw_d); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE else fprintf(stdout, " %29.20Le\n", hw_ld); #endif /* If the source is normalized values, print out error message; if it is * denormalized or special values, print out warning message.*/ if (++fails_all_tests >= max_fails) { if (run_test == TEST_NORMAL) puts(" maximum failures reached, aborting test..."); else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) puts(" maximum warnings reached, aborting test..."); puts(" (dst is library's conversion output. ans is compiler's conversion output.)"); goto done; } } if (!fails_all_tests) PASSED(); done: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); #ifdef HANDLE_SIGFPE if (run_test == TEST_NOOP || run_test == TEST_NORMAL) exit(MIN((int)fails_all_tests, 254)); else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) exit(EXIT_SUCCESS); assert(0 && "Should not reach this point!"); return 1; #else /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /* If the source is normalized values, treat the failures as error; * if it is denormalized or special values, treat the failure as warning.*/ if (run_test == TEST_NOOP || run_test == TEST_NORMAL) return (int)fails_all_tests; else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) return 0; #endif error: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); #ifdef HANDLE_SIGFPE if (run_test == TEST_NOOP || run_test == TEST_NORMAL) exit(MIN(MAX((int)fails_all_tests, 1), 254)); else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) exit(EXIT_FAILURE); assert(0 && "Should not reach this point!"); return 1; #else /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); if (run_test == TEST_NOOP || run_test == TEST_NORMAL) return MAX((int)fails_all_tests, 1); else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) return 1; #endif } /*------------------------------------------------------------------------- * Function: test_conv_int_fp * * Purpose: Test conversion between 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 * *------------------------------------------------------------------------- */ static int test_conv_int_fp(const char *name, int run_test, 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*/ bool except_set = false; /*whether user's exception handling is set*/ 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 sendian; /*source endianness */ int dendian; /*destination endianness */ size_t src_size, dst_size; /*type sizes */ unsigned char *buf = NULL; /*buffer for conversion */ unsigned char *saved = NULL; /*original values */ size_t 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) { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); H5_FAILED(); puts(" 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)) { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); H5_FAILED(); puts(" 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)) { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); H5_FAILED(); puts(" 2. Not a float-integer conversion."); 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) { snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); fflush(stdout); fails_this_test = 0; } else { if (run_test == TEST_NORMAL) snprintf(str, sizeof(str), "Testing %s normalized %s -> %s conversions", name, src_type_name, dst_type_name); else if (run_test == TEST_DENORM) snprintf(str, sizeof(str), "Testing %s denormalized %s -> %s conversions", name, src_type_name, dst_type_name); else snprintf(str, sizeof(str), "Testing %s special %s -> %s conversions", name, src_type_name, dst_type_name); printf("%-70s", str); fflush(stdout); fails_this_test = 0; } /* Some information about datatypes */ sendian = H5Tget_order(src); dendian = H5Tget_order(dst); 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 = calloc((size_t)1, 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. */ memset(&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_NORM 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_test == TEST_NORMAL) { INIT_FP_NORM(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); } else if (run_test == TEST_DENORM) { INIT_FP_DENORM(float, FLT_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); } else { INIT_FP_SPECIAL(src_size, src_nbits, sendian, FLT_MANT_DIG, dst_size, buf, saved, nelmts); } } else if (src_type == FLT_DOUBLE) { if (run_test == TEST_NORMAL) { INIT_FP_NORM(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); } else if (run_test == TEST_DENORM) { INIT_FP_DENORM(double, DBL_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); } else { INIT_FP_SPECIAL(src_size, src_nbits, sendian, DBL_MANT_DIG, dst_size, buf, saved, nelmts); } #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE } else if (src_type == FLT_LDOUBLE) { if (run_test == TEST_NORMAL) { INIT_FP_NORM(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, src_size, dst_size, buf, saved, nelmts); } else if (run_test == TEST_DENORM) { INIT_FP_DENORM(long double, LDBL_MANT_DIG, src_size, src_nbits, sendian, dst_size, buf, saved, nelmts); } else { INIT_FP_SPECIAL(src_size, src_nbits, sendian, LDBL_MANT_DIG, dst_size, buf, saved, nelmts); } #endif } else goto error; /* 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 == src_type || FLT_DOUBLE == src_type || FLT_LDOUBLE == src_type) if (my_isnan(src_type, saved + j * src_size)) continue; if (FLT_FLOAT == dst_type) { hw = (unsigned char *)&hw_float; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_float = (float)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_float = (float)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_float = (float)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_float = (float)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_float = (float)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_float = (float)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_float = (float)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_float = (float)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_float = (float)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_float = (float)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (FLT_DOUBLE == dst_type) { hw = (unsigned char *)&hw_double; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_double = (double)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_double = (double)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_double = (double)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_double = (double)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_double = (double)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_double = (double)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_double = (double)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_double = (double)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_double = (double)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_double = (double)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (FLT_LDOUBLE == dst_type) { hw = (unsigned char *)&hw_ldouble; switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); hw_ldouble = (long double)(*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); hw_ldouble = (long double)(*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); hw_ldouble = (long double)(*((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); hw_ldouble = (long double)(*((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); hw_ldouble = (long double)(*((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); hw_ldouble = (long double)(*((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); hw_ldouble = (long double)(*((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); hw_ldouble = (long double)(*((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); hw_ldouble = (long double)(*((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); hw_ldouble = (long double)(*((unsigned long long *)aligned)); break; case FLT_FLOAT: case FLT_DOUBLE: case FLT_LDOUBLE: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_SCHAR == dst_type) { hw = (unsigned char *)&hw_schar; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_schar = (signed char)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_schar = (signed char)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_schar = (signed char)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_UCHAR == dst_type) { hw = (unsigned char *)&hw_uchar; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_uchar = (unsigned char)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_uchar = (unsigned char)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_uchar = (unsigned char)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_SHORT == dst_type) { hw = (unsigned char *)&hw_short; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_short = (short)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_short = (short)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_short = (short)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_USHORT == dst_type) { hw = (unsigned char *)&hw_ushort; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_ushort = (unsigned short)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_ushort = (unsigned short)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_ushort = (unsigned short)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_INT == dst_type) { hw = (unsigned char *)&hw_int; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_int = (int)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_int = (int)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_int = (int)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_UINT == dst_type) { hw = (unsigned char *)&hw_uint; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_uint = (unsigned int)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_uint = (unsigned int)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_uint = (unsigned int)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_LONG == dst_type) { hw = (unsigned char *)&hw_long; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_long = (long)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_long = (long)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_long = (long)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_ULONG == dst_type) { hw = (unsigned char *)&hw_ulong; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_ulong = (unsigned long)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_ulong = (unsigned long)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_ulong = (unsigned long)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_LLONG == dst_type) { hw = (unsigned char *)&hw_llong; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_llong = (long long)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_llong = (long long)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_llong = (long long)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); break; } } else if (INT_ULLONG == dst_type) { hw = (unsigned char *)&hw_ullong; switch (src_type) { case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); hw_ullong = (unsigned long long)(*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); hw_ullong = (unsigned long long)(*((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); hw_ullong = (unsigned long long)(*((long double *)aligned)); break; case INT_SCHAR: case INT_UCHAR: case INT_SHORT: case INT_USHORT: case INT_INT: case INT_UINT: case INT_LONG: case INT_ULONG: case INT_LLONG: case INT_ULLONG: case OTHER: default: assert(0 && "Unknown type"); 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 (dendian == H5T_ORDER_LE && dst_type == FLT_LDOUBLE) { size_t 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*/ /* * 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. */ if ((FLT_FLOAT == src_type || FLT_DOUBLE == src_type) && sendian == H5T_ORDER_VAX) { for (k = 0; k < src_size; k += 2) { src_bits[k] = saved[j * src_size + (src_size - 2) - k]; src_bits[k + 1] = saved[j * src_size + (src_size - 1) - k]; } } else { for (k = 0; k < src_size; k++) src_bits[src_size - (k + 1)] = saved[j * src_size + ENDIAN(src_size, k, sendian)]; } for (k = 0; k < dst_size; k++) dst_bits[dst_size - (k + 1)] = buf[j * dst_size + ENDIAN(dst_size, k, dendian)]; /* 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, (size_t)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, (size_t)1) && H5T__bit_find(dst_bits, (size_t)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, (size_t)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, (size_t)1) && H5T__bit_find(dst_bits, (size_t)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, (size_t)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, (size_t)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, (size_t)0, dst_nbits, H5T_BIT_LSB, 0) < 0) continue; /*no error*/ } else { if (dst_bits[0] == fill_value) continue; /*no error*/ } } } /* Print errors */ if (0 == fails_this_test++) { if (run_test == TEST_NORMAL) { H5_FAILED(); } else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) { H5_WARNING(); } } printf(" elmt %u: \n", (unsigned)j); printf(" src = "); for (k = 0; k < src_size; k++) printf(" %02x", saved[j * src_size + ENDIAN(src_size, k, sendian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)dst_size - (ssize_t)src_size)), ""); switch (src_type) { case INT_SCHAR: memcpy(aligned, saved + j * sizeof(signed char), sizeof(signed char)); printf(" %29d\n", (int)*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, saved + j * sizeof(unsigned char), sizeof(unsigned char)); printf(" %29u\n", (unsigned)*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, saved + j * sizeof(short), sizeof(short)); printf(" %29hd\n", *((short *)aligned)); break; case INT_USHORT: memcpy(aligned, saved + j * sizeof(unsigned short), sizeof(unsigned short)); printf(" %29hu\n", *((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, saved + j * sizeof(int), sizeof(int)); printf(" %29d\n", *((int *)aligned)); break; case INT_UINT: memcpy(aligned, saved + j * sizeof(unsigned), sizeof(unsigned)); printf(" %29u\n", *((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, saved + j * sizeof(long), sizeof(long)); printf(" %29ld\n", *((long *)aligned)); break; case INT_ULONG: memcpy(aligned, saved + j * sizeof(unsigned long), sizeof(unsigned long)); printf(" %29lu\n", *((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, saved + j * sizeof(long long), sizeof(long long)); fprintf(stdout, " %29lld\n", *((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, saved + j * sizeof(unsigned long long), sizeof(unsigned long long)); fprintf(stdout, " %29llu\n", *((unsigned long long *)aligned)); break; case FLT_FLOAT: memcpy(aligned, saved + j * sizeof(float), sizeof(float)); printf(" %29f\n", (double)*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, saved + j * sizeof(double), sizeof(double)); printf(" %29f\n", *((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, saved + j * sizeof(long double), sizeof(long double)); printf(" %29Lf\n", *((long double *)aligned)); break; case OTHER: default: assert(0 && "Unknown type"); break; } printf(" dst = "); for (k = 0; k < dst_size; k++) printf(" %02x", buf[j * dst_size + ENDIAN(dst_size, k, dendian)]); printf("%*s", (int)(3 * MAX(0, (ssize_t)src_size - (ssize_t)dst_size)), ""); switch (dst_type) { case INT_SCHAR: memcpy(aligned, buf + j * sizeof(signed char), sizeof(signed char)); printf(" %29d\n", (int)*((signed char *)aligned)); break; case INT_UCHAR: memcpy(aligned, buf + j * sizeof(unsigned char), sizeof(unsigned char)); printf(" %29u\n", (unsigned)*((unsigned char *)aligned)); break; case INT_SHORT: memcpy(aligned, buf + j * sizeof(short), sizeof(short)); printf(" %29hd\n", *((short *)aligned)); break; case INT_USHORT: memcpy(aligned, buf + j * sizeof(unsigned short), sizeof(unsigned short)); printf(" %29hu\n", *((unsigned short *)aligned)); break; case INT_INT: memcpy(aligned, buf + j * sizeof(int), sizeof(int)); printf(" %29d\n", *((int *)aligned)); break; case INT_UINT: memcpy(aligned, buf + j * sizeof(unsigned), sizeof(unsigned)); printf(" %29u\n", *((unsigned *)aligned)); break; case INT_LONG: memcpy(aligned, buf + j * sizeof(long), sizeof(long)); printf(" %29ld\n", *((long *)aligned)); break; case INT_ULONG: memcpy(aligned, buf + j * sizeof(unsigned long), sizeof(unsigned long)); printf(" %29lu\n", *((unsigned long *)aligned)); break; case INT_LLONG: memcpy(aligned, buf + j * sizeof(long long), sizeof(long long)); fprintf(stdout, " %29lld\n", *((long long *)aligned)); break; case INT_ULLONG: memcpy(aligned, buf + j * sizeof(unsigned long long), sizeof(unsigned long long)); fprintf(stdout, " %29llu\n", *((unsigned long long *)aligned)); break; case FLT_FLOAT: memcpy(aligned, buf + j * sizeof(float), sizeof(float)); printf(" %29f\n", (double)*((float *)aligned)); break; case FLT_DOUBLE: memcpy(aligned, buf + j * sizeof(double), sizeof(double)); printf(" %29f\n", *((double *)aligned)); break; case FLT_LDOUBLE: memcpy(aligned, buf + j * sizeof(long double), sizeof(long double)); printf(" %29Lf\n", *((long double *)aligned)); break; case OTHER: default: assert(0 && "Unknown type"); break; } printf(" and = "); for (k = 0; k < dst_size; k++) printf(" %02x", hw[ENDIAN(dst_size, k, dendian)]); 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 *)((void *)hw))); break; case INT_UCHAR: printf(" %29u\n", (unsigned)*((unsigned char *)((void *)hw))); break; case INT_SHORT: printf(" %29hd\n", *((short *)((void *)hw))); break; case INT_USHORT: printf(" %29hu\n", *((unsigned short *)((void *)hw))); break; case INT_INT: printf(" %29d\n", *((int *)((void *)hw))); break; case INT_UINT: printf(" %29u\n", *((unsigned int *)((void *)hw))); break; case INT_LONG: printf(" %29ld\n", *((long *)((void *)hw))); break; case INT_ULONG: printf(" %29lu\n", *((unsigned long *)((void *)hw))); break; case INT_LLONG: fprintf(stdout, " %29lld\n", *((long long *)((void *)hw))); break; case INT_ULLONG: fprintf(stdout, " %29llu\n", *((unsigned long long *)((void *)hw))); break; case FLT_FLOAT: printf(" %29f\n", (double)*((float *)((void *)hw))); break; case FLT_DOUBLE: printf(" %29f\n", *((double *)((void *)hw))); break; case FLT_LDOUBLE: printf(" %29Lf\n", *((long double *)((void *)hw))); break; case OTHER: default: assert(0 && "Unknown type"); break; } /* If the source is normalized values, print out error message; if it is * denormalized or special values, print out warning message.*/ if (++fails_all_tests >= max_fails) { if (run_test == TEST_NORMAL) puts(" maximum failures reached, aborting test..."); else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) puts(" maximum warnings reached, aborting test..."); puts(" (dst is library's conversion output. ans is compiler's conversion output.)"); goto done; } } if (!fails_all_tests) PASSED(); done: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ /* If the source is normalized floating values, treat the failures as error; * if it is denormalized or special floating values, treat the failure as warning.*/ if (run_test == TEST_NORMAL) return (int)fails_all_tests; else if (run_test == TEST_DENORM || run_test == TEST_SPECIAL) return 0; error: if (buf) aligned_free(buf); if (saved) aligned_free(saved); if (aligned) free(aligned); fflush(stdout); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /*print statistics*/ if (run_test == TEST_NORMAL) return MAX((int)fails_all_tests, 1); else { assert(run_test == TEST_DENORM || run_test == TEST_SPECIAL); return 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 * *------------------------------------------------------------------------- */ static bool overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits) { bool 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]; memset(bits, 0, (size_t)32); memset(mant_bits, 0, (size_t)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); memcpy(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, (size_t)0, bits, (size_t)0, mant_digits); /* restore implicit bit if normalization is implied*/ if (norm == H5T_NORM_IMPLIED) { H5T__bit_inc(mant_bits, mant_digits, (size_t)1); mant_digits++; } /* shift significand */ H5T__bit_shift(mant_bits, (ssize_t)(expt - expt_digits), (size_t)0, (size_t)(32 * 8)); indx = H5T__bit_find(mant_bits, (size_t)0, (size_t)(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 * *------------------------------------------------------------------------- */ 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 * *------------------------------------------------------------------------- */ static int run_fp_tests(const char *name) { int nerrors = 0; if (!strcmp(name, "noop")) { nerrors += test_conv_flt_1("noop", TEST_NOOP, H5T_NATIVE_FLOAT, H5T_NATIVE_FLOAT); nerrors += test_conv_flt_1("noop", TEST_NOOP, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE); nerrors += test_conv_flt_1("noop", TEST_NOOP, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LDOUBLE); goto done; } /*Test normalized values. TEST_NORMAL indicates normalized values.*/ nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE); nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE); nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT); nerrors += test_conv_flt_1(name, TEST_NORMAL, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE); #endif /*Test denormalized values. TEST_DENORM indicates denormalized values.*/ nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE); nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE); #ifndef H5_DISABLE_SOME_LDOUBLE_CONV nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT); #else { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s denormalized %s -> %s conversions", name, "long double", "float"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to the conversion problem on IBM ppc64le cpu."); } #endif nerrors += test_conv_flt_1(name, TEST_DENORM, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE); #endif /*Test special values, +/-0, +/-infinity, +/-QNaN, +/-SNaN.*/ nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE); nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT); #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE); #ifndef H5_DISABLE_SOME_LDOUBLE_CONV nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT); nerrors += test_conv_flt_1(name, TEST_SPECIAL, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE); #else { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s special %s -> %s conversions", name, "long double", "float or double"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to the conversion problem on IBM ppc64le cpu."); } #endif #endif done: return nerrors; } /*------------------------------------------------------------------------- * Function: run_int_fp_conv * * Purpose: Runs all integer-float tests. * * Return: Number of errors * *------------------------------------------------------------------------- */ static int run_int_fp_conv(const char *name) { int nerrors = 0; nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SCHAR, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SCHAR, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UCHAR, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UCHAR, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SHORT, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SHORT, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_USHORT, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_USHORT, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_INT, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_INT, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UINT, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UINT, H5T_NATIVE_DOUBLE); #if H5_SIZEOF_LONG != H5_SIZEOF_INT nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LONG, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LONG, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULONG, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULONG, H5T_NATIVE_DOUBLE); #endif #if H5_SIZEOF_LONG_LONG != H5_SIZEOF_LONG nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LLONG, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LLONG, H5T_NATIVE_DOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULLONG, H5T_NATIVE_FLOAT); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULLONG, H5T_NATIVE_DOUBLE); #endif #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SCHAR, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UCHAR, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_SHORT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_USHORT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_INT, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_UINT, H5T_NATIVE_LDOUBLE); #if H5_SIZEOF_LONG != H5_SIZEOF_INT #if !defined(H5_LONG_TO_LDOUBLE_SPECIAL) && !defined(H5_DISABLE_SOME_LDOUBLE_CONV) nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LONG, H5T_NATIVE_LDOUBLE); nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULONG, H5T_NATIVE_LDOUBLE); #else { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "(unsigned) long", "long double"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to the special algorithm of hardware conversion."); } #endif #endif /* H5_SIZEOF_LONG!=H5_SIZEOF_INT */ #if H5_SIZEOF_LONG_LONG != H5_SIZEOF_LONG #if H5_LLONG_TO_LDOUBLE_CORRECT nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_LLONG, H5T_NATIVE_LDOUBLE); #else /* H5_LLONG_TO_LDOUBLE_CORRECT */ { char str[256]; /*hello string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "long long", "long double"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to compiler error in handling conversion."); } #endif /* H5_LLONG_TO_LDOUBLE_CORRECT */ #if H5_LLONG_TO_LDOUBLE_CORRECT nerrors += test_conv_int_fp(name, TEST_NORMAL, H5T_NATIVE_ULLONG, H5T_NATIVE_LDOUBLE); #else /* H5_LLONG_TO_LDOUBLE_CORRECT */ { char str[256]; /*hello string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "unsigned long long", "long double"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to compiler not handling conversion."); } #endif /* H5_LLONG_TO_LDOUBLE_CORRECT */ #endif #endif return nerrors; } /*------------------------------------------------------------------------- * Function: run_fp_int_conv * * Purpose: Runs all float-integer tests. * * Return: Number of errors * *------------------------------------------------------------------------- */ static int run_fp_int_conv(const char *name) { int nerrors = 0; int test_values; for (test_values = TEST_NORMAL; test_values <= TEST_SPECIAL; test_values++) { nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_SCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_SCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_UCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_UCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_SHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_SHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_USHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_USHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_INT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_INT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_UINT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_UINT); #if H5_SIZEOF_LONG != H5_SIZEOF_INT nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_LONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_LONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_ULONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULONG); #endif #if H5_SIZEOF_LONG_LONG != H5_SIZEOF_LONG if (!strcmp(name, "hw")) { /* Hardware conversion */ nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG); } else { /* Software conversion */ nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG); } nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_FLOAT, H5T_NATIVE_ULLONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULLONG); #endif #if H5_SIZEOF_LONG_DOUBLE != H5_SIZEOF_DOUBLE if (test_values != TEST_SPECIAL) { nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_USHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_INT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UINT); } else { #ifndef H5_DISABLE_SOME_LDOUBLE_CONV nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UCHAR); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_USHORT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_INT); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UINT); #else char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s special %s -> %s conversions", name, "long double", "signed and unsigned char, short, int, long"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to the conversion problem on IBM ppc64le cpu."); #endif } #if H5_SIZEOF_LONG != H5_SIZEOF_INT #ifndef H5_LDOUBLE_TO_LONG_SPECIAL if (test_values != TEST_SPECIAL && test_values != TEST_NORMAL) { nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULONG); } else { #ifndef H5_DISABLE_SOME_LDOUBLE_CONV nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LONG); nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULONG); #endif } #else { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "long double", "(unsigned) long"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to the special algorithm of hardware conversion."); } #endif #endif /*H5_SIZEOF_LONG!=H5_SIZEOF_INT */ #if H5_SIZEOF_LONG_LONG != H5_SIZEOF_LONG #ifdef H5_LDOUBLE_TO_LLONG_ACCURATE nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LLONG); #else /*H5_LDOUBLE_TO_LLONG_ACCURATE*/ { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "long double", "long long"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to hardware conversion error."); } #endif /*H5_LDOUBLE_TO_LLONG_ACCURATE*/ #if defined(H5_LDOUBLE_TO_LLONG_ACCURATE) nerrors += test_conv_int_fp(name, test_values, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULLONG); #else /*H5_LDOUBLE_TO_LLONG_ACCURATE*/ { char str[256]; /*string */ snprintf(str, sizeof(str), "Testing %s %s -> %s conversions", name, "long double", "unsigned long long"); printf("%-70s", str); SKIPPED(); puts(" Test skipped due to hardware conversion error."); } #endif /*H5_LDOUBLE_TO_LLONG_ACCURATE*/ #endif #endif } /* end for */ return nerrors; } /*------------------------------------------------------------------------- * Function: main * * Purpose: Test the data type(integer and floating-point number). * * Return: Success: * * Failure: * *------------------------------------------------------------------------- */ int main(void) { unsigned long nerrors = 0; /* Set the random # seed */ HDsrandom((unsigned)HDtime(NULL)); reset_hdf5(); if (ALIGNMENT) printf("Testing non-aligned conversions (ALIGNMENT=%d)....\n", ALIGNMENT); /* Do the tests */ /* Test H5Tcompiler_conv() for querying hard conversion. */ nerrors += (unsigned long)test_hard_query(); /* Test user-define, query functions and software conversion * for user-defined floating-point types */ nerrors += (unsigned long)test_derived_flt(); /* Test user-define, query functions and software conversion * for user-defined integer types */ nerrors += (unsigned long)test_derived_integer(); /* Test degenerate cases */ nerrors += (unsigned long)run_fp_tests("noop"); /* Test hardware floating-point conversion functions */ nerrors += (unsigned long)run_fp_tests("hard"); /* Test hardware integer conversion functions */ nerrors += (unsigned long)run_integer_tests("hard"); /* Test hardware integer-float conversion functions */ nerrors += (unsigned long)run_int_fp_conv("hard"); /* Test hardware float-integer conversion functions */ nerrors += (unsigned long)run_fp_int_conv("hard"); /* Test a few special values for hardware float-integer conversions */ nerrors += (unsigned long)test_particular_fp_integer(); /*---------------------------------------------------------------------- * Software tests *---------------------------------------------------------------------- */ without_hardware_g = true; /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /* Test software floating-point conversion functions */ nerrors += (unsigned long)run_fp_tests("soft"); /* Test software integer conversion functions */ nerrors += (unsigned long)test_conv_int_2(); nerrors += (unsigned long)run_integer_tests("soft"); /* Test software float-integer conversion functions */ nerrors += (unsigned long)run_fp_int_conv("soft"); /* Test software integer-float conversion functions */ nerrors += (unsigned long)run_int_fp_conv("soft"); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); reset_hdf5(); /* Restore the default error handler (set in h5_reset()) */ h5_restore_err(); if (nerrors) { printf("***** %lu FAILURE%s! *****\n", nerrors, 1 == nerrors ? "" : "S"); exit(EXIT_FAILURE); } printf("All data type tests passed.\n"); return 0; }