From 138abc32798ac465e600965cf63aba760d968f88 Mon Sep 17 00:00:00 2001 From: Raymond Lu Date: Tue, 5 Apr 2005 15:25:21 -0500 Subject: [svn-r10550] Purpose: Splitting dtypes.c Description: Yesterday, dtypes.c was split into two programs, dtypes.c and dt_atomic.c. After dt_atomic.c was added, dtypes.c remained the same as before waiting for the daily test to pass. While dt_atomic.c turned out fine today, the test in dt_atomic.c was taken out from dtypes.c. There's also some minor change in dt_atomic.c. Platforms tested: h5committest --- test/dt_atomic.c | 104 + test/dtypes.c | 6165 +++++++----------------------------------------------- 2 files changed, 878 insertions(+), 5391 deletions(-) diff --git a/test/dt_atomic.c b/test/dt_atomic.c index e1386f3..f60be7c 100644 --- a/test/dt_atomic.c +++ b/test/dt_atomic.c @@ -69,6 +69,20 @@ typedef enum dtype_t { FLT_LDOUBLE, OTHER } dtype_t; +/* Skip overflow tests if non-zero */ +static int skip_overflow_tests_g = 0; + +/* + * Although we check whether a floating point overflow generates a SIGFPE and + * turn off overflow tests in that case, it might still be possible for an + * overflow condition to occur. Once a SIGFPE is raised the program cannot + * be allowed to continue (cf. Posix signals) so in order to recover from a + * SIGFPE we run tests that might generate one in a child process. + */ +#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID) +#define HANDLE_SIGFPE +#endif + /* Don't use hardware conversions if set */ static int without_hardware_g = 0; @@ -314,6 +328,7 @@ static int without_hardware_g = 0; free(value); \ } +void some_dummy_func(float x); static hbool_t overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits); static int my_isnan(dtype_t type, void *val); @@ -433,6 +448,92 @@ except_func(int except_type, hid_t UNUSED src_id, hid_t UNUSED dst_id, void UNUS /*------------------------------------------------------------------------- + * Function: some_dummy_func + * + * Purpose: A dummy function to help check for overflow. + * + * Note: DO NOT DECLARE THIS FUNCTION STATIC OR THE COMPILER MIGHT + * PROMOTE ARGUMENT `x' TO DOUBLE AND DEFEAT THE OVERFLOW + * CHECKING. + * + * Return: void + * + * Programmer: Robb Matzke + * Tuesday, July 21, 1998 + * + * Modifications: + * + *------------------------------------------------------------------------- + */ +void +some_dummy_func(float x) +{ + char s[128]; + sprintf(s, "%g", x); +} + + +/*------------------------------------------------------------------------- + * Function: generates_sigfpe + * + * Purpose: Determines if SIGFPE is generated from overflows. We must be + * able to fork() and waitpid() in order for this test to work + * properly. Sets skip_overflow_tests_g to non-zero if they + * would generate SIGBUS, zero otherwise. + * + * Programmer: Robb Matzke + * Tuesday, July 21, 1998 + * + * Modifications: + * + *------------------------------------------------------------------------- + */ +static void +generates_sigfpe(void) +{ +#if defined(H5_HAVE_FORK) && defined(H5_HAVE_WAITPID) + pid_t pid; + int status; + size_t i, j; + double d; + unsigned char *dp = (unsigned char*)&d; + float f; + + HDfflush(stdout); + HDfflush(stderr); + if ((pid=fork())<0) { + HDperror("fork"); + HDexit(1); + } else if (0==pid) { + for (i=0; i<2000; i++) { + for (j=0; jinteger and integer->floating-point conversion test. - * This algorithm is mainly to avoid any casting and comparison between source and destination types - * for compiler, because we're testing conversions. */ -#define INIT_INTEGER(TYPE, SRC_MAX, SRC_MIN, SRC_SIZE, DST_SIZE, SRC_PREC, BUF, SAVED, NELMTS) \ -{ \ - unsigned char *buf_p, *saved_p; \ - int n; \ - TYPE value1 = 1; \ - TYPE value2 = 0; \ - \ - /* Allocate buffers */ \ - NELMTS=(SRC_PREC-1)*3+1; \ - BUF = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \ - SAVED = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \ - \ - buf_p = BUF; \ - saved_p = SAVED; \ - \ - /*positive values, ascending order. VALUE1 starts from 00000001, to 00000010, until 10000000*/ \ - /*VALUE2 ascends from 00000000, to 00000011, 00000111,..., until 11111111.*/ \ - for(n=0; n=SRC_MIN) { \ - memcpy(buf_p, &value1, SRC_SIZE); \ - memcpy(saved_p, &value1, SRC_SIZE); \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - } \ - if(value2<=SRC_MAX && value2>=SRC_MIN) { \ - memcpy(buf_p, &value2, SRC_SIZE); \ - memcpy(saved_p, &value2, SRC_SIZE); \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - } \ - \ - value1 <<= 1; \ - value2 = (value1 - 1) | value1; \ - } \ - \ - /* negative values for signed; descending positive values for unsigned */ \ - /* VALUE2 descends from 11111111 to 11111110, 11111100, ..., until 10000000. */ \ - for(n=0; n=SRC_MIN) { \ - memcpy(buf_p, &value2, SRC_SIZE); \ - memcpy(saved_p, &value2, SRC_SIZE); \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - } \ - value2 <<= 1; \ - } \ -} - -/* Change a buffer's byte order from big endian to little endian. It's mainly for library's - * bit operations which handle only little endian order. - */ -#define CHANGE_ORDER(EBUF, EORDER, ESIZE) \ -{ \ - if (H5T_ORDER_BE==EORDER) { \ - int m; \ - unsigned char mediator; \ - size_t half_size = ESIZE/2; \ - for (m=0; m=100 && SRC_MAX_10_EXP<400) { /*for double*/ \ - factor = 2; \ - multiply = 10000; \ - } else { /*for long double*/ \ - factor = 3; \ - multiply = 100000000; \ - } \ - \ - /*The number of values if multiplied by 10 for each step.*/ \ - num_norm = (SRC_MAX_10_EXP - SRC_MIN_10_EXP); \ - /*Reduce the number of values by 2^factor. MULTIPLY=10^(2^factor). Using this algorithm \ - *instead of arithmatic operation to avoid any conversion*/ \ - num_norm >>= factor; \ - \ - /*Total number of values*/ \ - NELMTS = 2 * /*both positive and negative*/ \ - (num_norm + /*number of normalized values*/ \ - 1 + /*maximal normalized value*/ \ - SRC_MANT_DIG - 1); /*number of denormalized values*/ \ - \ - /* Allocate buffers */ \ - BUF = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \ - SAVED = (unsigned char*)aligned_malloc(NELMTS*MAX(SRC_SIZE, DST_SIZE)); \ - tmp1 = (unsigned char*)malloc(SRC_SIZE); \ - tmp2 = (unsigned char*)malloc(SRC_SIZE); \ - \ - buf_p = BUF; \ - saved_p = SAVED; \ - \ - /*Normalized values*/ \ - value1 = SRC_MIN; \ - value2 = -SRC_MIN; \ - for(n=0; n-SRC_MAX) { /*negative*/ \ - memcpy(buf_p, &value2, SRC_SIZE); \ - memcpy(saved_p, &value2, SRC_SIZE); \ - value2 *= multiply; \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - } \ - } \ - \ - value1 = SRC_MAX; /*maximal value*/ \ - memcpy(buf_p, &value1, SRC_SIZE); \ - memcpy(saved_p, &value1, SRC_SIZE); \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - \ - value2 = -SRC_MAX; /*negative value*/ \ - memcpy(buf_p, &value2, SRC_SIZE); \ - memcpy(saved_p, &value2, SRC_SIZE); \ - buf_p += SRC_SIZE; \ - saved_p += SRC_SIZE; \ - \ - /*Denormalized values. Exponent is 0. Let mantissa starts from 00000001, 00000011, \ - *00000111,..., until 11111111.*/ \ - memset(tmp1, 0, SRC_SIZE); \ - memset(tmp2, 0, SRC_SIZE); \ - H5T_bit_set (tmp2, SRC_PREC-1, 1, TRUE); /*the negative value*/ \ - for(n=0; n set offset -> set precision -> set size. - * All these properties must be set before the type can function. Other - * properties can be set anytime. Derived type size cannot be expanded - * bigger than original size but can be decreased. There should be no - * holes among the significant bits. Exponent bias usually is set - * 2^(n-1)-1, where n is the exponent size. - *-----------------------------------------------------------------------*/ - if(H5Tset_fields(tid1, 44, 34, 10, 3, 31)<0) { - H5_FAILED(); - printf("Can't set fields\n"); - goto error; - } - if(H5Tset_offset(tid1, 3)<0) { - H5_FAILED(); - printf("Can't set offset\n"); - goto error; - } - if(H5Tset_precision(tid1, 42)<0) { - H5_FAILED(); - printf("Can't set precision 1\n"); - goto error; + /* Predefined types cannot be modified or closed */ + H5E_BEGIN_TRY { + status = H5Tset_precision (H5T_NATIVE_INT, 256); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Predefined types should not be modifiable!"); + goto error; } - if(H5Tset_size(tid1, 7)<0) { - H5_FAILED(); - printf("Can't set size\n"); - goto error; + H5E_BEGIN_TRY { + status = H5Tclose (H5T_NATIVE_INT); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Predefined types should not be closable!"); + goto error; } - if(H5Tset_ebias(tid1, 511)<0) { - H5_FAILED(); - printf("Can't set exponent bias\n"); - goto error; - } - if(H5Tset_pad(tid1, H5T_PAD_ZERO, H5T_PAD_ZERO)<0) { - H5_FAILED(); - printf("Can't set padding\n"); - goto error; - } + /* Copying a predefined type results in a modifiable copy */ + if ((type=H5Tcopy (H5T_NATIVE_INT))<0) goto error; + if (H5Tset_precision (type, 256)<0) goto error; - if(H5Tcommit(file, "new float type 1", tid1)<0) { - H5_FAILED(); - printf("Can't set inpad\n"); - goto error; - } - if(H5Tclose(tid1)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; + /* It should not be possible to create an attribute for a transient type */ + H5E_BEGIN_TRY { + status = H5Acreate (type, "attr1", H5T_NATIVE_INT, space, H5P_DEFAULT); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Attributes should not be allowed for transient types!"); + goto error; } - if((tid1 = H5Topen(file, "new float type 1"))<0) { - H5_FAILED(); - printf("Can't open datatype\n"); - goto error; - } - if(H5Tget_fields(tid1, &spos, &epos, &esize, &mpos, &msize)<0) { - H5_FAILED(); - printf("Can't get fields\n"); - goto error; - } - if(spos!=44 || epos!=34 || esize!=10 || mpos!=3 || msize!=31) { - H5_FAILED(); - printf("Wrong field values\n"); - goto error; - } + /* Create a dataset from a transient data type */ + if (H5Tclose (type)<0) goto error; + if ((type = H5Tcopy (H5T_NATIVE_INT))<0) goto error; + if ((dset=H5Dcreate (file, "dset1", type, space, H5P_DEFAULT))<0) + goto error; - if((precision = H5Tget_precision(tid1))!=42) { - H5_FAILED(); - printf("Can't get precision or wrong precision\n"); - goto error; - } - if((offset = H5Tget_offset(tid1))!=3) { - H5_FAILED(); - printf("Can't get offset or wrong offset\n"); - goto error; - } - if((size = H5Tget_size(tid1))!=7) { - H5_FAILED(); - printf("Can't get size or wrong size\n"); - goto error; - } - if((ebias = H5Tget_ebias(tid1))!=511) { - H5_FAILED(); - printf("Can't get exponent bias or wrong bias\n"); - goto error; + /* The type returned from a dataset should not be modifiable */ + if ((t2 = H5Dget_type (dset))<0) goto error; + H5E_BEGIN_TRY { + status = H5Tset_precision (t2, 256); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Dataset data types should not be modifiable!"); + goto error; } + if (H5Tclose (t2)<0) 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. + /* + * Close the dataset and reopen it, testing that it's type is still + * read-only. */ - src_size = H5Tget_size(H5T_NATIVE_INT); - endian = H5Tget_order(H5T_NATIVE_INT); - buf = (unsigned char*)malloc(nelmts*(MAX(src_size, size))); - saved_buf = (unsigned char*)malloc(nelmts*src_size); - aligned = (int*)malloc(src_size); - - for(i=0; i=0) { + H5_FAILED(); + HDputs (" Dataset data types should not be modifiable!"); + goto error; } + if (H5Tclose (t2)<0) goto error; - /* Are the values still the same?*/ - for(i=0; i derived floating-point conversions"); - printf("%-70s", str); - HDfflush(stdout); - H5_FAILED(); - } - printf(" test %u elmt %u: \n", 1, (unsigned)i); - - printf(" src = "); - for (j=0; j=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto error; - } + h5_fixname(FILENAME[1], fapl, filename, sizeof filename); + if ((file=H5Fcreate (filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl))<0) { + goto error; } + if ((space = H5Screate_simple (2, ds_size, ds_size))<0) goto error; - fails_this_test = 0; - if(buf) free(buf); - if(saved_buf) free(saved_buf); - if(aligned) free(aligned); - buf = NULL; - saved_buf = NULL; - aligned = NULL; - - /*-------------------------------------------------------------------------- - * 2nd floating-point type - * size=3 byte, precision=24 bits, offset=0 bits, mantissa size=16 bits, - * mantissa position=0, exponent size=7 bits, exponent position=16, exponent - * bias=63. It can be illustrated in little-endian order as - * - * 2 1 0 - * SEEEEEEE MMMMMMMM MMMMMMMM - *--------------------------------------------------------------------------*/ - if(H5Tset_fields(tid2, 23, 16, 7, 0, 16)<0) { - H5_FAILED(); - printf("Can't set fields\n"); - goto error; - } - if(H5Tset_offset(tid2, 0)<0) { - H5_FAILED(); - printf("Can't set offset\n"); - goto error; - } - if(H5Tset_precision(tid2, 24)<0) { - H5_FAILED(); - printf("Can't set precision 2\n"); - goto error; - } - if(H5Tset_size(tid2, 3)<0) { - H5_FAILED(); - printf("Can't set size\n"); - goto error; - } - if(H5Tset_ebias(tid2, 63)<0) { - H5_FAILED(); - printf("Can't set size\n"); - goto error; - } - if(H5Tset_pad(tid2, H5T_PAD_ZERO, H5T_PAD_ZERO)<0) { - H5_FAILED(); - printf("Can't set padding\n"); - goto error; + /* Predefined types cannot be committed */ + H5E_BEGIN_TRY { + status = H5Tcommit (file, "test_named_1 (should not exist)", + H5T_NATIVE_INT); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Predefined types should not be committable!"); + goto error; } - if(H5Tcommit(file, "new float type 2", tid2)<0) { - H5_FAILED(); - printf("Can't set inpad\n"); - goto error; - } - if(H5Tclose(tid2)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; + /* Copy a predefined data type and commit the copy */ + if ((type = H5Tcopy (H5T_NATIVE_INT))<0) goto error; + if (H5Tcommit (file, "native-int", type)<0) goto error; + if ((status=H5Tcommitted (type))<0) goto error; + if (0==status) { + H5_FAILED(); + HDputs (" H5Tcommitted() returned false!"); + goto error; } - if((tid2 = H5Topen(file, "new float type 2"))<0) { - H5_FAILED(); - printf("Can't open datatype\n"); - goto error; - } - if(H5Tget_fields(tid2, &spos, &epos, &esize, &mpos, &msize)<0) { - H5_FAILED(); - printf("Can't get fields\n"); - goto error; - } - if(spos!=23 || epos!=16 || esize!=7 || mpos!=0 || msize!=16) { - H5_FAILED(); - printf("Wrong field values\n"); - goto error; + /* We should not be able to modify a type after it has been committed. */ + H5E_BEGIN_TRY { + status = H5Tset_precision (type, 256); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Committed type is not constant!"); + goto error; } - if((precision = H5Tget_precision(tid2))!=24) { - H5_FAILED(); - printf("Can't get precision or wrong precision\n"); - goto error; - } - if((offset = H5Tget_offset(tid2))!=0) { - H5_FAILED(); - printf("Can't get offset or wrong offset\n"); - goto error; - } - if((size = H5Tget_size(tid2))!=3) { - H5_FAILED(); - printf("Can't get size or wrong size\n"); - goto error; + /* We should not be able to re-commit a committed type */ + H5E_BEGIN_TRY { + status = H5Tcommit(file, "test_named_2 (should not exist)", type); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + HDputs (" Committed types should not be recommitted!"); + goto error; } - if((ebias = H5Tget_ebias(tid2))!=63) { - H5_FAILED(); - printf("Can't get exponent bias or wrong bias\n"); - goto error; + + /* It should be possible to define an attribute for the named type */ + if ((attr1=H5Acreate (type, "attr1", H5T_NATIVE_UCHAR, space, + H5P_DEFAULT))<0) goto error; + for (i=0; i derived floating-point conversions"); - printf("%-70s", str); - HDfflush(stdout); - H5_FAILED(); - } - printf(" test %u elmt %u: \n", 1, (unsigned)i); - - printf(" src = "); - for (j=0; j=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto error; - } + /* Get the dataset's data type and make sure it's a named type */ + if ((t2 = H5Dget_type (dset))<0) goto error; + if ((status=H5Tcommitted (t2))<0) goto error; + if (!status) { + H5_FAILED(); + HDputs (" Dataset type should be a named type!"); + goto error; } - if (buf) free(buf); - if (saved_buf) free(saved_buf); + /* Close the dataset, then close its type, then reopen the dataset */ + if (H5Dclose (dset)<0) goto error; + if (H5Tclose (t2)<0) goto error; + if ((dset = H5Dopen (file, "dset1"))<0) goto error; - if(H5Tclose(tid1)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; + /* Get the dataset's type and make sure it's named */ + if ((t2 = H5Dget_type (dset))<0) goto error; + if ((status=H5Tcommitted (t2))<0) goto error; + if (!status) { + H5_FAILED(); + HDputs (" Dataset type should be a named type!"); + goto error; } - if(H5Tclose(tid2)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; + /* + * Close the dataset and create another with the type returned from the + * first dataset. + */ + if (H5Dclose (dset)<0) goto error; + if ((dset=H5Dcreate (file, "dset2", t2, space, H5P_DEFAULT))<0) { + goto error; } - if(H5Pclose(dxpl_id)<0) { - H5_FAILED(); - printf("Can't close property list\n"); - goto error; + /* Reopen the second dataset and make sure the type is shared */ + if (H5Tclose (t2)<0) goto error; + if (H5Dclose (dset)<0) goto error; + if ((dset = H5Dopen (file, "dset2"))<0) goto error; + if ((t2 = H5Dget_type (dset))<0) goto error; + if ((status=H5Tcommitted (t2))<0) goto error; + if (!status) { + H5_FAILED(); + HDputs (" Dataset type should be a named type!"); + goto error; } + if (H5Tclose (t2)<0) goto error; + + /* + * Get the dataset data type by applying H5Tcopy() to the dataset. The + * result should be modifiable. + */ + if ((t2=H5Tcopy (dset))<0) goto error; + if (H5Tset_precision (t2, 256)<0) goto error; + if (H5Tclose (t2)<0) goto error; - if(H5Fclose(file)<0) { - H5_FAILED(); - printf("Can't close file\n"); - goto error; - } /* end if */ - + /* Clean up */ + if (H5Dclose (dset)<0) goto error; + if (H5Tclose (type)<0) goto error; + if (H5Sclose (space)<0) goto error; + if (H5Fclose (file)<0) goto error; PASSED(); - reset_hdf5(); /*print statistics*/ - return 0; error: - if (buf) free(buf); - if (saved_buf) free(saved_buf); - if (aligned) free(aligned); - HDfflush(stdout); H5E_BEGIN_TRY { - H5Tclose (tid1); - H5Tclose (tid2); - H5Pclose (dxpl_id); - H5Fclose (file); + H5Tclose (t2); + H5Tclose (type); + H5Sclose (space); + H5Dclose (dset); + H5Fclose (file); } H5E_END_TRY; - reset_hdf5(); /*print statistics*/ - return MAX((int)fails_this_test, 1); + return 1; } /*------------------------------------------------------------------------- - * Function: test_derived_integer + * Function: mkstr * - * Purpose: Tests user-define and query functions of integer types. + * Purpose: Create a new string data type * - * Return: Success: 0 - * - * Failure: number of errors + * Return: Success: New type + * + * Failure: -1 + * + * Programmer: Robb Matzke + * Monday, August 10, 1998 * - * Programmer: Raymond Lu - * Saturday, Jan 29, 2005 - * * Modifications: * *------------------------------------------------------------------------- */ -static int -test_derived_integer(void) +static hid_t +mkstr(size_t len, H5T_str_t strpad) { - hid_t file=-1, tid1=-1, tid2=-1; - hid_t dxpl_id=-1; - char filename[1024]; - size_t precision, size; - int offset; - H5T_order_t order; - H5T_sign_t sign; - size_t src_size, dst_size; - unsigned char *buf=NULL, *saved_buf=NULL; - int *aligned=NULL; - int endian; /*endianess */ - size_t nelmts = NTESTELEM; - unsigned int fails_this_test = 0; - const size_t max_fails=40; /*max number of failures*/ - char str[256]; /*message string */ - unsigned int i, j; - - TESTING("user-define and query functions of integer types"); - - /* Create File */ - h5_fixname(FILENAME[7], H5P_DEFAULT, filename, sizeof filename); - if((file=H5Fcreate(filename, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT))<0) { - H5_FAILED(); - printf("Can't create file\n"); - goto error; - } - - if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0) { - H5_FAILED(); - printf("Can't create data transfer property list\n"); - goto error; - } - - if((tid1 = H5Tcopy(H5T_STD_I32LE))<0) { - H5_FAILED(); - printf("Can't copy data type\n"); - goto error; - } - - if((tid2 = H5Tcopy(H5T_STD_U64LE))<0) { - H5_FAILED(); - printf("Can't copy data type\n"); - goto error; - } - - /*-------------------------------------------------------------------------- - * 1st integer type - * size=3 byte, precision=24 bits, offset=0 bits, order=big endian. - * It can be illustrated in big-endian order as - * - * 0 1 2 - * SIIIIIII IIIIIIII IIIIIIII - * - * There's no specific order for these functions to define the attributes - * of a new integer type, H5Tset_precision, H5Tset_offset, H5Tset_size, - * H5Tset_order, H5Tset_pad, H5Tset_sign. - *--------------------------------------------------------------------------*/ - if(H5Tset_offset(tid1,0)<0) { - H5_FAILED(); - printf("Can't set offset\n"); - goto error; - } - - if(H5Tset_size(tid1, 3)<0) { - H5_FAILED(); - printf("Can't set size\n"); - goto error; - } - - if(H5Tset_precision(tid1,24)<0) { - H5_FAILED(); - printf("Can't set precision\n"); - goto error; - } - - if(H5Tset_order(tid1, H5T_ORDER_BE)<0) { - H5_FAILED(); - printf("Can't set order\n"); - goto error; - } - - if(H5Tcommit(file, "new integer type 1", tid1)<0) { - H5_FAILED(); - printf("Can't commit data type\n"); - goto error; - } - - if(H5Tclose(tid1)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; - } - - if((tid1 = H5Topen(file, "new integer type 1"))<0) { - H5_FAILED(); - printf("Can't open datatype\n"); - goto error; - } - - if((precision = H5Tget_precision(tid1))!=24) { - H5_FAILED(); - printf("Can't get precision or wrong precision\n"); - goto error; - } - if((offset = H5Tget_offset(tid1))!=0) { - H5_FAILED(); - printf("Can't get offset or wrong offset\n"); - goto error; - } - if((size = H5Tget_size(tid1))!=3) { - H5_FAILED(); - printf("Can't get size or wrong size\n"); - goto error; - } - if((order = H5Tget_order(tid1))!=H5T_ORDER_BE) { - H5_FAILED(); - printf("Can't get order or wrong order\n"); - goto error; - } - - /*-------------------------------------------------------------------------- - * 2nd integer type - * size=8 byte, precision=48 bits, offset=10 bits, order=little endian. - * It can be illustrated in little-endian order as - * - * 7 6 5 4 3 2 1 0 - * ??????SI IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIII?? ???????? - *--------------------------------------------------------------------------*/ - if(H5Tset_precision(tid2,48)<0) { - H5_FAILED(); - printf("Can't set precision\n"); - goto error; - } - - if(H5Tset_offset(tid2,10)<0) { - H5_FAILED(); - printf("Can't set offset\n"); - goto error; - } - - if(H5Tset_sign(tid2,H5T_SGN_2)<0) { - H5_FAILED(); - printf("Can't set offset\n"); - goto error; - } - - if(H5Tcommit(file, "new integer type 2", tid2)<0) { - H5_FAILED(); - printf("Can't commit data type\n"); - goto error; - } - - if(H5Tclose(tid2)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; - } - - if((tid2 = H5Topen(file, "new integer type 2"))<0) { - H5_FAILED(); - printf("Can't open datatype\n"); - goto error; - } - - if((precision = H5Tget_precision(tid2))!=48) { - H5_FAILED(); - printf("Can't get precision or wrong precision\n"); - goto error; - } - if((offset = H5Tget_offset(tid2))!=10) { - H5_FAILED(); - printf("Can't get offset or wrong offset\n"); - goto error; - } - if((size = H5Tget_size(tid2))!=8) { - H5_FAILED(); - printf("Can't get size or wrong size\n"); - goto error; - } - if((sign = H5Tget_sign(tid2))!=H5T_SGN_2) { - H5_FAILED(); - printf("Can't get sign or wrong sign\n"); - goto error; - } - - /* Convert data from the 1st to the 2nd derived integer type. - * Then convert data from the 2nd type back to the 1st type. - * Compare the final data with the original data. - */ - src_size = H5Tget_size(tid1); - dst_size = H5Tget_size(tid2); - endian = H5Tget_order(tid1); - buf = (unsigned char*)malloc(nelmts*(MAX(src_size, dst_size))); - saved_buf = (unsigned char*)malloc(nelmts*src_size); - - for(i=0; i derived integer conversions"); - printf("%-70s", str); - HDfflush(stdout); - H5_FAILED(); - } - printf(" test %u elmt %u: \n", 1, (unsigned)i); - - printf(" src = "); - for (j=0; j=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto error; - } - } - - if(H5Tclose(tid1)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; - } - - if(H5Tclose(tid2)<0) { - H5_FAILED(); - printf("Can't close datatype\n"); - goto error; - } - - if(H5Pclose(dxpl_id)<0) { - H5_FAILED(); - printf("Can't close property list\n"); - goto error; - } - - if(H5Fclose(file)<0) { - H5_FAILED(); - printf("Can't close file\n"); - goto error; - } /* end if */ - - PASSED(); - reset_hdf5(); /*print statistics*/ - - return 0; + hid_t t; - error: - if (buf) free(buf); - if (saved_buf) free(saved_buf); - if (aligned) free(aligned); - HDfflush(stdout); - H5E_BEGIN_TRY { - H5Tclose (tid1); - H5Tclose (tid2); - H5Pclose (dxpl_id); - H5Fclose (file); - } H5E_END_TRY; - reset_hdf5(); /*print statistics*/ - return MAX((int)fails_this_test, 1); + if ((t=H5Tcopy(H5T_C_S1))<0) return -1; + if (H5Tset_size(t, len)<0) return -1; + if (H5Tset_strpad(t, strpad)<0) return -1; + return t; } /*------------------------------------------------------------------------- - * Function: test_transient + * Function: test_conv_str_1 * - * Purpose: Tests transient data types. + * Purpose: Test string conversions * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke - * Thursday, June 4, 1998 + * Monday, August 10, 1998 * * Modifications: * *------------------------------------------------------------------------- */ static int -test_transient (hid_t fapl) +test_conv_str_1(void) { - static hsize_t ds_size[2] = {10, 20}; - hid_t file=-1, type=-1, space=-1, dset=-1, t2=-1; - char filename[1024]; - herr_t status; - - TESTING("transient data types"); - - h5_fixname(FILENAME[0], fapl, filename, sizeof filename); - if ((file=H5Fcreate (filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl))<0) { - goto error; - } - if ((space = H5Screate_simple (2, ds_size, ds_size))<0) goto error; - - /* Predefined types cannot be modified or closed */ - H5E_BEGIN_TRY { - status = H5Tset_precision (H5T_NATIVE_INT, 256); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Predefined types should not be modifiable!"); - goto error; - } - H5E_BEGIN_TRY { - status = H5Tclose (H5T_NATIVE_INT); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Predefined types should not be closable!"); - goto error; - } - - /* Copying a predefined type results in a modifiable copy */ - if ((type=H5Tcopy (H5T_NATIVE_INT))<0) goto error; - if (H5Tset_precision (type, 256)<0) goto error; - - /* It should not be possible to create an attribute for a transient type */ - H5E_BEGIN_TRY { - status = H5Acreate (type, "attr1", H5T_NATIVE_INT, space, H5P_DEFAULT); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Attributes should not be allowed for transient types!"); - goto error; - } - - /* Create a dataset from a transient data type */ - if (H5Tclose (type)<0) goto error; - if ((type = H5Tcopy (H5T_NATIVE_INT))<0) goto error; - if ((dset=H5Dcreate (file, "dset1", type, space, H5P_DEFAULT))<0) - goto error; - - /* The type returned from a dataset should not be modifiable */ - if ((t2 = H5Dget_type (dset))<0) goto error; - H5E_BEGIN_TRY { - status = H5Tset_precision (t2, 256); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Dataset data types should not be modifiable!"); - goto error; - } - if (H5Tclose (t2)<0) goto error; + char *buf=NULL; + hid_t src_type, dst_type; - /* - * Close the dataset and reopen it, testing that it's type is still - * read-only. - */ - if (H5Dclose (dset)<0) goto error; - if ((dset=H5Dopen (file, "dset1"))<0) goto error; - if ((t2 = H5Dget_type (dset))<0) goto error; - H5E_BEGIN_TRY { - status = H5Tset_precision (t2, 256); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Dataset data types should not be modifiable!"); - goto error; - } - if (H5Tclose (t2)<0) goto error; + TESTING("string conversions"); /* - * Get the dataset data type by applying H5Tcopy() to the dataset. The - * result should be modifiable. - */ - if ((t2=H5Tcopy (dset))<0) goto error; - if (H5Tset_precision (t2, 256)<0) goto error; - if (H5Tclose (t2)<0) goto error; - - - H5Dclose (dset); - H5Fclose (file); - H5Tclose (type); - H5Sclose (space); - PASSED(); - return 0; - - error: - H5E_BEGIN_TRY { - H5Tclose (t2); - H5Tclose (type); - H5Sclose (space); - H5Dclose (dset); - H5Fclose (file); - } H5E_END_TRY; - return 1; -} - - -/*------------------------------------------------------------------------- - * Function: test_named - * - * Purpose: Tests named data types. - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Monday, June 1, 1998 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_named (hid_t fapl) -{ - hid_t file=-1, type=-1, space=-1, dset=-1, t2=-1, attr1=-1; - herr_t status; - static hsize_t ds_size[2] = {10, 20}; - hsize_t i,j; - unsigned attr_data[10][20]; - char filename[1024]; - - TESTING("named data types"); - - h5_fixname(FILENAME[1], fapl, filename, sizeof filename); - if ((file=H5Fcreate (filename, H5F_ACC_TRUNC, H5P_DEFAULT, fapl))<0) { - goto error; - } - if ((space = H5Screate_simple (2, ds_size, ds_size))<0) goto error; - - /* Predefined types cannot be committed */ - H5E_BEGIN_TRY { - status = H5Tcommit (file, "test_named_1 (should not exist)", - H5T_NATIVE_INT); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Predefined types should not be committable!"); - goto error; - } - - /* Copy a predefined data type and commit the copy */ - if ((type = H5Tcopy (H5T_NATIVE_INT))<0) goto error; - if (H5Tcommit (file, "native-int", type)<0) goto error; - if ((status=H5Tcommitted (type))<0) goto error; - if (0==status) { - H5_FAILED(); - HDputs (" H5Tcommitted() returned false!"); - goto error; - } - - /* We should not be able to modify a type after it has been committed. */ - H5E_BEGIN_TRY { - status = H5Tset_precision (type, 256); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Committed type is not constant!"); - goto error; - } - - /* We should not be able to re-commit a committed type */ - H5E_BEGIN_TRY { - status = H5Tcommit(file, "test_named_2 (should not exist)", type); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - HDputs (" Committed types should not be recommitted!"); - goto error; - } - - /* It should be possible to define an attribute for the named type */ - if ((attr1=H5Acreate (type, "attr1", H5T_NATIVE_UCHAR, space, - H5P_DEFAULT))<0) goto error; - for (i=0; i1) { - sprintf(s, "Testing random string conversion speed (test %d/%d)", - (int)(i+1), (int)ntests); - } else { - sprintf(s, "Testing random string conversion speed"); - } - printf("%-70s", s); - HDfflush(stdout); - if (H5Tconvert(c_type, f_type, nelmts, buf, NULL, H5P_DEFAULT)<0) - goto error; - if (H5Tconvert(f_type, c_type, nelmts, buf, NULL, H5P_DEFAULT)<0) - goto error; - PASSED(); - } - ret_value = 0; - - error: - if (buf) HDfree(buf); - reset_hdf5(); - return ret_value; -} - - -/*------------------------------------------------------------------------- - * Function: test_conv_enum_1 - * - * Purpose: Test conversion speed for enum data types - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Tuesday, January 5, 1999 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_conv_enum_1(void) -{ - const size_t nelmts=NTESTELEM; - const int ntests=NTESTS; - int i, val, *buf=NULL; - hid_t t1, t2; - char s[80]; - int ret_value = 1; - - /* Build the data types */ - t1 = H5Tcreate(H5T_ENUM, sizeof(int)); - t2 = H5Tenum_create(H5T_NATIVE_INT); - s[1] = '\0'; - for (i=0; i<26; i++) { - s[0] = 'A'+i; - H5Tenum_insert(t1, s, &i); - H5Tenum_insert(t2, s, (val=i*1000+i, &val)); - } - - /* Initialize the buffer */ - buf = (int*)HDmalloc(nelmts*MAX(H5Tget_size(t1), H5Tget_size(t2))); - for (i=0; i<(int)nelmts; i++) - buf[i] = HDrand() % 26; - - /* Conversions */ - for (i=0; i1) { - sprintf(s, "Testing random enum conversion O(N) (test %d/%d)", - i+1, ntests); - } else { - sprintf(s, "Testing random enum conversion O(N)"); - } - printf("%-70s", s); - HDfflush(stdout); - if (H5Tconvert(t1, t2, nelmts, buf, NULL, H5P_DEFAULT)<0) goto error; - PASSED(); - } - - for (i=0; i1) { - sprintf(s, "Testing random enum conversion O(N log N) " - "(test %d/%d)", i+1, ntests); - } else { - sprintf(s, "Testing random enum conversion O(N log N)"); - } - printf("%-70s", s); - HDfflush(stdout); - if (H5Tconvert(t2, t1, nelmts, buf, NULL, H5P_DEFAULT)<0) goto error; - PASSED(); - } - ret_value = 0; - - error: - H5Tclose(t1); - H5Tclose(t2); - if (buf) HDfree(buf); - reset_hdf5(); - return ret_value; -} - - -/*------------------------------------------------------------------------- - * Function: test_conv_enum_2 - * - * Purpose: Tests enumeration conversions where source isn't a native type. - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke, LLNL, 2003-06-09 - * - * Modifications: - *------------------------------------------------------------------------- - */ -static int -test_conv_enum_2(void) -{ - hid_t srctype=-1, dsttype=-1, oddsize=-1; - int *data=NULL, i, nerrors=0; - const char *mname[] = { "RED", - "GREEN", - "BLUE", - "YELLOW", - "PINK", - "PURPLE", - "ORANGE", - "WHITE" }; - - TESTING("non-native enumeration type conversion"); - - /* Source enum type */ - oddsize = H5Tcopy(H5T_STD_I32BE); - H5Tset_size(oddsize, 3); /*reduce to 24 bits, not corresponding to any native size*/ - srctype = H5Tenum_create(oddsize); - for (i=7; i>=0; --i) { - char pattern[3]; - pattern[2] = i; - pattern[0] = pattern[1] = 0; - H5Tenum_insert(srctype, mname[i], pattern); - } - - /* Destination enum type */ - dsttype = H5Tenum_create(H5T_NATIVE_INT); - assert(H5Tget_size(dsttype)>H5Tget_size(srctype)); - for (i=0; i<8; i++) - H5Tenum_insert(dsttype, mname[i], &i); - - /* Source data */ - data = (int*)malloc(NTESTELEM*sizeof(int)); - for (i=0; icommand) num_opaque_conversions_g++; - return 0; -} - - -/*------------------------------------------------------------------------- - * Function: test_opaque - * - * Purpose: Driver function to test opaque datatypes - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Raymond Lu - * June 2, 2004 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_opaque(void) -{ - int num_errors = 0; - - TESTING("opaque datatypes"); - - /* Test opaque types with tags */ - num_errors += opaque_check(0); - /* Test opaque types without tag */ - num_errors += opaque_check(1); - - if(num_errors) - goto error; - - PASSED(); - return 0; - - error: - H5_FAILED(); - return num_errors; -} - - -/*------------------------------------------------------------------------- - * Function: opaque_check - * - * Purpose: Test opaque datatypes - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Thursday, May 20, 1999 - * - * Modifications: - * Raymond Lu - * June 2, 2004 - * Made tag for one opaque type optional. - * - *------------------------------------------------------------------------- - */ -static int -opaque_check(int tag_it) -{ -#define OPAQUE_NELMTS 1000 - hid_t st=-1, dt=-1; - herr_t status; - char buf[1]; /*not really used*/ - int saved; - - saved = num_opaque_conversions_g = 0; - - /* Build source and destination types */ - if ((st=H5Tcreate(H5T_OPAQUE, 4))<0) goto error; - if (H5Tset_tag(st, "opaque source type")<0) goto error; - - if ((dt=H5Tcreate(H5T_OPAQUE, 4))<0) goto error; - if (tag_it) { - if (H5Tset_tag(dt, "opaque destination type")<0) - goto error; - } - - /* Make sure that we can't convert between the types yet */ - H5E_BEGIN_TRY { - status = H5Tconvert(st, dt, OPAQUE_NELMTS, buf, NULL, H5P_DEFAULT); - } H5E_END_TRY; - if (status>=0) { - H5_FAILED(); - printf(" opaque conversion should have failed but succeeded\n"); - goto error; - } - - /* Register a conversion function */ - if (H5Tregister(H5T_PERS_HARD, "o_test", st, dt, convert_opaque)<0) - goto error; - - /* Try the conversion again, this time it should work */ - if (H5Tconvert(st, dt, OPAQUE_NELMTS, buf, NULL, H5P_DEFAULT)<0) goto error; - if (saved+1 != num_opaque_conversions_g) { - H5_FAILED(); - printf(" unexpected number of opaque conversions\n"); - goto error; - } - - /* Unregister conversion function */ - if (H5Tunregister(H5T_PERS_HARD, "o_test", st, dt, convert_opaque)<0) - goto error; - - H5Tclose(st); - H5Tclose(dt); - return 0; - - error: - if (st>0) H5Tclose(st); - if (dt>0) H5Tclose(dt); - H5_FAILED(); - return 1; -} - - -/*------------------------------------------------------------------------- - * Function: test_conv_int - * - * Purpose: Test atomic number conversions. - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Wednesday, June 10, 1998 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_conv_int (void) -{ - unsigned char byte[4]; - - /*--------------------------------------------------------------------- - * Test some specific overflow/underflow cases. - *--------------------------------------------------------------------- - */ - TESTING("integer overflow conversions"); - - /* (unsigned)0x80000000 -> (unsigned)0xffff */ - byte[0] = byte[1] = byte[2] = 0; - byte[3] = 0x80; - if (H5Tconvert (H5T_STD_U32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd\0abcd\0abcdefghi\0", 20)) { + H5_FAILED(); + HDputs(" Truncated C-string test failed"); goto error; } - if (byte[0]!=0xff || byte[1]!=0xff) { + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd\0\0\0\0\0\0abcd\0\0\0\0\0\0", 20)) { H5_FAILED(); - printf(" src: 0x80000000 unsigned\n"); - printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]); - printf(" ans: 0xffff unsigned\n"); + HDputs(" Extended C-string test failed"); goto error; } + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; - /* (unsigned)0xffffffff -> (signed)0x7fff */ - byte[0] = byte[1] = byte[2] = byte[3] = 0xff; - if (H5Tconvert (H5T_STD_U32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + /* + * Convert a null padded string to a shorter and then longer string. + */ + src_type = mkstr(10, H5T_STR_NULLPAD); + dst_type = mkstr(5, H5T_STR_NULLPAD); + buf = (char*)HDcalloc(2, 10); + HDmemcpy(buf, "abcdefghijabcdefghij", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdeabcdeabcdefghij", 20)) { + H5_FAILED(); + HDputs(" Truncated C buffer test failed"); goto error; } - if (byte[0]!=0xff || byte[1]!=0x7f) { + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcde\0\0\0\0\0abcde\0\0\0\0\0", 20)) { H5_FAILED(); - printf(" src: 0xffffffff unsigned\n"); - printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]); - printf(" ans: 0x7fff signed\n"); + HDputs(" Extended C buffer test failed"); goto error; } + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; - /* (signed)0xffffffff -> (unsigned)0x0000 */ - byte[0] = byte[1] = byte[2] = byte[3] = 0xff; - if (H5Tconvert (H5T_STD_I32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + /* + * Convert a space-padded string to a shorter and then longer string. + */ + src_type = mkstr(10, H5T_STR_SPACEPAD); + dst_type = mkstr(5, H5T_STR_SPACEPAD); + buf = (char*)HDcalloc(2, 10); + HDmemcpy(buf, "abcdefghijabcdefghij", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdeabcdeabcdefghij", 20)) { + H5_FAILED(); + HDputs(" Truncated Fortran-string test failed"); goto error; } - if (byte[0]!=0x00 || byte[1]!=0x00) { + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcde abcde ", 20)) { H5_FAILED(); - printf(" src: 0xffffffff signed\n"); - printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]); - printf(" ans: 0x0000 unsigned\n"); + HDputs(" Extended Fortran-string test failed"); goto error; } + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; - /* (signed)0x7fffffff -> (unsigned)0xffff */ - byte[0] = byte[1] = byte[2] = 0xff; - byte[3] = 0x7f; - if (H5Tconvert (H5T_STD_I32LE, H5T_STD_U16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + /* + * What happens if a null-terminated string is not null terminated? If + * the conversion is to an identical string then nothing happens but if + * the destination is a different size or type of string then the right + * thing should happen. + */ + src_type = mkstr(10, H5T_STR_NULLTERM); + dst_type = mkstr(10, H5T_STR_NULLTERM); + buf = (char*)HDcalloc(2, 10); + HDmemcpy(buf, "abcdefghijabcdefghij", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdefghijabcdefghij", 20)) { + H5_FAILED(); + HDputs(" Non-terminated string test 1"); goto error; } - if (byte[0]!=0xff || byte[1]!=0xff) { + H5Tclose(dst_type); + dst_type = mkstr(5, H5T_STR_NULLTERM); + HDmemcpy(buf, "abcdefghijabcdefghij", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd\0abcd\0abcdefghij", 20)) { H5_FAILED(); - printf(" src: 0x7fffffff signed\n"); - printf(" dst: 0x%02x%02x unsigned\n", byte[1], byte[0]); - printf(" ans: 0xffff unsigned\n"); + HDputs(" Non-terminated string test 2"); goto error; } - - /* (signed)0x7fffffff -> (signed)0x7fff */ - byte[0] = byte[1] = byte[2] = 0xff; - byte[3] = 0x7f; - if (H5Tconvert (H5T_STD_I32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + HDmemcpy(buf, "abcdeabcdexxxxxxxxxx", 20); + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcde\0\0\0\0\0abcde\0\0\0\0\0", 20)) { + H5_FAILED(); + HDputs(" Non-terminated string test 2"); goto error; } - if (byte[0]!=0xff || byte[1]!=0x7f) { + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; + + /* + * Test C string to Fortran and vice versa. + */ + src_type = mkstr(10, H5T_STR_NULLTERM); + dst_type = mkstr(10, H5T_STR_SPACEPAD); + buf = (char*)HDcalloc(2, 10); + HDmemcpy(buf, "abcdefghi\0abcdefghi\0", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdefghi abcdefghi ", 20)) { H5_FAILED(); - printf(" src: 0x7fffffff signed\n"); - printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]); - printf(" ans: 0x7fff signed\n"); + HDputs(" C string to Fortran test 1"); goto error; } - - /* (signed)0xbfffffff -> (signed)0x8000 */ - byte[0] = byte[1] = byte[2] = 0xff; - byte[3] = 0xbf; - if (H5Tconvert (H5T_STD_I32LE, H5T_STD_I16LE, 1, byte, NULL, H5P_DEFAULT)<0) { + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdefghi\0abcdefghi\0", 20)) { + H5_FAILED(); + HDputs(" Fortran to C string test 1"); goto error; } - if (byte[0]!=0x00 || byte[1]!=0x80) { + if (H5Tclose(dst_type)<0) goto error; + dst_type = mkstr(5, H5T_STR_SPACEPAD); + HDmemcpy(buf, "abcdefgh\0\0abcdefgh\0\0", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdeabcdeabcdefgh\0\0", 20)) { H5_FAILED(); - printf(" src: 0xbfffffff signed\n"); - printf(" dst: 0x%02x%02x signed\n", byte[1], byte[0]); - printf(" ans: 0x8000 signed\n"); + HDputs(" C string to Fortran test 2"); goto error; } - - PASSED(); - reset_hdf5(); - return 0; - - error: - reset_hdf5(); - return 1; -} - - -/*------------------------------------------------------------------------- - * Function: test_conv_int_1 - * - * Purpose: Test conversion of random integer values from SRC to DST. - * These types should be any combination of: - * - * H5T_NATIVE_SCHAR H5T_NATIVE_UCHAR - * H5T_NATIVE_SHORT H5T_NATIVE_USHORT - * H5T_NATIVE_INT H5T_NATIVE_UINT - * H5T_NATIVE_LONG H5T_NATIVE_ULONG - * H5T_NATIVE_LLONG H5T_NATIVE_ULLONG - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Monday, November 16, 1998 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_conv_int_1(const char *name, hid_t src, hid_t dst) -{ - const size_t ntests=NTESTS; /*number of tests */ - size_t nelmts=0; /*num values per test */ - const size_t max_fails=8; /*max number of failures*/ - size_t fails_all_tests=0; /*number of failures */ - size_t fails_this_test; /*fails for this test */ - char str[256]; /*hello string */ - dtype_t src_type, dst_type; /*data types */ - const char *src_type_name=NULL; /*source type name */ - const char *dst_type_name=NULL; /*destination type name */ - int endian; /*machine endianess */ - size_t src_size, dst_size; /*type sizes */ - unsigned char *buf=NULL; /*buffer for conversion */ - unsigned char *saved=NULL; /*original values */ - size_t i, j, k; /*counters */ - unsigned char *hw=NULL; /*hardware conv result */ - unsigned char src_bits[32]; /*src value in LE order */ - unsigned char dst_bits[32]; /*dest value in LE order*/ - size_t src_nbits; /*source length in bits */ - size_t dst_nbits; /*dst length in bits */ - H5T_sign_t src_sign; /*source sign type */ - H5T_sign_t dst_sign; /*dst sign type */ - void *aligned=NULL; /*aligned temp buffer */ - signed char hw_char; - unsigned char hw_uchar; - short hw_short; - unsigned short hw_ushort; - int hw_int; - unsigned hw_uint; - long hw_long; - unsigned long hw_ulong; - long_long hw_llong; - unsigned long_long hw_ullong; - - /* What are the names of the source and destination types */ - if (H5Tequal(src, H5T_NATIVE_SCHAR)) { - src_type_name = "signed char"; - src_type = INT_SCHAR; - } else if (H5Tequal(src, H5T_NATIVE_UCHAR)) { - src_type_name = "unsigned char"; - src_type = INT_UCHAR; - } else if (H5Tequal(src, H5T_NATIVE_SHORT)) { - src_type_name = "short"; - src_type = INT_SHORT; - } else if (H5Tequal(src, H5T_NATIVE_USHORT)) { - src_type_name = "unsigned short"; - src_type = INT_USHORT; - } else if (H5Tequal(src, H5T_NATIVE_INT)) { - src_type_name = "int"; - src_type = INT_INT; - } else if (H5Tequal(src, H5T_NATIVE_UINT)) { - src_type_name = "unsigned int"; - src_type = INT_UINT; - } else if (H5Tequal(src, H5T_NATIVE_LONG)) { - src_type_name = "long"; - src_type = INT_LONG; - } else if (H5Tequal(src, H5T_NATIVE_ULONG)) { - src_type_name = "unsigned long"; - src_type = INT_ULONG; - } else if (H5Tequal(src, H5T_NATIVE_LLONG)) { - src_type_name = "long long"; - src_type = INT_LLONG; - } else if (H5Tequal(src, H5T_NATIVE_ULLONG)) { - src_type_name = "unsigned long long"; - src_type = INT_ULLONG; - } else { - src_type_name = "UNKNOWN"; - src_type = OTHER; + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcde\0\0\0\0\0abcde\0\0\0\0\0", 20)) { + H5_FAILED(); + HDputs(" Fortran to C string test 2"); + goto error; } - - 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; + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; + src_type = mkstr(5, H5T_STR_NULLTERM); + dst_type = mkstr(10, H5T_STR_SPACEPAD); + HDmemcpy(buf, "abcd\0abcd\0xxxxxxxxxx", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd abcd ", 20)) { + H5_FAILED(); + HDputs(" C string to Fortran test 3"); + goto error; } - - /* Sanity checks */ - if (OTHER==src_type || OTHER==dst_type) { - sprintf(str, "Testing random %s %s -> %s conversions", - name, src_type_name, dst_type_name); - printf("%-70s", str); + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd\0abcd\0abcd ", 20)) { H5_FAILED(); - HDputs(" Unknown data type."); + HDputs(" Fortran to C string test 3"); goto error; } - + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; - /* Some information about datatypes */ - endian = H5Tget_order(H5T_NATIVE_INT); - src_size = H5Tget_size(src); - dst_size = H5Tget_size(dst); - src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */ - dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */ - src_sign = H5Tget_sign(src); - dst_sign = H5Tget_sign(dst); - aligned = HDmalloc(sizeof(long_long)); -#ifdef SHOW_OVERFLOWS - noverflows_g = 0; -#endif - - /* Allocate and initialize the source buffer through macro INIT_INTEGER. The BUF - * will be used for the conversion while the SAVED buffer will be - * used for the comparison later. + /* + * Test C buffer to Fortran and vice versa. */ - if(src_type == INT_SCHAR) { - INIT_INTEGER(signed char, SCHAR_MAX, SCHAR_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_UCHAR) { - INIT_INTEGER(unsigned char, UCHAR_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_SHORT) { - INIT_INTEGER(short, SHRT_MAX, SHRT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_USHORT) { - INIT_INTEGER(unsigned short, USHRT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_INT) { - INIT_INTEGER(int, INT_MAX, INT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_UINT) { - INIT_INTEGER(unsigned int, UINT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_LONG) { - INIT_INTEGER(long, LONG_MAX, LONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_ULONG) { - INIT_INTEGER(unsigned long, ULONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_LLONG) { - INIT_INTEGER(long_long, LLONG_MAX, LLONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_ULLONG) { - INIT_INTEGER(unsigned long_long, ULLONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else - goto error; - - /* The tests */ - for (i=0; i1) { - sprintf(str, "Testing random %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); - } else { - sprintf(str, "Testing random %s %s -> %s conversions", - name, src_type_name, dst_type_name); - } - printf("%-70s", str); - HDfflush(stdout); - fails_this_test=0; - - /* Perform the conversion */ - if (H5Tconvert(src, dst, nelmts, buf, NULL, H5P_DEFAULT)<0) - goto error; - - /* Check the results from the library against hardware */ - for (j=0; jdst_nbits) { - if(0==H5T_bit_get_d(src_bits, src_nbits-1, 1) && - H5T_bit_find(src_bits, dst_nbits-1, (src_nbits-dst_nbits), - H5T_BIT_MSB, 1)>=0) { - /* - * Source is positive and the magnitude is too large for - * the destination. The destination should be set to the - * maximum possible value: 0x7f...f - */ - if (0==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) && - H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 0)<0) - continue; /*no error*/ - } else if (1==H5T_bit_get_d(src_bits, src_nbits-1, 1) && - H5T_bit_find(src_bits, 0, src_nbits-1, H5T_BIT_MSB, - 0)+1>=(ssize_t)dst_nbits) { - /* - * Source is negative but the magnitude is too large for - * the destination. The destination should be set to the - * smallest possible value: 0x80...0 - */ - if (1==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) && - H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 1)<0) - continue; /*no error*/ - } - } else if(src_nbitsdst_nbits && - H5T_bit_find(src_bits, dst_nbits-1, - src_nbits-dst_nbits, H5T_BIT_LSB, 1)>=0) { - /* - * The source is a value with a magnitude too large for - * the destination. The destination should be the - * largest possible value: 0xff...f - */ - if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 0)<0) - continue; /*no error*/ - } - } else if (H5T_SGN_NONE==src_sign && H5T_SGN_2==dst_sign) { - if (src_nbits>=dst_nbits && - H5T_bit_find(src_bits, dst_nbits-1, (src_nbits-dst_nbits)+1, - H5T_BIT_LSB, 1)>=0) { - /* - * The source value has a magnitude that is larger than - * the destination can handle. The destination should be - * set to the largest possible positive value: 0x7f...f - */ - if (0==H5T_bit_get_d(dst_bits, dst_nbits-1, 1) && - H5T_bit_find(dst_bits, 0, dst_nbits-1, H5T_BIT_LSB, 0)<0) - continue; /*no error*/ - } - } else { - if (src_nbits>dst_nbits && - H5T_bit_find(src_bits, dst_nbits, src_nbits-dst_nbits, - H5T_BIT_LSB, 1)>=0) { - /* - * The unsigned source has a value which is too large for - * the unsigned destination. The destination should be - * set to the largest possible value: 0xff...f - */ - if (H5T_bit_find(dst_bits, 0, dst_nbits, H5T_BIT_LSB, 0)<0) - continue; /*no error*/ - } - } - - /* Print errors */ - if (0==fails_this_test++) - H5_FAILED(); - printf(" test %u elmt %u\n", (unsigned)i+1, (unsigned)j); - - printf(" src = "); - for (k=0; k=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto done; - } - } - PASSED(); + src_type = mkstr(10, H5T_STR_NULLPAD); + dst_type = mkstr(10, H5T_STR_SPACEPAD); + buf = (char*)HDcalloc(2, 10); + HDmemcpy(buf, "abcdefghijabcdefghij", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdefghijabcdefghij", 20)) { + H5_FAILED(); + HDputs(" C buffer to Fortran test 1"); + goto error; } -#ifdef SHOW_OVERFLOWS - if (noverflows_g>0) { - printf(" %d overflow%s in previous test\n", - noverflows_g, 1==noverflows_g?"":"s"); + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdefghijabcdefghij", 20)) { + H5_FAILED(); + HDputs(" Fortran to C buffer test 1"); + goto error; } -#endif - -done: - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); - reset_hdf5(); /*print statistics*/ - return (int)fails_all_tests; - -error: - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); - reset_hdf5(); /*print statistics*/ - return MAX((int)fails_all_tests, 1); -} - - -/*------------------------------------------------------------------------- - * Function: test_conv_int_2 - * - * Purpose: Tests overlap calculates in H5T_conv_i_i(), which should be - * the same as for H5T_conv_f_f() and H5T_conv_s_s(). - * - * Return: Success: 0 - * - * Failure: number of errors - * - * Programmer: Robb Matzke - * Friday, April 30, 1999 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -test_conv_int_2(void) -{ - int i, j; - hid_t src_type, dst_type; - char buf[32*100]; - - printf("%-70s", "Testing overlap calculations"); - HDfflush(stdout); - - HDmemset(buf, 0, sizeof buf); - for (i=1; i<=32; i++) { - for (j=1; j<=32; j++) { - - /* Source type */ - src_type = H5Tcopy(H5T_NATIVE_CHAR); - H5Tset_size(src_type, (size_t)i); - - /* Destination type */ - dst_type = H5Tcopy(H5T_NATIVE_CHAR); - H5Tset_size(dst_type, (size_t)j); - - /* - * Conversion. If overlap calculations aren't right then an - * assertion will fail in H5T_conv_i_i() - */ - H5Tconvert(src_type, dst_type, 100, buf, NULL, H5P_DEFAULT); - H5Tclose(src_type); - H5Tclose(dst_type); - } + if (H5Tclose(dst_type)<0) goto error; + dst_type = mkstr(5, H5T_STR_SPACEPAD); + HDmemcpy(buf, "abcdefgh\0\0abcdefgh\0\0", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcdeabcdeabcdefgh\0\0", 20)) { + H5_FAILED(); + HDputs(" C buffer to Fortran test 2"); + goto error; } - PASSED(); - return 0; -} - - -/*------------------------------------------------------------------------- - * Function: my_isnan - * - * Purpose: Determines whether VAL points to NaN. - * - * Return: TRUE or FALSE - * - * Programmer: Robb Matzke - * Monday, July 6, 1998 - * - * Modifications: - * - *------------------------------------------------------------------------- - */ -static int -my_isnan(dtype_t type, void *val) -{ - int retval; - char s[256]; - - if (FLT_FLOAT==type) { - float x; - HDmemcpy(&x, val, sizeof(float)); - retval = (x!=x); - } else if (FLT_DOUBLE==type) { - double x; - HDmemcpy(&x, val, sizeof(double)); - retval = (x!=x); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (FLT_LDOUBLE==type) { - long double x; - HDmemcpy(&x, val, sizeof(long double)); - retval = (x!=x); -#endif - } else { - return 0; + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcde\0\0\0\0\0abcde\0\0\0\0\0", 20)) { + H5_FAILED(); + HDputs(" Fortran to C buffer test 2"); + goto error; } - - /* - * Sometimes NaN==NaN (e.g., DEC Alpha) so we try to print it and see if - * the result contains a NaN string. - */ - if (!retval) { - if (FLT_FLOAT==type) { - float x; - HDmemcpy(&x, val, sizeof(float)); - sprintf(s, "%g", x); - } else if (FLT_DOUBLE==type) { - double x; - HDmemcpy(&x, val, sizeof(double)); - sprintf(s, "%g", x); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (FLT_LDOUBLE==type) { - long double x; - HDmemcpy(&x, val, sizeof(long double)); - sprintf(s, "%Lg", x); -#endif - } else { - return 0; - } - if (HDstrstr(s, "NaN") || HDstrstr(s, "NAN") || HDstrstr(s, "nan")) - retval = 1; + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; + src_type = mkstr(5, H5T_STR_NULLPAD); + dst_type = mkstr(10, H5T_STR_SPACEPAD); + HDmemcpy(buf, "abcd\0abcd\0xxxxxxxxxx", 20); + if (H5Tconvert(src_type, dst_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd abcd ", 20)) { + H5_FAILED(); + HDputs(" C buffer to Fortran test 3"); + goto error; + } + if (H5Tconvert(dst_type, src_type, 2, buf, NULL, H5P_DEFAULT)<0) goto error; + if (HDmemcmp(buf, "abcd\0abcd\0abcd ", 20)) { + H5_FAILED(); + HDputs(" Fortran to C buffer test 3"); + goto error; } + HDfree(buf); + if (H5Tclose(src_type)<0) goto error; + if (H5Tclose(dst_type)<0) goto error; + + PASSED(); + reset_hdf5(); + return 0; - return retval; + error: + reset_hdf5(); + return 1; } /*------------------------------------------------------------------------- - * Function: test_conv_flt_1 + * Function: test_conv_str_2 * - * Purpose: Test conversion of random floating point values from SRC to - * DST. These types should be H5T_NATIVE_FLOAT, - * H5T_NATIVE_DOUBLE, or H5T_NATIVE_LDOUBLE. + * Purpose: Tests C-to-Fortran and Fortran-to-C string conversion speed. * * Return: Success: 0 * * Failure: number of errors * * Programmer: Robb Matzke - * Tuesday, June 23, 1998 + * Monday, August 10, 1998 * * Modifications: - * Albert Cheng, Apr 16, 2004 - * Check for underflow condition. If the src number is - * smaller than the dst MIN float number, consider it okay - * if the converted sw and hw dst are both less than or - * equal to the dst MIN float number. * *------------------------------------------------------------------------- */ static int -test_conv_flt_1 (const char *name, hbool_t run_special, hid_t src, hid_t dst) +test_conv_str_2(void) { - dtype_t src_type, dst_type; /*data types */ - const size_t ntests=NTESTS; /*number of tests */ - size_t nelmts=0; /*num values per test */ - const size_t max_fails=8; /*max number of failures*/ - size_t fails_all_tests=0; /*number of failures */ - size_t fails_this_test; /*fails for this test */ - const char *src_type_name = NULL; /*source type name */ - const char *dst_type_name = NULL; /*destination type name */ - size_t src_size, dst_size; /*type sizes */ - unsigned char *buf = NULL; /*buffer for conversion */ - unsigned char *saved = NULL; /*original values */ - char str[256]; /*hello string */ - float hw_f; /*hardware-converted */ - double hw_d; /*hardware-converted */ - void *aligned=NULL; /*aligned buffer */ -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - long double hw_ld; /*hardware-converted */ -#endif - unsigned char *hw=NULL; /*ptr to hardware-conv'd*/ - int underflow; /*underflow occurred */ - int uflow=0; /*underflow debug counters*/ - size_t i, j, k; /*counters */ - int endian; /*machine endianess */ - size_t dst_ebias; /* Destination type's exponent bias */ - size_t src_epos; /* Source type's exponent position */ - size_t src_esize; /* Source type's exponent size */ - size_t dst_epos; /* Destination type's exponent position */ - size_t dst_esize; /* Destination type's exponent size */ - size_t dst_msize; /* Destination type's mantissa size */ - size_t src_nbits; /* source length in bits */ - size_t dst_nbits; /* dst length in bits */ - -#ifdef HANDLE_SIGFPE - pid_t child_pid; /*process ID of child */ - int status; /*child exit status */ - - /* - * Some systems generage SIGFPE during floating point overflow and we - * cannot assume that we can continue from such a signal. Therefore, we - * fork here and let the child run the test and return the number of - * failures with the exit status. - */ - HDfflush(stdout); - HDfflush(stderr); - if ((child_pid=fork())<0) { - HDperror("fork"); - return 1; - } else if (child_pid>0) { - while (child_pid!=waitpid(child_pid, &status, 0)) /*void*/; - if (WIFEXITED(status) && 255==WEXITSTATUS(status)) { - return 0; /*child exit after catching SIGFPE*/ - } else if (WIFEXITED(status)) { - return WEXITSTATUS(status); - } else { - HDputs(" Child didn't exit normally."); - return 1; - } - } -#endif + char *buf=NULL, s[80]; + hid_t c_type, f_type; + const size_t nelmts = NTESTELEM, ntests=NTESTS; + size_t i, j, nchars; + int ret_value = 1; /* - * The remainder of this function is executed only by the child if - * HANDLE_SIGFPE is defined. - */ - HDsignal(SIGFPE,fpe_handler); - - /* What are the names of the source and destination types */ - if (H5Tequal(src, H5T_NATIVE_FLOAT)) { - src_type_name = "float"; - src_type = FLT_FLOAT; - } else if (H5Tequal(src, H5T_NATIVE_DOUBLE)) { - src_type_name = "double"; - src_type = FLT_DOUBLE; -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (H5Tequal(src, H5T_NATIVE_LDOUBLE)) { - src_type_name = "long double"; - src_type = FLT_LDOUBLE; -#endif - } else { - src_type_name = "UNKNOWN"; - src_type = OTHER; - } - - if (H5Tequal(dst, H5T_NATIVE_FLOAT)) { - dst_type_name = "float"; - dst_type = FLT_FLOAT; - } else if (H5Tequal(dst, H5T_NATIVE_DOUBLE)) { - dst_type_name = "double"; - dst_type = FLT_DOUBLE; -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (H5Tequal(dst, H5T_NATIVE_LDOUBLE)) { - dst_type_name = "long double"; - dst_type = FLT_LDOUBLE; -#endif - } else { - dst_type_name = "UNKNOWN"; - dst_type = OTHER; - } - - /* Sanity checks */ - if(sizeof(float)==sizeof(double)) - HDputs("Sizeof(float)==sizeof(double) - some tests may not be sensible."); - if (OTHER==src_type || OTHER==dst_type) { - if(!strcmp(name, "noop")) - sprintf(str, "Testing %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else if(run_special) - sprintf(str, "Testing special %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else - sprintf(str, "Testing regular %s %s -> %s conversions", - name, src_type_name, dst_type_name); - - printf("%-70s", str); - H5_FAILED(); - HDputs(" Unknown data type."); - goto error; - } - - /* Get "interesting" values */ - src_size = H5Tget_size(src); - dst_size = H5Tget_size(dst); - src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */ - dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */ - dst_ebias=H5Tget_ebias(dst); - H5Tget_fields(src,NULL,&src_epos,&src_esize,NULL,NULL); - H5Tget_fields(dst,NULL,&dst_epos,&dst_esize,NULL,&dst_msize); - endian = H5Tget_order(H5T_NATIVE_FLOAT); - - /* Allocate buffers */ - aligned = HDmalloc(MAX(sizeof(long double), sizeof(double))); -#ifdef SHOW_OVERFLOWS - noverflows_g = 0; -#endif - - /* Allocate and initialize the source buffer through macro INIT_FP or INIT_FP_SPECIAL. - * The BUF will be used for the conversion while the SAVED buffer will be used for - * the comparison later. INIT_FP will fill in the buffer with regular values like - * normalized and denormalized values; INIT_FP_SPECIAL will fill with special values - * like infinity, NaN. + * Initialize types and buffer. */ - switch (run_special) { - case FALSE: - if(src_type == FLT_FLOAT) { - INIT_FP(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, FLT_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); - } else if(src_type == FLT_DOUBLE) { - INIT_FP(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, DBL_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if(src_type == FLT_LDOUBLE) { - INIT_FP(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, LDBL_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); -#endif - } else - goto error; - - break; - case TRUE: - if(src_type == FLT_FLOAT) { - INIT_FP_SPECIAL(src_size, src_nbits, endian, FLT_MANT_DIG, dst_size, - buf, saved, nelmts); - } else if(src_type == FLT_DOUBLE) { - INIT_FP_SPECIAL(src_size, src_nbits, endian, DBL_MANT_DIG, dst_size, - buf, saved, nelmts); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if(src_type == FLT_LDOUBLE) { - INIT_FP_SPECIAL(src_size, src_nbits, endian, LDBL_MANT_DIG, dst_size, - buf, saved, nelmts); -#endif - } else - goto error; - - break; - default: - goto error; + c_type = mkstr(8, H5T_STR_NULLPAD); + f_type = mkstr(8, H5T_STR_SPACEPAD); + buf = (char*)HDcalloc(nelmts, 8); + for (i=0; i1) { - if(!strcmp(name, "noop")) - sprintf(str, "Testing %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); - else if(run_special) - sprintf(str, "Testing special %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); - else - sprintf(str, "Testing regular %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); + sprintf(s, "Testing random string conversion speed (test %d/%d)", + (int)(i+1), (int)ntests); } else { - if(!strcmp(name, "noop")) - sprintf(str, "Testing %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); - else if(run_special) - sprintf(str, "Testing special %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else - sprintf(str, "Testing regular %s %s -> %s conversions", - name, src_type_name, dst_type_name); + sprintf(s, "Testing random string conversion speed"); } - printf("%-70s", str); + printf("%-70s", s); HDfflush(stdout); - fails_this_test = 0; - - /* Perform the conversion in software */ - if (H5Tconvert(src, dst, nelmts, buf, NULL, H5P_DEFAULT)<0) + if (H5Tconvert(c_type, f_type, nelmts, buf, NULL, H5P_DEFAULT)<0) + goto error; + if (H5Tconvert(f_type, c_type, nelmts, buf, NULL, H5P_DEFAULT)<0) goto error; - - /* Check the software results against the hardware */ - for (j=0; j=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto done; - } - } PASSED(); } -#ifdef SHOW_OVERFLOWS - if (noverflows_g>0) - printf(" %d overflow%s in previous test\n", - noverflows_g, 1==noverflows_g?"":"s"); -#endif - - done: -#ifdef AKCDEBUG - printf("uflow=%d, fails_all_tests=%d\n", uflow, fails_all_tests); -#endif - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); -#ifdef HANDLE_SIGFPE - HDexit(MIN((int)fails_all_tests, 254)); -#else - reset_hdf5(); - return (int)fails_all_tests; -#endif + ret_value = 0; error: - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); -#ifdef HANDLE_SIGFPE - HDexit(MIN(MAX((int)fails_all_tests, 1), 254)); -#else + if (buf) HDfree(buf); reset_hdf5(); - return MAX((int)fails_all_tests, 1); -#endif + return ret_value; } /*------------------------------------------------------------------------- - * Function: test_conv_int_float - * - * Purpose: Test conversion between random integer and float values - * from SRC to DST. These types should be any combination of: + * Function: test_conv_enum_1 * - * H5T_NATIVE_SCHAR H5T_NATIVE_FLOAT - * H5T_NATIVE_SHORT H5T_NATIVE_DOUBLE - * H5T_NATIVE_INT H5T_NATIVE_LDOUBLE - * H5T_NATIVE_LONG - * H5T_NATIVE_LLONG + * Purpose: Test conversion speed for enum data types * * Return: Success: 0 * * Failure: number of errors * - * Programmer: Raymond Lu - * Thursday, November 6, 2003 + * Programmer: Robb Matzke + * Tuesday, January 5, 1999 * * Modifications: * *------------------------------------------------------------------------- */ static int -test_conv_int_float(const char *name, hbool_t run_special, hid_t src, hid_t dst) +test_conv_enum_1(void) { - hid_t dxpl_id; /*dataset transfer property list*/ - int fill_value=9; /*fill value for conversion exception*/ - H5T_conv_except_func_t op; /*returned callback function for conversion exception*/ - void *user_data; /*returned pointer to user data passed in to the callback*/ - hbool_t except_set = FALSE; /*whether user's exception handling is set*/ - const size_t ntests=NTESTS; /*number of tests */ - size_t nelmts=0; /*num values per test */ - const size_t max_fails=40; /*max number of failures*/ - size_t fails_all_tests=0; /*number of failures */ - size_t fails_this_test; /*fails for this test */ - char str[256]; /*hello string */ - dtype_t src_type; /*data types */ - dtype_t dst_type; /*data types */ - const char *src_type_name=NULL; /*source type name */ - const char *dst_type_name=NULL; /*destination type name */ - int endian; /*machine endianess */ - size_t src_size, dst_size; /*type sizes */ - unsigned char *buf=NULL; /*buffer for conversion */ - unsigned char *saved=NULL; /*original values */ - size_t i, j, k; /*counters */ - unsigned char *hw=NULL; /*hardware conv result */ - unsigned char src_bits[32]; /*src value in LE order */ - unsigned char dst_bits[32]; /*dest value in LE order*/ - size_t src_nbits; /*source length in bits */ - size_t dst_nbits; /*dst length in bits */ - void *aligned=NULL; /*aligned temp buffer */ - float hw_float=0; - double hw_double=0; - long double hw_ldouble=0; - signed char hw_schar=0; - unsigned char hw_uchar=0; - short hw_short=0; - unsigned short hw_ushort=0; - int hw_int=0; - unsigned hw_uint=0; - long hw_long=0; - unsigned long hw_ulong=0; - long_long hw_llong=0; - unsigned long_long hw_ullong=0; - - /* What is the name of the source type */ - if (H5Tequal(src, H5T_NATIVE_SCHAR)) { - src_type_name = "signed char"; - src_type = INT_SCHAR; - } else if (H5Tequal(src, H5T_NATIVE_UCHAR)) { - src_type_name = "unsigned char"; - src_type = INT_UCHAR; - } else if (H5Tequal(src, H5T_NATIVE_SHORT)) { - src_type_name = "short"; - src_type = INT_SHORT; - } else if (H5Tequal(src, H5T_NATIVE_USHORT)) { - src_type_name = "unsigned short"; - src_type = INT_USHORT; - } else if (H5Tequal(src, H5T_NATIVE_INT)) { - src_type_name = "int"; - src_type = INT_INT; - } else if (H5Tequal(src, H5T_NATIVE_UINT)) { - src_type_name = "unsigned int"; - src_type = INT_UINT; - } else if (H5Tequal(src, H5T_NATIVE_LONG)) { - src_type_name = "long"; - src_type = INT_LONG; - } else if (H5Tequal(src, H5T_NATIVE_ULONG)) { - src_type_name = "unsigned long"; - src_type = INT_ULONG; - } else if (H5Tequal(src, H5T_NATIVE_LLONG)) { - src_type_name = "long long"; - src_type = INT_LLONG; - } else if (H5Tequal(src, H5T_NATIVE_ULLONG)) { - src_type_name = "unsigned long long"; - src_type = INT_ULLONG; - } else if (H5Tequal(src, H5T_NATIVE_FLOAT)) { - src_type_name = "float"; - src_type = FLT_FLOAT; - } else if (H5Tequal(src, H5T_NATIVE_DOUBLE)) { - src_type_name = "double"; - src_type = FLT_DOUBLE; -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (H5Tequal(src, H5T_NATIVE_LDOUBLE)) { - src_type_name = "long double"; - src_type = FLT_LDOUBLE; -#endif - } else { - src_type_name = "UNKNOWN"; - src_type = OTHER; - } - - /* What is the name of the destination type */ - if (H5Tequal(dst, H5T_NATIVE_SCHAR)) { - dst_type_name = "signed char"; - dst_type = INT_SCHAR; - } else if (H5Tequal(dst, H5T_NATIVE_UCHAR)) { - dst_type_name = "unsigned char"; - dst_type = INT_UCHAR; - } else if (H5Tequal(dst, H5T_NATIVE_SHORT)) { - dst_type_name = "short"; - dst_type = INT_SHORT; - } else if (H5Tequal(dst, H5T_NATIVE_USHORT)) { - dst_type_name = "unsigned short"; - dst_type = INT_USHORT; - } else if (H5Tequal(dst, H5T_NATIVE_INT)) { - dst_type_name = "int"; - dst_type = INT_INT; - } else if (H5Tequal(dst, H5T_NATIVE_UINT)) { - dst_type_name = "unsigned int"; - dst_type = INT_UINT; - } else if (H5Tequal(dst, H5T_NATIVE_LONG)) { - dst_type_name = "long"; - dst_type = INT_LONG; - } else if (H5Tequal(dst, H5T_NATIVE_ULONG)) { - dst_type_name = "unsigned long"; - dst_type = INT_ULONG; - } else if (H5Tequal(dst, H5T_NATIVE_LLONG)) { - dst_type_name = "long long"; - dst_type = INT_LLONG; - } else if (H5Tequal(dst, H5T_NATIVE_ULLONG)) { - dst_type_name = "unsigned long long"; - dst_type = INT_ULLONG; - } else if (H5Tequal(dst, H5T_NATIVE_FLOAT)) { - dst_type_name = "float"; - dst_type = FLT_FLOAT; - } else if (H5Tequal(dst, H5T_NATIVE_DOUBLE)) { - dst_type_name = "double"; - dst_type = FLT_DOUBLE; -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if (H5Tequal(dst, H5T_NATIVE_LDOUBLE)) { - dst_type_name = "long double"; - dst_type = FLT_LDOUBLE; -#endif - } else { - dst_type_name = "UNKNOWN"; - dst_type = OTHER; - } - - /* Sanity checks */ - if (OTHER==src_type || OTHER==dst_type) { - sprintf(str, "Testing random %s %s -> %s conversions", - name, src_type_name, dst_type_name); - printf("%-70s", str); - H5_FAILED(); - HDputs(" Unknown data type."); - goto error; - } - - if ((INT_SCHAR==src_type || INT_UCHAR==src_type || INT_SHORT==src_type || - INT_USHORT==src_type || INT_INT==src_type || INT_UINT==src_type || - INT_LONG==src_type || INT_ULONG==src_type || INT_LLONG==src_type || - INT_ULLONG==src_type) && - (FLT_FLOAT!=dst_type && FLT_DOUBLE!=dst_type && FLT_LDOUBLE!=dst_type)) { - if(!run_special) - sprintf(str, "Testing regular %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else - sprintf(str, "Testing special %s %s -> %s conversions", - name, src_type_name, dst_type_name); - printf("%-70s", str); - H5_FAILED(); - HDputs(" 1. Not an integer-float conversion."); - goto error; - } + const size_t nelmts=NTESTELEM; + const int ntests=NTESTS; + int i, val, *buf=NULL; + hid_t t1, t2; + char s[80]; + int ret_value = 1; - if ((FLT_FLOAT==src_type || FLT_DOUBLE==src_type || FLT_LDOUBLE==src_type) - && (INT_SCHAR!=dst_type && INT_UCHAR!=dst_type && INT_SHORT!=dst_type - && INT_USHORT!=dst_type && INT_INT!=dst_type && INT_UINT!=dst_type - && INT_LONG!=dst_type && INT_ULONG!=dst_type && INT_LLONG!=dst_type - && INT_ULLONG!=dst_type)) { - if(!run_special) - sprintf(str, "Testing regular %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else - sprintf(str, "Testing special %s %s -> %s conversions", - name, src_type_name, dst_type_name); - printf("%-70s", str); - H5_FAILED(); - HDputs(" 2. Not a float-integer conversion."); - goto error; + /* Build the data types */ + t1 = H5Tcreate(H5T_ENUM, sizeof(int)); + t2 = H5Tenum_create(H5T_NATIVE_INT); + s[1] = '\0'; + for (i=0; i<26; i++) { + s[0] = 'A'+i; + H5Tenum_insert(t1, s, &i); + H5Tenum_insert(t2, s, (val=i*1000+i, &val)); } - - /* Some information about datatypes */ - endian = H5Tget_order(H5T_NATIVE_INT); - src_size = H5Tget_size(src); - dst_size = H5Tget_size(dst); - src_nbits = H5Tget_precision(src); /* not 8*src_size, esp on J90 - QAK */ - dst_nbits = H5Tget_precision(dst); /* not 8*dst_size, esp on J90 - QAK */ - aligned = HDmalloc(MAX(sizeof(long double), sizeof(long_long))); -#ifdef SHOW_OVERFLOWS - noverflows_g = 0; -#endif - - /* This is for some Linux systems where long double has the size - * 12 bytes but precision is 10 bytes. The 2 unused bytes may - * have garbage causing wrong value comparison. - */ - HDmemset(&hw_ldouble, 0, sizeof(long double)); - - /* Create a dataset transfer property list and datatype conversion - * exception handler function and pass in fill value. This is mainly - * for NetCDF compatibility, which requests fill in fill value when - * conversion exception happens. We only test (unsigned) int - float - * and float - (unsigned) int conversions, which should cover more cases. - */ - if((dxpl_id = H5Pcreate(H5P_DATASET_XFER))<0) - goto error; - - if((src_type == INT_INT && dst_type == FLT_FLOAT) || - (src_type == INT_UINT && dst_type == FLT_FLOAT) || - (src_type == FLT_FLOAT && dst_type == INT_UINT) || - (src_type == FLT_FLOAT && dst_type == INT_INT)) { - if(H5Pset_type_conv_cb(dxpl_id, except_func, &fill_value)<0) - goto error; - else - except_set = TRUE; - - if(H5Pget_type_conv_cb(dxpl_id, &op, &user_data)<0) - goto error; - if(op != except_func || *(int*)user_data != fill_value) - goto error; - } + /* Initialize the buffer */ + buf = (int*)HDmalloc(nelmts*MAX(H5Tget_size(t1), H5Tget_size(t2))); + for (i=0; i<(int)nelmts; i++) + buf[i] = HDrand() % 26; - /* Allocate and initialize the source buffer through macro INIT_INTEGER if the source is integer, - * INIT_FP if floating-point. The BUF will be used for the conversion while the SAVED buffer will be - * used for the comparison later. - */ - if(src_type == INT_SCHAR) { - INIT_INTEGER(signed char, SCHAR_MAX, SCHAR_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_UCHAR) { - INIT_INTEGER(unsigned char, UCHAR_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_SHORT) { - INIT_INTEGER(short, SHRT_MAX, SHRT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_USHORT) { - INIT_INTEGER(unsigned short, USHRT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_INT) { - INIT_INTEGER(int, INT_MAX, INT_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_UINT) { - INIT_INTEGER(unsigned int, UINT_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_LONG) { - INIT_INTEGER(long, LONG_MAX, LONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_ULONG) { - INIT_INTEGER(unsigned long, ULONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_LLONG) { - INIT_INTEGER(long_long, LLONG_MAX, LLONG_MIN, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == INT_ULLONG) { - INIT_INTEGER(unsigned long_long, ULLONG_MAX, 0, src_size, dst_size, src_nbits, buf, saved, nelmts); - } else if(src_type == FLT_FLOAT) { - if(!run_special) { - INIT_FP(float, FLT_MAX, FLT_MIN, FLT_MAX_10_EXP, FLT_MIN_10_EXP, FLT_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); - } else { - INIT_FP_SPECIAL(src_size, src_nbits, endian, FLT_MANT_DIG, dst_size, buf, saved, nelmts); - } - } else if(src_type == FLT_DOUBLE) { - if(!run_special) { - INIT_FP(double, DBL_MAX, DBL_MIN, DBL_MAX_10_EXP, DBL_MIN_10_EXP, DBL_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); - } else { - INIT_FP_SPECIAL(src_size, src_nbits, endian, DBL_MANT_DIG, dst_size, buf, saved, nelmts); - } -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - } else if(src_type == FLT_LDOUBLE) { - if(!run_special) { - INIT_FP(long double, LDBL_MAX, LDBL_MIN, LDBL_MAX_10_EXP, LDBL_MIN_10_EXP, LDBL_MANT_DIG, - src_size, src_nbits, endian, dst_size, buf, saved, nelmts); - } else { - INIT_FP_SPECIAL(src_size, src_nbits, endian, LDBL_MANT_DIG, dst_size, buf, saved, nelmts); - } -#endif - } else - goto error; - - /* The tests */ + /* Conversions */ for (i=0; i1) { - if(!run_special) - sprintf(str, "Testing regular %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); - else - sprintf(str, "Testing special %s %s -> %s conversions (test %d/%d)", - name, src_type_name, dst_type_name, (int)i+1, (int)ntests); + sprintf(s, "Testing random enum conversion O(N) (test %d/%d)", + i+1, ntests); } else { - if(!run_special) - sprintf(str, "Testing regular %s %s -> %s conversions", - name, src_type_name, dst_type_name); - else - sprintf(str, "Testing special %s %s -> %s conversions", - name, src_type_name, dst_type_name); + sprintf(s, "Testing random enum conversion O(N)"); } - printf("%-70s", str); + printf("%-70s", s); HDfflush(stdout); - fails_this_test=0; - - /* Perform the conversion */ - if (H5Tconvert(src, dst, nelmts, buf, NULL, dxpl_id)<0) - goto error; - - /* Check the results from the library against hardware */ - for (j=0; j=max_fails) { - HDputs(" maximum failures reached, aborting test..."); - goto done; - } + for (i=0; i1) { + sprintf(s, "Testing random enum conversion O(N log N) " + "(test %d/%d)", i+1, ntests); + } else { + sprintf(s, "Testing random enum conversion O(N log N)"); } + printf("%-70s", s); + HDfflush(stdout); + if (H5Tconvert(t2, t1, nelmts, buf, NULL, H5P_DEFAULT)<0) goto error; PASSED(); } -#ifdef SHOW_OVERFLOWS - if (noverflows_g>0) { - printf(" %d overflow%s in previous test\n", - noverflows_g, 1==noverflows_g?"":"s"); - } -#endif - - done: - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); - reset_hdf5(); /*print statistics*/ - return (int)fails_all_tests; + ret_value = 0; error: - if (buf) aligned_free(buf); - if (saved) aligned_free(saved); - if (aligned) HDfree(aligned); - HDfflush(stdout); - reset_hdf5(); /*print statistics*/ - return MAX((int)fails_all_tests, 1); + H5Tclose(t1); + H5Tclose(t2); + if (buf) HDfree(buf); + reset_hdf5(); + return ret_value; } /*------------------------------------------------------------------------- - * Function: overflows - * - * Purpose: When convert from float or double to any integer type, - * check if overflow occurs. - * - * - * Return: TRUE: overflow happens + * Function: test_conv_enum_2 * - * FALSE: no overflow + * Purpose: Tests enumeration conversions where source isn't a native type. * - * Programmer: Raymond Lu - * Monday, Nov 17, 2003 + * Return: Success: 0 + * + * Failure: number of errors * + * Programmer: Robb Matzke, LLNL, 2003-06-09 + * * Modifications: - * *------------------------------------------------------------------------- */ -static hbool_t -overflows(unsigned char *origin_bits, hid_t src_id, size_t dst_num_bits) +static int +test_conv_enum_2(void) { - hbool_t ret_value=FALSE; - hsize_t expt; - size_t mant_digits=0, expt_digits=0, bias=0; - size_t epos, mpos; - size_t src_prec=0; /*source type precision in bits*/ - H5T_norm_t norm; - ssize_t indx; - unsigned char bits[32], mant_bits[32]; - - HDmemset(bits, 0, 32); - HDmemset(mant_bits, 0, 32); - - /* - * Sometimes, type size isn't equal to the precision like Linux's "long - * double", where size is 96 bits and precision is 80 bits. - */ - - src_prec = H5Tget_precision(src_id); - H5Tget_fields(src_id, NULL, &epos, &expt_digits, &mpos, &mant_digits); - bias = H5Tget_ebias(src_id); - norm = H5Tget_norm(src_id); - - HDmemcpy(bits, origin_bits, src_prec/8+1); - - /*Check for special cases: +Inf, -Inf*/ - if (H5T_bit_find (bits, mpos, mant_digits, H5T_BIT_LSB, TRUE)<0) { - if (H5T_bit_find (bits, epos, expt_digits, H5T_BIT_LSB, FALSE)<0) { - ret_value=TRUE; - goto done; - } - } else if (H5T_NORM_NONE==norm && H5T_bit_find (bits, mpos, mant_digits-1, - H5T_BIT_LSB, TRUE)<0 && H5T_bit_find (bits, epos, expt_digits, - H5T_BIT_LSB, FALSE)<0) { - /*This is a special case for the source of no implied mantissa bit. - *If the exponent bits are all 1s and only the 1st bit of mantissa - *is set to 1. It's infinity. The Intel-Linux "long double" is this case.*/ - ret_value=TRUE; - goto done; - } + hid_t srctype=-1, dsttype=-1, oddsize=-1; + int *data=NULL, i, nerrors=0; + const char *mname[] = { "RED", + "GREEN", + "BLUE", + "YELLOW", + "PINK", + "PURPLE", + "ORANGE", + "WHITE" }; - /* get exponent */ - expt = H5T_bit_get_d(bits, mant_digits, expt_digits) - bias; - - if(expt>=(dst_num_bits-1)) { - ret_value=TRUE; - goto done; - } + TESTING("non-native enumeration type conversion"); + + /* Source enum type */ + oddsize = H5Tcopy(H5T_STD_I32BE); + H5Tset_size(oddsize, 3); /*reduce to 24 bits, not corresponding to any native size*/ + srctype = H5Tenum_create(oddsize); + for (i=7; i>=0; --i) { + char pattern[3]; + pattern[2] = i; + pattern[0] = pattern[1] = 0; + H5Tenum_insert(srctype, mname[i], pattern); + } - /* get significand */ - H5T_bit_copy (mant_bits, 0, bits, 0, mant_digits); + /* Destination enum type */ + dsttype = H5Tenum_create(H5T_NATIVE_INT); + assert(H5Tget_size(dsttype)>H5Tget_size(srctype)); + for (i=0; i<8; i++) + H5Tenum_insert(dsttype, mname[i], &i); - - /* restore implicit bit if normalization is implied*/ - if(norm == H5T_NORM_IMPLIED) { - H5T_bit_inc(mant_bits, mant_digits, 1); - mant_digits++; + /* Source data */ + data = (int*)malloc(NTESTELEM*sizeof(int)); + for (i=0; i=dst_num_bits) - ret_value=TRUE; + /* Cleanup */ + free(data); + H5Tclose(srctype); + H5Tclose(dsttype); + H5Tclose(oddsize); -done: - return ret_value; + /* Failure */ + if (nerrors) { + printf("total of %d conversion errors out of %d elements for enums\n", + nerrors, NTESTELEM); + return 1; + } + + PASSED(); + return 0; } /*------------------------------------------------------------------------- - * Function: run_integer_tests + * Function: test_conv_bitfield + * + * Purpose: Test bitfield conversions. * - * Purpose: Runs all integer tests. + * Return: Success: 0 * - * Return: Number of errors + * Failure: number of errors * * Programmer: Robb Matzke - * Tuesday, November 24, 1998 + * Thursday, May 20, 1999 * * Modifications: * *------------------------------------------------------------------------- */ static int -run_integer_tests(const char *name) +test_conv_bitfield(void) { - 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 + unsigned char buf[4]; + hid_t st=-1, dt=-1; - 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 + TESTING("bitfield conversions"); -#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 + /* + * First test a simple bitfield conversion: + * 1010101010101010 + * ________________1010101010101010 + */ + st = H5Tcopy(H5T_STD_B16LE); + dt = H5Tcopy(H5T_STD_B32LE); + buf[0] = buf[1] = 0xAA; + buf[2] = buf[3] = 0x55; /*irrelevant*/ + if (H5Tconvert(st, dt, 1, buf, NULL, H5P_DEFAULT)<0) goto error; + if (buf[0]!=0xAA || buf[1]!=0xAA || buf[2]!=0 || buf[3]!=0) { + H5_FAILED(); + printf(" s=0xaaaa, d=0x%02x%02x%02x%02x (test 1)\n", + buf[3], buf[2], buf[1], buf[0]); + goto error; + } -#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 + /* + * Test2: Offset a 12-byte value in the middle of a 16 and 32 byte + * field. + * __10 1010 1010 10__ + * ____ ____ __10 1010 1010 10__ ____ ____ + */ + H5Tset_precision(st, 12); + H5Tset_offset(st, 2); + H5Tset_precision(dt, 12); + H5Tset_offset(dt, 10); + buf[0] = 0xA8; buf[1] = 0x2A; buf[2] = buf[3] = 0; + if (H5Tconvert(st, dt, 1, buf, NULL, H5P_DEFAULT)<0) goto error; + if (buf[0]!=0 || buf[1]!=0xA8 || buf[2]!=0x2A || buf[3]!=0) { + H5_FAILED(); + printf(" s=0x2AA8 d=0x%02x%02x%02x%02x (test 2)\n", + buf[3], buf[2], buf[1], buf[0]); + goto error; + } -#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 + /* + * Same as previous test except unused bits of the destination will + * be filled with ones. + */ + H5Tset_pad(dt, H5T_PAD_ONE, H5T_PAD_ONE); + buf[0] = 0xA8; buf[1] = 0x2A; buf[2] = buf[3] = 0; + if (H5Tconvert(st, dt, 1, buf, NULL, H5P_DEFAULT)<0) goto error; + if (buf[0]!=0xff || buf[1]!=0xAB || buf[2]!=0xEA || buf[3]!=0xff) { + H5_FAILED(); + printf(" s=0x2AA8 d=0x%02x%02x%02x%02x (test 3)\n", + buf[3], buf[2], buf[1], buf[0]); + goto error; + } -#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 + H5Tclose(st); + H5Tclose(dt); + PASSED(); + reset_hdf5(); + return 0; - return nerrors; + error: + H5Tclose(st); + H5Tclose(dt); + reset_hdf5(); + return 1; } /*------------------------------------------------------------------------- - * Function: run_fp_tests + * Function: convert_opaque * - * Purpose: Runs all floating-point tests. + * Purpose: A fake opaque conversion functions * - * Return: Number of errors + * Return: Success: 0 * - * Programmer: Raymond Lu - * Tuesday, March 22, 2005 + * Failure: -1 + * + * Programmer: Robb Matzke + * Friday, June 4, 1999 * * Modifications: * *------------------------------------------------------------------------- */ -static int -run_fp_tests(const char *name) +static herr_t +convert_opaque(hid_t UNUSED st, hid_t UNUSED dt, H5T_cdata_t *cdata, + size_t UNUSED nelmts, size_t UNUSED buf_stride, + size_t UNUSED bkg_stride, void UNUSED *_buf, + void UNUSED *bkg, hid_t UNUSED dset_xfer_plid) { - int nerrors = 0; - - if(!strcmp(name, "noop")) { - nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_FLOAT); - nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE); - nerrors += test_conv_flt_1("noop", FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LDOUBLE); - goto done; - } - - /*Test normalized and denormalized values. FALSE indicates non-special values.*/ - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE); - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT); - nerrors += test_conv_flt_1(name, FALSE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE); -#endif - - /*Test special values, +/-0, +/-infinity, +/-QNaN, +/-SNaN.*/ - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_FLOAT, H5T_NATIVE_DOUBLE); - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_DOUBLE, H5T_NATIVE_FLOAT); -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_FLOAT, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_DOUBLE, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_FLOAT); - nerrors += test_conv_flt_1(name, TRUE, H5T_NATIVE_LDOUBLE, H5T_NATIVE_DOUBLE); -#endif - -done: - return nerrors; + if (H5T_CONV_CONV==cdata->command) num_opaque_conversions_g++; + return 0; } /*------------------------------------------------------------------------- - * Function: run_int_fp_conv + * Function: test_opaque * - * Purpose: Runs all integer-float tests. + * Purpose: Driver function to test opaque datatypes + * + * Return: Success: 0 * - * Return: Number of errors + * Failure: number of errors * * Programmer: Raymond Lu - * Monday, November 10, 2003 + * June 2, 2004 * * Modifications: * *------------------------------------------------------------------------- */ static int -run_int_fp_conv(const char *name) +test_opaque(void) { - int nerrors = 0; - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_DOUBLE); + int num_errors = 0; -#if H5_SIZEOF_LONG!=H5_SIZEOF_INT - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_DOUBLE); - - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_DOUBLE); -#endif - -#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_DOUBLE); - -#ifdef H5_ULLONG_TO_FP_CAST_WORKS - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_FLOAT); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_DOUBLE); -#else /* H5_ULLONG_TO_FP_CAST_WORKS */ - { - char str[256]; /*hello string */ - - sprintf(str, "Testing random %s %s -> %s conversions", - name, "unsigned long long", "float"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to compiler not handling conversion."); - - sprintf(str, "Testing random %s %s -> %s conversions", - name, "unsigned long long", "double"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to compiler not handling conversion."); - } -#endif /* H5_ULLONG_TO_FP_CAST_WORKS */ -#endif + TESTING("opaque datatypes"); -#if H5_SW_INTEGER_TO_LDOUBLE_WORKS -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SCHAR, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UCHAR, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_SHORT, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_USHORT, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_INT, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_UINT, H5T_NATIVE_LDOUBLE); -#if H5_SIZEOF_LONG!=H5_SIZEOF_INT - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LONG, H5T_NATIVE_LDOUBLE); - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULONG, H5T_NATIVE_LDOUBLE); -#endif -#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_LLONG, H5T_NATIVE_LDOUBLE); -#if H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS - nerrors += test_conv_int_float(name, FALSE, H5T_NATIVE_ULLONG, H5T_NATIVE_LDOUBLE); -#else /* H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS */ - { - char str[256]; /*hello string */ + /* Test opaque types with tags */ + num_errors += opaque_check(0); + /* Test opaque types without tag */ + num_errors += opaque_check(1); - sprintf(str, "Testing random %s %s -> %s conversions", - name, "unsigned long long", "long double"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to compiler not handling conversion."); - } -#endif /* H5_ULLONG_TO_FP_CAST_WORKS && H5_ULLONG_TO_LDOUBLE_PRECISION_WORKS */ -#endif -#endif -#else /*H5_SW_INTEGER_TO_LDOUBLE_WORKS*/ - { - char str[256]; /*string */ + if(num_errors) + goto error; - sprintf(str, "Testing random %s %s -> %s conversions", - name, "all integers", "long double"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - } -#endif /*H5_SW_INTEGER_TO_LDOUBLE_WORKS*/ + PASSED(); + return 0; - return nerrors; + error: + H5_FAILED(); + return num_errors; } /*------------------------------------------------------------------------- - * Function: run_fp_int_conv + * Function: opaque_check * - * Purpose: Runs all float-integer tests. + * Purpose: Test opaque datatypes * - * Return: Number of errors + * Return: Success: 0 * - * Programmer: Raymond Lu - * Monday, November 10, 2003 + * Failure: number of errors + * + * Programmer: Robb Matzke + * Thursday, May 20, 1999 * * Modifications: + * Raymond Lu + * June 2, 2004 + * Made tag for one opaque type optional. * *------------------------------------------------------------------------- */ static int -run_fp_int_conv(const char *name) +opaque_check(int tag_it) { - int nerrors = 0; - hbool_t special; - int i; - - for(i=0; i<2; i++) { - if(i==0) - special = FALSE; - else - special = TRUE; +#define OPAQUE_NELMTS 1000 + hid_t st=-1, dt=-1; + herr_t status; + char buf[1]; /*not really used*/ + int saved; - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_SCHAR); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_SCHAR); - - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_UCHAR); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_UCHAR); + saved = num_opaque_conversions_g = 0; - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_SHORT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_SHORT); - - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_USHORT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_USHORT); + /* Build source and destination types */ + if ((st=H5Tcreate(H5T_OPAQUE, 4))<0) goto error; + if (H5Tset_tag(st, "opaque source type")<0) goto error; + + if ((dt=H5Tcreate(H5T_OPAQUE, 4))<0) goto error; + if (tag_it) { + if (H5Tset_tag(dt, "opaque destination type")<0) + goto error; + } - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_INT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_INT); - - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_UINT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_UINT); - -#if H5_SIZEOF_LONG!=H5_SIZEOF_INT - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LONG); - - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_ULONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULONG); -#endif - -#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG - if(!strcmp(name, "hw")) { /* Hardware conversion */ - /* Windows .NET 2003 doesn't work for hardware conversion of this case. - * .NET should define this macro H5_HW_FP_TO_LLONG_NOT_WORKS. */ -#ifndef H5_HW_FP_TO_LLONG_NOT_WORKS - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG); -#endif /*H5_HW_FP_TO_LLONG_NOT_WORKS*/ - } else { /* Software conversion */ - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_LLONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_LLONG); - } -#ifdef H5_FP_TO_ULLONG_RIGHT_MAXIMUM - nerrors += test_conv_int_float(name, special, H5T_NATIVE_FLOAT, H5T_NATIVE_ULLONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_DOUBLE, H5T_NATIVE_ULLONG); -#else /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/ - { - char str[256]; /*hello string */ - - sprintf(str, "Testing random %s %s -> %s conversions", - name, "float", "unsigned long long"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - - sprintf(str, "Testing random %s %s -> %s conversions", - name, "double", "unsigned long long"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - } -#endif /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/ -#endif + /* Make sure that we can't convert between the types yet */ + H5E_BEGIN_TRY { + status = H5Tconvert(st, dt, OPAQUE_NELMTS, buf, NULL, H5P_DEFAULT); + } H5E_END_TRY; + if (status>=0) { + H5_FAILED(); + printf(" opaque conversion should have failed but succeeded\n"); + goto error; + } -#if H5_SW_LDOUBLE_TO_INTEGER_WORKS -#if H5_SIZEOF_LONG_DOUBLE!=H5_SIZEOF_DOUBLE - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SCHAR); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UCHAR); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_SHORT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_USHORT); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_INT); -#if H5_CV_LDOUBLE_TO_UINT_WORKS - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_UINT); -#else /*H5_CV_LDOUBLE_TO_UINT_WORKS*/ - { - char str[256]; /*string */ - - sprintf(str, "Testing %s %s -> %s conversions", - name, "long double", "unsigned int"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - } -#endif /*H5_CV_LDOUBLE_TO_UINT_WORKS*/ -#if H5_SIZEOF_LONG!=H5_SIZEOF_INT - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LONG); - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULONG); -#endif + /* Register a conversion function */ + if (H5Tregister(H5T_PERS_HARD, "o_test", st, dt, convert_opaque)<0) + goto error; -#if H5_SIZEOF_LONG_LONG!=H5_SIZEOF_LONG - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_LLONG); -#ifdef H5_FP_TO_ULLONG_RIGHT_MAXIMUM - nerrors += test_conv_int_float(name, special, H5T_NATIVE_LDOUBLE, H5T_NATIVE_ULLONG); -#else /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/ - { - char str[256]; /*string */ - - sprintf(str, "Testing %s %s -> %s conversions", - name, "long double", "unsigned long long"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - } -#endif /*H5_FP_TO_ULLONG_RIGHT_MAXIMUM*/ -#endif -#endif -#else /*H5_SW_LDOUBLE_TO_INTEGER_WORKS*/ - { - char str[256]; /*hello string */ - - sprintf(str, "Testing %s %s -> %s conversions", - name, "long double", "all integers"); - printf("%-70s", str); - SKIPPED(); - HDputs(" Test skipped due to hardware conversion error."); - } -#endif /*H5_SW_LDOUBLE_TO_INTEGER_WORKS*/ + /* Try the conversion again, this time it should work */ + if (H5Tconvert(st, dt, OPAQUE_NELMTS, buf, NULL, H5P_DEFAULT)<0) goto error; + if (saved+1 != num_opaque_conversions_g) { + H5_FAILED(); + printf(" unexpected number of opaque conversions\n"); + goto error; } + + /* Unregister conversion function */ + if (H5Tunregister(H5T_PERS_HARD, "o_test", st, dt, convert_opaque)<0) + goto error; + + H5Tclose(st); + H5Tclose(dt); + return 0; - return nerrors; + error: + if (st>0) H5Tclose(st); + if (dt>0) H5Tclose(dt); + H5_FAILED(); + return 1; } @@ -8379,60 +3811,11 @@ main(void) nerrors += test_compound_10(); nerrors += test_compound_11(); nerrors += test_compound_12(); - nerrors += test_conv_int (); nerrors += test_conv_enum_1(); nerrors += test_conv_enum_2(); nerrors += test_conv_bitfield(); nerrors += test_opaque(); - /* Test user-define, query functions and software conversion - * for user-defined floating-point types */ - nerrors += test_derived_flt(); - - /* Test user-define, query functions and software conversion - * for user-defined integer types */ - nerrors += test_derived_integer(); - - /* Does floating point overflow generate a SIGFPE? */ - generates_sigfpe(); - - /* Test degenerate cases */ - nerrors += run_fp_tests("noop"); - - /* Test hardware floating-point conversion functions */ - nerrors += run_fp_tests("hw"); - - /* Test hardware integer conversion functions */ - nerrors += run_integer_tests("hw"); - - /* Test hardware integer-float conversion functions */ - nerrors += run_int_fp_conv("hw"); - - /* Test hardware float-integer conversion functions */ - nerrors += run_fp_int_conv("hw"); - - /*---------------------------------------------------------------------- - * Software tests - *---------------------------------------------------------------------- - */ - without_hardware_g = TRUE; - reset_hdf5(); - - /* Test software floating-point conversion functions */ - nerrors += run_fp_tests("sw"); - - /* Test software integer conversion functions */ - nerrors += test_conv_int_2(); - nerrors += run_integer_tests("sw"); - - /* Test software float-integer conversion functions */ - nerrors += run_fp_int_conv("sw"); - - /* Test software integer-float conversion functions */ - nerrors += run_int_fp_conv("sw"); - - reset_hdf5(); - if (nerrors) { printf("***** %lu FAILURE%s! *****\n", nerrors, 1==nerrors?"":"S"); -- cgit v0.12