diff options
author | Leon Arber <larber@ncsa.uiuc.edu> | 2004-05-06 22:44:38 (GMT) |
---|---|---|
committer | Leon Arber <larber@ncsa.uiuc.edu> | 2004-05-06 22:44:38 (GMT) |
commit | 1867c05efb2527e435776808720572e31ab95bb7 (patch) | |
tree | bd196861d102803a7c170872ff7ed0920942e95c /src | |
parent | 28eb73c97268c016728bce5e3e8d95ceb1f66b39 (diff) | |
download | hdf5-1867c05efb2527e435776808720572e31ab95bb7.zip hdf5-1867c05efb2527e435776808720572e31ab95bb7.tar.gz hdf5-1867c05efb2527e435776808720572e31ab95bb7.tar.bz2 |
[svn-r8488] Purpose:
Bug fix.
Description:
Fixed handling of trivial data transform expressions (like 5/3 + 3) and some
data coversion fixes. Also added more tests to dtransform.c
Solution:
Added some more checks in the H5Z_xform_reduce_tree function to see if perhaps
the transform expression is complicated and is a non-trivial reduction.
Added tests for data conversion to dtransform as well as tests for a trivial
data transform expression.
Platforms tested:
h5committest'ed, except used arabica instead of sol and didn't do on copper
b/c no logon there. Problem noted with mtime test...doesn't appear to be
related to anything having to do with data transforms.
Diffstat (limited to 'src')
-rw-r--r-- | src/H5D.c | 1 | ||||
-rw-r--r-- | src/H5Dio.c | 53 | ||||
-rw-r--r-- | src/H5P.c | 1 | ||||
-rw-r--r-- | src/H5Zprivate.h | 3 | ||||
-rw-r--r-- | src/H5Ztrans.c | 833 |
5 files changed, 417 insertions, 474 deletions
@@ -453,7 +453,6 @@ H5D_term_interface(void) FUNC_LEAVE_NOAPI(n) } - /*------------------------------------------------------------------------- * Function: H5D_xfer_xform_del * diff --git a/src/H5Dio.c b/src/H5Dio.c index 8fc2982..4bcc9aa 100644 --- a/src/H5Dio.c +++ b/src/H5Dio.c @@ -1207,18 +1207,22 @@ H5D_contig_read(hsize_t nelmts, H5D_t *dataset, const H5T_t *mem_type, HGOTO_ERROR (H5E_IO, H5E_READERROR, FAIL, "mem gather failed") } /* end if */ - /* + /* * Perform data type conversion. */ if (H5T_convert(tpath, src_id, dst_id, smine_nelmts, 0, 0, tconv_buf, bkg_buf, dxpl_id)<0) HGOTO_ERROR(H5E_DATASET, H5E_CANTINIT, FAIL, "data type conversion failed") + - - /* Added by LA to do data transforms */ + /* Do the data transform after the conversion (since we're using type mem_type) */ if(!H5Z_xform_noop(dxpl_cache->data_xform_prop)) - H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, mem_type); + { + if( H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, mem_type) < 0) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Error performing data transform") + } /* end of LA additions */ + /* * Scatter the data into memory. */ @@ -1445,12 +1449,23 @@ H5D_contig_write(hsize_t nelmts, H5D_t *dataset, const H5T_t *mem_type, const H5 HGOTO_ERROR (H5E_IO, H5E_WRITEERROR, FAIL, "file gather failed") } /* end if */ - /* + + + /* * Perform data type conversion. */ if (H5T_convert(tpath, src_id, dst_id, smine_nelmts, 0, 0, tconv_buf, bkg_buf, dxpl_id)<0) HGOTO_ERROR(H5E_DATASET, H5E_CANTINIT, FAIL, "data type conversion failed") + /* Do the data transform after the type conversion (since we're using dataset->type). */ + if(!H5Z_xform_noop(dxpl_cache->data_xform_prop)) + { + if( H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, dataset->type) < 0) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Error performing data transform") + } + /* End of LA modifications */ + + /* * Scatter the data out to the file. */ @@ -1458,15 +1473,11 @@ H5D_contig_write(hsize_t nelmts, H5D_t *dataset, const H5T_t *mem_type, const H5 H5_timer_begin(&timer); #endif - /* LA additions for data transforms */ - if(!H5Z_xform_noop(dxpl_cache->data_xform_prop)) - H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, dataset->type); - /* End of LA modifications */ - - status = H5S_select_fscat(dataset->ent.file, &(dataset->layout), + status = H5S_select_fscat(dataset->ent.file, &(dataset->layout), &dataset->dcpl_cache, (H5D_storage_t *)&(dataset->efl), file_space, &file_iter, smine_nelmts, dxpl_cache, dxpl_id, tconv_buf); + #ifdef H5S_DEBUG H5_timer_end(&(sconv->stats[0].scat_timer), &timer); sconv->stats[0].scat_nbytes += smine_nelmts * dst_type_size; @@ -1734,9 +1745,12 @@ UNUSED tconv_buf, bkg_buf, dxpl_id)<0) HGOTO_ERROR(H5E_DATASET, H5E_CANTINIT, FAIL, "data type conversion failed") - /* LA additions for data transforms */ + /* Do the data transform after the conversion (since we're using type mem_type) */ if(!H5Z_xform_noop(dxpl_cache->data_xform_prop)) - H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, mem_type); + { + if( H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, mem_type) < 0) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Error performing data transform") + } /* * Scatter the data into memory. @@ -2070,17 +2084,22 @@ nelmts, H5D_t *dataset, const H5T_t *mem_type, const H5S_t *mem_space, if (n!=smine_nelmts) HGOTO_ERROR (H5E_IO, H5E_WRITEERROR, FAIL, "file gather failed") } /* end if */ - - /* + + + + /* * Perform data type conversion. */ if (H5T_convert(tpath, src_id, dst_id, smine_nelmts, 0, 0, tconv_buf, bkg_buf, dxpl_id)<0) HGOTO_ERROR(H5E_DATASET, H5E_CANTINIT, FAIL, "data type conversion failed") - /* LA additions for data transforms */ + /* Do the data transform after the type conversion (since we're using dataset->type) */ if(!H5Z_xform_noop(dxpl_cache->data_xform_prop)) - H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, dataset->type); + { + if( H5Z_xform_eval(dxpl_cache->data_xform_prop, tconv_buf, smine_nelmts, dataset->type) < 0) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Error performing data transform") + } /* * Scatter the data out to the file. @@ -4455,6 +4455,7 @@ H5P_get(H5P_genplist_t *plist, const char *name, void *value) assert(plist); assert(name); + assert(value); /* Check if the property has been deleted */ diff --git a/src/H5Zprivate.h b/src/H5Zprivate.h index 33b7da9..0c43bf7 100644 --- a/src/H5Zprivate.h +++ b/src/H5Zprivate.h @@ -66,10 +66,11 @@ H5_DLL herr_t H5Z_delete(struct H5O_pline_t *pline, H5Z_filter_t filter); /* Data Transform Functions */ typedef struct H5Z_data_xform_t H5Z_data_xform_t; /* Defined in H5Ztrans.c */ + H5_DLL H5Z_data_xform_t *H5Z_xform_create(const char *expr); H5_DLL herr_t H5Z_xform_copy(H5Z_data_xform_t **data_xform_prop); H5_DLL herr_t H5Z_xform_destroy(H5Z_data_xform_t *data_xform_prop); -H5_DLL void H5Z_xform_eval(const H5Z_data_xform_t *data_xform_prop, void* array, hsize_t array_size, const H5T_t *buf_type); +H5_DLL herr_t H5Z_xform_eval(const H5Z_data_xform_t *data_xform_prop, void* array, hsize_t array_size, const H5T_t *buf_type); H5_DLL hbool_t H5Z_xform_noop(const H5Z_data_xform_t *data_xform_prop); #endif diff --git a/src/H5Ztrans.c b/src/H5Ztrans.c index 04f3a52..ff16bee 100644 --- a/src/H5Ztrans.c +++ b/src/H5Ztrans.c @@ -720,154 +720,103 @@ done: FUNC_LEAVE_NOAPI(new_node); } -#ifdef H5Z_XFORM_DEBUG - -/*------------------------------------------------------------------------- - * Function: H5Z_XFORM_DEBUG - * Purpose: Print out the expression in a nice format which displays - * the precedences explicitly with parentheses. - * Return: Nothing - * Programmer: Bill Wendling - * 26. August 2003 - * Modifications: - * - *------------------------------------------------------------------------- - */ - -static void -H5Z_XFORM_DEBUG(H5Z_node *tree) -{ - int i; - printf("Expression: "); - H5Z_result res; - res = H5Z_eval(tree); - if(res.type == H5Z_INTEGER) - printf("H5Z_result is: %d",res.value.int_val ); - else if(res.type == H5Z_FLOAT) - printf("H5Z_result is: %f",res.value.float_val ); - else - { - printf("H5Z_result is "); - for(i=0; i<5; i++) - printf("%d ", array[i]); - } - printf("\n"); -} - - /*------------------------------------------------------------------------- - * Function: H5Z_print - * Purpose: Print out the expression in a nice format which displays - * the precedences explicitly with parentheses. - * Return: Nothing - * Programmer: Bill Wendling - * 26. August 2003 + * Function: H5Z_eval_full + * Purpose: If the transform is trivial, this function applies it. + * Otherwise, it calls H5Z_xform_eval_full to do the full + * transform. + * Return: SUCCEED if transform applied succesfully, FAIL otherwise + * Programmer: Leon Arber + * 5/1/04 * Modifications: * *------------------------------------------------------------------------- */ -static void -H5Z_print(H5Z_node *tree, FILE *stream) -{ - /* check args */ - assert(stream); - - if (!tree) - return; - - if (tree->type == H5Z_INTEGER) { - fprintf(stream, "%ld", tree->value.int_val); - } else if (tree->type == H5Z_FLOAT) { - fprintf(stream, "%f", tree->value.float_val); - } else if (tree->type == SYMBOL) { - fprintf(stream, "%s", tree->value.sym_val); - } else { - fprintf(stream, "("); - H5Z_print(tree->lchild, stream); - - switch (tree->type) { - case PLUS: fprintf(stream, "+"); break; - case MINUS: fprintf(stream, "-"); break; - case MULT: fprintf(stream, "*"); break; - case DIVIDE: fprintf(stream, "/"); break; - default: fprintf(stream, "Invalid expression tree\n"); - return; - } - - H5Z_print(tree->rchild, stream); - fprintf(stream, ")"); - } -} -#endif /* H5Z_XFORM_DEBUG */ - - -void H5Z_xform_eval(const H5Z_data_xform_t *data_xform_prop, void* array, hsize_t array_size, const H5T_t *buf_type) +herr_t H5Z_xform_eval(const H5Z_data_xform_t *data_xform_prop, void* array, hsize_t array_size, const H5T_t *buf_type) { H5Z_node *tree; hid_t array_type; unsigned int i; int n; float f; - - FUNC_ENTER_NOAPI_NOINIT_NOFUNC(H5Z_xform_eval) - - assert(data_xform_prop); + H5Z_result res; + herr_t ret_value = SUCCEED; + + FUNC_ENTER_NOAPI(H5Z_xform_eval, FAIL) + + assert(data_xform_prop); tree=data_xform_prop->parse_root; /* Get the datatype ID for the buffer's type */ - array_type = H5Z_xform_find_type(buf_type); - - if( tree->type == H5Z_INTEGER) - { - if(array_type == H5T_NATIVE_INT) - { - n = tree->value.int_val; - for(i=0; i<array_size; i++) - *((int*)array + i) = n; - - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = (float)tree->value.int_val; - for(i=0; i<array_size; i++) - *((float*)array + i) = f; - } - } - else if (tree->type == H5Z_FLOAT) - { - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = (int)tree->value.float_val; - for(i=0; i<array_size; i++) - *((int*)array + i) = n; - - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = tree->value.float_val; - for(i=0; i<array_size; i++) - *((float*)array + i) = f; - } - } - } - else - H5Z_eval_full(tree, array, array_size, array_type); - - FUNC_LEAVE_NOAPI_VOID + if( (array_type = H5Z_xform_find_type(buf_type)) < 0) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Could not find matching data type for data transform.") + + /* If it's a trivial data transform, perform it */ + if( tree->type == H5Z_INTEGER) + { + if(array_type == H5T_NATIVE_INT) + { + n = tree->value.int_val; + for(i=0; i<array_size; i++) + *((int*)array + i) = n; + + } + else if(array_type == H5T_NATIVE_FLOAT) + { + f = (float)tree->value.int_val; + for(i=0; i<array_size; i++) + *((float*)array + i) = f; + } + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Unsupported data type for data transform") + + + } + else if (tree->type == H5Z_FLOAT) + { + if(array_type == H5T_NATIVE_INT) + { + n = (int)tree->value.float_val; + for(i=0; i<array_size; i++) + *((int*)array + i) = n; + + } + else if(array_type == H5T_NATIVE_FLOAT) + { + f = tree->value.float_val; + for(i=0; i<array_size; i++) + *((float*)array + i) = f; + } + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "Unsupported data type for data transform") + + } + /* Otherwise, do the full data transform */ + else + { + res = H5Z_eval_full(tree, array, array_size, array_type); + if(res.type == ERROR) + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, FAIL, "error while performing data transform") + } + + +done: + FUNC_LEAVE_NOAPI(ret_value); } /*------------------------------------------------------------------------- - * Function: H5Z_eval + * Function: H5Z_eval_full + * Purpose: Does a full evaluation of the parse tree contained in tree + * and applies this transform to array. * Return: Nothing * Programmer: Leon Arber + * 5/1/04 * Modifications: * *------------------------------------------------------------------------- @@ -876,20 +825,13 @@ static H5Z_result H5Z_eval_full(H5Z_node *tree, void* array, hsize_t array_size, hid_t array_type) { - H5Z_result res, resl, resr; - int n; + H5Z_result res, resl, resr, ret_value, error; unsigned int i; - float f; - H5Z_result ret_value; - ret_value.type = ERROR; + + error.type = ERROR; FUNC_ENTER_NOAPI(H5Z_eval_full, ret_value); - /*if (!tree) - { - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, NULL, "tree was NULL") - }*/ - /* check args */ assert(tree); @@ -922,161 +864,155 @@ H5Z_eval_full(H5Z_node *tree, void* array, hsize_t array_size, hid_t array_type if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) { res.type = SYMBOL; - for(i=0; i<array_size; i++) - { if(array_type == H5T_NATIVE_INT) { - n = resr.value.int_val + *((int*)array + i); - *((int*)array + i) = n; + for(i=0; i<array_size; i++) + *((int*)array + i) = resr.value.int_val + *((int*)array + i); } else if(array_type == H5T_NATIVE_FLOAT) { - f = resr.value.int_val +*((float*)array + i); - *((float*)array + i) = f; + for(i=0; i<array_size; i++) + *((float*)array + i) = resr.value.int_val + *((float*)array + i); } - } - HGOTO_DONE(res) + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + HGOTO_DONE(res) } else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) { res.type = SYMBOL; - for(i=0; i<array_size; i++) - { if(array_type == H5T_NATIVE_INT) { - n = resr.value.float_val + *((int*)array + i); - *((int*)array + i) = n; + for(i=0; i<array_size; i++) + *((int*)array + i) += resr.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) += resr.value.float_val; } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resr.value.float_val +*((float*)array + i); - *((float*)array + i) = f; - } - } + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + HGOTO_DONE(res) } else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) { res.type = SYMBOL; - for(i=0; i<array_size; i++) - { if(array_type == H5T_NATIVE_INT) { - n = resl.value.int_val + *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.int_val +*((float*)array + i); - *((float*)array + i) = f; - } - } + for(i=0; i<array_size; i++) + *((int*)array + i) += resl.value.int_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) += resl.value.int_val ; + } + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") HGOTO_DONE(res) } - else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) - { - - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.float_val + *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.float_val +*((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) + else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) + { + + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) += resl.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) += resl.value.float_val ; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) } else - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, ret_value, "Unexpected type conversion operation") + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Unexpected type conversion operation") HGOTO_DONE(res) break; - case MINUS: - - if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = *((int*)array + i) - resr.value.int_val; - *((int*)array + i) = n; - /* array[i] += resr.value.int_val;*/ - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f =*((float*)array + i) - resr.value.int_val; - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) /*we can't upgrade an array w/o allocating more memory, so we downgrade the float_val to an int.*/ - { + case MINUS: + + if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) -= resr.value.int_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) -= resr.value.int_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } + else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) /*we can't upgrade an array w/o allocating more memory, so we downgrade the float_val to an int.*/ + { res.type = SYMBOL; - for(i=0; i<array_size; i++) - { if(array_type == H5T_NATIVE_INT) { - n = *((int*)array + i) - resr.value.float_val; - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f =*((float*)array + i) - resr.value.float_val; - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.int_val - *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.int_val -*((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - - } - else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) + for(i=0; i<array_size; i++) + *((int*)array + i) -= resr.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) -= resr.value.float_val; + } + else + { + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + } + HGOTO_DONE(res) + } + else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) = resl.value.int_val - *((int*)array + i); + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) = resl.value.int_val - *((float*)array + i); + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + + } + else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.float_val - *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.float_val -*((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) = resl.value.float_val - *((int*)array + i); + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) = resl.value.float_val - *((float*)array + i); + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } else - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, ret_value, "Unexpected type conversion operation") + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Unexpected type conversion operation") HGOTO_DONE(res) @@ -1085,141 +1021,136 @@ H5Z_eval_full(H5Z_node *tree, void* array, hsize_t array_size, hid_t array_type case MULT: - if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resr.value.int_val * *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resr.value.int_val **((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resr.value.float_val * *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resr.value.float_val **((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.int_val * *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.int_val **((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - - } - else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) - { - - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.float_val * *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.float_val **((float*)array + i); - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, ret_value, "Unexpected type operation") - - HGOTO_DONE(res) - break; - - - case DIVIDE: + if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) *= resr.value.int_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) *= resr.value.int_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } + else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) *= resr.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) *= resr.value.float_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } + else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) *= resl.value.int_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) *= resl.value.int_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + + } + else if( (resr.type == SYMBOL) && (resl.type==H5Z_FLOAT)) + { + + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) *= resl.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) *= resl.value.float_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } + else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Unexpected type operation") + + HGOTO_DONE(res) + break; + + + case DIVIDE: + + if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) /= resr.value.int_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) /= resr.value.int_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + HGOTO_DONE(res) + } + else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) + { + res.type = SYMBOL; + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) /= resr.value.float_val; + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) /= resr.value.float_val; + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") - if( (resl.type == SYMBOL) && (resr.type==H5Z_INTEGER)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = *((int*)array + i) / resr.value.int_val; - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f =*((float*)array + i) / resr.value.int_val ; - *((float*)array + i) = f; - } - } - HGOTO_DONE(res) - } - else if( (resl.type == SYMBOL) && (resr.type==H5Z_FLOAT)) - { - res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = *((int*)array + i) / resr.value.float_val; - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f =*((float*)array + i) / resr.value.float_val; - *((float*)array + i) = f; - } - } HGOTO_DONE(res) } else if( (resr.type == SYMBOL) && (resl.type==H5Z_INTEGER)) { res.type = SYMBOL; - for(i=0; i<array_size; i++) { if(array_type == H5T_NATIVE_INT) - { - n = resl.value.int_val / *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.int_val / *((float*)array + i); - *((float*)array + i) = f; - } - } + { + for(i=0; i<array_size; i++) + *((int*)array + i) = resl.value.int_val / *((int*)array + i); + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) = resl.value.int_val / *((float*)array + i); + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + + } HGOTO_DONE(res) } @@ -1227,23 +1158,24 @@ H5Z_eval_full(H5Z_node *tree, void* array, hsize_t array_size, hid_t array_type { res.type = SYMBOL; - for(i=0; i<array_size; i++) - { - if(array_type == H5T_NATIVE_INT) - { - n = resl.value.float_val / *((int*)array + i); - *((int*)array + i) = n; - } - else if(array_type == H5T_NATIVE_FLOAT) - { - f = resl.value.float_val / *((float*)array + i); - *((float*)array + i) = f; - } + { + if(array_type == H5T_NATIVE_INT) + { + for(i=0; i<array_size; i++) + *((int*)array + i) = resl.value.float_val / *((int*)array + i); + } + else if(array_type == H5T_NATIVE_FLOAT) + { + for(i=0; i<array_size; i++) + *((float*)array + i) = resl.value.float_val / *((float*)array + i); + }else + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Cannot perform a data transform on this type") + } HGOTO_DONE(res) } else - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, ret_value, "Unexpected type operation") + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Unexpected type operation") HGOTO_DONE(res) break; @@ -1251,12 +1183,14 @@ H5Z_eval_full(H5Z_node *tree, void* array, hsize_t array_size, hid_t array_type default: - HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, ret_value, "Invalid expression tree") + HGOTO_ERROR(H5E_ARGS, H5E_BADVALUE, error, "Invalid expression tree") } } done: - FUNC_LEAVE_NOAPI(res) + + FUNC_LEAVE_NOAPI(ret_value) + } @@ -1273,80 +1207,63 @@ done: static hid_t H5Z_xform_find_type(const H5T_t* type) { hid_t ret_value = SUCCEED; - + FUNC_ENTER_NOAPI_NOINIT(H5Z_xform_find_type); assert(type); -/* H5T_NATIVE_CHAR - H5T_NATIVE_SHORT - H5T_NATIVE_INT - H5T_NATIVE_LONG - H5T_NATIVE_LLONG - - H5T_NATIVE_UCHAR - H5T_NATIVE_USHORT - H5T_NATIVE_UINT - H5T_NATIVE_ULONG - H5T_NATIVE_ULLONG - - H5T_NATIVE_FLOAT - H5T_NATIVE_DOUBLE - H5T_NATIVE_LDOUBLE -*/ - /* Check for SHORT type */ if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_SHORT,H5I_DATATYPE) ))==0) - HGOTO_DONE(H5T_NATIVE_SHORT) + HGOTO_DONE(H5T_NATIVE_SHORT) - /* Check for INT type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_INT,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_INT) + /* Check for INT type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_INT,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_INT) - /* Check for LONG type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LONG,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_LONG) + /* Check for LONG type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LONG,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_LONG) - /* Check for LONGLONG type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LLONG,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_LLONG) + /* Check for LONGLONG type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LLONG,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_LLONG) - /* Check for UCHAR type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_UCHAR,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_UCHAR) + /* Check for UCHAR type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_UCHAR,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_UCHAR) - /* Check for USHORT type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_USHORT,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_USHORT) + /* Check for USHORT type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_USHORT,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_USHORT) - /* Check for UINT type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_UINT,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_UINT) + /* Check for UINT type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_UINT,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_UINT) - /* Check for ULONG type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_ULONG,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_ULONG) + /* Check for ULONG type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_ULONG,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_ULONG) - /* Check for ULONGLONG type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_ULLONG,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_ULLONG) + /* Check for ULONGLONG type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_ULLONG,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_ULLONG) - /* Check for FLOAT type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_FLOAT,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_FLOAT) + /* Check for FLOAT type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_FLOAT,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_FLOAT) - /* Check for DOUBLE type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_DOUBLE,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_DOUBLE) + /* Check for DOUBLE type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_DOUBLE,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_DOUBLE) - /* Check for LONGDOUBLE type */ - else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LDOUBLE,H5I_DATATYPE)))==0) - HGOTO_DONE(H5T_NATIVE_LDOUBLE) + /* Check for LONGDOUBLE type */ + else if((H5T_cmp(type, H5I_object_verify(H5T_NATIVE_LDOUBLE,H5I_DATATYPE)))==0) + HGOTO_DONE(H5T_NATIVE_LDOUBLE) else - HGOTO_ERROR(H5E_ARGS, H5E_BADTYPE, FAIL, "could not find matching type"); + HGOTO_ERROR(H5E_ARGS, H5E_BADTYPE, FAIL, "could not find matching type"); + - done: FUNC_LEAVE_NOAPI(ret_value) @@ -1473,11 +1390,12 @@ void* H5Z_xform_copy_tree(H5Z_node* tree) */ void H5Z_xform_reduce_tree(H5Z_node* tree) { - + hid_t ret_value = SUCCEED; + FUNC_ENTER_NOAPI_NOINIT_NOFUNC(H5Z_xform_reduce_tree) if(!tree) - goto done; + HGOTO_DONE(SUCCEED) if((tree->type == PLUS) || (tree->type == DIVIDE) ||(tree->type == MULT) ||(tree->type == MINUS)) { @@ -1487,14 +1405,19 @@ void H5Z_xform_reduce_tree(H5Z_node* tree) { H5Z_xform_reduce_tree(tree->lchild); if(((tree->lchild->type == H5Z_INTEGER) || (tree->lchild->type == H5Z_FLOAT)) && ((tree->rchild->type == H5Z_INTEGER) || (tree->rchild->type == H5Z_FLOAT))) - H5Z_do_op(tree); + { + H5Z_do_op(tree); + HGOTO_DONE(SUCCEED) + } H5Z_xform_reduce_tree(tree->rchild); if(((tree->lchild->type == H5Z_INTEGER) || (tree->lchild->type == H5Z_FLOAT)) && ((tree->rchild->type == H5Z_INTEGER) || (tree->rchild->type == H5Z_FLOAT))) - H5Z_do_op(tree); + { + H5Z_do_op(tree); + HGOTO_DONE(SUCCEED) + } } } - goto done; done: FUNC_LEAVE_NOAPI_VOID; |