diff options
| -rw-r--r-- | hl/src/H5DS.c | 3245 | ||||
| -rw-r--r-- | hl/test/test_ds.c | 4942 |
2 files changed, 4151 insertions, 4036 deletions
diff --git a/hl/src/H5DS.c b/hl/src/H5DS.c index f0c83f8..83bb1ed 100644 --- a/hl/src/H5DS.c +++ b/hl/src/H5DS.c @@ -1,17 +1,17 @@ /* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * - * Copyright by The HDF Group. * - * Copyright by the Board of Trustees of the University of Illinois. * - * All rights reserved. * - * * - * This file is part of HDF5. The full HDF5 copyright notice, including * - * terms governing use, modification, and redistribution, is contained in * - * the files COPYING and Copyright.html. COPYING can be found at the root * - * of the source code distribution tree; Copyright.html can be found at the * - * root level of an installed copy of the electronic HDF5 document set and * - * is linked from the top-level documents page. It can also be found at * - * http://hdfgroup.org/HDF5/doc/Copyright.html. If you do not have * - * access to either file, you may request a copy from help@hdfgroup.org. * - * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ +* Copyright by The HDF Group. * +* Copyright by the Board of Trustees of the University of Illinois. * +* All rights reserved. * +* * +* This file is part of HDF5. The full HDF5 copyright notice, including * +* terms governing use, modification, and redistribution, is contained in * +* the files COPYING and Copyright.html. COPYING can be found at the root * +* of the source code distribution tree; Copyright.html can be found at the * +* root level of an installed copy of the electronic HDF5 document set and * +* is linked from the top-level documents page. It can also be found at * +* http://hdfgroup.org/HDF5/doc/Copyright.html. If you do not have * +* access to either file, you may request a copy from help@hdfgroup.org. * +* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */ #include <assert.h> #include <stdlib.h> @@ -24,1203 +24,1277 @@ /*------------------------------------------------------------------------- - * Function: H5DSset_scale - * - * Purpose: The dataset DSID is converted to a Dimension Scale dataset. - * Creates the CLASS attribute, set to the value "DIMENSION_SCALE" - * and an empty REFERENCE_LIST attribute. - * If DIMNAME is specified, then an attribute called NAME is created, - * with the value DIMNAME. - * - * Return: Success: SUCCEED, Failure: FAIL - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 04, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSset_scale +* +* Purpose: The dataset DSID is converted to a Dimension Scale dataset. +* Creates the CLASS attribute, set to the value "DIMENSION_SCALE" +* and an empty REFERENCE_LIST attribute. +* If DIMNAME is specified, then an attribute called NAME is created, +* with the value DIMNAME. +* +* Return: Success: SUCCEED, Failure: FAIL +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 04, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DSset_scale(hid_t dsid, const char *dimname) { - int has_dimlist; - H5I_type_t it; + int has_dimlist; + H5I_type_t it; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* get ID type */ - if ((it = H5Iget_type(dsid)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* get ID type */ + if ((it = H5Iget_type(dsid)) < 0) + return FAIL; - if (H5I_DATASET!=it) - return FAIL; + if (H5I_DATASET!=it) + return FAIL; -/*------------------------------------------------------------------------- - * check if the dataset is a dataset wich has references to dimension scales - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * check if the dataset is a dataset wich has references to dimension scales + *------------------------------------------------------------------------- + */ - /* try to find the attribute "DIMENSION_LIST" */ - if ((has_dimlist = H5LT_find_attribute(dsid,DIMENSION_LIST)) < 0) - return FAIL; + /* try to find the attribute "DIMENSION_LIST" */ + if ((has_dimlist = H5LT_find_attribute(dsid,DIMENSION_LIST)) < 0) + return FAIL; - if (has_dimlist == 1) - return FAIL; + if (has_dimlist == 1) + return FAIL; -/*------------------------------------------------------------------------- - * write the standard attributes for a Dimension Scale dataset - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * write the standard attributes for a Dimension Scale dataset + *------------------------------------------------------------------------- + */ - if (H5LT_set_attribute_string(dsid,"CLASS",DIMENSION_SCALE_CLASS) < 0) - return FAIL; + if (H5LT_set_attribute_string(dsid,"CLASS",DIMENSION_SCALE_CLASS) < 0) + return FAIL; - if (dimname!=NULL) - { - if (H5LT_set_attribute_string(dsid,"NAME",dimname) < 0) - return FAIL; - } + if (dimname!=NULL) + { + if (H5LT_set_attribute_string(dsid,"NAME",dimname) < 0) + return FAIL; + } - return SUCCEED; + return SUCCEED; } /*------------------------------------------------------------------------- - * Function: H5DSattach_scale - * - * Purpose: Define Dimension Scale DSID to be associated with dimension IDX - * of Dataset DID. Entries are created in the DIMENSION_LIST and - * REFERENCE_LIST attributes. - * - * Return: - * Success: SUCCEED - * Failure: FAIL - * - * Fails if: Bad arguments - * If DSID is not a Dimension Scale - * If DID is a Dimension Scale (A Dimension Scale cannot have scales) - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: December 20, 2004 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSattach_scale +* +* Purpose: Define Dimension Scale DSID to be associated with dimension IDX +* of Dataset DID. Entries are created in the DIMENSION_LIST and +* REFERENCE_LIST attributes. +* +* Return: +* Success: SUCCEED +* Failure: FAIL +* +* Fails if: Bad arguments +* If DSID is not a Dimension Scale +* If DID is a Dimension Scale (A Dimension Scale cannot have scales) +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: December 20, 2004 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DSattach_scale(hid_t did, hid_t dsid, unsigned int idx) { - int has_dimlist; - int has_reflist; - int is_ds; - hssize_t nelmts; - hid_t sid; /* space ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t aid = -1; /* attribute ID */ - int rank; /* rank of dataset */ - hsize_t *dims; /* dimension of the "REFERENCE_LIST" array */ - ds_list_t dsl; /* attribute data in the DS pointing to the dataset */ - ds_list_t *dsbuf; /* array of attribute data in the DS pointing to the dataset */ - hobj_ref_t ref_to_ds; /* reference to the DS */ - hobj_ref_t ref_j; /* iterator reference */ - hvl_t *buf; /* VL buffer to store in the attribute */ - hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ - H5O_info_t oi1, oi2; - H5I_type_t it1, it2; - int i, len; - int found_ds=0; + int has_dimlist; + int has_reflist; + int is_ds; + hssize_t nelmts; + hid_t sid; /* space ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t aid = -1; /* attribute ID */ + int rank; /* rank of dataset */ + hsize_t *dims=NULL; /* dimension of the "REFERENCE_LIST" array */ + ds_list_t dsl; /* attribute data in the DS pointing to the dataset */ + ds_list_t *dsbuf=NULL; /* array of attribute data in the DS pointing to the dataset */ + hobj_ref_t ref_to_ds; /* reference to the DS */ + hobj_ref_t ref_j; /* iterator reference */ + hvl_t *buf=NULL; /* VL buffer to store in the attribute */ + hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ + H5O_info_t oi1, oi2; + H5I_type_t it1, it2; + int i, len; + int found_ds=0; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* the dataset cannot be a DS dataset */ - if ((H5DSis_scale(did))==1) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* the dataset cannot be a DS dataset */ + if ((H5DSis_scale(did))==1) + return FAIL; - /* get info for the dataset in the parameter list */ - if(H5Oget_info(did, &oi1) < 0) - return FAIL; + /* get info for the dataset in the parameter list */ + if(H5Oget_info(did, &oi1) < 0) + return FAIL; - /* get info for the scale in the parameter list */ - if(H5Oget_info(dsid, &oi2) < 0) - return FAIL; + /* get info for the scale in the parameter list */ + if(H5Oget_info(dsid, &oi2) < 0) + return FAIL; - /* same object, not valid */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) - return FAIL; + /* same object, not valid */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) + return FAIL; - /* get ID type */ - if ((it1 = H5Iget_type(did)) < 0) - return FAIL; - if ((it2 = H5Iget_type(dsid)) < 0) - return FAIL; + /* get ID type */ + if ((it1 = H5Iget_type(did)) < 0) + return FAIL; + if ((it2 = H5Iget_type(dsid)) < 0) + return FAIL; - if (H5I_DATASET!=it1 || H5I_DATASET!=it2) - return FAIL; + if (H5I_DATASET!=it1 || H5I_DATASET!=it2) + return FAIL; - /* the DS dataset cannot have dimension scales */ - if (H5LT_find_attribute(dsid,DIMENSION_LIST)==1) - return FAIL; + /* the DS dataset cannot have dimension scales */ + if (H5LT_find_attribute(dsid,DIMENSION_LIST)==1) + return FAIL; - /* check if the dataset is a "reserved" dataset (image, table) */ - if (H5DS_is_reserved(did)==1) - return FAIL; + /* check if the dataset is a "reserved" dataset (image, table) */ + if (H5DS_is_reserved(did)==1) + return FAIL; -/*------------------------------------------------------------------------- - * The dataset may or may not have the associated DS attribute - * First we try to open to see if it is already there; if not, it is created. - * If it exists, the array of references is extended to hold the reference - * to the new DS - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * The dataset may or may not have the associated DS attribute + * First we try to open to see if it is already there; if not, it is created. + * If it exists, the array of references is extended to hold the reference + * to the new DS + *------------------------------------------------------------------------- + */ - /* get dataset space */ - if ((sid = H5Dget_space(did)) < 0) - return FAIL; + /* get dataset space */ + if ((sid = H5Dget_space(did)) < 0) + return FAIL; - /* get rank */ - if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; + /* get rank */ + if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; - /* scalar rank */ - if (rank==0) - rank=1; + /* scalar rank */ + if (rank==0) + rank=1; - /* close dataset space */ - if (H5Sclose(sid) < 0) - return FAIL; + /* close dataset space */ + if (H5Sclose(sid) < 0) + return FAIL; - /* parameter range checking */ - if (idx>(unsigned)rank-1) - return FAIL; + /* parameter range checking */ + if (idx>(unsigned)rank-1) + return FAIL; -/*------------------------------------------------------------------------- - * two references are created: one to the DS, saved in "DIMENSION_LIST" - * and one to the dataset, saved in "REFERENCE_LIST" - *------------------------------------------------------------------------- - */ - /* create a reference for the >>DS<< dataset */ - if (H5Rcreate(&ref_to_ds,dsid,".",H5R_OBJECT,-1) < 0) - return FAIL; - - /* create a reference for the >>data<< dataset */ - if (H5Rcreate(&dsl.ref,did,".",H5R_OBJECT,-1) < 0) - return FAIL; - - /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ - if ((has_dimlist = H5LT_find_attribute(did,DIMENSION_LIST)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * two references are created: one to the DS, saved in "DIMENSION_LIST" + * and one to the dataset, saved in "REFERENCE_LIST" + *------------------------------------------------------------------------- + */ + /* create a reference for the >>DS<< dataset */ + if (H5Rcreate(&ref_to_ds,dsid,".",H5R_OBJECT,-1) < 0) + return FAIL; -/*------------------------------------------------------------------------- - * it does not exist. we create the attribute and its reference data - *------------------------------------------------------------------------- - */ - if (has_dimlist == 0) - { - /* create one entry array */ - dims = (hsize_t*) malloc (1 * sizeof (hsize_t)); - - if(dims == NULL) - return FAIL; - - dims[0] = rank; - - /* space for the attribute */ - if((sid = H5Screate_simple(1, dims, NULL)) < 0) - return FAIL; - - /* create the type for the attribute "DIMENSION_LIST" */ - if((tid = H5Tvlen_create(H5T_STD_REF_OBJ)) < 0) - goto out; - - /* create the attribute */ - if((aid = H5Acreate2(did, DIMENSION_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - - if(buf == NULL) - goto out; - - for(i = 0; i < rank; i++) { - buf[i].len = 0; - buf[i].p = NULL; - } - - /* store the REF information in the index of the dataset that has the DS */ - buf[idx].len = 1; - buf[idx].p = malloc( 1 * sizeof(hobj_ref_t)); - ((hobj_ref_t *)buf[idx].p)[0] = ref_to_ds; - - /* write the attribute with the reference */ - if (H5Awrite(aid,tid,buf) < 0) - goto out; - - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - - if (dims) - free(dims); - if (buf) - free(buf); - - } + /* create a reference for the >>data<< dataset */ + if (H5Rcreate(&dsl.ref,did,".",H5R_OBJECT,-1) < 0) + return FAIL; -/*------------------------------------------------------------------------- - * the attribute already exists, open it, extend the buffer, - * and insert the new reference - *------------------------------------------------------------------------- - */ - - else if ( has_dimlist == 1 ) - { - if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) - goto out; - - if((tid = H5Aget_type(aid)) < 0) - goto out; - - if((sid = H5Aget_space(aid)) < 0) - goto out; - - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - - if(buf == NULL) - goto out; - - /* read */ - if(H5Aread(aid, tid, buf) < 0) - goto out; - - /* check to avoid inserting duplicates. it is not FAIL, just do nothing */ - /* iterate all the REFs in this dimension IDX */ - for (i=0; i<(int)buf[idx].len; i++) - { - /* get the reference */ - ref_j = ((hobj_ref_t *)buf[idx].p)[i]; - - /* get the scale id for this REF */ - if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref_j)) < 0) - goto out; - - /* get info for DS in the parameter list */ - if(H5Oget_info(dsid, &oi1) < 0) - goto out; - - /* get info for this DS */ - if(H5Oget_info(dsid_j, &oi2) < 0) - goto out; - - /* same object, so this DS scale is already in this DIM IDX */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) - found_ds = 1; - - /* close the dereferenced dataset */ - if (H5Dclose(dsid_j) < 0) - goto out; - } - - if (found_ds==0) - { - /* we are adding one more DS to this dimension */ - if ( buf[idx].len > 0 ) - { - buf[idx].len++; - len = buf[idx].len; - buf[idx].p = realloc( buf[idx].p, len * sizeof(hobj_ref_t)); - ((hobj_ref_t *)buf[idx].p)[ len-1 ] = ref_to_ds; - } - else - { - /* store the REF information in the index of the dataset that has the DS */ - buf[idx].len = 1; - buf[idx].p = malloc( 1 * sizeof(hobj_ref_t)); - ((hobj_ref_t *)buf[idx].p)[0] = ref_to_ds; - } - } - - /* write the attribute with the new references */ - if (H5Awrite(aid,tid,buf) < 0) - goto out; - - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (buf) - free(buf); - - } /* has_dimlist */ + /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ + if ((has_dimlist = H5LT_find_attribute(did,DIMENSION_LIST)) < 0) + return FAIL; -/*------------------------------------------------------------------------- - * save DS info on the >>DS<< dataset - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * it does not exist. we create the attribute and its reference data + *------------------------------------------------------------------------- + */ + if (has_dimlist == 0) + { + /* create one entry array */ + dims = (hsize_t*) malloc (1 * sizeof (hsize_t)); - /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ - if ((has_reflist = H5LT_find_attribute(dsid,REFERENCE_LIST)) < 0) - goto out; + if(dims == NULL) + return FAIL; -/*------------------------------------------------------------------------- - * it does not exist. we create the attribute and its reference data - *------------------------------------------------------------------------- - */ - if (has_reflist == 0) - { - /* create one entry array */ - dims = (hsize_t*) malloc (1 * sizeof (hsize_t)); - - if (dims == NULL) - goto out; - - dims[0] = 1; - - /* space for the attribute */ - if ((sid = H5Screate_simple(1,dims,NULL)) < 0) - goto out; - - /* create the compound datatype for the attribute "REFERENCE_LIST" */ - if((tid = H5Tcreate(H5T_COMPOUND, sizeof(ds_list_t))) < 0) - goto out; - - /* insert reference field */ - if(H5Tinsert(tid, "dataset", HOFFSET(ds_list_t,ref), H5T_STD_REF_OBJ) < 0) - goto out; - - /* insert dimension idx of the dataset field */ - if(H5Tinsert(tid, "dimension", HOFFSET(ds_list_t, dim_idx), H5T_NATIVE_INT) < 0) - goto out; - - /* create the attribute */ - if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - - /* store the IDX information */ - dsl.dim_idx = idx; - - /* write the attribute with the reference */ - if(H5Awrite(aid, tid, &dsl) < 0) - goto out; - - /* close */ - if(H5Sclose(sid) < 0) - goto out; - if(H5Tclose(tid) < 0) - goto out; - if(H5Aclose(aid) < 0) - goto out; - - if(dims) - free(dims); - } + dims[0] = rank; -/*------------------------------------------------------------------------- - * the "REFERENCE_LIST" array already exists, open it and extend it - *------------------------------------------------------------------------- - */ + /* space for the attribute */ + if((sid = H5Screate_simple(1, dims, NULL)) < 0) + return FAIL; - else if(has_reflist == 1) - { - if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) - goto out; + /* create the type for the attribute "DIMENSION_LIST" */ + if((tid = H5Tvlen_create(H5T_STD_REF_OBJ)) < 0) + goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; + /* create the attribute */ + if((aid = H5Acreate2(did, DIMENSION_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; - /* get and save the old reference(s) */ - if((sid = H5Aget_space(aid)) < 0) - goto out; + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) - goto out; + if(buf == NULL) + goto out; - nelmts++; + for(i = 0; i < rank; i++) { + buf[i].len = 0; + buf[i].p = NULL; + } - dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); + /* store the REF information in the index of the dataset that has the DS */ + buf[idx].len = 1; + buf[idx].p = malloc( 1 * sizeof(hobj_ref_t)); + ((hobj_ref_t *)buf[idx].p)[0] = ref_to_ds; - if (dsbuf == NULL) - goto out; + /* write the attribute with the reference */ + if (H5Awrite(aid,tid,buf) < 0) + goto out; - if (H5Aread(aid,tid,dsbuf) < 0) - goto out; + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; - /* close */ - if (H5Sclose(sid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; + if (dims) + { + free(dims); + dims = NULL; + } + if (buf) + { + free(buf); + buf = NULL; + } -/*------------------------------------------------------------------------- - * create a new attribute - *------------------------------------------------------------------------- - */ + } + + /*------------------------------------------------------------------------- + * the attribute already exists, open it, extend the buffer, + * and insert the new reference + *------------------------------------------------------------------------- + */ + + else if ( has_dimlist == 1 ) + { + if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) + goto out; - /* the attribute must be deleted, in order to the new one can reflect the changes*/ - if(H5Adelete(dsid, REFERENCE_LIST) < 0) - goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; - /* store the IDX information (index of the dataset that has the DS) */ - dsl.dim_idx = idx; - dsbuf[nelmts-1] = dsl; + if((sid = H5Aget_space(aid)) < 0) + goto out; - /* create a new data space for the new references array */ - dims = (hsize_t *)malloc((size_t)nelmts * sizeof(hsize_t)); - if(dims == NULL) - goto out; - dims[0] = nelmts; + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - if((sid = H5Screate_simple(1, dims, NULL)) < 0) - goto out; + if(buf == NULL) + goto out; + + /* read */ + if(H5Aread(aid, tid, buf) < 0) + goto out; - if(dims) - free(dims); + /* check to avoid inserting duplicates. it is not FAIL, just do nothing */ + /* iterate all the REFs in this dimension IDX */ + for (i=0; i<(int)buf[idx].len; i++) + { + /* get the reference */ + ref_j = ((hobj_ref_t *)buf[idx].p)[i]; - /* create the attribute again with the changes of space */ - if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; + /* get the scale id for this REF */ + if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref_j)) < 0) + goto out; - /* write the attribute with the new references */ - if(H5Awrite(aid, tid, dsbuf) < 0) - goto out; + /* get info for DS in the parameter list */ + if(H5Oget_info(dsid, &oi1) < 0) + goto out; - /* close */ - if(H5Sclose(sid) < 0) - goto out; - if(H5Tclose(tid) < 0) - goto out; - if(H5Aclose(aid) < 0) - goto out; + /* get info for this DS */ + if(H5Oget_info(dsid_j, &oi2) < 0) + goto out; - if (dsbuf) - free(dsbuf); + /* same object, so this DS scale is already in this DIM IDX */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) + found_ds = 1; - } /* has_reflist */ + /* close the dereferenced dataset */ + if (H5Dclose(dsid_j) < 0) + goto out; + } -/*------------------------------------------------------------------------- - * write the standard attributes for a Dimension Scale dataset - *------------------------------------------------------------------------- - */ + if (found_ds==0) + { + /* we are adding one more DS to this dimension */ + if ( buf[idx].len > 0 ) + { + buf[idx].len++; + len = buf[idx].len; + buf[idx].p = realloc( buf[idx].p, len * sizeof(hobj_ref_t)); + ((hobj_ref_t *)buf[idx].p)[ len-1 ] = ref_to_ds; + } + else + { + /* store the REF information in the index of the dataset that has the DS */ + buf[idx].len = 1; + buf[idx].p = malloc( 1 * sizeof(hobj_ref_t)); + ((hobj_ref_t *)buf[idx].p)[0] = ref_to_ds; + } + } + + /* write the attribute with the new references */ + if (H5Awrite(aid,tid,buf) < 0) + goto out; + + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (buf) + { + free(buf); + buf = NULL; + } + + } /* has_dimlist */ + + /*------------------------------------------------------------------------- + * save DS info on the >>DS<< dataset + *------------------------------------------------------------------------- + */ + + /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ + if ((has_reflist = H5LT_find_attribute(dsid,REFERENCE_LIST)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * it does not exist. we create the attribute and its reference data + *------------------------------------------------------------------------- + */ + if (has_reflist == 0) + { + /* create one entry array */ + dims = (hsize_t*) malloc (1 * sizeof (hsize_t)); + + if (dims == NULL) + goto out; + + dims[0] = 1; + + /* space for the attribute */ + if ((sid = H5Screate_simple(1,dims,NULL)) < 0) + goto out; + + /* create the compound datatype for the attribute "REFERENCE_LIST" */ + if((tid = H5Tcreate(H5T_COMPOUND, sizeof(ds_list_t))) < 0) + goto out; + + /* insert reference field */ + if(H5Tinsert(tid, "dataset", HOFFSET(ds_list_t,ref), H5T_STD_REF_OBJ) < 0) + goto out; + + /* insert dimension idx of the dataset field */ + if(H5Tinsert(tid, "dimension", HOFFSET(ds_list_t, dim_idx), H5T_NATIVE_INT) < 0) + goto out; + + /* create the attribute */ + if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + + /* store the IDX information */ + dsl.dim_idx = idx; + + /* write the attribute with the reference */ + if(H5Awrite(aid, tid, &dsl) < 0) + goto out; + + /* close */ + if(H5Sclose(sid) < 0) + goto out; + if(H5Tclose(tid) < 0) + goto out; + if(H5Aclose(aid) < 0) + goto out; + + if(dims) + { + free(dims); + dims = NULL; + } + } + + /*------------------------------------------------------------------------- + * the "REFERENCE_LIST" array already exists, open it and extend it + *------------------------------------------------------------------------- + */ + + else if(has_reflist == 1) + { + if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) + goto out; + + if((tid = H5Aget_type(aid)) < 0) + goto out; + + /* get and save the old reference(s) */ + if((sid = H5Aget_space(aid)) < 0) + goto out; + + if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) + goto out; + + nelmts++; + + dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); + + if (dsbuf == NULL) + goto out; + + if (H5Aread(aid,tid,dsbuf) < 0) + goto out; - if ((is_ds=H5DSis_scale(dsid)) < 0) - return FAIL; + /* close */ + if (H5Sclose(sid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + + /*------------------------------------------------------------------------- + * create a new attribute + *------------------------------------------------------------------------- + */ + + /* the attribute must be deleted, in order to the new one can reflect the changes*/ + if(H5Adelete(dsid, REFERENCE_LIST) < 0) + goto out; + + /* store the IDX information (index of the dataset that has the DS) */ + dsl.dim_idx = idx; + dsbuf[nelmts-1] = dsl; + + /* create a new data space for the new references array */ + dims = (hsize_t *)malloc((size_t)nelmts * sizeof(hsize_t)); + if(dims == NULL) + goto out; + dims[0] = nelmts; + + if((sid = H5Screate_simple(1, dims, NULL)) < 0) + goto out; + + if (dims) + { + free(dims); + dims = NULL; + } - if (is_ds == 0 ) - { - if (H5LT_set_attribute_string(dsid,"CLASS",DIMENSION_SCALE_CLASS) < 0) - return FAIL; - } + /* create the attribute again with the changes of space */ + if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; - return SUCCEED; + /* write the attribute with the new references */ + if(H5Awrite(aid, tid, dsbuf) < 0) + goto out; -/* error zone, gracefully close */ + /* close */ + if(H5Sclose(sid) < 0) + goto out; + if(H5Tclose(tid) < 0) + goto out; + if(H5Aclose(aid) < 0) + goto out; + + if (dsbuf) + { + free(dsbuf); + dsbuf = NULL; + } + + } /* has_reflist */ + + /*------------------------------------------------------------------------- + * write the standard attributes for a Dimension Scale dataset + *------------------------------------------------------------------------- + */ + + if ((is_ds=H5DSis_scale(dsid)) < 0) + return FAIL; + + if (is_ds == 0 ) + { + if (H5LT_set_attribute_string(dsid,"CLASS",DIMENSION_SCALE_CLASS) < 0) + return FAIL; + } + + return SUCCEED; + + /* error zone */ out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - } H5E_END_TRY; - return FAIL; + + if (dims) + { + free(dims); + dims = NULL; + } + if (buf) + { + free(buf); + buf = NULL; + } + if (dsbuf) + { + free(dsbuf); + dsbuf = NULL; + } + + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + } H5E_END_TRY; + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSdetach_scale - * - * Purpose: If possible, deletes association of Dimension Scale DSID with - * dimension IDX of Dataset DID. This deletes the entries in the - * DIMENSION_LIST and REFERENCE_LIST attributes. - * - * Return: - * Success: SUCCEED - * Failure: FAIL - * - * Fails if: Bad arguments - * The dataset DID or DSID do not exist. - * The DSID is not a Dimension Scale - * DSID is not attached to DID. - * Note that a scale may be associated with more than dimension of the same dataset. - * If so, the detach operation only deletes one of the associations, for DID. - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: December 20, 2004 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSdetach_scale +* +* Purpose: If possible, deletes association of Dimension Scale DSID with +* dimension IDX of Dataset DID. This deletes the entries in the +* DIMENSION_LIST and REFERENCE_LIST attributes. +* +* Return: +* Success: SUCCEED +* Failure: FAIL +* +* Fails if: Bad arguments +* The dataset DID or DSID do not exist. +* The DSID is not a Dimension Scale +* DSID is not attached to DID. +* Note that a scale may be associated with more than dimension of the same dataset. +* If so, the detach operation only deletes one of the associations, for DID. +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: December 20, 2004 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DSdetach_scale(hid_t did, hid_t dsid, unsigned int idx) { - int has_dimlist; - int has_reflist; - hssize_t nelmts; - hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ + int has_dimlist; + int has_reflist; + hssize_t nelmts; + hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ hid_t did_i; /* dataset ID in REFERENCE_LIST */ - hid_t sid; /* space ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t aid = -1; /* attribute ID */ - int rank; /* rank of dataset */ - ds_list_t *dsbuf = NULL; /* array of attribute data in the DS pointing to the dataset */ - ds_list_t *dsbufn = NULL; /* array of attribute data in the DS pointing to the dataset */ - hsize_t *dims = NULL; /* dimension of the "REFERENCE_LIST" array */ - hobj_ref_t ref; /* reference to the DS */ - hvl_t *buf = NULL; /* VL buffer to store in the attribute */ - unsigned i, j, jj; - H5O_info_t oi1, oi2, oi3, oi4; - int found_dset = 0, found_ds = 0; - H5I_type_t it1, it2; + hid_t sid; /* space ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t aid = -1; /* attribute ID */ + int rank; /* rank of dataset */ + ds_list_t *dsbuf = NULL; /* array of attribute data in the DS pointing to the dataset */ + ds_list_t *dsbufn = NULL; /* array of attribute data in the DS pointing to the dataset */ + hsize_t *dims = NULL; /* dimension of the "REFERENCE_LIST" array */ + hobj_ref_t ref; /* reference to the DS */ + hvl_t *buf = NULL; /* VL buffer to store in the attribute */ + unsigned i, j, jj; + H5O_info_t oi1, oi2, oi3, oi4; + int found_dset = 0, found_ds = 0; + H5I_type_t it1, it2; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* the dataset cannot be a DS dataset */ - if ((H5DSis_scale(did))==1) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* the dataset cannot be a DS dataset */ + if ((H5DSis_scale(did))==1) + return FAIL; - /* get info for the dataset in the parameter list */ - if(H5Oget_info(did, &oi1) < 0) - return FAIL; + /* get info for the dataset in the parameter list */ + if(H5Oget_info(did, &oi1) < 0) + return FAIL; - /* get info for the scale in the parameter list */ - if(H5Oget_info(dsid, &oi2) < 0) - return FAIL; + /* get info for the scale in the parameter list */ + if(H5Oget_info(dsid, &oi2) < 0) + return FAIL; - /* same object, not valid */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) - return FAIL; + /* same object, not valid */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) + return FAIL; - /* get ID type */ - if ((it1 = H5Iget_type(did)) < 0) - return FAIL; - if ((it2 = H5Iget_type(dsid)) < 0) - return FAIL; + /* get ID type */ + if ((it1 = H5Iget_type(did)) < 0) + return FAIL; + if ((it2 = H5Iget_type(dsid)) < 0) + return FAIL; - if (H5I_DATASET!=it1 || H5I_DATASET!=it2) - return FAIL; + if (H5I_DATASET!=it1 || H5I_DATASET!=it2) + return FAIL; -/*------------------------------------------------------------------------- - * Find "DIMENSION_LIST" - *------------------------------------------------------------------------- - */ - /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ - if ((has_dimlist = H5LT_find_attribute(did,DIMENSION_LIST)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * Find "DIMENSION_LIST" + *------------------------------------------------------------------------- + */ + /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ + if ((has_dimlist = H5LT_find_attribute(did,DIMENSION_LIST)) < 0) + return FAIL; - if (has_dimlist == 0) - return FAIL; + if (has_dimlist == 0) + return FAIL; - /* get dataset space */ - if ((sid = H5Dget_space(did)) < 0) - return FAIL; + /* get dataset space */ + if ((sid = H5Dget_space(did)) < 0) + return FAIL; - /* get rank */ - if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; + /* get rank */ + if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; - /* close dataset space */ - if (H5Sclose(sid) < 0) - return FAIL; + /* close dataset space */ + if (H5Sclose(sid) < 0) + return FAIL; -/*------------------------------------------------------------------------- - * find "REFERENCE_LIST" - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * find "REFERENCE_LIST" + *------------------------------------------------------------------------- + */ - /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ - if((has_reflist = H5LT_find_attribute(dsid, REFERENCE_LIST)) < 0) - return FAIL; + /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ + if((has_reflist = H5LT_find_attribute(dsid, REFERENCE_LIST)) < 0) + return FAIL; - if(has_reflist == 0) - return FAIL; + if(has_reflist == 0) + return FAIL; -/*------------------------------------------------------------------------- - * open "DIMENSION_LIST", and delete the reference - *------------------------------------------------------------------------- - */ - - if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) - return FAIL; - - if((tid = H5Aget_type(aid)) < 0) - goto out; - - if((sid = H5Aget_space(aid)) < 0) - goto out; - - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - - if(buf == NULL) - goto out; - - /* read */ - if (H5Aread(aid,tid,buf) < 0) - goto out; - - /* reset */ - if ( buf[idx].len > 0 ) - { - for (j=0; j<buf[idx].len; j++) - { - /* get the reference */ - ref = ((hobj_ref_t *)buf[idx].p)[j]; - - /* get the DS id */ - if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) - goto out; - - /* get info for DS in the parameter list */ - if(H5Oget_info(dsid, &oi1) < 0) - goto out; - - /* get info for this DS */ - if(H5Oget_info(dsid_j, &oi2) < 0) - goto out; - - /* same object, reset */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) { - for(jj=j; jj<buf[idx].len-1; jj++) - ((hobj_ref_t *)buf[idx].p)[jj] = ((hobj_ref_t *)buf[idx].p)[jj+1]; - buf[idx].len--; - - found_ds = 1; - - /* close the dereferenced dataset and break */ - if (H5Dclose(dsid_j) < 0) - goto out; - break; - } - - /* close the dereferenced dataset */ - if (H5Dclose(dsid_j) < 0) - goto out; - } /* j */ - } /* if */ - - /* write the attribute */ - if (H5Awrite(aid,tid,buf) < 0) - goto out; - - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if(H5Aclose(aid) < 0) - goto out; - if(buf) - free(buf); - - /* the scale must be present */ - if(found_ds == 0) - goto out; + /*------------------------------------------------------------------------- + * open "DIMENSION_LIST", and delete the reference + *------------------------------------------------------------------------- + */ -/*------------------------------------------------------------------------- - * the "REFERENCE_LIST" array exists, update - *------------------------------------------------------------------------- - */ + if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) + return FAIL; - if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) - goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; + if((sid = H5Aget_space(aid)) < 0) + goto out; - /* get and save the old reference(s) */ - if((sid = H5Aget_space(aid)) < 0) - goto out; + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) - goto out; + if(buf == NULL) + goto out; - dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); - if(dsbuf == NULL) - goto out; + /* read */ + if (H5Aread(aid,tid,buf) < 0) + goto out; - if (H5Aread(aid,tid,dsbuf) < 0) - goto out; + /* reset */ + if ( buf[idx].len > 0 ) + { + for (j=0; j<buf[idx].len; j++) + { + /* get the reference */ + ref = ((hobj_ref_t *)buf[idx].p)[j]; + + /* get the DS id */ + if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) + goto out; + + /* get info for DS in the parameter list */ + if(H5Oget_info(dsid, &oi1) < 0) + goto out; + + /* get info for this DS */ + if(H5Oget_info(dsid_j, &oi2) < 0) + goto out; + + /* same object, reset */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) { + for(jj=j; jj<buf[idx].len-1; jj++) + ((hobj_ref_t *)buf[idx].p)[jj] = ((hobj_ref_t *)buf[idx].p)[jj+1]; + buf[idx].len--; + + found_ds = 1; + + /* close the dereferenced dataset and break */ + if (H5Dclose(dsid_j) < 0) + goto out; + break; + } + + /* close the dereferenced dataset */ + if (H5Dclose(dsid_j) < 0) + goto out; + } /* j */ + } /* if */ + + /* write the attribute */ + if (H5Awrite(aid,tid,buf) < 0) + goto out; - for(i=0; i<nelmts; i++) - { - /* get the reference to the dataset */ - ref = dsbuf[i].ref; + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; - /* get the dataset id */ - if ((did_i = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) - goto out; + if (buf) + { + free(buf); + buf = NULL; + } - /* get info for dataset in the parameter list */ - if(H5Oget_info(did, &oi3) < 0) - goto out; + /* the scale must be present */ + if(found_ds == 0) + goto out; - /* get info for this dataset */ - if(H5Oget_info(did_i, &oi4) < 0) - goto out; + /*------------------------------------------------------------------------- + * the "REFERENCE_LIST" array exists, update + *------------------------------------------------------------------------- + */ - /* same object, reset. we want to detach only for this DIM */ - if(oi3.fileno == oi4.fileno && oi3.addr == oi4.addr && (int)idx == dsbuf[i].dim_idx) { - for(jj=i; jj<nelmts-1; jj++) - dsbuf[jj] = dsbuf[jj+1]; - nelmts--; - found_dset=1; + if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) + goto out; - /* close the dereferenced dataset and break */ - if (H5Dclose(did_i) < 0) - goto out; - break; - } /* if */ + if((tid = H5Aget_type(aid)) < 0) + goto out; - /* close the dereferenced dataset */ - if (H5Dclose(did_i) < 0) - goto out; + /* get and save the old reference(s) */ + if((sid = H5Aget_space(aid)) < 0) + goto out; - } /* i */ + if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) + goto out; - /* close space and attribute */ - if (H5Sclose(sid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; + dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); + if(dsbuf == NULL) + goto out; -/*------------------------------------------------------------------------- - * check if we found the pointed dataset - *------------------------------------------------------------------------- - */ + if (H5Aread(aid,tid,dsbuf) < 0) + goto out; - /* the pointed dataset must exist */ - if (found_dset == 0) - goto out; + for(i=0; i<nelmts; i++) + { + /* get the reference to the dataset */ + ref = dsbuf[i].ref; -/*------------------------------------------------------------------------- - * create a new attribute - *------------------------------------------------------------------------- - */ - - /* the attribute must be deleted, in order to the new one can reflect the changes*/ - if(H5Adelete(dsid, REFERENCE_LIST) < 0) - goto out; - - /* don't do anything for an empty array */ - if(nelmts) - { - /* create a new data space for the new references array */ - dims = (hsize_t*)malloc((size_t)nelmts * sizeof (hsize_t)); - if(dims == NULL) - goto out; - dims[0] = nelmts; - - dsbufn = malloc((size_t)nelmts * sizeof(ds_list_t)); - if(dsbufn == NULL) - goto out; - - /* store the new information */ - for(i = 0; i < nelmts; i++) - dsbufn[i] = dsbuf[i]; - - if((sid = H5Screate_simple(1, dims, NULL)) < 0) - goto out; - - /* create the attribute again with the changes of space */ - if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - - /* write the new attribute with the new references */ - if(H5Awrite(aid, tid, dsbufn) < 0) - goto out; - - /* close space and attribute */ - if(H5Sclose(sid) < 0) - goto out; - if(H5Aclose(aid) < 0) - goto out; - } /* nelmts */ - - /* close type */ - if (H5Tclose(tid) < 0) - goto out; - - if (dsbuf) { - free(dsbuf); - dsbuf=NULL; - } - if (dsbufn) { - free(dsbufn); - dsbufn=NULL; - } - if (dims) { - free(dims); - dims=NULL; - } - - - return SUCCEED; - -/* error zone, gracefully close */ + /* get the dataset id */ + if ((did_i = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) + goto out; + + /* get info for dataset in the parameter list */ + if(H5Oget_info(did, &oi3) < 0) + goto out; + + /* get info for this dataset */ + if(H5Oget_info(did_i, &oi4) < 0) + goto out; + + /* same object, reset. we want to detach only for this DIM */ + if(oi3.fileno == oi4.fileno && oi3.addr == oi4.addr && (int)idx == dsbuf[i].dim_idx) { + for(jj=i; jj<nelmts-1; jj++) + dsbuf[jj] = dsbuf[jj+1]; + nelmts--; + found_dset=1; + + /* close the dereferenced dataset and break */ + if (H5Dclose(did_i) < 0) + goto out; + break; + } /* if */ + + /* close the dereferenced dataset */ + if (H5Dclose(did_i) < 0) + goto out; + + } /* i */ + + /* close space and attribute */ + if (H5Sclose(sid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + + /*------------------------------------------------------------------------- + * check if we found the pointed dataset + *------------------------------------------------------------------------- + */ + + /* the pointed dataset must exist */ + if (found_dset == 0) + goto out; + + /*------------------------------------------------------------------------- + * create a new attribute + *------------------------------------------------------------------------- + */ + + /* the attribute must be deleted, in order to the new one can reflect the changes*/ + if (H5Adelete(dsid, REFERENCE_LIST) < 0) + goto out; + + /* don't do anything for an empty array */ + if(nelmts) + { + /* create a new data space for the new references array */ + dims = (hsize_t*)malloc((size_t)nelmts * sizeof (hsize_t)); + if(dims == NULL) + goto out; + dims[0] = nelmts; + + dsbufn = malloc((size_t)nelmts * sizeof(ds_list_t)); + if(dsbufn == NULL) + goto out; + + /* store the new information */ + for(i = 0; i < nelmts; i++) + dsbufn[i] = dsbuf[i]; + + if((sid = H5Screate_simple(1, dims, NULL)) < 0) + goto out; + + /* create the attribute again with the changes of space */ + if((aid = H5Acreate2(dsid, REFERENCE_LIST, tid, sid, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + + /* write the new attribute with the new references */ + if(H5Awrite(aid, tid, dsbufn) < 0) + goto out; + + /* close space and attribute */ + if(H5Sclose(sid) < 0) + goto out; + if(H5Aclose(aid) < 0) + goto out; + } /* nelmts */ + + /* close type */ + if (H5Tclose(tid) < 0) + goto out; + + if (dsbuf) + { + free(dsbuf); + dsbuf=NULL; + } + if (dsbufn) + { + free(dsbufn); + dsbufn=NULL; + } + if (dims) + { + free(dims); + dims=NULL; + } + + + return SUCCEED; + + /* error zone */ out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - if (dsbuf) - free(dsbuf); - if (dsbufn) - free(dsbufn); - if (dims) - free(dims); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + + if (dsbuf) + { + free(dsbuf); + dsbuf=NULL; + } + if (dsbufn) + { + free(dsbufn); + dsbufn=NULL; + } + if (dims) + { + free(dims); + dims=NULL; + } + if (buf) + { + free(buf); + buf=NULL; + } + } H5E_END_TRY; + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSis_attached - * - * Purpose: Report if dimension scale DSID is currently attached to - * dimension IDX of dataset DID by checking if DID has a pointer in the REFERENCE_LIST - * attribute and DSID (scale ) has a pointer in the DIMENSION_LIST attribute - * - * Return: - * 1: both the DS and the dataset pointers match - * 0: one of them or both do not match - * FAIL (-1): error - * - * Fails if: Bad arguments - * If DSID is not a Dimension Scale - * If DID is a Dimension Scale (A Dimension Scale cannot have scales) - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: February 18, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSis_attached +* +* Purpose: Report if dimension scale DSID is currently attached to +* dimension IDX of dataset DID by checking if DID has a pointer in the REFERENCE_LIST +* attribute and DSID (scale ) has a pointer in the DIMENSION_LIST attribute +* +* Return: +* 1: both the DS and the dataset pointers match +* 0: one of them or both do not match +* FAIL (-1): error +* +* Fails if: Bad arguments +* If DSID is not a Dimension Scale +* If DID is a Dimension Scale (A Dimension Scale cannot have scales) +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: February 18, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ htri_t H5DSis_attached(hid_t did, hid_t dsid, unsigned int idx) { - int has_dimlist; - int has_reflist; - hssize_t nelmts; - hid_t sid; /* space ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t aid = -1; /* attribute ID */ - int rank; /* rank of dataset */ - ds_list_t *dsbuf; /* array of attribute data in the DS pointing to the dataset */ - hobj_ref_t ref; /* reference to the DS */ - hvl_t *buf; /* VL buffer to store in the attribute */ - hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ + int has_dimlist; + int has_reflist; + hssize_t nelmts; + hid_t sid; /* space ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t aid = -1; /* attribute ID */ + int rank; /* rank of dataset */ + ds_list_t *dsbuf; /* array of attribute data in the DS pointing to the dataset */ + hobj_ref_t ref; /* reference to the DS */ + hvl_t *buf; /* VL buffer to store in the attribute */ + hid_t dsid_j; /* DS dataset ID in DIMENSION_LIST */ hid_t did_i; /* dataset ID in REFERENCE_LIST */ - H5O_info_t oi1, oi2, oi3, oi4; - H5I_type_t it1, it2; - int i; - int found_dset=0, found_ds=0; + H5O_info_t oi1, oi2, oi3, oi4; + H5I_type_t it1, it2; + int i; + int found_dset=0, found_ds=0; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* the dataset cannot be a DS dataset */ - if ((H5DSis_scale(did))==1) - return FAIL; - - /* get info for the dataset in the parameter list */ - if(H5Oget_info(did, &oi1) < 0) - return FAIL; - - /* get info for the scale in the parameter list */ - if(H5Oget_info(dsid, &oi2) < 0) - return FAIL; - - /* same object, not valid */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) - return FAIL; - - /* get ID type */ - if ((it1 = H5Iget_type(did)) < 0) - return FAIL; - if ((it2 = H5Iget_type(dsid)) < 0) - return FAIL; - - if (H5I_DATASET!=it1 || H5I_DATASET!=it2) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* the dataset cannot be a DS dataset */ + if ((H5DSis_scale(did))==1) + return FAIL; -/*------------------------------------------------------------------------- - * get space - *------------------------------------------------------------------------- - */ + /* get info for the dataset in the parameter list */ + if(H5Oget_info(did, &oi1) < 0) + return FAIL; - /* get dataset space */ - if ((sid = H5Dget_space(did)) < 0) - return FAIL; + /* get info for the scale in the parameter list */ + if(H5Oget_info(dsid, &oi2) < 0) + return FAIL; - /* get rank */ - if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; + /* same object, not valid */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) + return FAIL; - /* close dataset space */ - if (H5Sclose(sid) < 0) - goto out; + /* get ID type */ + if ((it1 = H5Iget_type(did)) < 0) + return FAIL; + if ((it2 = H5Iget_type(dsid)) < 0) + return FAIL; - /* parameter range checking */ - if(idx > ((unsigned)rank - 1)) - goto out; + if (H5I_DATASET!=it1 || H5I_DATASET!=it2) + return FAIL; - /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ - if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * get space + *------------------------------------------------------------------------- + */ -/*------------------------------------------------------------------------- - * open "DIMENSION_LIST" - *------------------------------------------------------------------------- - */ + /* get dataset space */ + if ((sid = H5Dget_space(did)) < 0) + return FAIL; + + /* get rank */ + if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; + + /* close dataset space */ + if (H5Sclose(sid) < 0) + goto out; - if(has_dimlist == 1) - { - if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) - goto out; + /* parameter range checking */ + if(idx > ((unsigned)rank - 1)) + goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; + /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ + if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) + return FAIL; - if((sid = H5Aget_space(aid)) < 0) - goto out; + /*------------------------------------------------------------------------- + * open "DIMENSION_LIST" + *------------------------------------------------------------------------- + */ - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); + if(has_dimlist == 1) + { + if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) + goto out; - if(buf == NULL) - goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; - /* read */ - if (H5Aread(aid,tid,buf) < 0) - goto out; + if((sid = H5Aget_space(aid)) < 0) + goto out; - /* iterate all the REFs in this dimension IDX */ - for (i=0; i<(int)buf[idx].len; i++) - { - /* get the reference */ - ref = ((hobj_ref_t *)buf[idx].p)[i]; + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - /* get the scale id for this REF */ - if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) - goto out; + if(buf == NULL) + goto out; - /* get info for DS in the parameter list */ - if(H5Oget_info(dsid, &oi1) < 0) - goto out; + /* read */ + if (H5Aread(aid,tid,buf) < 0) + goto out; - /* get info for this DS */ - if(H5Oget_info(dsid_j, &oi2) < 0) - goto out; + /* iterate all the REFs in this dimension IDX */ + for (i=0; i<(int)buf[idx].len; i++) + { + /* get the reference */ + ref = ((hobj_ref_t *)buf[idx].p)[i]; - /* same object */ - if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) - found_ds = 1; + /* get the scale id for this REF */ + if ((dsid_j = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) + goto out; - /* close the dereferenced dataset */ - if (H5Dclose(dsid_j) < 0) - goto out; + /* get info for DS in the parameter list */ + if(H5Oget_info(dsid, &oi1) < 0) + goto out; - } + /* get info for this DS */ + if(H5Oget_info(dsid_j, &oi2) < 0) + goto out; + /* same object */ + if(oi1.fileno == oi2.fileno && oi1.addr == oi2.addr) + found_ds = 1; - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (buf) - free(buf); - } /* has_dimlist */ + /* close the dereferenced dataset */ + if (H5Dclose(dsid_j) < 0) + goto out; -/*------------------------------------------------------------------------- - * info on the >>DS<< dataset - *------------------------------------------------------------------------- - */ + } - /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ - if((has_reflist = H5LT_find_attribute(dsid, REFERENCE_LIST)) < 0) - goto out; -/*------------------------------------------------------------------------- - * open "REFERENCE_LIST" - *------------------------------------------------------------------------- - */ + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (buf) + { + free(buf); + buf = NULL; + } + } /* has_dimlist */ - if(has_reflist == 1) - { - if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) - goto out; + /*------------------------------------------------------------------------- + * info on the >>DS<< dataset + *------------------------------------------------------------------------- + */ - if((tid = H5Aget_type(aid)) < 0) - goto out; + /* try to find the attribute "REFERENCE_LIST" on the >>DS<< dataset */ + if((has_reflist = H5LT_find_attribute(dsid, REFERENCE_LIST)) < 0) + goto out; - /* get and save the old reference(s) */ - if((sid = H5Aget_space(aid)) < 0) - goto out; + /*------------------------------------------------------------------------- + * open "REFERENCE_LIST" + *------------------------------------------------------------------------- + */ - if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) - goto out; + if(has_reflist == 1) + { + if((aid = H5Aopen(dsid, REFERENCE_LIST, H5P_DEFAULT)) < 0) + goto out; - dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); + if((tid = H5Aget_type(aid)) < 0) + goto out; - if (dsbuf == NULL) - goto out; + /* get and save the old reference(s) */ + if((sid = H5Aget_space(aid)) < 0) + goto out; - if (H5Aread(aid,tid,dsbuf) < 0) - goto out; + if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) + goto out; -/*------------------------------------------------------------------------- - * iterate - *------------------------------------------------------------------------- - */ - - for(i=0; i<nelmts; i++) - { - /* get the reference */ - ref = dsbuf[i].ref; - - /* the reference was not deleted */ - if (ref) - { - /* get the dataset id */ - if ((did_i = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) - goto out; - - /* get info for dataset in the parameter list */ - if(H5Oget_info(did, &oi3) < 0) - goto out; - - /* get info for this dataset */ - if(H5Oget_info(did_i, &oi4) < 0) - goto out; - - /* same object */ - if(oi3.fileno == oi4.fileno && oi3.addr == oi4.addr && (int)idx==dsbuf[i].dim_idx) - found_dset=1; - - /* close the dereferenced dataset */ - if (H5Dclose(did_i) < 0) - goto out; - } /* if */ - } /* i */ - - - /* close */ - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (dsbuf) - free(dsbuf); - } /* has_reflist */ - - if (found_ds && found_dset) - return 1; - else - return 0; - -/* error zone, gracefully close */ + dsbuf = malloc((size_t)nelmts * sizeof(ds_list_t)); + + if (dsbuf == NULL) + goto out; + + if (H5Aread(aid,tid,dsbuf) < 0) + goto out; + + /*------------------------------------------------------------------------- + * iterate + *------------------------------------------------------------------------- + */ + + for(i=0; i<nelmts; i++) + { + /* get the reference */ + ref = dsbuf[i].ref; + + /* the reference was not deleted */ + if (ref) + { + /* get the dataset id */ + if ((did_i = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) + goto out; + + /* get info for dataset in the parameter list */ + if(H5Oget_info(did, &oi3) < 0) + goto out; + + /* get info for this dataset */ + if(H5Oget_info(did_i, &oi4) < 0) + goto out; + + /* same object */ + if(oi3.fileno == oi4.fileno && oi3.addr == oi4.addr && (int)idx==dsbuf[i].dim_idx) + found_dset=1; + + /* close the dereferenced dataset */ + if (H5Dclose(did_i) < 0) + goto out; + } /* if */ + } /* i */ + + + /* close */ + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (dsbuf) + { + free(dsbuf); + dsbuf = NULL; + } + } /* has_reflist */ + + if (found_ds && found_dset) + return 1; + else + return 0; + + /* error zone */ out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + } H5E_END_TRY; + + if (buf) + { + free(buf); + buf = NULL; + } + if (dsbuf) + { + free(dsbuf); + dsbuf = NULL; + } + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSiterate_scales - * - * Purpose: H5DSiterate_scales iterates over the scales attached to dimension DIM - * of dataset DID. For each scale in the list, the visitor_data and some - * additional information, specified below, are passed to the visitor function. - * The iteration begins with the IDX object in the group and the next element - * to be processed by the operator is returned in IDX. If IDX is NULL, then the - * iterator starts at zero. - * - * Parameters: - * - * hid_t DID; IN: the dataset - * unsigned int dim; IN: the dimension of the dataset - * int *idx; IN/OUT: input the index to start iterating, output the - * next index to visit. If NULL, start at the first position. - * H5DS_iterate_t visitor; IN: the visitor function - * void *visitor_data; IN: arbitrary data to pass to the visitor function. - * - * Iterate over all scales of DIM, calling an application callback - * with the item, key and any operator data. - * - * The operator callback receives a pointer to the item , - * and the pointer to the operator data passed - * in to H5SL_iterate ('op_data'). The return values from an operator are: - * A. Zero causes the iterator to continue, returning zero when all - * nodes of that type have been processed. - * B. Positive causes the iterator to immediately return that positive - * value, indicating short-circuit success. - * C. Negative causes the iterator to immediately return that value, - * indicating failure. - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 31, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSiterate_scales +* +* Purpose: H5DSiterate_scales iterates over the scales attached to dimension DIM +* of dataset DID. For each scale in the list, the visitor_data and some +* additional information, specified below, are passed to the visitor function. +* The iteration begins with the IDX object in the group and the next element +* to be processed by the operator is returned in IDX. If IDX is NULL, then the +* iterator starts at zero. +* +* Parameters: +* +* hid_t DID; IN: the dataset +* unsigned int dim; IN: the dimension of the dataset +* int *idx; IN/OUT: input the index to start iterating, output the +* next index to visit. If NULL, start at the first position. +* H5DS_iterate_t visitor; IN: the visitor function +* void *visitor_data; IN: arbitrary data to pass to the visitor function. +* +* Iterate over all scales of DIM, calling an application callback +* with the item, key and any operator data. +* +* The operator callback receives a pointer to the item , +* and the pointer to the operator data passed +* in to H5SL_iterate ('op_data'). The return values from an operator are: +* A. Zero causes the iterator to continue, returning zero when all +* nodes of that type have been processed. +* B. Positive causes the iterator to immediately return that positive +* value, indicating short-circuit success. +* C. Negative causes the iterator to immediately return that value, +* indicating failure. +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 31, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DSiterate_scales(hid_t did, unsigned int dim, @@ -1228,162 +1302,165 @@ herr_t H5DSiterate_scales(hid_t did, H5DS_iterate_t visitor, void *visitor_data ) { - hid_t scale_id; - int rank; - hobj_ref_t ref; /* reference to the DS */ - hid_t sid; /* space ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t aid = -1; /* attribute ID */ - hvl_t *buf=NULL; /* VL buffer to store in the attribute */ - H5I_type_t it; /* ID type */ - herr_t ret_value=0; - int j_idx; - int nscales; - int has_dimlist; - int i; + hid_t scale_id; + int rank; + hobj_ref_t ref; /* reference to the DS */ + hid_t sid; /* space ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t aid = -1; /* attribute ID */ + hvl_t *buf=NULL; /* VL buffer to store in the attribute */ + H5I_type_t it; /* ID type */ + herr_t ret_value=0; + int j_idx; + int nscales; + int has_dimlist; + int i; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* get ID type */ - if ((it = H5Iget_type(did)) < 0) - return FAIL; - - if (H5I_DATASET!=it) - return FAIL; - - /* get the number of scales assotiated with this DIM */ - if ((nscales = H5DSget_num_scales(did,dim)) < 0) - return FAIL; - - /* get dataset space */ - if ((sid = H5Dget_space(did)) < 0) - return FAIL; - - /* get rank */ - if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; - - /* close dataset space */ - if(H5Sclose(sid) < 0) - goto out; - - /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ - if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) - return FAIL; - - if(has_dimlist == 0) - return SUCCEED; - - else if(has_dimlist == 1) - { - if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) - goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; - if((sid = H5Aget_space(aid)) < 0) - goto out; - - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - - if(buf == NULL) - goto out; - - /* read */ - if(H5Aread(aid, tid, buf) < 0) - goto out; - - if ( buf[dim].len > 0 ) - { - if (idx!=NULL) - j_idx = *idx; - else - j_idx=0; - - /* iterate */ - for(i=j_idx; i<nscales; i++) - { - /* get the reference */ - ref = ((hobj_ref_t *)buf[dim].p)[ i ]; - - /* disable error reporting, the ID might refer to a deleted dataset */ - H5E_BEGIN_TRY { - /* get the DS id */ - if ((scale_id = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) - goto out; - } H5E_END_TRY; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* get ID type */ + if ((it = H5Iget_type(did)) < 0) + return FAIL; + + if (H5I_DATASET!=it) + return FAIL; + + /* get the number of scales assotiated with this DIM */ + if ((nscales = H5DSget_num_scales(did,dim)) < 0) + return FAIL; + + /* get dataset space */ + if ((sid = H5Dget_space(did)) < 0) + return FAIL; + + /* get rank */ + if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; - if((ret_value=(visitor)(did,dim,scale_id,visitor_data))!=0) + /* close dataset space */ + if(H5Sclose(sid) < 0) + goto out; + + /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ + if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) + return FAIL; + + if(has_dimlist == 0) + return SUCCEED; + + else if(has_dimlist == 1) { - /* set the return IDX OUT value at current scale index and break */ - if (idx!=NULL) - { - *idx = i; - } + if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) + goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; + if((sid = H5Aget_space(aid)) < 0) + goto out; - /* close the DS id */ - if (H5Dclose(scale_id) < 0) - goto out; + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - break; - } + if(buf == NULL) + goto out; - /* close the DS id */ - if (H5Dclose(scale_id) < 0) - goto out; + /* read */ + if(H5Aread(aid, tid, buf) < 0) + goto out; - } /* i */ - } /* if */ + if ( buf[dim].len > 0 ) + { + if (idx!=NULL) + j_idx = *idx; + else + j_idx=0; + + /* iterate */ + for(i=j_idx; i<nscales; i++) + { + /* get the reference */ + ref = ((hobj_ref_t *)buf[dim].p)[ i ]; + + /* disable error reporting, the ID might refer to a deleted dataset */ + H5E_BEGIN_TRY { + /* get the DS id */ + if ((scale_id = H5Rdereference(did,H5R_OBJECT,&ref)) < 0) + goto out; + } H5E_END_TRY; + + if((ret_value=(visitor)(did,dim,scale_id,visitor_data))!=0) + { + /* set the return IDX OUT value at current scale index and break */ + if (idx!=NULL) + { + *idx = i; + } + + /* close the DS id */ + if (H5Dclose(scale_id) < 0) + goto out; + + break; + } + + /* close the DS id */ + if (H5Dclose(scale_id) < 0) + goto out; + + } /* i */ + } /* if */ - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (buf) - free(buf); + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (buf) + { + free(buf); + buf = NULL; + } - } /* if has_dimlist */ + } /* if has_dimlist */ - return ret_value; + return ret_value; out: - H5E_BEGIN_TRY { - if (buf) - { - H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf); - free(buf); - } - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + if (buf) + { + H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf); + free(buf); + } + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + } H5E_END_TRY; + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSset_label - * - * Purpose: Set label for the dimension IDX of dataset DID to the value LABEL - * - * Return: Success: SUCCEED, Failure: FAIL - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 11, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSset_label +* +* Purpose: Set label for the dimension IDX of dataset DID to the value LABEL +* +* Return: Success: SUCCEED, Failure: FAIL +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 11, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { int has_labels; @@ -1397,9 +1474,9 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { unsigned int i; /*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ + * parameter checking + *------------------------------------------------------------------------- + */ /* get ID type */ if ((it = H5Iget_type(did)) < 0) return FAIL; @@ -1408,9 +1485,9 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { return FAIL; /*------------------------------------------------------------------------- - * attribute "DIMENSION_LABELS" - *------------------------------------------------------------------------- - */ + * attribute "DIMENSION_LABELS" + *------------------------------------------------------------------------- + */ /* try to find the attribute "DIMENSION_LABELS" on the >>data<< dataset */ if ((has_labels = H5LT_find_attribute(did, DIMENSION_LABELS)) < 0) @@ -1429,9 +1506,9 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { goto out; /*------------------------------------------------------------------------- - * make the attribute and insert label - *------------------------------------------------------------------------- - */ + * make the attribute and insert label + *------------------------------------------------------------------------- + */ if (has_labels == 0) { dims[0] = rank; @@ -1448,7 +1525,7 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { /* create the attribute */ if ((aid = H5Acreate2(did, DIMENSION_LABELS, tid, sid, - H5P_DEFAULT, H5P_DEFAULT)) < 0) + H5P_DEFAULT, H5P_DEFAULT)) < 0) goto out; /* allocate and initialize */ @@ -1475,13 +1552,16 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { if (H5Aclose(aid) < 0) goto out; if (buf) + { free((void *) buf); + buf = NULL; + } } /*------------------------------------------------------------------------- - * just insert label - *------------------------------------------------------------------------- - */ + * just insert label + *------------------------------------------------------------------------- + */ else { if ((aid = H5Aopen(did, DIMENSION_LABELS, H5P_DEFAULT)) < 0) @@ -1526,13 +1606,16 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { if (H5Aclose(aid) < 0) goto out; if (buf) + { free((void *) buf); + buf = NULL; + } } return SUCCEED; - /* error zone, gracefully close */ - out: + /* error zone */ +out: if (buf) { if (buf[idx]) /* check if we errored during H5Awrite */ buf[idx] = NULL; /* don't free label */ @@ -1544,35 +1627,35 @@ herr_t H5DSset_label(hid_t did, unsigned int idx, const char *label) { free(buf); } H5E_BEGIN_TRY - { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - }H5E_END_TRY; + { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + }H5E_END_TRY; return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSget_label - * - * Purpose: Read the label LABEL for dimension IDX of dataset DID - * Up to 'size' characters are stored in 'label' followed by a '\0' string - * terminator. If the label is longer than 'size'-1, - * the string terminator is stored in the last position of the buffer to - * properly terminate the string. - * - * Return: 0 if no label found, size of label if found, Failure: FAIL - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 11, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSget_label +* +* Purpose: Read the label LABEL for dimension IDX of dataset DID +* Up to 'size' characters are stored in 'label' followed by a '\0' string +* terminator. If the label is longer than 'size'-1, +* the string terminator is stored in the last position of the buffer to +* properly terminate the string. +* +* Return: 0 if no label found, size of label if found, Failure: FAIL +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 11, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { int has_labels; hid_t sid; /* space ID */ @@ -1586,9 +1669,9 @@ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { int i; /*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ + * parameter checking + *------------------------------------------------------------------------- + */ /* get ID type */ if ((it = H5Iget_type(did)) < 0) return FAIL; @@ -1597,20 +1680,19 @@ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { return FAIL; /*------------------------------------------------------------------------- - * attribute "DIMENSION_LABELS" - *------------------------------------------------------------------------- - */ + * attribute "DIMENSION_LABELS" + *------------------------------------------------------------------------- + */ /* try to find the attribute "DIMENSION_LABELS" on the >>data<< dataset */ if ((has_labels = H5LT_find_attribute(did, DIMENSION_LABELS)) < 0) return FAIL; /* return 0 and NULL for label if no label found */ - if (has_labels == 0) { + if (has_labels == 0) + { if (label) -/* label = NULL; - */ - label[0] = 0; + label[0] = 0; return 0; } @@ -1627,9 +1709,9 @@ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { goto out; /*------------------------------------------------------------------------- - * open the attribute and read label - *------------------------------------------------------------------------- - */ + * open the attribute and read label + *------------------------------------------------------------------------- + */ assert (has_labels == 1); if ((aid = H5Aopen(did, DIMENSION_LABELS, H5P_DEFAULT)) < 0) @@ -1674,12 +1756,15 @@ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { if (H5Aclose(aid) < 0) goto out; if (buf) + { free(buf); + buf = NULL; + } return (ssize_t) nbytes; - /* error zone, gracefully close */ - out: + /* error zone */ +out: if (buf) { /* free all the ptr's from the H5Aread() */ for (i = 0; i < (unsigned int) rank; i++) { @@ -1689,411 +1774,423 @@ ssize_t H5DSget_label(hid_t did, unsigned int idx, char *label, size_t size) { free(buf); } H5E_BEGIN_TRY - { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - }H5E_END_TRY; + { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + }H5E_END_TRY; return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSget_scale_name - * - * Purpose: Read the name of dataset scale DID into buffer NAME - * Up to 'size' characters are stored in 'name' followed by a '\0' string - * terminator. If the name is longer than 'size'-1, - * the string terminator is stored in the last position of the buffer to - * properly terminate the string. - * - * Return: size of name if found, zero if not found, Failure: FAIL - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 04, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSget_scale_name +* +* Purpose: Read the name of dataset scale DID into buffer NAME +* Up to 'size' characters are stored in 'name' followed by a '\0' string +* terminator. If the name is longer than 'size'-1, +* the string terminator is stored in the last position of the buffer to +* properly terminate the string. +* +* Return: size of name if found, zero if not found, Failure: FAIL +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 04, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ ssize_t H5DSget_scale_name(hid_t did, char *name, size_t size) { - hid_t aid; /* attribute ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t sid; /* space ID */ - H5I_type_t it; /* ID type */ - size_t nbytes; - size_t copy_len; - int has_name; - char *buf=NULL; + hid_t aid; /* attribute ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t sid; /* space ID */ + H5I_type_t it; /* ID type */ + size_t nbytes; + size_t copy_len; + int has_name; + char *buf=NULL; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* get ID type */ - if ((it = H5Iget_type(did)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* get ID type */ + if ((it = H5Iget_type(did)) < 0) + return FAIL; - if (H5I_DATASET!=it) - return FAIL; + if (H5I_DATASET!=it) + return FAIL; - if ((H5DSis_scale(did))<=0) - return FAIL; + if ((H5DSis_scale(did))<=0) + return FAIL; -/*------------------------------------------------------------------------- - * check if the DS has a name - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * check if the DS has a name + *------------------------------------------------------------------------- + */ + + /* try to find the attribute "NAME" on the >>DS<< dataset */ + if ((has_name = H5LT_find_attribute(did, "NAME")) < 0) + return FAIL; - /* try to find the attribute "NAME" on the >>DS<< dataset */ - if((has_name = H5LT_find_attribute(did, "NAME")) < 0) - return FAIL; + if (has_name == 0) + return 0; - if(has_name == 0) - return 0; + /*------------------------------------------------------------------------- + * open the attribute + *------------------------------------------------------------------------- + */ -/*------------------------------------------------------------------------- - * open the attribute - *------------------------------------------------------------------------- - */ - - if((aid = H5Aopen(did, "NAME", H5P_DEFAULT)) < 0) - return FAIL; - - /* get space */ - if((sid = H5Aget_space(aid)) < 0) - goto out; - - /* get type */ - if((tid = H5Aget_type(aid)) < 0) - goto out; - - /* get the size */ - if((nbytes = H5Tget_size(tid)) == 0) - goto out; - - /* allocate a temporary buffer */ - buf = (char*)malloc(nbytes * sizeof(char)); - if (buf == NULL) - goto out; - - /* read */ - if (H5Aread(aid,tid,buf) < 0) - goto out; - - /* compute the string length which will fit into the user's buffer */ - copy_len = MIN(size-1, nbytes); - - /* copy all/some of the name */ - if (name) { - memcpy(name, buf, copy_len); - - /* terminate the string */ - name[copy_len]='\0'; - } - - /* close */ - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (buf) - { - free(buf); - buf=NULL; - } - - return (ssize_t) nbytes; - - /* error zone, gracefully close */ + if((aid = H5Aopen(did, "NAME", H5P_DEFAULT)) < 0) + return FAIL; + + /* get space */ + if((sid = H5Aget_space(aid)) < 0) + goto out; + + /* get type */ + if((tid = H5Aget_type(aid)) < 0) + goto out; + + /* get the size */ + if((nbytes = H5Tget_size(tid)) == 0) + goto out; + + /* allocate a temporary buffer */ + buf = (char*)malloc(nbytes * sizeof(char)); + if (buf == NULL) + goto out; + + /* read */ + if (H5Aread(aid,tid,buf) < 0) + goto out; + + /* compute the string length which will fit into the user's buffer */ + copy_len = MIN(size-1, nbytes); + + /* copy all/some of the name */ + if (name) { + memcpy(name, buf, copy_len); + + /* terminate the string */ + name[copy_len]='\0'; + } + + /* close */ + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (buf) + { + free(buf); + buf=NULL; + } + + return (ssize_t) nbytes; + + /* error zone */ out: - H5E_BEGIN_TRY { - H5Aclose(aid); - H5Tclose(tid); - H5Sclose(sid); - } H5E_END_TRY; - if (buf) - free(buf); - return FAIL; + H5E_BEGIN_TRY { + H5Aclose(aid); + H5Tclose(tid); + H5Sclose(sid); + } H5E_END_TRY; + if (buf) + { + free(buf); + buf=NULL; + } + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSis_scale - * - * Purpose: check if the dataset DID is a dimension scale - * - * Return: 1, is, 0, not, FAIL, error - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 04, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSis_scale +* +* Purpose: check if the dataset DID is a dimension scale +* +* Return: 1, is, 0, not, FAIL, error +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 04, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ htri_t H5DSis_scale(hid_t did) { - hid_t tid = -1; /* attribute type ID */ - hid_t aid; /* attribute ID */ - herr_t has_class; /* has the "CLASS" attribute */ - htri_t is_ds; /* boolean return value */ - H5I_type_t it; /* ID type */ - char buf[20]; + hid_t tid = -1; /* attribute type ID */ + hid_t aid; /* attribute ID */ + herr_t has_class; /* has the "CLASS" attribute */ + htri_t is_ds; /* boolean return value */ + H5I_type_t it; /* ID type */ + char buf[20]; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* get ID type */ - if ((it = H5Iget_type(did)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* get ID type */ + if ((it = H5Iget_type(did)) < 0) + return FAIL; - if(H5I_DATASET != it) - return FAIL; + if(H5I_DATASET != it) + return FAIL; - /* try to find the attribute "CLASS" on the dataset */ - if((has_class = H5LT_find_attribute(did, "CLASS")) < 0) - return FAIL; + /* try to find the attribute "CLASS" on the dataset */ + if((has_class = H5LT_find_attribute(did, "CLASS")) < 0) + return FAIL; - if(has_class == 0) - is_ds = 0; + if(has_class == 0) + is_ds = 0; - else - { - if((aid = H5Aopen(did, "CLASS", H5P_DEFAULT)) < 0) - goto out; + else + { + if((aid = H5Aopen(did, "CLASS", H5P_DEFAULT)) < 0) + goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; - if(H5Aread(aid, tid, buf) < 0) - goto out; + if(H5Aread(aid, tid, buf) < 0) + goto out; - if(strcmp(buf, DIMENSION_SCALE_CLASS)==0) - is_ds = 1; - else - is_ds = 0; + if(strcmp(buf, DIMENSION_SCALE_CLASS)==0) + is_ds = 1; + else + is_ds = 0; - if(H5Tclose(tid) < 0) - goto out; + if(H5Tclose(tid) < 0) + goto out; - if (H5Aclose(aid) < 0) - goto out; - } + if (H5Aclose(aid) < 0) + goto out; + } - return is_ds; + return is_ds; -/* error zone, gracefully close */ + /* error zone */ out: - H5E_BEGIN_TRY { - H5Aclose(aid); - H5Tclose(tid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Aclose(aid); + H5Tclose(tid); + } H5E_END_TRY; + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DSget_num_scales - * - * Purpose: get the number of scales linked to the IDX dimension of dataset DID - * - * Return: - * Success: number of scales - * Failure: FAIL - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: January 13, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DSget_num_scales +* +* Purpose: get the number of scales linked to the IDX dimension of dataset DID +* +* Return: +* Success: number of scales +* Failure: FAIL +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: January 13, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ int H5DSget_num_scales(hid_t did, unsigned int idx) { - int has_dimlist; - hid_t sid; /* space ID */ - hid_t tid = -1; /* attribute type ID */ - hid_t aid = -1; /* attribute ID */ - int rank; /* rank of dataset */ - hvl_t *buf; /* VL buffer to store in the attribute */ - H5I_type_t it; /* ID type */ - int nscales; + int has_dimlist; + hid_t sid; /* space ID */ + hid_t tid = -1; /* attribute type ID */ + hid_t aid = -1; /* attribute ID */ + int rank; /* rank of dataset */ + hvl_t *buf; /* VL buffer to store in the attribute */ + H5I_type_t it; /* ID type */ + int nscales; -/*------------------------------------------------------------------------- - * parameter checking - *------------------------------------------------------------------------- - */ - /* get ID type */ - if ((it = H5Iget_type(did)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * parameter checking + *------------------------------------------------------------------------- + */ + /* get ID type */ + if ((it = H5Iget_type(did)) < 0) + return FAIL; - if (H5I_DATASET!=it) - return FAIL; + if (H5I_DATASET!=it) + return FAIL; -/*------------------------------------------------------------------------- - * the attribute "DIMENSION_LIST" on the >>data<< dataset must exist - *------------------------------------------------------------------------- - */ - /* get dataset space */ - if ((sid = H5Dget_space(did)) < 0) - return FAIL; + /*------------------------------------------------------------------------- + * the attribute "DIMENSION_LIST" on the >>data<< dataset must exist + *------------------------------------------------------------------------- + */ + /* get dataset space */ + if ((sid = H5Dget_space(did)) < 0) + return FAIL; - /* get rank */ - if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; + /* get rank */ + if ((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; - /* close dataset space */ - if (H5Sclose(sid) < 0) - goto out; + /* close dataset space */ + if (H5Sclose(sid) < 0) + goto out; - /* dimemsion index IDX range checking */ - if (idx>=(unsigned int )rank) - return FAIL; + /* dimemsion index IDX range checking */ + if (idx>=(unsigned int )rank) + return FAIL; - /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ - if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) - return FAIL; + /* try to find the attribute "DIMENSION_LIST" on the >>data<< dataset */ + if((has_dimlist = H5LT_find_attribute(did, DIMENSION_LIST)) < 0) + return FAIL; - /* it does not exist */ - if(has_dimlist == 0) - return 0; + /* it does not exist */ + if(has_dimlist == 0) + return 0; -/*------------------------------------------------------------------------- - * the attribute exists, open it - *------------------------------------------------------------------------- - */ - - else - { - if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) - goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; - if((sid = H5Aget_space(aid)) < 0) - goto out; - - /* allocate and initialize the VL */ - buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); - - if(buf == NULL) - goto out; - - /* read */ - if(H5Aread(aid, tid, buf) < 0) - goto out; - - nscales=(int)buf[idx].len; - - /* close */ - if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) - goto out; - if (H5Sclose(sid) < 0) - goto out; - if (H5Tclose(tid) < 0) - goto out; - if (H5Aclose(aid) < 0) - goto out; - if (buf) - free(buf); - - } /* has_dimlist */ - - return nscales; - -/* error zone, gracefully close */ + /*------------------------------------------------------------------------- + * the attribute exists, open it + *------------------------------------------------------------------------- + */ + + else + { + if((aid = H5Aopen(did, DIMENSION_LIST, H5P_DEFAULT)) < 0) + goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; + if((sid = H5Aget_space(aid)) < 0) + goto out; + + /* allocate and initialize the VL */ + buf = (hvl_t*)malloc((size_t)rank * sizeof(hvl_t)); + + if(buf == NULL) + goto out; + + /* read */ + if(H5Aread(aid, tid, buf) < 0) + goto out; + + nscales=(int)buf[idx].len; + + /* close */ + if (H5Dvlen_reclaim(tid,sid,H5P_DEFAULT,buf) < 0) + goto out; + if (H5Sclose(sid) < 0) + goto out; + if (H5Tclose(tid) < 0) + goto out; + if (H5Aclose(aid) < 0) + goto out; + if (buf) + { + free(buf); + buf = NULL; + } + + } /* has_dimlist */ + + return nscales; + + /* error zone */ out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Aclose(aid); - H5Tclose(tid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Aclose(aid); + H5Tclose(tid); + } H5E_END_TRY; + + if (buf) + { + free(buf); + buf = NULL; + } + return FAIL; } /*------------------------------------------------------------------------- - * Function: H5DS_is_reserved - * - * Purpose: Verify that a dataset's CLASS is either an image, palette or table - * - * Return: true, false, fail - * - * Programmer: pvn@ncsa.uiuc.edu - * - * Date: March 19, 2005 - * - * Comments: - * - * Modifications: - * - *------------------------------------------------------------------------- - */ +* Function: H5DS_is_reserved +* +* Purpose: Verify that a dataset's CLASS is either an image, palette or table +* +* Return: true, false, fail +* +* Programmer: pvn@ncsa.uiuc.edu +* +* Date: March 19, 2005 +* +* Comments: +* +* Modifications: +* +*------------------------------------------------------------------------- +*/ herr_t H5DS_is_reserved(hid_t did) { - int has_class; - hid_t tid = -1; - hid_t aid = -1; - char buf[40]; - herr_t ret; + int has_class; + hid_t tid = -1; + hid_t aid = -1; + char buf[40]; + herr_t ret; - /* try to find the attribute "CLASS" on the dataset */ - if((has_class = H5LT_find_attribute(did, "CLASS")) < 0) - return -1; + /* try to find the attribute "CLASS" on the dataset */ + if((has_class = H5LT_find_attribute(did, "CLASS")) < 0) + return -1; - if(has_class == 0) - return 0; + if(has_class == 0) + return 0; - assert(has_class == 1); - if((aid = H5Aopen(did, "CLASS", H5P_DEFAULT)) < 0) - goto out; + assert(has_class == 1); + if((aid = H5Aopen(did, "CLASS", H5P_DEFAULT)) < 0) + goto out; - if((tid = H5Aget_type(aid)) < 0) - goto out; + if((tid = H5Aget_type(aid)) < 0) + goto out; - if(H5Aread(aid, tid, buf) < 0) - goto out; + if(H5Aread(aid, tid, buf) < 0) + goto out; - if(strcmp(buf, IMAGE_CLASS) == 0 || - strcmp(buf, PALETTE_CLASS) == 0 || - strcmp(buf, TABLE_CLASS) == 0 ) - ret = 1; - else - ret = 0; + if(strcmp(buf, IMAGE_CLASS) == 0 || + strcmp(buf, PALETTE_CLASS) == 0 || + strcmp(buf, TABLE_CLASS) == 0 ) + ret = 1; + else + ret = 0; - if (H5Tclose(tid) < 0) - goto out; + if (H5Tclose(tid) < 0) + goto out; - if (H5Aclose(aid) < 0) - goto out; + if (H5Aclose(aid) < 0) + goto out; - return ret; + return ret; -/* error zone, gracefully close */ + /* error zone */ out: - H5E_BEGIN_TRY { - H5Tclose(tid); - H5Aclose(aid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Tclose(tid); + H5Aclose(aid); + } H5E_END_TRY; + return FAIL; } diff --git a/hl/test/test_ds.c b/hl/test/test_ds.c index 8fb8604..e4d4260 100644 --- a/hl/test/test_ds.c +++ b/hl/test/test_ds.c @@ -76,22 +76,22 @@ static int read_data( const char* fname, int ndims, hsize_t *dims, float **buf ) */ int main(void) { - int nerrors=0; + int nerrors=0; - nerrors += test_simple() < 0 ?1:0; - nerrors += test_errors() < 0 ?1:0; - nerrors += test_rank() < 0 ?1:0; - nerrors += test_iterators() < 0 ?1:0; - nerrors += test_types() < 0 ?1:0; - nerrors += test_data() < 0 ?1:0; + nerrors += test_simple() < 0 ?1:0; + nerrors += test_errors() < 0 ?1:0; + nerrors += test_rank() < 0 ?1:0; + nerrors += test_iterators() < 0 ?1:0; + nerrors += test_types() < 0 ?1:0; + nerrors += test_data() < 0 ?1:0; - if(nerrors) goto error; - printf("All dimension scales tests passed.\n"); - return 0; + if(nerrors) goto error; + printf("All dimension scales tests passed.\n"); + return 0; error: - printf("***** %d DIMENSION SCALES TEST%s FAILED! *****\n",nerrors, 1 == nerrors ? "" : "S"); - return 1; + printf("***** %d DIMENSION SCALES TEST%s FAILED! *****\n",nerrors, 1 == nerrors ? "" : "S"); + return 1; } @@ -116,1332 +116,1332 @@ error: static int test_simple(void) { - hid_t fid; /* file ID */ - hid_t did = -1; /* dataset ID */ - hid_t dsid; /* DS dataset ID */ - hid_t sid; /* space ID */ - hid_t gid; /* group ID */ - int rank = RANK; /* rank of data dataset */ - int rankds = 1; /* rank of DS dataset */ - hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ - int buf[DIM_DATA] = {1,2,3,4,5,6,7,8,9,10,11,12}; /* data of data dataset */ - hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ - hsize_t s2_dim[1] = {DIM2_SIZE}; /* size of DS 2 dataset */ - char sname[30]; /* scale name buffer */ - char dname[30]; /* dataset name */ - int s1_wbuf[DIM1_SIZE] = {10,20,30}; /* data of DS 1 dataset */ - int s11_wbuf[DIM1_SIZE] = {10,100,300}; /* data of DS 1 dataset */ - int s2_wbuf[DIM2_SIZE] = {100,200,300,400}; /* data of DS 2 dataset */ - int s21_wbuf[DIM2_SIZE] = {10,20,30,40}; /* data of DS 2 dataset */ - int s22_wbuf[DIM2_SIZE] = {5,10,50,300}; /* data of DS 2 dataset */ - char dim0_label[16]; /* read label for DIM 0 */ - char dim1_label[16]; /* read label for DIM 1 */ - char *dim0_labeld; /* read label for DIM 0 */ - char *dim1_labeld; /* read label for DIM 1 */ - char dim0_labels[3]; /* read label for DIM 0 */ - char dim1_labels[3]; /* read label for DIM 1 */ - ssize_t dim0_label_size; /* lenght of label buffer */ - ssize_t dim1_label_size; /* lenght of label buffer */ - unsigned int dim; /* dataset dimension index */ - int scale_idx; /* scale index */ - int nscales; /* number of scales in DIM */ - ssize_t name_len; /* lenght of name buffer */ - char *name_out=NULL; /* scale name buffer */ - char snames[3]; /* scale name buffer */ - int i, j; - - - printf("Testing API functions\n"); + hid_t fid; /* file ID */ + hid_t did = -1; /* dataset ID */ + hid_t dsid; /* DS dataset ID */ + hid_t sid; /* space ID */ + hid_t gid; /* group ID */ + int rank = RANK; /* rank of data dataset */ + int rankds = 1; /* rank of DS dataset */ + hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ + int buf[DIM_DATA] = {1,2,3,4,5,6,7,8,9,10,11,12}; /* data of data dataset */ + hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ + hsize_t s2_dim[1] = {DIM2_SIZE}; /* size of DS 2 dataset */ + char sname[30]; /* scale name buffer */ + char dname[30]; /* dataset name */ + int s1_wbuf[DIM1_SIZE] = {10,20,30}; /* data of DS 1 dataset */ + int s11_wbuf[DIM1_SIZE] = {10,100,300}; /* data of DS 1 dataset */ + int s2_wbuf[DIM2_SIZE] = {100,200,300,400}; /* data of DS 2 dataset */ + int s21_wbuf[DIM2_SIZE] = {10,20,30,40}; /* data of DS 2 dataset */ + int s22_wbuf[DIM2_SIZE] = {5,10,50,300}; /* data of DS 2 dataset */ + char dim0_label[16]; /* read label for DIM 0 */ + char dim1_label[16]; /* read label for DIM 1 */ + char *dim0_labeld; /* read label for DIM 0 */ + char *dim1_labeld; /* read label for DIM 1 */ + char dim0_labels[3]; /* read label for DIM 0 */ + char dim1_labels[3]; /* read label for DIM 1 */ + ssize_t dim0_label_size; /* lenght of label buffer */ + ssize_t dim1_label_size; /* lenght of label buffer */ + unsigned int dim; /* dataset dimension index */ + int scale_idx; /* scale index */ + int nscales; /* number of scales in DIM */ + ssize_t name_len; /* lenght of name buffer */ + char *name_out=NULL; /* scale name buffer */ + char snames[3]; /* scale name buffer */ + int i, j; + + + printf("Testing API functions\n"); + + /*------------------------------------------------------------------------- + * create a file for the test + *------------------------------------------------------------------------- + */ + + /* create a file using default properties */ + if((fid=H5Fcreate(FILE1,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * create datasets: 1 "data" dataset and 4 dimension scales + *------------------------------------------------------------------------- + */ + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,buf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_int(fid,DS_1_NAME,rankds,s1_dim,s1_wbuf) < 0) + goto out; + + /* make a DS dataset with an alternate scale for the 2nd dimension */ + if(H5LTmake_dataset_int(fid,DS_11_NAME,rankds,s1_dim,s11_wbuf) < 0) + goto out; + + /* make a DS dataset for the second dimension */ + if(H5LTmake_dataset_int(fid,DS_2_NAME,rankds,s2_dim,s2_wbuf) < 0) + goto out; + + /* make a DS dataset with an alternate scale for the 2nd dimension */ + if(H5LTmake_dataset_int(fid,DS_21_NAME,rankds,s2_dim,s21_wbuf) < 0) + goto out; + + /* make a DS dataset with an alternate scale for the 2nd dimension */ + if(H5LTmake_dataset_int(fid,DS_22_NAME,rankds,s2_dim,s22_wbuf) < 0) + goto out; + + + /*------------------------------------------------------------------------- + * test 1: attach scale + *------------------------------------------------------------------------- + */ + + TESTING2("attach scales"); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach the DS_1_NAME dimension scale to "dset_a" + *------------------------------------------------------------------------- + */ + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) + goto out; + + /* attach the DS_1_NAME dimension scale to "dset_a" at dimension 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach the DS_11_NAME dimension scale to "dset_a" + *------------------------------------------------------------------------- + */ + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_11_NAME, H5P_DEFAULT)) < 0) + goto out; + + /* attach the DS_11_NAME dimension scale to "dset_a" at dimension 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach the DS_2_NAME dimension scale to "dset_a" + *------------------------------------------------------------------------- + */ -/*------------------------------------------------------------------------- - * create a file for the test - *------------------------------------------------------------------------- - */ - - /* create a file using default properties */ - if((fid=H5Fcreate(FILE1,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) - goto out; - -/*------------------------------------------------------------------------- - * create datasets: 1 "data" dataset and 4 dimension scales - *------------------------------------------------------------------------- - */ - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,buf) < 0) - goto out; - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_int(fid,DS_1_NAME,rankds,s1_dim,s1_wbuf) < 0) - goto out; - - /* make a DS dataset with an alternate scale for the 2nd dimension */ - if(H5LTmake_dataset_int(fid,DS_11_NAME,rankds,s1_dim,s11_wbuf) < 0) - goto out; - - /* make a DS dataset for the second dimension */ - if(H5LTmake_dataset_int(fid,DS_2_NAME,rankds,s2_dim,s2_wbuf) < 0) - goto out; - - /* make a DS dataset with an alternate scale for the 2nd dimension */ - if(H5LTmake_dataset_int(fid,DS_21_NAME,rankds,s2_dim,s21_wbuf) < 0) - goto out; - - /* make a DS dataset with an alternate scale for the 2nd dimension */ - if(H5LTmake_dataset_int(fid,DS_22_NAME,rankds,s2_dim,s22_wbuf) < 0) - goto out; - - -/*------------------------------------------------------------------------- - * test 1: attach scale - *------------------------------------------------------------------------- - */ - - TESTING2("attach scales"); - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the DS_1_NAME dimension scale to "dset_a" - *------------------------------------------------------------------------- - */ - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) - goto out; - - /* attach the DS_1_NAME dimension scale to "dset_a" at dimension 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the DS_11_NAME dimension scale to "dset_a" - *------------------------------------------------------------------------- - */ - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_11_NAME, H5P_DEFAULT)) < 0) - goto out; - - /* attach the DS_11_NAME dimension scale to "dset_a" at dimension 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the DS_2_NAME dimension scale to "dset_a" - *------------------------------------------------------------------------- - */ - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) - goto out; - - /* attach the "ds2" dimension scale to "dset_a" as the 2nd dimension */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the DS_21_NAME dimension scale to "dset_a" - *------------------------------------------------------------------------- - */ - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) - goto out; - - /* attach the DS_21_NAME dimension scale to "dset_a" as the 2nd dimension */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the DS_22_NAME dimension scale to "dset_a" - *------------------------------------------------------------------------- - */ - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) - goto out; - - /* attach the "ds22" dimension scale to "dset_a" as the 2nd dimension */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - - -/*------------------------------------------------------------------------- - * verify attachment - *------------------------------------------------------------------------- - */ - - if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM0)<=0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM1)<=0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM1)<=0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM1)<=0) - goto out; - if(H5Dclose(dsid)) - goto out; - - - /* close dataset ID of "dset_a" */ - if(H5Dclose(did) < 0) - goto out; - - - PASSED(); - - -/*------------------------------------------------------------------------- - * test 2: get number of scales - *------------------------------------------------------------------------- - */ - - TESTING2("get number of scales"); - -/*------------------------------------------------------------------------- - * verify that "dset_a" has dimension scales - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* verify that "dset_a" has 1 dimension scale at DIM 0 */ - if((nscales = H5DSget_num_scales(did,0)) < 0) - goto out; - if(nscales!=2) - goto out; - - /* verify that "dset_a" has 3 dimension scales at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=3) - goto out; - - /* close dataset ID of "dset_a" */ - if(H5Dclose(did) < 0) - goto out; - - -/*------------------------------------------------------------------------- - * create datasets: 1 "data" dataset and 1 dimension scale - *------------------------------------------------------------------------- - */ - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_b",rank,dims,buf) < 0) - goto out; - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_int(fid,"ds_b_1",rankds,s1_dim,s1_wbuf) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach the scale to "dset_b" - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) - goto out; - if((dsid = H5Dopen2(fid,"ds_b_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * verify if "dset_b" has dimension scales - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_b" */ - if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) - goto out; - - /* verify that "dset_b" has 1 dimension scale at DIM 0 */ - if((nscales = H5DSget_num_scales(did,0)) < 0) - goto out; - if(nscales!=1) - goto out; - - /* verify that "dset_b" has 0 dimension scales at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=0) - goto out; - - /* close dataset ID of "dset_b" */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * test 3: detach scales - *------------------------------------------------------------------------- - */ - - TESTING2("detach scales "); - - -/*------------------------------------------------------------------------- - * create datasets: one "data" dataset and 4 dimension scales - *------------------------------------------------------------------------- - */ - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_c",rank,dims,buf) < 0) - goto out; - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_int(fid,"ds_c_1",rankds,s1_dim,s1_wbuf) < 0) - goto out; - - /* make a DS dataset for the second dimension */ - if(H5LTmake_dataset_int(fid,"ds_c_2",rankds,s2_dim,s2_wbuf) < 0) - goto out; - - /* make a DS dataset with an alternate scale for the 2nd dimension */ - if(H5LTmake_dataset_int(fid,"ds_c_21",rankds,s2_dim,s2_wbuf) < 0) - goto out; - - /* make a DS dataset with an alternate scale for the 2nd dimension */ - if(H5LTmake_dataset_int(fid,"ds_c_22",rankds,s2_dim,s2_wbuf) < 0) - goto out; - - -/*------------------------------------------------------------------------- - * attach the scales to "dset_c" - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_c_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_c_2", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,1) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_c_21", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,1) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_c_22", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,1) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * verify if "dset_c" has dimension scales - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - /* verify that "dset_c" has 1 dimension scale at DIM 0 */ - if((nscales = H5DSget_num_scales(did,0)) < 0) - goto out; - if(nscales!=1) - goto out; - /* verify that "dset_c" has 3 dimension scales at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=3) - goto out; - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * detach the "ds_c_21" dimension scale to "dset_c" - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_c" */ - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_c_21", H5P_DEFAULT)) < 0) - goto out; - - /* detach the "ds_c_21" dimension scale to "dset_c" in DIM 1 */ - if(H5DSdetach_scale(did,dsid,1) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_c" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * "dset_c" must have now 2 dimension scales at DIM 1 - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - /* verify that "dset_c" has 2 dimension scales at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=2) - goto out; - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * detach the "ds_c_22" dimension scale to "dset_c" - *------------------------------------------------------------------------- - */ + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) + goto out; - /* get the dataset id for "dset_c" */ - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; + /* attach the "ds2" dimension scale to "dset_a" as the 2nd dimension */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_c_22", H5P_DEFAULT)) < 0) - goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; - /* detach the "ds_c_22" dimension scale to "dset_c" in DIM 1 */ - if(H5DSdetach_scale(did,dsid,1) < 0) - goto out; + /*------------------------------------------------------------------------- + * attach the DS_21_NAME dimension scale to "dset_a" + *------------------------------------------------------------------------- + */ + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) + goto out; + + /* attach the DS_21_NAME dimension scale to "dset_a" as the 2nd dimension */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; - /* close dataset ID of "dset_c" */ - if(H5Dclose(did) < 0) - goto out; + /*------------------------------------------------------------------------- + * attach the DS_22_NAME dimension scale to "dset_a" + *------------------------------------------------------------------------- + */ -/*------------------------------------------------------------------------- - * "dset_c" must have now 1 dimension scale at DIM 1 - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - /* verify that "dset_c" has 1 dimension scale at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=1) - goto out; - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * detach the "ds_c_2" dimension scale to "dset_c" - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_c" */ - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_c_2", H5P_DEFAULT)) < 0) - goto out; - - /* detach the "ds_c_2" dimension scale to "dset_c" in DIM 1 */ - if(H5DSdetach_scale(did,dsid,1) < 0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_c" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * "dset_c" must have now 0 dimension scales at DIM 1 - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) - goto out; - /* verify that "dset_c" has 1 dimension scale at DIM 1 */ - if((nscales = H5DSget_num_scales(did,1)) < 0) - goto out; - if(nscales!=0) - goto out; - if(H5Dclose(did) < 0) - goto out; - - -/*------------------------------------------------------------------------- - * create 3 datasets: 1 "data" dataset and 2 dimension scales - *------------------------------------------------------------------------- - */ - if(H5LTmake_dataset_int(fid,"dset_d",rank,dims,NULL) < 0) - goto out; - if(H5LTmake_dataset_int(fid,"ds_d_1",rankds,s1_dim,NULL) < 0) - goto out; - if(H5LTmake_dataset_int(fid,"ds_d_2",rankds,s2_dim,NULL) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach them - *------------------------------------------------------------------------- - */ - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,1) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * verify - *------------------------------------------------------------------------- - */ - - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM0)<=0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM1)<=0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if(H5Dclose(did) < 0) - goto out; - - -/*------------------------------------------------------------------------- - * detach - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_d" */ - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) - goto out; - - /* detach the dimension scale to "dset_d" in DIM 0 */ - if(H5DSdetach_scale(did,dsid,DIM0) < 0) - goto out; - - /* verify attach, it must return 0 for no attach */ - if(H5DSis_attached(did,dsid,DIM0)!=0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_d" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach again - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_d" */ - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) - goto out; - - /* attach "ds_d_1" again in DIM 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - - /* verify attach, it must return 1 for attach */ - if(H5DSis_attached(did,dsid,DIM0)!=1) - goto out; - - /* verify that "ds_d_1" has only 1 scale at DIM0 */ - if((nscales = H5DSget_num_scales(did,DIM0)) < 0) - goto out; - if(nscales!=1) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_d" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * detach/detach - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_d" */ - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) - goto out; - - /* detach the "ds_d_2" dimension scale to "dset_d" in DIM 1 */ - if(H5DSdetach_scale(did,dsid,DIM1) < 0) - goto out; - - /* detach again, it should fail */ - if(H5DSdetach_scale(did,dsid,DIM1)==SUCCEED) - goto out; - - /* verify attach, it must return 0 for no attach */ - if(H5DSis_attached(did,dsid,DIM1)!=0) - goto out; - - /* verify that "ds_d_1" has no scale at DIM1 */ - if((nscales = H5DSget_num_scales(did,DIM1)) < 0) - goto out; - if(nscales!=0) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_d" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach twice - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "dset_d" */ - if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) - goto out; + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) + goto out; - /* attach "ds_d_2" in DIM 1 */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; + /* attach the "ds22" dimension scale to "dset_a" as the 2nd dimension */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; - /* verify attach, it must return 1 for attach */ - if(H5DSis_attached(did,dsid,DIM1)!=1) - goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + + + /*------------------------------------------------------------------------- + * verify attachment + *------------------------------------------------------------------------- + */ + + if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM0)<=0) + goto out; + if(H5Dclose(dsid)) + goto out; + + if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM1)<=0) + goto out; + if(H5Dclose(dsid)) + goto out; + + if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM1)<=0) + goto out; + if(H5Dclose(dsid)) + goto out; + + if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM1)<=0) + goto out; + if(H5Dclose(dsid)) + goto out; + + + /* close dataset ID of "dset_a" */ + if(H5Dclose(did) < 0) + goto out; + + + PASSED(); + + + /*------------------------------------------------------------------------- + * test 2: get number of scales + *------------------------------------------------------------------------- + */ + + TESTING2("get number of scales"); + + /*------------------------------------------------------------------------- + * verify that "dset_a" has dimension scales + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* verify that "dset_a" has 1 dimension scale at DIM 0 */ + if((nscales = H5DSget_num_scales(did,0)) < 0) + goto out; + if(nscales!=2) + goto out; + + /* verify that "dset_a" has 3 dimension scales at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=3) + goto out; + + /* close dataset ID of "dset_a" */ + if(H5Dclose(did) < 0) + goto out; - /* verify that "ds_d_2" has only 1 scale at DIM1 */ - if((nscales = H5DSget_num_scales(did,DIM0)) < 0) - goto out; - if(nscales!=1) - goto out; - /* attach "ds_d_2" again in DIM 1 */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; + /*------------------------------------------------------------------------- + * create datasets: 1 "data" dataset and 1 dimension scale + *------------------------------------------------------------------------- + */ + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_b",rank,dims,buf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_int(fid,"ds_b_1",rankds,s1_dim,s1_wbuf) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach the scale to "dset_b" + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) + goto out; + if((dsid = H5Dopen2(fid,"ds_b_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * verify if "dset_b" has dimension scales + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_b" */ + if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) + goto out; + + /* verify that "dset_b" has 1 dimension scale at DIM 0 */ + if((nscales = H5DSget_num_scales(did,0)) < 0) + goto out; + if(nscales!=1) + goto out; + + /* verify that "dset_b" has 0 dimension scales at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=0) + goto out; + + /* close dataset ID of "dset_b" */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + + /*------------------------------------------------------------------------- + * test 3: detach scales + *------------------------------------------------------------------------- + */ + + TESTING2("detach scales "); + + + /*------------------------------------------------------------------------- + * create datasets: one "data" dataset and 4 dimension scales + *------------------------------------------------------------------------- + */ + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_c",rank,dims,buf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_int(fid,"ds_c_1",rankds,s1_dim,s1_wbuf) < 0) + goto out; + + /* make a DS dataset for the second dimension */ + if(H5LTmake_dataset_int(fid,"ds_c_2",rankds,s2_dim,s2_wbuf) < 0) + goto out; + + /* make a DS dataset with an alternate scale for the 2nd dimension */ + if(H5LTmake_dataset_int(fid,"ds_c_21",rankds,s2_dim,s2_wbuf) < 0) + goto out; + + /* make a DS dataset with an alternate scale for the 2nd dimension */ + if(H5LTmake_dataset_int(fid,"ds_c_22",rankds,s2_dim,s2_wbuf) < 0) + goto out; + + + /*------------------------------------------------------------------------- + * attach the scales to "dset_c" + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_c_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_c_2", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,1) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_c_21", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,1) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_c_22", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,1) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * verify if "dset_c" has dimension scales + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + /* verify that "dset_c" has 1 dimension scale at DIM 0 */ + if((nscales = H5DSget_num_scales(did,0)) < 0) + goto out; + if(nscales!=1) + goto out; + /* verify that "dset_c" has 3 dimension scales at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=3) + goto out; + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * detach the "ds_c_21" dimension scale to "dset_c" + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_c" */ + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_c_21", H5P_DEFAULT)) < 0) + goto out; + + /* detach the "ds_c_21" dimension scale to "dset_c" in DIM 1 */ + if(H5DSdetach_scale(did,dsid,1) < 0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_c" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * "dset_c" must have now 2 dimension scales at DIM 1 + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + /* verify that "dset_c" has 2 dimension scales at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=2) + goto out; + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * detach the "ds_c_22" dimension scale to "dset_c" + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_c" */ + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_c_22", H5P_DEFAULT)) < 0) + goto out; + + /* detach the "ds_c_22" dimension scale to "dset_c" in DIM 1 */ + if(H5DSdetach_scale(did,dsid,1) < 0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_c" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * "dset_c" must have now 1 dimension scale at DIM 1 + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + /* verify that "dset_c" has 1 dimension scale at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=1) + goto out; + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * detach the "ds_c_2" dimension scale to "dset_c" + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_c" */ + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_c_2", H5P_DEFAULT)) < 0) + goto out; + + /* detach the "ds_c_2" dimension scale to "dset_c" in DIM 1 */ + if(H5DSdetach_scale(did,dsid,1) < 0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_c" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * "dset_c" must have now 0 dimension scales at DIM 1 + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_c", H5P_DEFAULT)) < 0) + goto out; + /* verify that "dset_c" has 1 dimension scale at DIM 1 */ + if((nscales = H5DSget_num_scales(did,1)) < 0) + goto out; + if(nscales!=0) + goto out; + if(H5Dclose(did) < 0) + goto out; + + + /*------------------------------------------------------------------------- + * create 3 datasets: 1 "data" dataset and 2 dimension scales + *------------------------------------------------------------------------- + */ + if(H5LTmake_dataset_int(fid,"dset_d",rank,dims,NULL) < 0) + goto out; + if(H5LTmake_dataset_int(fid,"ds_d_1",rankds,s1_dim,NULL) < 0) + goto out; + if(H5LTmake_dataset_int(fid,"ds_d_2",rankds,s2_dim,NULL) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach them + *------------------------------------------------------------------------- + */ + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,1) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * verify + *------------------------------------------------------------------------- + */ + + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM0)<=0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM1)<=0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if(H5Dclose(did) < 0) + goto out; + + + /*------------------------------------------------------------------------- + * detach + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_d" */ + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) + goto out; + + /* detach the dimension scale to "dset_d" in DIM 0 */ + if(H5DSdetach_scale(did,dsid,DIM0) < 0) + goto out; + + /* verify attach, it must return 0 for no attach */ + if(H5DSis_attached(did,dsid,DIM0)!=0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_d" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach again + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_d" */ + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_d_1", H5P_DEFAULT)) < 0) + goto out; + + /* attach "ds_d_1" again in DIM 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + + /* verify attach, it must return 1 for attach */ + if(H5DSis_attached(did,dsid,DIM0)!=1) + goto out; + + /* verify that "ds_d_1" has only 1 scale at DIM0 */ + if((nscales = H5DSget_num_scales(did,DIM0)) < 0) + goto out; + if(nscales!=1) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_d" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * detach/detach + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_d" */ + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) + goto out; + + /* detach the "ds_d_2" dimension scale to "dset_d" in DIM 1 */ + if(H5DSdetach_scale(did,dsid,DIM1) < 0) + goto out; + + /* detach again, it should fail */ + if(H5DSdetach_scale(did,dsid,DIM1)==SUCCEED) + goto out; + + /* verify attach, it must return 0 for no attach */ + if(H5DSis_attached(did,dsid,DIM1)!=0) + goto out; + + /* verify that "ds_d_1" has no scale at DIM1 */ + if((nscales = H5DSget_num_scales(did,DIM1)) < 0) + goto out; + if(nscales!=0) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_d" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach twice + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "dset_d" */ + if((did = H5Dopen2(fid,"dset_d", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_d_2", H5P_DEFAULT)) < 0) + goto out; + + /* attach "ds_d_2" in DIM 1 */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + + /* verify attach, it must return 1 for attach */ + if(H5DSis_attached(did,dsid,DIM1)!=1) + goto out; + + /* verify that "ds_d_2" has only 1 scale at DIM1 */ + if((nscales = H5DSget_num_scales(did,DIM0)) < 0) + goto out; + if(nscales!=1) + goto out; + + /* attach "ds_d_2" again in DIM 1 */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + + /* verify attach, it must return 1 for attach */ + if(H5DSis_attached(did,dsid,DIM1)!=1) + goto out; + + /* verify that "ds_d_2" has only 1 scale at DIM1 */ + if((nscales = H5DSget_num_scales(did,DIM0)) < 0) + goto out; + if(nscales != 1) + goto out; + + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + + /* close dataset ID of "dset_d" */ + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * create 10 datasets: 5 "data" dataset and 5 dimension scales + *------------------------------------------------------------------------- + */ + + /* create a group */ + if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + + /* create the data space for the dataset */ + if((sid = H5Screate_simple(rank,dims,NULL)) < 0) + goto out; + + for(i = 0; i < 5; i++) { + sprintf(dname,"dset_%d",i); + if((did = H5Dcreate2(gid, dname, H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + sprintf(sname,"ds_%d",i); + if((dsid = H5Dcreate2(gid, sname, H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,"scale") < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + } + + /*------------------------------------------------------------------------- + * attach for DIM 0 + *------------------------------------------------------------------------- + */ + + for(i = 0; i < 5; i++) { + sprintf(dname, "dset_%d", i); + if((did = H5Dopen2(gid, dname, H5P_DEFAULT)) < 0) + goto out; + for(j = 0; j < 5; j++) { + sprintf(sname, "ds_%d", j); + if((dsid = H5Dopen2(gid, sname, H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did, dsid, DIM0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + } + if(H5Dclose(did) < 0) + goto out; + } + + /*------------------------------------------------------------------------- + * dettach for DIM0 + *------------------------------------------------------------------------- + */ + + for(i = 0; i < 5; i++) { + sprintf(dname, "dset_%d", i); + if((did = H5Dopen2(gid, dname, H5P_DEFAULT)) < 0) + goto out; + for(j = 0; j < 5; j++) { + sprintf(sname, "ds_%d", j); + if((dsid = H5Dopen2(gid, sname, H5P_DEFAULT)) < 0) + goto out; + if(H5DSdetach_scale(did, dsid, DIM0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + } + if(H5Dclose(did) < 0) + goto out; + } + + + /*------------------------------------------------------------------------- + * attach again for DIM0 + *------------------------------------------------------------------------- + */ + + for(i=0; i<5; i++) + { + sprintf(dname,"dset_%d",i); + if((did = H5Dopen2(gid,dname, H5P_DEFAULT)) < 0) + goto out; + for(j=0; j<5; j++) + { + sprintf(sname,"ds_%d",j); + if((dsid = H5Dopen2(gid,sname, H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + } + if(H5Dclose(did) < 0) + goto out; + } + + /* close */ + if(H5Sclose(sid) < 0) + goto out; + if(H5Gclose(gid) < 0) + goto out; + + + PASSED(); + + + /*------------------------------------------------------------------------- + * create a dataset and attach only to 1 dimension + *------------------------------------------------------------------------- + */ + + TESTING2("attach only to 1 dimension"); + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_e",rank,dims,NULL) < 0) + goto out; + + /* make a scale */ + if(H5LTmake_dataset_int(fid,"ds_e_1",rankds,s1_dim,NULL) < 0) + goto out; + + /* attach the DS to dimension 1 */ + if((did = H5Dopen2(fid,"dset_e", H5P_DEFAULT)) < 0) + goto out; + if((dsid = H5Dopen2(fid,"ds_e_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + if(H5DSis_attached(did,dsid,DIM1)<=0) + goto out; + + + /* try to detach all dimensions. for dimensions 0 and 2, it is an error */ + for(i=0; i<rank; i++) + { + if( i==1 ) + { + if(H5DSdetach_scale(did,dsid,(unsigned)i) < 0) + goto out; + } + else + { + if(H5DSdetach_scale(did,dsid,(unsigned)i)!=FAIL) + goto out; + } + } + + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + + /*------------------------------------------------------------------------- + * test 4: set/get label + *------------------------------------------------------------------------- + */ + + TESTING2("set/get label"); + + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * set label + *------------------------------------------------------------------------- + */ + + if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) + goto out; + if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) + goto out; + + /*------------------------------------------------------------------------- + * get the scale name using a static buffer + *------------------------------------------------------------------------- + */ + + if(H5DSget_label(did,DIM0,dim0_label,sizeof(dim0_label)) < 0) + goto out; + if(H5DSget_label(did,DIM1,dim1_label,sizeof(dim1_label)) < 0) + goto out; + + if(strcmp(DIM0_LABEL,dim0_label)!=0) + goto out; + if(strcmp(DIM1_LABEL,dim1_label)!=0) + goto out; + + /*------------------------------------------------------------------------- + * get the scale name using a dynamic buffer + *------------------------------------------------------------------------- + */ + + if((dim0_label_size=H5DSget_label(did,DIM0,NULL,(size_t)0)) < 0) + goto out; + if((dim1_label_size=H5DSget_label(did,DIM1,NULL,(size_t)0)) < 0) + goto out; + + /* allocate */ + dim0_labeld = (char*)malloc(dim0_label_size * sizeof(char)); + dim1_labeld = (char*)malloc(dim1_label_size * sizeof(char)); + if( dim0_labeld==NULL || dim1_labeld==NULL) + goto out; + + if(H5DSget_label(did,DIM0,dim0_labeld,(size_t)dim0_label_size) < 0) + goto out; + if(H5DSget_label(did,DIM1,dim1_labeld,(size_t)dim1_label_size) < 0) + goto out; + + if(strncmp(DIM0_LABEL,dim0_labeld,(size_t)(dim0_label_size-1))!=0) + goto out; + if(strncmp(DIM1_LABEL,dim1_labeld,(size_t)(dim1_label_size-1))!=0) + goto out; + + if(dim0_labeld) + { + free(dim0_labeld); + dim0_labeld=NULL; + } + if(dim1_labeld) + { + free(dim1_labeld); + dim1_labeld=NULL; + } + + + /*------------------------------------------------------------------------- + * get the label using a static buffer smaller than the string lenght + *------------------------------------------------------------------------- + */ + + if(H5DSget_label(did,DIM0,dim0_labels,sizeof(dim0_labels)) < 0) + goto out; + if(H5DSget_label(did,DIM1,dim1_labels,sizeof(dim1_labels)) < 0) + goto out; + + if(strncmp(DIM0_LABEL,dim0_label,sizeof(dim0_labels)-1)!=0) + goto out; + if(strncmp(DIM1_LABEL,dim1_label,sizeof(dim1_labels)-1)!=0) + goto out; + + + + if(H5Dclose(did)) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * test 5: set scale/get scale name + *------------------------------------------------------------------------- + */ + TESTING2("set scale/get scale name"); + + if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) + goto out; + + if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) + goto out; + + /* verify that DS_1_NAME is a dimension scale dataset */ + if((H5DSis_scale(dsid))==0) + goto out; + + /*------------------------------------------------------------------------- + * get the scale name using a dynamic buffer + *------------------------------------------------------------------------- + */ + + /* get the lenght of the scale name (pass NULL in name) */ + if((name_len=H5DSget_scale_name(dsid,NULL,(size_t)0)) < 0) + goto out; + + /* allocate a buffer */ + name_out = (char*)malloc(name_len * sizeof(char)); + if(name_out == NULL) + goto out; + + /* get the scale name using this buffer */ + if(H5DSget_scale_name(dsid,name_out,(size_t)name_len) < 0) + goto out; + + if(strcmp(SCALE_1_NAME,name_out)!=0) + goto out; + + if(name_out) + { + free(name_out); + name_out=NULL; + } + + /*------------------------------------------------------------------------- + * get the scale name using a static buffer + *------------------------------------------------------------------------- + */ + + /* get the scale name using this buffer */ + if(H5DSget_scale_name(dsid,sname,sizeof(sname)) < 0) + goto out; + + if(strcmp(SCALE_1_NAME,sname)!=0) + goto out; + + /*------------------------------------------------------------------------- + * get the scale name using a static buffer smaller than the string lenght + *------------------------------------------------------------------------- + */ + + /* get the scale name using this buffer */ + if(H5DSget_scale_name(dsid,snames,sizeof(snames)) < 0) + goto out; + + if(strncmp(SCALE_1_NAME,snames,sizeof(snames)-1)!=0) + goto out; + + if(H5Dclose(dsid)) + goto out; + + /*------------------------------------------------------------------------- + * add scale names + *------------------------------------------------------------------------- + */ - /* verify attach, it must return 1 for attach */ - if(H5DSis_attached(did,dsid,DIM1)!=1) - goto out; - - /* verify that "ds_d_2" has only 1 scale at DIM1 */ - if((nscales = H5DSget_num_scales(did,DIM0)) < 0) - goto out; - if(nscales != 1) - goto out; - - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - - /* close dataset ID of "dset_d" */ - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * create 10 datasets: 5 "data" dataset and 5 dimension scales - *------------------------------------------------------------------------- - */ - - /* create a group */ - if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - - /* create the data space for the dataset */ - if((sid = H5Screate_simple(rank,dims,NULL)) < 0) - goto out; - - for(i = 0; i < 5; i++) { - sprintf(dname,"dset_%d",i); - if((did = H5Dcreate2(gid, dname, H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - sprintf(sname,"ds_%d",i); - if((dsid = H5Dcreate2(gid, sname, H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,"scale") < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - } - -/*------------------------------------------------------------------------- - * attach for DIM 0 - *------------------------------------------------------------------------- - */ - - for(i = 0; i < 5; i++) { - sprintf(dname, "dset_%d", i); - if((did = H5Dopen2(gid, dname, H5P_DEFAULT)) < 0) - goto out; - for(j = 0; j < 5; j++) { - sprintf(sname, "ds_%d", j); - if((dsid = H5Dopen2(gid, sname, H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did, dsid, DIM0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - } - if(H5Dclose(did) < 0) - goto out; - } - -/*------------------------------------------------------------------------- - * dettach for DIM0 - *------------------------------------------------------------------------- - */ - - for(i = 0; i < 5; i++) { - sprintf(dname, "dset_%d", i); - if((did = H5Dopen2(gid, dname, H5P_DEFAULT)) < 0) - goto out; - for(j = 0; j < 5; j++) { - sprintf(sname, "ds_%d", j); - if((dsid = H5Dopen2(gid, sname, H5P_DEFAULT)) < 0) - goto out; - if(H5DSdetach_scale(did, dsid, DIM0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - } - if(H5Dclose(did) < 0) - goto out; - } - - -/*------------------------------------------------------------------------- - * attach again for DIM0 - *------------------------------------------------------------------------- - */ - - for(i=0; i<5; i++) - { - sprintf(dname,"dset_%d",i); - if((did = H5Dopen2(gid,dname, H5P_DEFAULT)) < 0) - goto out; - for(j=0; j<5; j++) - { - sprintf(sname,"ds_%d",j); - if((dsid = H5Dopen2(gid,sname, H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - } - if(H5Dclose(did) < 0) - goto out; - } - - /* close */ - if(H5Sclose(sid) < 0) - goto out; - if(H5Gclose(gid) < 0) - goto out; - - - PASSED(); - - -/*------------------------------------------------------------------------- - * create a dataset and attach only to 1 dimension - *------------------------------------------------------------------------- - */ - - TESTING2("attach only to 1 dimension"); - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_e",rank,dims,NULL) < 0) - goto out; - - /* make a scale */ - if(H5LTmake_dataset_int(fid,"ds_e_1",rankds,s1_dim,NULL) < 0) - goto out; - - /* attach the DS to dimension 1 */ - if((did = H5Dopen2(fid,"dset_e", H5P_DEFAULT)) < 0) - goto out; - if((dsid = H5Dopen2(fid,"ds_e_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - if(H5DSis_attached(did,dsid,DIM1)<=0) - goto out; - - - /* try to detach all dimensions. for dimensions 0 and 2, it is an error */ - for(i=0; i<rank; i++) - { - if( i==1 ) - { - if(H5DSdetach_scale(did,dsid,(unsigned)i) < 0) - goto out; - } - else - { - if(H5DSdetach_scale(did,dsid,(unsigned)i)!=FAIL) - goto out; - } - } - - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * test 4: set/get label - *------------------------------------------------------------------------- - */ - - TESTING2("set/get label"); - - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - -/*------------------------------------------------------------------------- - * set label - *------------------------------------------------------------------------- - */ - - if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) - goto out; - if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) - goto out; - -/*------------------------------------------------------------------------- - * get the scale name using a static buffer - *------------------------------------------------------------------------- - */ - - if(H5DSget_label(did,DIM0,dim0_label,sizeof(dim0_label)) < 0) - goto out; - if(H5DSget_label(did,DIM1,dim1_label,sizeof(dim1_label)) < 0) - goto out; - - if(strcmp(DIM0_LABEL,dim0_label)!=0) - goto out; - if(strcmp(DIM1_LABEL,dim1_label)!=0) - goto out; - -/*------------------------------------------------------------------------- - * get the scale name using a dynamic buffer - *------------------------------------------------------------------------- - */ - - if((dim0_label_size=H5DSget_label(did,DIM0,NULL,(size_t)0)) < 0) - goto out; - if((dim1_label_size=H5DSget_label(did,DIM1,NULL,(size_t)0)) < 0) - goto out; - - /* allocate */ - dim0_labeld = (char*)malloc(dim0_label_size * sizeof(char)); - dim1_labeld = (char*)malloc(dim1_label_size * sizeof(char)); - if( dim0_labeld==NULL || dim1_labeld==NULL) - goto out; - - if(H5DSget_label(did,DIM0,dim0_labeld,(size_t)dim0_label_size) < 0) - goto out; - if(H5DSget_label(did,DIM1,dim1_labeld,(size_t)dim1_label_size) < 0) - goto out; - - if(strncmp(DIM0_LABEL,dim0_labeld,(size_t)(dim0_label_size-1))!=0) - goto out; - if(strncmp(DIM1_LABEL,dim1_labeld,(size_t)(dim1_label_size-1))!=0) - goto out; - - if(dim0_labeld) - { - free(dim0_labeld); - dim0_labeld=NULL; - } - if(dim1_labeld) - { - free(dim1_labeld); - dim1_labeld=NULL; - } - - -/*------------------------------------------------------------------------- - * get the label using a static buffer smaller than the string lenght - *------------------------------------------------------------------------- - */ - - if(H5DSget_label(did,DIM0,dim0_labels,sizeof(dim0_labels)) < 0) - goto out; - if(H5DSget_label(did,DIM1,dim1_labels,sizeof(dim1_labels)) < 0) - goto out; - - if(strncmp(DIM0_LABEL,dim0_label,sizeof(dim0_labels)-1)!=0) - goto out; - if(strncmp(DIM1_LABEL,dim1_label,sizeof(dim1_labels)-1)!=0) - goto out; - - - - if(H5Dclose(did)) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * test 5: set scale/get scale name - *------------------------------------------------------------------------- - */ - TESTING2("set scale/get scale name"); - - if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) - goto out; - - if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) - goto out; - - /* verify that DS_1_NAME is a dimension scale dataset */ - if((H5DSis_scale(dsid))==0) - goto out; - -/*------------------------------------------------------------------------- - * get the scale name using a dynamic buffer - *------------------------------------------------------------------------- - */ + if((dsid = H5Dopen2(fid,DS_11_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,SCALE_11_NAME) < 0) + goto out; + if(H5Dclose(dsid)) + goto out; - /* get the lenght of the scale name (pass NULL in name) */ - if((name_len=H5DSget_scale_name(dsid,NULL,(size_t)0)) < 0) - goto out; + if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) + goto out; + if(H5Dclose(dsid)) + goto out; - /* allocate a buffer */ - name_out = (char*)malloc(name_len * sizeof(char)); - if(name_out == NULL) - goto out; - - /* get the scale name using this buffer */ - if(H5DSget_scale_name(dsid,name_out,(size_t)name_len) < 0) - goto out; - - if(strcmp(SCALE_1_NAME,name_out)!=0) - goto out; - - if(name_out) - { - free(name_out); - name_out=NULL; - } - -/*------------------------------------------------------------------------- - * get the scale name using a static buffer - *------------------------------------------------------------------------- - */ - - /* get the scale name using this buffer */ - if(H5DSget_scale_name(dsid,sname,sizeof(sname)) < 0) - goto out; - - if(strcmp(SCALE_1_NAME,sname)!=0) - goto out; - -/*------------------------------------------------------------------------- - * get the scale name using a static buffer smaller than the string lenght - *------------------------------------------------------------------------- - */ - - /* get the scale name using this buffer */ - if(H5DSget_scale_name(dsid,snames,sizeof(snames)) < 0) - goto out; - - if(strncmp(SCALE_1_NAME,snames,sizeof(snames)-1)!=0) - goto out; - - if(H5Dclose(dsid)) - goto out; - -/*------------------------------------------------------------------------- - * add scale names - *------------------------------------------------------------------------- - */ - - if((dsid = H5Dopen2(fid,DS_11_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,SCALE_11_NAME) < 0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,SCALE_21_NAME) < 0) - goto out; - if(H5Dclose(dsid)) - goto out; - - if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,SCALE_22_NAME) < 0) - goto out; - if(H5Dclose(dsid)) - goto out; - - - PASSED(); - -/*------------------------------------------------------------------------- - * test 6: test iterate scales with a function verify_scale - *------------------------------------------------------------------------- - */ - TESTING2("iterate scales (verify scale)"); - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - dim = 0; - - /* iterate trough the 1st dimension of "dset_a" and verify that its DS is valid */ - if(H5DSiterate_scales(did,dim,NULL,verify_scale,NULL) < 0) - goto out; - - /* iterate trough the 2nd dimension of "dset_a" and verify that its DS is valid + if((dsid = H5Dopen2(fid,DS_21_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,SCALE_21_NAME) < 0) + goto out; + if(H5Dclose(dsid)) + goto out; + + if((dsid = H5Dopen2(fid,DS_22_NAME, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,SCALE_22_NAME) < 0) + goto out; + if(H5Dclose(dsid)) + goto out; + + + PASSED(); + + /*------------------------------------------------------------------------- + * test 6: test iterate scales with a function verify_scale + *------------------------------------------------------------------------- + */ + TESTING2("iterate scales (verify scale)"); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + dim = 0; + + /* iterate trough the 1st dimension of "dset_a" and verify that its DS is valid */ + if(H5DSiterate_scales(did,dim,NULL,verify_scale,NULL) < 0) + goto out; + + /* iterate trough the 2nd dimension of "dset_a" and verify that its DS is valid start at DS index 2 */ - dim = 1; - scale_idx = 2; + dim = 1; + scale_idx = 2; - if(H5DSiterate_scales(did,dim,&scale_idx,verify_scale,NULL) < 0) - goto out; + if(H5DSiterate_scales(did,dim,&scale_idx,verify_scale,NULL) < 0) + goto out; - /* close dataset ID of "dset_a" */ - if(H5Dclose(did) < 0) - goto out; + /* close dataset ID of "dset_a" */ + if(H5Dclose(did) < 0) + goto out; - PASSED(); + PASSED(); -/*------------------------------------------------------------------------- - * test 7: test iterate scales with a function read_scale - *------------------------------------------------------------------------- - */ - TESTING2("iterate scales (read scale values)"); + /*------------------------------------------------------------------------- + * test 7: test iterate scales with a function read_scale + *------------------------------------------------------------------------- + */ + TESTING2("iterate scales (read scale values)"); - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; - dim = 0; + dim = 0; - /* iterate trough the 1st dimension of "dset_a" and read the DS */ - if(H5DSiterate_scales(did,dim,NULL,read_scale,s1_wbuf) < 0) - goto out; + /* iterate trough the 1st dimension of "dset_a" and read the DS */ + if(H5DSiterate_scales(did,dim,NULL,read_scale,s1_wbuf) < 0) + goto out; - /* iterate trough the 2nd dimension of "dset_a" and read the DS + /* iterate trough the 2nd dimension of "dset_a" and read the DS start at DS index 2 */ - dim = 1; - scale_idx = 2; - - if(H5DSiterate_scales(did,dim,&scale_idx,read_scale,s22_wbuf) < 0) - goto out; - - /* close dataset ID of "dset_a" */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * test 8: test iterate scales with a function match_dim_scale - *------------------------------------------------------------------------- - */ - TESTING2("iterate scales (verify the scale sizes match)"); - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* get dataset space */ - if((sid = H5Dget_space(did)) < 0) - goto out; - - /* get rank */ - if((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; - - /* get dimensions of dataset */ - if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) - goto out; - - { - int match_size; /* does this scale size matches the dataset DIM size */ - int idx=0; /* scale index to start iterating, on return, index where iterator stoped */ - - /* iterate trough all the dimensions */ - for(dim=0; dim<(unsigned)rank; dim++) - { - if((match_size=H5DSiterate_scales(did,dim,&idx,match_dim_scale,NULL)) < 0) - goto out; - - /* "dset_a" was defined with all dimension scales size matching the size of its dimensions */ - if(match_size==0) - goto out; - - /* both DS_1_NAME and DS_2_NAME are the on the first index */ - if(idx!=0) - goto out; - } - } - - - /* close */ - if(H5Dclose(did) < 0) - goto out; - if(H5Sclose(sid) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * test 9: test iterate scales with a function match_dim_scale - *------------------------------------------------------------------------- - */ - TESTING2("iterate scales (verify the scale sizes do not match)"); - -/*------------------------------------------------------------------------- - * create 3 datasets: 1 "data" dataset and dimension scales (some are empty) - *------------------------------------------------------------------------- - */ - if(H5LTmake_dataset_int(fid,"dset_f",rank,dims,buf) < 0) - goto out; - if(H5LTmake_dataset_int(fid,"ds_f_1",rankds,s1_dim,NULL) < 0) - goto out; - if(H5LTmake_dataset_int(fid,"ds_f_11",rankds,s1_dim,s1_wbuf) < 0) - goto out; - if(H5LTmake_dataset_int(fid,"ds_f_2",rankds,s2_dim,NULL) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach them - *------------------------------------------------------------------------- - */ - if((did = H5Dopen2(fid,"dset_f", H5P_DEFAULT)) < 0) - goto out; - - if((dsid = H5Dopen2(fid,"ds_f_1", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if((dsid = H5Dopen2(fid,"ds_f_11", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if((dsid = H5Dopen2(fid,"ds_f_2", H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - if(H5Dclose(did) < 0) - goto out; - -/*------------------------------------------------------------------------- - * verify match - *------------------------------------------------------------------------- - */ - /* get the dataset id for "dset_f" */ - if((did = H5Dopen2(fid,"dset_f", H5P_DEFAULT)) < 0) - goto out; - - /* get dataset space */ - if((sid = H5Dget_space(did)) < 0) - goto out; - - /* get rank */ - if((rank=H5Sget_simple_extent_ndims(sid)) < 0) - goto out; - - /* get dimensions of dataset */ - if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) - goto out; - - { - int match_size; /* does this scale size matches the dataset DIM size */ - int idx; /* scale index to start iterating, on return, index where iterator stoped */ - - /* iterate trough all the dimensions */ - for(dim=0; dim<(unsigned)rank; dim++) - { - /* always start at 1st scale */ - idx=0; - - if((match_size=H5DSiterate_scales(did,dim,&idx,match_dim_scale,NULL)) < 0) - goto out; - - /* "dset_e" was defined with : - dim 0: 2 scales, first is empty - dim 1: 1 scale, empty */ - switch(dim) - { - case 0: /* for DIM 0, we get a valid scale at IDX 1 */ - if(match_size!=1 && idx!=1) - goto out; - break; - case 1: /* for DIM 1, we get no valid scales */ - if(match_size!=0 && idx!=0) - goto out; - }/*switch*/ - }/*for*/ - } - - /* close */ - if(H5Dclose(did) < 0) - goto out; - if(H5Sclose(sid) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * end - *------------------------------------------------------------------------- - */ - - /* close */ - H5Fclose(fid); - - return 0; - - /* error zone, gracefully close */ + dim = 1; + scale_idx = 2; + + if(H5DSiterate_scales(did,dim,&scale_idx,read_scale,s22_wbuf) < 0) + goto out; + + /* close dataset ID of "dset_a" */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * test 8: test iterate scales with a function match_dim_scale + *------------------------------------------------------------------------- + */ + TESTING2("iterate scales (verify the scale sizes match)"); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* get dataset space */ + if((sid = H5Dget_space(did)) < 0) + goto out; + + /* get rank */ + if((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; + + /* get dimensions of dataset */ + if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) + goto out; + + { + int match_size; /* does this scale size matches the dataset DIM size */ + int idx=0; /* scale index to start iterating, on return, index where iterator stoped */ + + /* iterate trough all the dimensions */ + for(dim=0; dim<(unsigned)rank; dim++) + { + if((match_size=H5DSiterate_scales(did,dim,&idx,match_dim_scale,NULL)) < 0) + goto out; + + /* "dset_a" was defined with all dimension scales size matching the size of its dimensions */ + if(match_size==0) + goto out; + + /* both DS_1_NAME and DS_2_NAME are the on the first index */ + if(idx!=0) + goto out; + } + } + + + /* close */ + if(H5Dclose(did) < 0) + goto out; + if(H5Sclose(sid) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * test 9: test iterate scales with a function match_dim_scale + *------------------------------------------------------------------------- + */ + TESTING2("iterate scales (verify the scale sizes do not match)"); + + /*------------------------------------------------------------------------- + * create 3 datasets: 1 "data" dataset and dimension scales (some are empty) + *------------------------------------------------------------------------- + */ + if(H5LTmake_dataset_int(fid,"dset_f",rank,dims,buf) < 0) + goto out; + if(H5LTmake_dataset_int(fid,"ds_f_1",rankds,s1_dim,NULL) < 0) + goto out; + if(H5LTmake_dataset_int(fid,"ds_f_11",rankds,s1_dim,s1_wbuf) < 0) + goto out; + if(H5LTmake_dataset_int(fid,"ds_f_2",rankds,s2_dim,NULL) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach them + *------------------------------------------------------------------------- + */ + if((did = H5Dopen2(fid,"dset_f", H5P_DEFAULT)) < 0) + goto out; + + if((dsid = H5Dopen2(fid,"ds_f_1", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if((dsid = H5Dopen2(fid,"ds_f_11", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if((dsid = H5Dopen2(fid,"ds_f_2", H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + if(H5Dclose(did) < 0) + goto out; + + /*------------------------------------------------------------------------- + * verify match + *------------------------------------------------------------------------- + */ + /* get the dataset id for "dset_f" */ + if((did = H5Dopen2(fid,"dset_f", H5P_DEFAULT)) < 0) + goto out; + + /* get dataset space */ + if((sid = H5Dget_space(did)) < 0) + goto out; + + /* get rank */ + if((rank=H5Sget_simple_extent_ndims(sid)) < 0) + goto out; + + /* get dimensions of dataset */ + if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) + goto out; + + { + int match_size; /* does this scale size matches the dataset DIM size */ + int idx; /* scale index to start iterating, on return, index where iterator stoped */ + + /* iterate trough all the dimensions */ + for(dim=0; dim<(unsigned)rank; dim++) + { + /* always start at 1st scale */ + idx=0; + + if((match_size=H5DSiterate_scales(did,dim,&idx,match_dim_scale,NULL)) < 0) + goto out; + + /* "dset_e" was defined with : + dim 0: 2 scales, first is empty + dim 1: 1 scale, empty */ + switch(dim) + { + case 0: /* for DIM 0, we get a valid scale at IDX 1 */ + if(match_size!=1 && idx!=1) + goto out; + break; + case 1: /* for DIM 1, we get no valid scales */ + if(match_size!=0 && idx!=0) + goto out; + }/*switch*/ + }/*for*/ + } + + /* close */ + if(H5Dclose(did) < 0) + goto out; + if(H5Sclose(sid) < 0) + goto out; + + PASSED(); + + + /*------------------------------------------------------------------------- + * end + *------------------------------------------------------------------------- + */ + + /* close */ + H5Fclose(fid); + + return 0; + + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Dclose(did); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Dclose(did); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + return FAIL; } @@ -1465,25 +1465,25 @@ out: static herr_t verify_scale(hid_t dset, unsigned dim, hid_t scale_id, void *visitor_data) { - /* define a default zero value for return. This will cause the iterator to continue */ - int ret = 0; + /* define a default zero value for return. This will cause the iterator to continue */ + int ret = 0; - /* unused */ - dset=dset; - dim=dim; - visitor_data=visitor_data; + /* unused */ + dset=dset; + dim=dim; + visitor_data=visitor_data; - /* define a positive value for return value. This will cause the iterator to + /* define a positive value for return value. This will cause the iterator to immediately return that positive value, indicating short-circuit success - */ + */ - /* the parameter DS dataset must be a valid DS dataset */ - if((H5DSis_scale(scale_id))==1) - { - ret = 1; - } + /* the parameter DS dataset must be a valid DS dataset */ + if((H5DSis_scale(scale_id))==1) + { + ret = 1; + } - return ret; + return ret; } @@ -1507,78 +1507,78 @@ static herr_t verify_scale(hid_t dset, unsigned dim, hid_t scale_id, void *visit static herr_t read_scale(hid_t dset, unsigned dim, hid_t scale_id, void *visitor_data) { - int ret = 0; /* define a default zero value for return. This will cause the iterator to continue */ - hid_t sid; /* space ID */ - hid_t tid = -1; /* file type ID */ - hid_t mtid = -1; /* memory type ID */ - hssize_t nelmts; /* number of data elements */ - char *buf=NULL; /* data buffer */ - size_t size; - int i; - char *data=visitor_data; - - /* unused */ - dset=dset; - dim=dim; - - /* get space */ - if((sid = H5Dget_space(scale_id)) < 0) - goto out; - /* get type */ - if((tid = H5Dget_type(scale_id)) < 0) - goto out; - /* get size of the DS array */ - if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) - goto out; - /* get type */ - if((mtid=H5Tget_native_type(tid,H5T_DIR_DEFAULT)) < 0) - goto out; - /* get type size */ - if((size=H5Tget_size(mtid))==0) - goto out; - - if(nelmts) - { - buf=(char *) malloc((size_t)(nelmts*size)); - if( buf==NULL) - goto out; - if(H5Dread(scale_id,mtid,H5S_ALL,H5S_ALL,H5P_DEFAULT,buf) < 0) - goto out; - - for(i=0; i<nelmts; i++) - { - if(buf[i] != data[i]) - { - printf("read and write buffers differ\n"); - goto out; - } - } - - } /* if */ - - if(H5Sclose(sid) < 0) - goto out; - if(H5Tclose(tid) < 0) - goto out; - if(H5Tclose(mtid) < 0) - goto out; - if(buf) - free(buf); - - - return ret; - - /* error zone, gracefully close */ - out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Tclose(tid); - H5Tclose(mtid); - if(buf) - free(buf); - } H5E_END_TRY; - - return FAIL; + int ret = 0; /* define a default zero value for return. This will cause the iterator to continue */ + hid_t sid; /* space ID */ + hid_t tid = -1; /* file type ID */ + hid_t mtid = -1; /* memory type ID */ + hssize_t nelmts; /* number of data elements */ + char *buf=NULL; /* data buffer */ + size_t size; + int i; + char *data=visitor_data; + + /* unused */ + dset=dset; + dim=dim; + + /* get space */ + if((sid = H5Dget_space(scale_id)) < 0) + goto out; + /* get type */ + if((tid = H5Dget_type(scale_id)) < 0) + goto out; + /* get size of the DS array */ + if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) + goto out; + /* get type */ + if((mtid=H5Tget_native_type(tid,H5T_DIR_DEFAULT)) < 0) + goto out; + /* get type size */ + if((size=H5Tget_size(mtid))==0) + goto out; + + if(nelmts) + { + buf=(char *) malloc((size_t)(nelmts*size)); + if( buf==NULL) + goto out; + if(H5Dread(scale_id,mtid,H5S_ALL,H5S_ALL,H5P_DEFAULT,buf) < 0) + goto out; + + for(i=0; i<nelmts; i++) + { + if(buf[i] != data[i]) + { + printf("read and write buffers differ\n"); + goto out; + } + } + + } /* if */ + + if(H5Sclose(sid) < 0) + goto out; + if(H5Tclose(tid) < 0) + goto out; + if(H5Tclose(mtid) < 0) + goto out; + if(buf) + free(buf); + + + return ret; + + /* error zone, gracefully close */ +out: + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Tclose(tid); + H5Tclose(mtid); + if(buf) + free(buf); + } H5E_END_TRY; + + return FAIL; } @@ -1601,66 +1601,66 @@ static herr_t read_scale(hid_t dset, unsigned dim, hid_t scale_id, void *visitor static herr_t match_dim_scale(hid_t did, unsigned dim, hid_t dsid, void *visitor_data) { - int ret = 0; /* define a default zero value for return. This will cause the iterator to continue */ - hid_t sid; /* space ID */ - hssize_t nelmts; /* size of a dimension scale array */ - hsize_t dims[H5S_MAX_RANK]; /* dimensions of dataset */ - hsize_t storage_size; + int ret = 0; /* define a default zero value for return. This will cause the iterator to continue */ + hid_t sid; /* space ID */ + hssize_t nelmts; /* size of a dimension scale array */ + hsize_t dims[H5S_MAX_RANK]; /* dimensions of dataset */ + hsize_t storage_size; /* Stop compiler from whining about "unused parameters" */ visitor_data = visitor_data; -/*------------------------------------------------------------------------- - * get DID (dataset) space info - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * get DID (dataset) space info + *------------------------------------------------------------------------- + */ - /* get dataset space */ - if((sid = H5Dget_space(did)) < 0) - goto out; + /* get dataset space */ + if((sid = H5Dget_space(did)) < 0) + goto out; - /* get dimensions of dataset */ - if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) - goto out; + /* get dimensions of dataset */ + if(H5Sget_simple_extent_dims(sid,dims,NULL) < 0) + goto out; - /* close the dataspace id */ - if(H5Sclose(sid) < 0) - goto out; + /* close the dataspace id */ + if(H5Sclose(sid) < 0) + goto out; -/*------------------------------------------------------------------------- - * get DSID (scale) space info - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * get DSID (scale) space info + *------------------------------------------------------------------------- + */ - /* get the space for the scale */ - if((sid = H5Dget_space(dsid)) < 0) - goto out; + /* get the space for the scale */ + if((sid = H5Dget_space(dsid)) < 0) + goto out; - /* get size of the DS array */ - if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) - goto out; + /* get size of the DS array */ + if((nelmts = H5Sget_simple_extent_npoints(sid)) < 0) + goto out; - /* close */ - if(H5Sclose(sid) < 0) - goto out; + /* close */ + if(H5Sclose(sid) < 0) + goto out; - /* the size of the DS array must match the dimension of the dataset */ - if(nelmts == (hssize_t)dims[dim]) - ret = 1; + /* the size of the DS array must match the dimension of the dataset */ + if(nelmts == (hssize_t)dims[dim]) + ret = 1; - /* if the scale is empty assume it cannot be used */ - storage_size=H5Dget_storage_size(dsid); + /* if the scale is empty assume it cannot be used */ + storage_size=H5Dget_storage_size(dsid); - if(storage_size==0) - ret = 0; + if(storage_size==0) + ret = 0; - return ret; + return ret; out: - H5E_BEGIN_TRY { - H5Sclose(sid); - } H5E_END_TRY; - return FAIL; + H5E_BEGIN_TRY { + H5Sclose(sid); + } H5E_END_TRY; + return FAIL; } @@ -1688,16 +1688,12 @@ static herr_t op_bogus(hid_t dset, unsigned dim, hid_t scale_id, void *visitor_d scale_id = scale_id; visitor_data = visitor_data; - /* define a default zero value for return. This will cause the iterator to continue */ - return 0; + /* define a default zero value for return. This will cause the iterator to continue */ + return 0; } - - - - /*------------------------------------------------------------------------- * test error conditions *------------------------------------------------------------------------- @@ -1705,374 +1701,374 @@ static herr_t op_bogus(hid_t dset, unsigned dim, hid_t scale_id, void *visitor_d static int test_errors(void) { - hid_t fid; /* file ID */ - int rank = RANK; /* rank of data dataset */ - int rankds = 1; /* rank of DS dataset */ - hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ - hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ - hid_t did = -1; /* dataset ID */ - hid_t dsid = -1; /* scale ID */ - hid_t gid = -1; /* group ID */ - hid_t sid = -1; /* space ID */ - hid_t sidds = -1; /* space ID */ - hsize_t pal_dims[] = {9,3}; - - printf("Testing error conditions\n"); + hid_t fid; /* file ID */ + int rank = RANK; /* rank of data dataset */ + int rankds = 1; /* rank of DS dataset */ + hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ + hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ + hid_t did = -1; /* dataset ID */ + hid_t dsid = -1; /* scale ID */ + hid_t gid = -1; /* group ID */ + hid_t sid = -1; /* space ID */ + hid_t sidds = -1; /* space ID */ + hsize_t pal_dims[] = {9,3}; + + printf("Testing error conditions\n"); + + /*------------------------------------------------------------------------- + * create a file, spaces, dataset and group ids + *------------------------------------------------------------------------- + */ + + /* create a file using default properties */ + if((fid = H5Fcreate(FILE2, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* create a group */ + if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* create the data space for the dataset */ + if((sid = H5Screate_simple(rank, dims, NULL)) < 0) + goto out; + /* create the data space for the scale */ + if((sidds = H5Screate_simple(rankds, s1_dim, NULL)) < 0) + goto out; + /* create a dataset */ + if((did = H5Dcreate2(fid, "dset_a", H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* create a dataset for the scale */ + if((dsid = H5Dcreate2(fid, "ds_a", H5T_NATIVE_INT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attempt to attach a dataset to itself, it should fail + *------------------------------------------------------------------------- + */ + + TESTING2("attach a dataset to itself"); + + if(H5DSattach_scale(did, did, 0) == SUCCEED) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * attempt to attach a group with a dataset, it should fail + *------------------------------------------------------------------------- + */ + TESTING2("attach a group with a dataset"); + + if(H5DSattach_scale(gid,dsid,0)==SUCCEED) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * attempt to attach a dataset with a group, it should fail + *------------------------------------------------------------------------- + */ + TESTING2("attach a dataset with a group"); + + if(H5DSattach_scale(did,gid,0)==SUCCEED) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * attempt to set scale for a group, it should fail + *------------------------------------------------------------------------- + */ + TESTING2("set scale for a group"); + + if(H5DSset_scale(gid,"scale 1")==SUCCEED) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * close IDs for this set + *------------------------------------------------------------------------- + */ + + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + if(H5Sclose(sid) < 0) + goto out; + if(H5Sclose(sidds) < 0) + goto out; + if(H5Gclose(gid) < 0) + goto out; + + + /*------------------------------------------------------------------------- + * try to attach a scale that has scales + *------------------------------------------------------------------------- + */ + + TESTING2("attach a scale that has scales"); + + /* create the data space for the scale */ + if((sidds = H5Screate_simple(rankds, s1_dim, NULL)) < 0) + goto out; + + /* create a dataset "ds_b" for the scale */ + if((dsid = H5Dcreate2(fid, "ds_b", H5T_NATIVE_INT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + + /* open the previous written "ds_a" */ + if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) + goto out; + + /* attach "ds_b" to "ds_a", valid */ + if(H5DSattach_scale(did, dsid, 0) < 0) + goto out; + + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + if(H5Sclose(sidds) < 0) + goto out; + + /* open the previous written "dset_a" */ + if((did = H5Dopen2(fid, "dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* open the previous written "ds_a" */ + if((dsid = H5Dopen2(fid, "ds_a", H5P_DEFAULT)) < 0) + goto out; -/*------------------------------------------------------------------------- - * create a file, spaces, dataset and group ids - *------------------------------------------------------------------------- - */ + /* try to attach "ds_a" to "dset_a", not valid */ + if(H5DSattach_scale(did,dsid,0)==SUCCEED) + goto out; - /* create a file using default properties */ - if((fid = H5Fcreate(FILE2, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* create a group */ - if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* create the data space for the dataset */ - if((sid = H5Screate_simple(rank, dims, NULL)) < 0) - goto out; - /* create the data space for the scale */ - if((sidds = H5Screate_simple(rankds, s1_dim, NULL)) < 0) - goto out; - /* create a dataset */ - if((did = H5Dcreate2(fid, "dset_a", H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* create a dataset for the scale */ - if((dsid = H5Dcreate2(fid, "ds_a", H5T_NATIVE_INT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + + /* open the previous written "ds_a" */ + if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) + goto out; + + /* open the previous written "ds_b" */ + if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) + goto out; + + /* detach "ds_b" to "ds_a" */ + if(H5DSdetach_scale(did,dsid,0) < 0) + goto out; + + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; -/*------------------------------------------------------------------------- - * attempt to attach a dataset to itself, it should fail - *------------------------------------------------------------------------- - */ + PASSED(); - TESTING2("attach a dataset to itself"); + /*------------------------------------------------------------------------- + * try to attach a dataset that is a scale + *------------------------------------------------------------------------- + */ - if(H5DSattach_scale(did, did, 0) == SUCCEED) - goto out; + TESTING2("attach to a dataset that is a scale"); - PASSED(); + /* open the previous written "ds_b", that is a scale */ + if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) + goto out; -/*------------------------------------------------------------------------- - * attempt to attach a group with a dataset, it should fail - *------------------------------------------------------------------------- - */ - TESTING2("attach a group with a dataset"); + /* open the previous written "ds_a" */ + if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) + goto out; - if(H5DSattach_scale(gid,dsid,0)==SUCCEED) - goto out; + /* try to attach "ds_a" to "ds_b", not valid */ + if(H5DSattach_scale(dsid,did,0)==SUCCEED) + goto out; + + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; - PASSED(); + PASSED(); -/*------------------------------------------------------------------------- - * attempt to attach a dataset with a group, it should fail - *------------------------------------------------------------------------- - */ - TESTING2("attach a dataset with a group"); + /*------------------------------------------------------------------------- + * try to attach a scale to an image, pallete or table + *------------------------------------------------------------------------- + */ - if(H5DSattach_scale(did,gid,0)==SUCCEED) - goto out; + TESTING2("attach to a dataset that is a reserved class dataset"); - PASSED(); + /* make an image */ + if(H5IMmake_image_8bit(fid,"image",(hsize_t)100,(hsize_t)50,NULL) < 0) + goto out; -/*------------------------------------------------------------------------- - * attempt to set scale for a group, it should fail - *------------------------------------------------------------------------- - */ - TESTING2("set scale for a group"); + /* make a palette */ + if(H5IMmake_palette(fid,"pallete",pal_dims,NULL) < 0) + goto out; - if(H5DSset_scale(gid,"scale 1")==SUCCEED) - goto out; + /* open the previous written "ds_b" */ + if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) + goto out; - PASSED(); + /* open the image dataset */ + if((did = H5Dopen2(fid,"image", H5P_DEFAULT)) < 0) + goto out; -/*------------------------------------------------------------------------- - * close IDs for this set - *------------------------------------------------------------------------- - */ + /* try to attach "ds_a" to the image, not valid */ + if(H5DSattach_scale(did,dsid,0)==SUCCEED) + goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - if(H5Sclose(sid) < 0) - goto out; - if(H5Sclose(sidds) < 0) - goto out; - if(H5Gclose(gid) < 0) - goto out; + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + PASSED(); -/*------------------------------------------------------------------------- - * try to attach a scale that has scales - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * is scale + *------------------------------------------------------------------------- + */ - TESTING2("attach a scale that has scales"); + TESTING2("is scale"); - /* create the data space for the scale */ - if((sidds = H5Screate_simple(rankds, s1_dim, NULL)) < 0) - goto out; + /* open a non scale dataset */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; - /* create a dataset "ds_b" for the scale */ - if((dsid = H5Dcreate2(fid, "ds_b", H5T_NATIVE_INT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; + /* verify that it is not a dimension scale dataset */ + if((H5DSis_scale(did))==1) + goto out; - /* open the previous written "ds_a" */ - if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) - goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; - /* attach "ds_b" to "ds_a", valid */ - if(H5DSattach_scale(did, dsid, 0) < 0) - goto out; + /* open the group. */ + if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) + goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - if(H5Sclose(sidds) < 0) - goto out; + /* verify that it is not a dimension scale dataset */ + if((H5DSis_scale(gid))==1) + goto out; - /* open the previous written "dset_a" */ - if((did = H5Dopen2(fid, "dset_a", H5P_DEFAULT)) < 0) - goto out; + /* close */ + if(H5Gclose(gid) < 0) + goto out; - /* open the previous written "ds_a" */ - if((dsid = H5Dopen2(fid, "ds_a", H5P_DEFAULT)) < 0) - goto out; + PASSED(); - /* try to attach "ds_a" to "dset_a", not valid */ - if(H5DSattach_scale(did,dsid,0)==SUCCEED) - goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; + /*------------------------------------------------------------------------- + * detach + *------------------------------------------------------------------------- + */ - /* open the previous written "ds_a" */ - if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) - goto out; + TESTING2("detach scale from dataset it is not attached to"); - /* open the previous written "ds_b" */ - if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) - goto out; + /* open the previous written "ds_a" */ + if((dsid = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) + goto out; - /* detach "ds_b" to "ds_a" */ - if(H5DSdetach_scale(did,dsid,0) < 0) - goto out; + /* open the previous written "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; + /* try to detach "ds_a" from "dset_a" */ + if(H5DSdetach_scale(did,dsid,0)==SUCCEED) + goto out; - PASSED(); + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; -/*------------------------------------------------------------------------- - * try to attach a dataset that is a scale - *------------------------------------------------------------------------- - */ + PASSED(); - TESTING2("attach to a dataset that is a scale"); - /* open the previous written "ds_b", that is a scale */ - if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) - goto out; + /*------------------------------------------------------------------------- + * detach + *------------------------------------------------------------------------- + */ - /* open the previous written "ds_a" */ - if((did = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) - goto out; + TESTING2("detach scale from group"); - /* try to attach "ds_a" to "ds_b", not valid */ - if(H5DSattach_scale(dsid,did,0)==SUCCEED) - goto out; + /* open the previous written "ds_a" */ + if((dsid = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) + goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; + /* open the group. */ + if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) + goto out; - PASSED(); + /* try to detach "ds_a" from "grp" */ + if(H5DSdetach_scale(gid,dsid,0)==SUCCEED) + goto out; -/*------------------------------------------------------------------------- - * try to attach a scale to an image, pallete or table - *------------------------------------------------------------------------- - */ + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + if(H5Gclose(gid) < 0) + goto out; - TESTING2("attach to a dataset that is a reserved class dataset"); + PASSED(); - /* make an image */ - if(H5IMmake_image_8bit(fid,"image",(hsize_t)100,(hsize_t)50,NULL) < 0) - goto out; - /* make a palette */ - if(H5IMmake_palette(fid,"pallete",pal_dims,NULL) < 0) - goto out; + /*------------------------------------------------------------------------- + * detach + *------------------------------------------------------------------------- + */ - /* open the previous written "ds_b" */ - if((dsid = H5Dopen2(fid,"ds_b", H5P_DEFAULT)) < 0) - goto out; + TESTING2("detach scale when scale is group"); - /* open the image dataset */ - if((did = H5Dopen2(fid,"image", H5P_DEFAULT)) < 0) - goto out; + /* open the previous written "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; - /* try to attach "ds_a" to the image, not valid */ - if(H5DSattach_scale(did,dsid,0)==SUCCEED) - goto out; + /* open the group. */ + if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) + goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; + /* try to detach "grp" from "dset_a" */ + if(H5DSdetach_scale(did,gid,0)==SUCCEED) + goto out; - PASSED(); + /* close */ + if(H5Dclose(did) < 0) + goto out; + if(H5Gclose(gid) < 0) + goto out; -/*------------------------------------------------------------------------- - * is scale - *------------------------------------------------------------------------- - */ - - TESTING2("is scale"); - - /* open a non scale dataset */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* verify that it is not a dimension scale dataset */ - if((H5DSis_scale(did))==1) - goto out; - - /* close */ - if(H5Dclose(did) < 0) - goto out; - - /* open the group. */ - if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) - goto out; - - /* verify that it is not a dimension scale dataset */ - if((H5DSis_scale(gid))==1) - goto out; - - /* close */ - if(H5Gclose(gid) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * detach - *------------------------------------------------------------------------- - */ - - TESTING2("detach scale from dataset it is not attached to"); - - /* open the previous written "ds_a" */ - if((dsid = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) - goto out; - - /* open the previous written "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* try to detach "ds_a" from "dset_a" */ - if(H5DSdetach_scale(did,dsid,0)==SUCCEED) - goto out; - - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * detach - *------------------------------------------------------------------------- - */ - - TESTING2("detach scale from group"); - - /* open the previous written "ds_a" */ - if((dsid = H5Dopen2(fid,"ds_a", H5P_DEFAULT)) < 0) - goto out; - - /* open the group. */ - if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) - goto out; - - /* try to detach "ds_a" from "grp" */ - if(H5DSdetach_scale(gid,dsid,0)==SUCCEED) - goto out; - - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - if(H5Gclose(gid) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * detach - *------------------------------------------------------------------------- - */ + PASSED(); - TESTING2("detach scale when scale is group"); - /* open the previous written "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; + /* close */ + if(H5Fclose(fid) < 0) + goto out; - /* open the group. */ - if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) - goto out; + return 0; - /* try to detach "grp" from "dset_a" */ - if(H5DSdetach_scale(did,gid,0)==SUCCEED) - goto out; - - /* close */ - if(H5Dclose(did) < 0) - goto out; - if(H5Gclose(gid) < 0) - goto out; - - PASSED(); - - - /* close */ - if(H5Fclose(fid) < 0) - goto out; - - return 0; - - /* error zone, gracefully close */ + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Sclose(sid); - H5Sclose(sidds); - H5Dclose(did); - H5Dclose(dsid); - H5Gclose(gid); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Sclose(sid); + H5Sclose(sidds); + H5Dclose(did); + H5Dclose(dsid); + H5Gclose(gid); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + return FAIL; } @@ -2084,182 +2080,182 @@ out: static int test_iterators(void) { - hid_t fid; /* file ID */ - int rank = RANK; /* rank of data dataset */ - int rankds = 1; /* rank of DS dataset */ - hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ - hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ - hid_t gid = -1; /* group ID */ - hid_t did; /* dataset ID */ - hid_t dsid; /* scale ID */ - char dname[30]; /* dataset name */ - int i; - - printf("Testing iterators\n"); - -/*------------------------------------------------------------------------- - * create a file, spaces, dataset and group ids - *------------------------------------------------------------------------- - */ - - /* create a file using default properties */ - if((fid=H5Fcreate(FILE3,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) - goto out; - /* create a group */ - if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* close */ - if(H5Gclose(gid) < 0) - goto out; - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,NULL) < 0) - goto out; - /* make a DS dataset */ - if(H5LTmake_dataset_int(fid,"ds_a",rankds,s1_dim,NULL) < 0) - goto out; - -/*------------------------------------------------------------------------- - * iterate when the dataset has no scales - *------------------------------------------------------------------------- - */ - - TESTING2("iterate when the dataset has no scales "); - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* try to iterate trough the 1st dimension of "dset_a", return error */ - if(H5DSiterate_scales(did,0,NULL,verify_scale,NULL) < 0) - goto out; - - /* close */ - if(H5Dclose(did) < 0) - goto out; + hid_t fid; /* file ID */ + int rank = RANK; /* rank of data dataset */ + int rankds = 1; /* rank of DS dataset */ + hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ + hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ + hid_t gid = -1; /* group ID */ + hid_t did; /* dataset ID */ + hid_t dsid; /* scale ID */ + char dname[30]; /* dataset name */ + int i; + + printf("Testing iterators\n"); + + /*------------------------------------------------------------------------- + * create a file, spaces, dataset and group ids + *------------------------------------------------------------------------- + */ + + /* create a file using default properties */ + if((fid=H5Fcreate(FILE3,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) + goto out; + /* create a group */ + if((gid = H5Gcreate2(fid, "grp", H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* close */ + if(H5Gclose(gid) < 0) + goto out; + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,NULL) < 0) + goto out; + /* make a DS dataset */ + if(H5LTmake_dataset_int(fid,"ds_a",rankds,s1_dim,NULL) < 0) + goto out; + + /*------------------------------------------------------------------------- + * iterate when the dataset has no scales + *------------------------------------------------------------------------- + */ + + TESTING2("iterate when the dataset has no scales "); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* try to iterate trough the 1st dimension of "dset_a", return error */ + if(H5DSiterate_scales(did,0,NULL,verify_scale,NULL) < 0) + goto out; + + /* close */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + + /*------------------------------------------------------------------------- + * iterate on dimension that is outside the rank + *------------------------------------------------------------------------- + */ + + TESTING2("iterate on dimension that is outside the rank "); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* try to iterate trough the 3rd dimension of "dset_a", return error */ + if(H5DSiterate_scales(did,3,NULL,verify_scale,NULL)==SUCCEED) + goto out; + + /* close */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * iterate for dimension with many scales + *------------------------------------------------------------------------- + */ + + TESTING2("iterate for dimension with many scales "); + + /* open the previously written "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + for(i=0; i<100; i++) + { + /* make a DS */ + sprintf(dname,"ds_%d",i); + if(H5LTmake_dataset_int(fid,dname,rankds,s1_dim,NULL) < 0) + goto out; + /* open */ + if((dsid = H5Dopen2(fid,dname, H5P_DEFAULT)) < 0) + goto out; + /* attach */ + if(H5DSattach_scale(did,dsid,0) < 0) + goto out; + /* close */ + if(H5Dclose(dsid) < 0) + goto out; + } - PASSED(); + /* iterate trough the 1st dimension of "dset_a" */ + if(H5DSiterate_scales(did,0,NULL,op_bogus,NULL) < 0) + goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; -/*------------------------------------------------------------------------- - * iterate on dimension that is outside the rank - *------------------------------------------------------------------------- - */ + PASSED(); - TESTING2("iterate on dimension that is outside the rank "); + /*------------------------------------------------------------------------- + * iterate on group + *------------------------------------------------------------------------- + */ - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; + TESTING2("iterate on group "); - /* try to iterate trough the 3rd dimension of "dset_a", return error */ - if(H5DSiterate_scales(did,3,NULL,verify_scale,NULL)==SUCCEED) - goto out; + /* open */ + if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) + goto out; - /* close */ - if(H5Dclose(did) < 0) - goto out; + /* try to iterate, return error */ + if(H5DSiterate_scales(gid,0,NULL,verify_scale,NULL)==SUCCEED) + goto out; - PASSED(); - -/*------------------------------------------------------------------------- - * iterate for dimension with many scales - *------------------------------------------------------------------------- - */ + /* close */ + if(H5Gclose(gid) < 0) + goto out; - TESTING2("iterate for dimension with many scales "); - - /* open the previously written "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - for(i=0; i<100; i++) - { - /* make a DS */ - sprintf(dname,"ds_%d",i); - if(H5LTmake_dataset_int(fid,dname,rankds,s1_dim,NULL) < 0) - goto out; - /* open */ - if((dsid = H5Dopen2(fid,dname, H5P_DEFAULT)) < 0) - goto out; - /* attach */ - if(H5DSattach_scale(did,dsid,0) < 0) - goto out; - /* close */ - if(H5Dclose(dsid) < 0) - goto out; - } - - /* iterate trough the 1st dimension of "dset_a" */ - if(H5DSiterate_scales(did,0,NULL,op_bogus,NULL) < 0) - goto out; - - /* close */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); + PASSED(); -/*------------------------------------------------------------------------- - * iterate on group - *------------------------------------------------------------------------- - */ - TESTING2("iterate on group "); - - /* open */ - if((gid = H5Gopen2(fid, "grp", H5P_DEFAULT)) < 0) - goto out; - - /* try to iterate, return error */ - if(H5DSiterate_scales(gid,0,NULL,verify_scale,NULL)==SUCCEED) - goto out; - - /* close */ - if(H5Gclose(gid) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * iterate in deleted scales - *------------------------------------------------------------------------- - */ + /*------------------------------------------------------------------------- + * iterate in deleted scales + *------------------------------------------------------------------------- + */ - TESTING2("iterate in deleted scales "); + TESTING2("iterate in deleted scales "); - if(H5Ldelete(fid, "ds_0", H5P_DEFAULT) < 0) - goto out; + if(H5Ldelete(fid, "ds_0", H5P_DEFAULT) < 0) + goto out; - /* open the previously written "dset_a" */ - if((did = H5Dopen2(fid, "dset_a", H5P_DEFAULT)) < 0) - goto out; + /* open the previously written "dset_a" */ + if((did = H5Dopen2(fid, "dset_a", H5P_DEFAULT)) < 0) + goto out; - /* iterate */ - if(H5DSiterate_scales(did, 0, NULL, op_bogus, NULL) == SUCCEED) - goto out; + /* iterate */ + if(H5DSiterate_scales(did, 0, NULL, op_bogus, NULL) == SUCCEED) + goto out; - /* close */ - if(H5Dclose(did) < 0) - goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; - PASSED(); + PASSED(); - /* close */ - if(H5Fclose(fid) < 0) - goto out; + /* close */ + if(H5Fclose(fid) < 0) + goto out; - return 0; + return 0; - /* error zone, gracefully close */ + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Gclose(gid); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Gclose(gid); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + return FAIL; } @@ -2270,204 +2266,204 @@ out: static int test_rank(void) { - hid_t fid; /* file ID */ - hid_t did = -1; /* dataset ID */ - hid_t dsid = -1; /* scale ID */ - hid_t sid; /* space ID */ - hid_t sidds; /* space ID */ - hsize_t dims1[1] = {DIM1_SIZE}; /* size of data dataset */ - hsize_t dims2[2] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ - hsize_t dims3[3] = {DIM1_SIZE,DIM2_SIZE,DIM3_SIZE}; /* size of data dataset */ - hsize_t dimss[2] = {1,1}; /* size of data dataset */ - char name[30]; /* dataset name buffer */ - char names[30]; /* dataset scale name buffer */ - char namel[30]; /* dataset label name buffer */ - int bufi[1]={2}; - float buff[1]={1}; - int i; - - printf("Testing ranks\n"); - -/*------------------------------------------------------------------------- - * create a file, a dataset, scales - *------------------------------------------------------------------------- - */ - - /* create a file using default properties */ - if((fid=H5Fcreate(FILE4,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) - goto out; - - /* make a dataset a 3D data dataset */ - if(H5LTmake_dataset_int(fid,"dset_a",3,dims3,NULL) < 0) - goto out; - - /* make a 1D scale dataset */ - if(H5LTmake_dataset_int(fid,"ds_a_0",1,dims1,NULL) < 0) - goto out; - - /* make a 2D scale dataset */ - if(H5LTmake_dataset_int(fid,"ds_a_1",2,dims2,NULL) < 0) - goto out; - - /* make a 3D scale dataset */ - if(H5LTmake_dataset_int(fid,"ds_a_2",3,dims3,NULL) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach - *------------------------------------------------------------------------- - */ - - TESTING2("attach"); - - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - for(i=0; i<3; i++) - { - sprintf(name,"ds_a_%d",i); - if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) - goto out; - if(H5DSattach_scale(did,dsid,(unsigned)i) < 0) - goto out; - if(H5DSis_attached(did,dsid,(unsigned)i)<=0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - } - - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - - -/*------------------------------------------------------------------------- - * detach - *------------------------------------------------------------------------- - */ - - TESTING2("detach"); - - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - for(i=0; i<3; i++) - { - sprintf(name,"ds_a_%d",i); - if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) - goto out; - if(H5DSdetach_scale(did,dsid,(unsigned)i) < 0) - goto out; - if(H5DSis_attached(did,dsid,(unsigned)i)!=0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - } - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * attach, set, get names, labels - *------------------------------------------------------------------------- - */ - - TESTING2("attach, set, get names, labels"); - - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - for(i=0; i<3; i++) - { - sprintf(name,"ds_a_%d",i); - if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) - goto out; - if(H5DSset_scale(dsid,name) < 0) - goto out; - if(H5DSattach_scale(did,dsid,(unsigned)i) < 0) - goto out; - if(H5DSis_attached(did,dsid,(unsigned)i)<=0) - goto out; - if(H5DSget_scale_name(dsid,names,sizeof(names)) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - if(H5DSset_label(did,(unsigned)i,name) < 0) - goto out; - if(H5DSget_label(did,(unsigned)i,namel,sizeof(namel)) < 0) - goto out; - if(strcmp(name,names)!=0) - goto out; - if(strcmp(name,namel)!=0) - goto out; - } - - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * attach a scalar scale - *------------------------------------------------------------------------- - */ - - TESTING2("attach a scalar scale"); - - /* create the data space for the dataset */ - if((sid = H5Screate_simple(2, dimss, NULL)) < 0) - goto out; - /* create a dataset of rank 2 */ - if((did = H5Dcreate2(fid, "dset_b", H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* create a scalar space */ - if((sidds = H5Screate(H5S_SCALAR)) < 0) - goto out; - /* create a dataset of scalar rank for the scale */ - if((dsid = H5Dcreate2(fid, "ds_b_1", H5T_NATIVE_FLOAT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) - goto out; - /* write */ - if(H5Dwrite(did, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, bufi) < 0) - goto out; - if(H5Dwrite(dsid, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buff) < 0) - goto out; - /* attach */ - if(H5DSattach_scale(did, dsid, 0) < 0) - goto out; - if(H5DSattach_scale(did, dsid, 1) < 0) - goto out; - /* close */ - if(H5Sclose(sid) < 0) - goto out; - if(H5Sclose(sidds) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - if(H5Dclose(dsid) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * close - *------------------------------------------------------------------------- - */ - if(H5Fclose(fid) < 0) - goto out; - - return 0; - - /* error zone, gracefully close */ + hid_t fid; /* file ID */ + hid_t did = -1; /* dataset ID */ + hid_t dsid = -1; /* scale ID */ + hid_t sid; /* space ID */ + hid_t sidds; /* space ID */ + hsize_t dims1[1] = {DIM1_SIZE}; /* size of data dataset */ + hsize_t dims2[2] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ + hsize_t dims3[3] = {DIM1_SIZE,DIM2_SIZE,DIM3_SIZE}; /* size of data dataset */ + hsize_t dimss[2] = {1,1}; /* size of data dataset */ + char name[30]; /* dataset name buffer */ + char names[30]; /* dataset scale name buffer */ + char namel[30]; /* dataset label name buffer */ + int bufi[1]={2}; + float buff[1]={1}; + int i; + + printf("Testing ranks\n"); + + /*------------------------------------------------------------------------- + * create a file, a dataset, scales + *------------------------------------------------------------------------- + */ + + /* create a file using default properties */ + if((fid=H5Fcreate(FILE4,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) + goto out; + + /* make a dataset a 3D data dataset */ + if(H5LTmake_dataset_int(fid,"dset_a",3,dims3,NULL) < 0) + goto out; + + /* make a 1D scale dataset */ + if(H5LTmake_dataset_int(fid,"ds_a_0",1,dims1,NULL) < 0) + goto out; + + /* make a 2D scale dataset */ + if(H5LTmake_dataset_int(fid,"ds_a_1",2,dims2,NULL) < 0) + goto out; + + /* make a 3D scale dataset */ + if(H5LTmake_dataset_int(fid,"ds_a_2",3,dims3,NULL) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach + *------------------------------------------------------------------------- + */ + + TESTING2("attach"); + + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + for(i=0; i<3; i++) + { + sprintf(name,"ds_a_%d",i); + if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) + goto out; + if(H5DSattach_scale(did,dsid,(unsigned)i) < 0) + goto out; + if(H5DSis_attached(did,dsid,(unsigned)i)<=0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + } + + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + + /*------------------------------------------------------------------------- + * detach + *------------------------------------------------------------------------- + */ + + TESTING2("detach"); + + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + for(i=0; i<3; i++) + { + sprintf(name,"ds_a_%d",i); + if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) + goto out; + if(H5DSdetach_scale(did,dsid,(unsigned)i) < 0) + goto out; + if(H5DSis_attached(did,dsid,(unsigned)i)!=0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + } + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * attach, set, get names, labels + *------------------------------------------------------------------------- + */ + + TESTING2("attach, set, get names, labels"); + + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + for(i=0; i<3; i++) + { + sprintf(name,"ds_a_%d",i); + if((dsid = H5Dopen2(fid,name, H5P_DEFAULT)) < 0) + goto out; + if(H5DSset_scale(dsid,name) < 0) + goto out; + if(H5DSattach_scale(did,dsid,(unsigned)i) < 0) + goto out; + if(H5DSis_attached(did,dsid,(unsigned)i)<=0) + goto out; + if(H5DSget_scale_name(dsid,names,sizeof(names)) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + if(H5DSset_label(did,(unsigned)i,name) < 0) + goto out; + if(H5DSget_label(did,(unsigned)i,namel,sizeof(namel)) < 0) + goto out; + if(strcmp(name,names)!=0) + goto out; + if(strcmp(name,namel)!=0) + goto out; + } + + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * attach a scalar scale + *------------------------------------------------------------------------- + */ + + TESTING2("attach a scalar scale"); + + /* create the data space for the dataset */ + if((sid = H5Screate_simple(2, dimss, NULL)) < 0) + goto out; + /* create a dataset of rank 2 */ + if((did = H5Dcreate2(fid, "dset_b", H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* create a scalar space */ + if((sidds = H5Screate(H5S_SCALAR)) < 0) + goto out; + /* create a dataset of scalar rank for the scale */ + if((dsid = H5Dcreate2(fid, "ds_b_1", H5T_NATIVE_FLOAT, sidds, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT)) < 0) + goto out; + /* write */ + if(H5Dwrite(did, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, bufi) < 0) + goto out; + if(H5Dwrite(dsid, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buff) < 0) + goto out; + /* attach */ + if(H5DSattach_scale(did, dsid, 0) < 0) + goto out; + if(H5DSattach_scale(did, dsid, 1) < 0) + goto out; + /* close */ + if(H5Sclose(sid) < 0) + goto out; + if(H5Sclose(sidds) < 0) + goto out; + if(H5Dclose(did) < 0) + goto out; + if(H5Dclose(dsid) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * close + *------------------------------------------------------------------------- + */ + if(H5Fclose(fid) < 0) + goto out; + + return 0; + + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Dclose(did); - H5Dclose(dsid); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Dclose(did); + H5Dclose(dsid); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + return FAIL; } @@ -2478,173 +2474,173 @@ out: static int test_types(void) { - hid_t fid; /* file ID */ - hid_t did = -1; /* dataset ID */ - hid_t dsid = -1; /* DS dataset ID */ - int rank = RANK; /* rank of data dataset */ - int rankds = 1; /* rank of DS dataset */ - hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ - int buf[DIM_DATA] = {1,2,3,4,5,6,7,8,9,10,11,12}; /* data of data dataset */ - hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ - hsize_t s2_dim[1] = {DIM2_SIZE}; /* size of DS 2 dataset */ - float s1_float[DIM1_SIZE] = {10,20,30}; /* data of DS 1 dataset */ - unsigned short s2_ushort[DIM2_SIZE] = {10,20,30,40}; /* data of DS 2 dataset */ - const char *s1_str = "ABC"; - const char *s2_str = "ABCD"; - - printf("Testing scales with several datatypes\n"); - -/*------------------------------------------------------------------------- - * create a file for the test - *------------------------------------------------------------------------- - */ - /* create a file using default properties */ - if((fid=H5Fcreate(FILE5,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) - goto out; - -/*------------------------------------------------------------------------- - * create datasets: 1 "data" dataset and 2 dimension scales - *------------------------------------------------------------------------- - */ - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,buf) < 0) - goto out; - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_float(fid,DS_1_NAME,rankds,s1_dim,s1_float) < 0) - goto out; - - /* make a DS dataset for the second dimension */ - if(H5LTmake_dataset(fid,DS_2_NAME,rankds,s2_dim,H5T_NATIVE_USHORT,s2_ushort) < 0) - goto out; - -/*------------------------------------------------------------------------- - * floating point and short scales - *------------------------------------------------------------------------- - */ - - TESTING2("floating point and short scales"); - - /* get the dataset id for "dset_a" */ - if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_1_NAME dimension scale to "dset_a" at dimension 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_2_NAME dimension scale to "dset_a" at dimension 1 */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* set a label */ - if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) - goto out; - if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) - goto out; - /* close */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * create datasets: 1 "data" dataset and 2 dimension scales - *------------------------------------------------------------------------- - */ - - /* make a dataset */ - if(H5LTmake_dataset_int(fid,"dset_b",rank,dims,buf) < 0) - goto out; - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_string(fid,"ds_b_1",s1_str) < 0) - goto out; - - /* make a DS dataset for the second dimension */ - if(H5LTmake_dataset_string(fid,"ds_b_2",s2_str) < 0) - goto out; - -/*------------------------------------------------------------------------- - * floating point and short scales - *------------------------------------------------------------------------- - */ - - TESTING2("string scales"); - - /* get the dataset id for "dset_b" */ - if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) - goto out; - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_b_1", H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_1_NAME dimension scale to "dset_b" at dimension 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"ds_b_2", H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_2_NAME dimension scale to "dset_b" at dimension 1 */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* set a label */ - if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) - goto out; - if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) - goto out; - /* close */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - -/*------------------------------------------------------------------------- - * close - *------------------------------------------------------------------------- - */ - if(H5Fclose(fid) < 0) - goto out; - - return 0; - - /* error zone, gracefully close */ + hid_t fid; /* file ID */ + hid_t did = -1; /* dataset ID */ + hid_t dsid = -1; /* DS dataset ID */ + int rank = RANK; /* rank of data dataset */ + int rankds = 1; /* rank of DS dataset */ + hsize_t dims[RANK] = {DIM1_SIZE,DIM2_SIZE}; /* size of data dataset */ + int buf[DIM_DATA] = {1,2,3,4,5,6,7,8,9,10,11,12}; /* data of data dataset */ + hsize_t s1_dim[1] = {DIM1_SIZE}; /* size of DS 1 dataset */ + hsize_t s2_dim[1] = {DIM2_SIZE}; /* size of DS 2 dataset */ + float s1_float[DIM1_SIZE] = {10,20,30}; /* data of DS 1 dataset */ + unsigned short s2_ushort[DIM2_SIZE] = {10,20,30,40}; /* data of DS 2 dataset */ + const char *s1_str = "ABC"; + const char *s2_str = "ABCD"; + + printf("Testing scales with several datatypes\n"); + + /*------------------------------------------------------------------------- + * create a file for the test + *------------------------------------------------------------------------- + */ + /* create a file using default properties */ + if((fid=H5Fcreate(FILE5,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * create datasets: 1 "data" dataset and 2 dimension scales + *------------------------------------------------------------------------- + */ + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_a",rank,dims,buf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_float(fid,DS_1_NAME,rankds,s1_dim,s1_float) < 0) + goto out; + + /* make a DS dataset for the second dimension */ + if(H5LTmake_dataset(fid,DS_2_NAME,rankds,s2_dim,H5T_NATIVE_USHORT,s2_ushort) < 0) + goto out; + + /*------------------------------------------------------------------------- + * floating point and short scales + *------------------------------------------------------------------------- + */ + + TESTING2("floating point and short scales"); + + /* get the dataset id for "dset_a" */ + if((did = H5Dopen2(fid,"dset_a", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_1_NAME, H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_1_NAME dimension scale to "dset_a" at dimension 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,DS_2_NAME, H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_2_NAME dimension scale to "dset_a" at dimension 1 */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* set a label */ + if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) + goto out; + if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) + goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * create datasets: 1 "data" dataset and 2 dimension scales + *------------------------------------------------------------------------- + */ + + /* make a dataset */ + if(H5LTmake_dataset_int(fid,"dset_b",rank,dims,buf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_string(fid,"ds_b_1",s1_str) < 0) + goto out; + + /* make a DS dataset for the second dimension */ + if(H5LTmake_dataset_string(fid,"ds_b_2",s2_str) < 0) + goto out; + + /*------------------------------------------------------------------------- + * floating point and short scales + *------------------------------------------------------------------------- + */ + + TESTING2("string scales"); + + /* get the dataset id for "dset_b" */ + if((did = H5Dopen2(fid,"dset_b", H5P_DEFAULT)) < 0) + goto out; + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_b_1", H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_1_NAME dimension scale to "dset_b" at dimension 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"ds_b_2", H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_2_NAME dimension scale to "dset_b" at dimension 1 */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* set a label */ + if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) + goto out; + if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) + goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + /*------------------------------------------------------------------------- + * close + *------------------------------------------------------------------------- + */ + if(H5Fclose(fid) < 0) + goto out; + + return 0; + + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Dclose(did); - H5Dclose(dsid); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Dclose(did); + H5Dclose(dsid); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + return FAIL; } /*------------------------------------------------------------------------- @@ -2654,149 +2650,162 @@ out: static int test_data(void) { - hid_t fid; /* file ID */ - hid_t did = -1; /* dataset ID */ - hid_t dsid = -1; /* DS dataset ID */ - hid_t dcpl; /* dataset creation property list */ - hid_t sid; /* dataspace ID */ - float *vals=NULL; /* array to hold data values */ - float *latbuf=NULL; /* array to hold the latitude values */ - float *lonbuf=NULL; /* array to hold the longitude values */ - hsize_t dims[2]; /* array to hold dimensions */ - hsize_t latdims[1]; /* array to hold dimensions */ - hsize_t londims[1]; /* array to hold dimensions */ - float fill=-99; /* fill value */ - - - printf("Testing reading ASCII data and generate HDF5 data with scales\n"); - -/*------------------------------------------------------------------------- - * create a file for the test - *------------------------------------------------------------------------- - */ - /* create a file using default properties */ - if((fid=H5Fcreate(FILE6,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) - goto out; - -/*------------------------------------------------------------------------- - * read data - *------------------------------------------------------------------------- - */ - - /* read ASCII bathymetry data */ - if(read_data("dsdata.txt",2,dims,&vals) < 0) - goto out; - - /* read the latitude */ - if(read_data("dslat.txt",1,latdims,&latbuf) < 0) - goto out; - - /* read the longitude */ - if(read_data("dslon.txt",1,londims,&lonbuf) < 0) - goto out; - -/*------------------------------------------------------------------------- - * generating scales - *------------------------------------------------------------------------- - */ - - TESTING2("generating scales"); - -/*------------------------------------------------------------------------- - * create datasets: 1 "data" dataset and 2 dimension scales - *------------------------------------------------------------------------- - */ - - /* make a DS dataset for the first dimension */ - if(H5LTmake_dataset_float(fid, "lat", 1, latdims, latbuf) < 0) - goto out; - - /* make a DS dataset for the second dimension */ - if(H5LTmake_dataset_float(fid, "lon", 1, londims, lonbuf) < 0) - goto out; - - /* make a dataset for the data. a fill value is set */ - if((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) - goto out; - if(H5Pset_fill_value(dcpl, H5T_NATIVE_FLOAT, &fill) < 0) - goto out; - if((sid = H5Screate_simple(2, dims, NULL)) < 0) - goto out; - if((did = H5Dcreate2(fid, "data", H5T_NATIVE_FLOAT, sid, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) - goto out; - if(H5Dwrite(did, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, vals) < 0) - goto out; - if(H5Dclose(did) < 0) - goto out; - if(H5Pclose(dcpl) < 0) - goto out; - if(H5Sclose(sid) < 0) - goto out; - -/*------------------------------------------------------------------------- - * attach - *------------------------------------------------------------------------- - */ - - /* get the dataset id for "data" */ - if((did = H5Dopen2(fid,"data", H5P_DEFAULT)) < 0) - goto out; - - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"lat", H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_1_NAME dimension scale to "data" at dimension 0 */ - if(H5DSattach_scale(did,dsid,DIM0) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* get the DS dataset id */ - if((dsid = H5Dopen2(fid,"lon", H5P_DEFAULT)) < 0) - goto out; - /* attach the DS_2_NAME dimension scale to "data" at dimension 1 */ - if(H5DSattach_scale(did,dsid,DIM1) < 0) - goto out; - /* set name */ - if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) - goto out; - /* close DS id */ - if(H5Dclose(dsid) < 0) - goto out; - /* set a label */ - if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) - goto out; - if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) - goto out; - /* close */ - if(H5Dclose(did) < 0) - goto out; - - PASSED(); - - - -/*------------------------------------------------------------------------- - * close - *------------------------------------------------------------------------- - */ - if(H5Fclose(fid) < 0) - goto out; - - return 0; - - /* error zone, gracefully close */ + hid_t fid; /* file ID */ + hid_t did = -1; /* dataset ID */ + hid_t dsid = -1; /* DS dataset ID */ + hid_t dcpl; /* dataset creation property list */ + hid_t sid; /* dataspace ID */ + float *vals=NULL; /* array to hold data values */ + float *latbuf=NULL; /* array to hold the latitude values */ + float *lonbuf=NULL; /* array to hold the longitude values */ + hsize_t dims[2]; /* array to hold dimensions */ + hsize_t latdims[1]; /* array to hold dimensions */ + hsize_t londims[1]; /* array to hold dimensions */ + float fill=-99; /* fill value */ + + + printf("Testing reading ASCII data and generate HDF5 data with scales\n"); + + /*------------------------------------------------------------------------- + * create a file for the test + *------------------------------------------------------------------------- + */ + /* create a file using default properties */ + if((fid=H5Fcreate(FILE6,H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT)) < 0) + goto out; + + /*------------------------------------------------------------------------- + * generating scales + *------------------------------------------------------------------------- + */ + + TESTING2("generating scales"); + + /*------------------------------------------------------------------------- + * create datasets: 1 "data" dataset and 2 dimension scales + *------------------------------------------------------------------------- + */ + + /* read the latitude */ + if(read_data("dslat.txt",1,latdims,&latbuf) < 0) + goto out; + + /* make a DS dataset for the first dimension */ + if(H5LTmake_dataset_float(fid, "lat", 1, latdims, latbuf) < 0) + goto out; + + free( latbuf ); + latbuf = NULL; + + /* read the longitude */ + if(read_data("dslon.txt",1,londims,&lonbuf) < 0) + goto out; + + /* make a DS dataset for the second dimension */ + if(H5LTmake_dataset_float(fid, "lon", 1, londims, lonbuf) < 0) + goto out; + + free( lonbuf ); + lonbuf = NULL; + + /* make a dataset for the data. a fill value is set */ + if((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0) + goto out; + if(H5Pset_fill_value(dcpl, H5T_NATIVE_FLOAT, &fill) < 0) + goto out; + + /* read ASCII bathymetry data and dimensions to create dataset */ + if(read_data("dsdata.txt",2,dims,&vals) < 0) + goto out; + + if((sid = H5Screate_simple(2, dims, NULL)) < 0) + goto out; + if((did = H5Dcreate2(fid, "data", H5T_NATIVE_FLOAT, sid, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0) + goto out; + if(H5Dwrite(did, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, vals) < 0) + goto out; + + free ( vals ); + vals = NULL; + + if(H5Dclose(did) < 0) + goto out; + if(H5Pclose(dcpl) < 0) + goto out; + if(H5Sclose(sid) < 0) + goto out; + + /*------------------------------------------------------------------------- + * attach + *------------------------------------------------------------------------- + */ + + /* get the dataset id for "data" */ + if((did = H5Dopen2(fid,"data", H5P_DEFAULT)) < 0) + goto out; + + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"lat", H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_1_NAME dimension scale to "data" at dimension 0 */ + if(H5DSattach_scale(did,dsid,DIM0) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_1_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* get the DS dataset id */ + if((dsid = H5Dopen2(fid,"lon", H5P_DEFAULT)) < 0) + goto out; + /* attach the DS_2_NAME dimension scale to "data" at dimension 1 */ + if(H5DSattach_scale(did,dsid,DIM1) < 0) + goto out; + /* set name */ + if(H5DSset_scale(dsid,SCALE_2_NAME) < 0) + goto out; + /* close DS id */ + if(H5Dclose(dsid) < 0) + goto out; + /* set a label */ + if(H5DSset_label(did,DIM0,DIM0_LABEL) < 0) + goto out; + if(H5DSset_label(did,DIM1,DIM1_LABEL) < 0) + goto out; + /* close */ + if(H5Dclose(did) < 0) + goto out; + + PASSED(); + + + + /*------------------------------------------------------------------------- + * close + *------------------------------------------------------------------------- + */ + if(H5Fclose(fid) < 0) + goto out; + + return 0; + + /* error zone, gracefully close */ out: - H5E_BEGIN_TRY { - H5Dclose(did); - H5Dclose(dsid); - H5Fclose(fid); - } H5E_END_TRY; - H5_FAILED(); - return FAIL; + H5E_BEGIN_TRY { + H5Dclose(did); + H5Dclose(dsid); + H5Fclose(fid); + } H5E_END_TRY; + H5_FAILED(); + + if (latbuf) + free( latbuf ); + if (lonbuf) + free( lonbuf ); + if (vals) + free( vals ); + return FAIL; } @@ -2814,52 +2823,61 @@ out: *------------------------------------------------------------------------- */ -static int read_data( const char* fname, int ndims, hsize_t *dims, float **buf ) +static int read_data( const char* fname, + int ndims, + hsize_t *dims, + float **buf ) { - int i, n; - unsigned j; - char str[20]; - size_t nelms; - FILE *f; - float val; - char *srcdir = getenv("srcdir"); /* the source directory */ - char data_file[512]; /* buffer to hold name of existing data file */ - - strcpy(data_file, ""); - /* compose the name of the file to open, using the srcdir, if appropriate */ - if(srcdir) - { - strcpy(data_file, srcdir); - strcat(data_file, "/"); - } - /* read first data file */ - strcat(data_file,fname); - - f = fopen(data_file, "r"); - if( f == NULL ) - { - printf( "Could not open file %s\n", data_file ); - return -1; - } - - for(i=0, nelms=1; i < ndims; i++) - { - fscanf( f, "%s %u", str, &j); - fscanf( f, "%d",&n ); - dims[i] = n; - nelms *= n; - } - - *buf = (float*) malloc (nelms * sizeof( float )); - - for(j = 0; j < nelms; j++) - { - fscanf( f, "%f",&val ); - (*buf)[j] = val; - } - fclose(f); - - return 1; + int i, n; + unsigned j; + char str[20]; + size_t nelms; + FILE *f; + float val; + char *srcdir = getenv("srcdir"); /* the source directory */ + char data_file[512]; /* buffer to hold name of existing data file */ + + strcpy(data_file, ""); + /* compose the name of the file to open, using the srcdir, if appropriate */ + if(srcdir) + { + strcpy(data_file, srcdir); + strcat(data_file, "/"); + } + /* read first data file */ + strcat(data_file,fname); + + f = fopen(data_file, "r"); + if( f == NULL ) + { + printf( "Could not open file %s\n", data_file ); + return -1; + } + + for(i=0, nelms=1; i < ndims; i++) + { + fscanf( f, "%s %u", str, &j); + fscanf( f, "%d",&n ); + dims[i] = n; + nelms *= n; + } + + *buf = (float*) malloc (nelms * sizeof( float )); + + if ( *buf == NULL ) + { + printf( "memory allocation failed\n" ); + return -1; + } + + for(j = 0; j < nelms; j++) + { + fscanf( f, "%f",&val ); + (*buf)[j] = val; + } + fclose(f); + + return 1; } |