path: root/hl/test/test_ds.c

/****************************************************************************
 * NCSA HDF                                                                 *
 * Scientific Data Technologies                                             *
 * National Center for Supercomputing Applications                          *
 * University of Illinois at Urbana-Champaign                               *
 * 605 E. Springfield, Champaign IL 61820                                   *
 *                                                                          *
 * For conditions of distribution and use, see the accompanying             *
 * hdf/COPYING f.                                                           *
 *                                                                          *
 ****************************************************************************/



#include "H5DS.h"
#include "H5LT.h"
#include <stdlib.h>
#include <string.h>
#include "H5IM.h"


/* operator functions */
static herr_t verifiy_scale(hid_t dset, unsigned dim, hid_t scale, void *visitor_data);
static herr_t read_scale(hid_t dset, unsigned dim, hid_t scale, void *visitor_data);
static herr_t match_dim_scale(hid_t did, unsigned dim, hid_t dsid, void *visitor_data);


/* prototypes */
static int test_simple(void);
static int test_errors(void);
static int test_rank(void);



#define RANK          2
#define DIM_DATA      12
#define DIM1_SIZE     3
#define DIM2_SIZE     4
#define DIM3_SIZE     2
#define DIM0          0
#define DIM1          1
#define DIM2          2

#define DS_1_NAME   "ds_a_1"
#define DS_11_NAME  "ds_a_11"
#define DS_2_NAME   "ds_a_2"
#define DS_21_NAME  "ds_a_21"
#define DS_22_NAME  "ds_a_22"
#define DS_3_NAME   "ds_a_3"


#define SCALE_1_NAME   "Latitude set 0"
#define SCALE_11_NAME  "Latitude set 1"
#define SCALE_2_NAME   "Longitude set 0"
#define SCALE_21_NAME  "Longitude set 1"
#define SCALE_22_NAME  "Longitude set 2"

#define DIM0_LABEL   "Latitude"
#define DIM1_LABEL   "Longitude"




/*-------------------------------------------------------------------------
 * the main program
 *-------------------------------------------------------------------------
 */
int main(void)
{
 int	nerrors=0;

 nerrors += test_simple()<0 	?1:0;
 nerrors += test_errors()<0 	?1:0;
 nerrors += test_rank()<0 	?1: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;
}


/*-------------------------------------------------------------------------
 * DS API test
 *
 * Functions tested:
 *
 * H5DSattach_scale
 * H5DSdetach_scale
 * H5DSset_label
 * H5DSget_label
 * H5DSset_scale
 * H5DSget_scale_name
 * H5DSis_scale
 *
 *-------------------------------------------------------------------------
 */


static int test_simple(void)
{
 hid_t   fid;                                              /* file ID */
 hid_t   did;                                              /* 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("test_ds1.h5",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 = H5Dopen(fid,"dset_a"))<0)
  goto out;

/*-------------------------------------------------------------------------
 * attach the DS_1_NAME dimension scale to "dset_a"
 *-------------------------------------------------------------------------
 */  

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,DS_1_NAME))<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 = H5Dopen(fid,DS_11_NAME))<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 = H5Dopen(fid,DS_2_NAME))<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 = H5Dopen(fid,DS_21_NAME))<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 = H5Dopen(fid,DS_22_NAME))<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 = H5Dopen(fid,DS_1_NAME))<0)
  goto out;
 if (H5DSis_attached(did,dsid,DIM0)<=0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_2_NAME))<0)
  goto out;
 if (H5DSis_attached(did,dsid,DIM1)<=0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_21_NAME))<0)
  goto out;
 if (H5DSis_attached(did,dsid,DIM1)<=0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_22_NAME))<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 = H5Dopen(fid,"dset_a"))<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 = H5Dopen(fid,"dset_b"))<0)
  goto out;
 if ((dsid = H5Dopen(fid,"ds_b_1"))<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 = H5Dopen(fid,"dset_b"))<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 = H5Dopen(fid,"dset_c"))<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_c_1"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,0)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_c_2"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,1)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_c_21"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,1)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_c_22"))<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 = H5Dopen(fid,"dset_c"))<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 = H5Dopen(fid,"dset_c"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_c_21"))<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 = H5Dopen(fid,"dset_c"))<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 = H5Dopen(fid,"dset_c"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_c_22"))<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 = H5Dopen(fid,"dset_c"))<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 = H5Dopen(fid,"dset_c"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_c_2"))<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 = H5Dopen(fid,"dset_c"))<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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_d_1"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,0)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;
 if ((dsid = H5Dopen(fid,"ds_d_2"))<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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_d_1"))<0)
  goto out;
 if (H5DSis_attached(did,dsid,DIM0)<=0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_d_2"))<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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_d_1"))<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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_d_1"))<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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_d_2"))<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)==SUCCESS) 
  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 = H5Dopen(fid,"dset_d"))<0)
  goto out;

 /* get the DS dataset id */
 if ((dsid = H5Dopen(fid,"ds_d_2"))<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=H5Gcreate(fid,"grp",0))<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 = H5Dcreate(gid,dname,H5T_NATIVE_INT,sid,H5P_DEFAULT))<0)
   goto out;
  sprintf(sname,"ds_%d",i);
  if((dsid = H5Dcreate(gid,sname,H5T_NATIVE_INT,sid,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 = H5Dopen(gid,dname))<0)
   goto out;
  for (j=0; j<5; j++)  
  {
   sprintf(sname,"ds_%d",j);
   if((dsid = H5Dopen(gid,sname))<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 = H5Dopen(gid,dname))<0)
   goto out;
  for (j=0; j<5; j++)  
  {
   sprintf(sname,"ds_%d",j);
   if((dsid = H5Dopen(gid,sname))<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 = H5Dopen(gid,dname))<0)
   goto out;
  for (j=0; j<5; j++)  
  {
   sprintf(sname,"ds_%d",j);
   if((dsid = H5Dopen(gid,sname))<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();




/*-------------------------------------------------------------------------
 * test 4: set/get label
 *-------------------------------------------------------------------------
 */ 
 
 TESTING2("set/get label");

 if ((did = H5Dopen(fid,"dset_a"))<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,0))<0)
  goto out;
 if ((dim1_label_size=H5DSget_label(did,DIM1,NULL,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,dim0_label_size)<0)
  goto out;
 if (H5DSget_label(did,DIM1,dim1_labeld,dim1_label_size)<0)
  goto out;
 
 if (strncmp(DIM0_LABEL,dim0_labeld,dim0_label_size-1)!=0)
  goto out;
 if (strncmp(DIM1_LABEL,dim1_labeld,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 = H5Dopen(fid,DS_1_NAME))<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,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,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 = H5Dopen(fid,DS_11_NAME))<0)
  goto out;
 if (H5DSset_scale(dsid,SCALE_11_NAME)<0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_2_NAME))<0)
  goto out;
 if (H5DSset_scale(dsid,SCALE_2_NAME)<0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_21_NAME))<0)
  goto out;
 if (H5DSset_scale(dsid,SCALE_21_NAME)<0)
  goto out;
 if (H5Dclose(dsid))
  goto out;

 if ((dsid = H5Dopen(fid,DS_22_NAME))<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 verifiy_scale
 *-------------------------------------------------------------------------
 */ 
 TESTING2("iterate scales (verify scale)");

 /* get the dataset id for "dset_a" */
 if ((did = H5Dopen(fid,"dset_a"))<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,verifiy_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;

 if (H5DSiterate_scales(did,dim,&scale_idx,verifiy_scale,NULL)<0)
  goto out;

 /* close dataset ID of "dset_a" */
 if (H5Dclose(did)<0)
  goto out;

 PASSED();


/*-------------------------------------------------------------------------
 * 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 = H5Dopen(fid,"dset_a"))<0)
  goto out;

 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 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 = H5Dopen(fid,"dset_a"))<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_e",rank,dims,buf)<0)
  goto out;
 if (H5LTmake_dataset_int(fid,"ds_e_1",rankds,s1_dim,NULL)<0)
  goto out;
 if (H5LTmake_dataset_int(fid,"ds_e_11",rankds,s1_dim,s1_wbuf)<0)
  goto out;
 if (H5LTmake_dataset_int(fid,"ds_e_2",rankds,s2_dim,NULL)<0)
  goto out;

/*-------------------------------------------------------------------------
 * attach them
 *-------------------------------------------------------------------------
 */  
 if ((did = H5Dopen(fid,"dset_e"))<0)
  goto out;

 if ((dsid = H5Dopen(fid,"ds_e_1"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,DIM0)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;
 if ((dsid = H5Dopen(fid,"ds_e_11"))<0)
  goto out;
 if (H5DSattach_scale(did,dsid,DIM0)<0)
  goto out;
 if (H5Dclose(dsid)<0)
  goto out;
 if ((dsid = H5Dopen(fid,"ds_e_2"))<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_e" */
 if ((did = H5Dopen(fid,"dset_e"))<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;
}



/*-------------------------------------------------------------------------
 * Function: verifiy_scale
 *
 * Purpose: example operator function used by H5DSiterate_scales, used
 *  to verify that SCALE_ID refers to a valid DS dataset
 *
 * Return: 
 * The return values from an operator are: 
 * Zero causes the iterator to continue, returning zero when all group members have been processed. 
 * Positive causes the iterator to immediately return that positive value, indicating 
 *  short-circuit success. The iterator can be restarted at the next group member. 
 * Negative causes the iterator to immediately return that value, indicating failure. 
 *  The iterator can be restarted at the next group member. 
 *
 *-------------------------------------------------------------------------
 */

static herr_t verifiy_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;  
  
 /* unused */
 dset=dset;
 dim=dim;
 visitor_data=visitor_data;

 /* 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;
 }
 
 return ret;
} 


/*-------------------------------------------------------------------------
 * Function: read_scale
 *
 * Purpose: example operator function used by H5DSiterate_scales, used
 *  to read data from DS dataset. compare read and write buffers
 *
 * Return: 
 * The return values from an operator are: 
 * Zero causes the iterator to continue, returning zero when all group members have been processed. 
 * Positive causes the iterator to immediately return that positive value, indicating 
 *  short-circuit success. The iterator can be restarted at the next group member. 
 * Negative causes the iterator to immediately return that value, indicating failure. 
 *  The iterator can be restarted at the next group member 
 *
 *
 *-------------------------------------------------------------------------
 */

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;       /* file type ID */
 hid_t    mtid;      /* 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((unsigned)(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;
} 


/*-------------------------------------------------------------------------
 * Function: match_dim_scale
 *
 * Purpose: example operator function used by H5DSiterate_scales, used
 *  to verify the the DSID scale size matches the dataset DIM size
 *
 * Return: 
 * The return values from an operator are: 
 * Zero causes the iterator to continue, returning zero when all group members have been processed. 
 * Positive causes the iterator to immediately return that positive value, indicating 
 *  short-circuit success. The iterator can be restarted at the next group member. 
 * Negative causes the iterator to immediately return that value, indicating failure. 
 *  The iterator can be restarted at the next group member. 
 *
 *-------------------------------------------------------------------------
 */

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 */
 int       rank;                 /* rank of dataset */
 hsize_t   dims[H5S_MAX_RANK];   /* dimensions of dataset */
 hsize_t   storage_size;

/*-------------------------------------------------------------------------
 * get DID (dataset) space info
 *-------------------------------------------------------------------------
 */ 

 /* 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;

 /* close the dataspace id */
 if (H5Sclose(sid)<0)
  goto out;
 
/*-------------------------------------------------------------------------
 * get DSID (scale) space info
 *-------------------------------------------------------------------------
 */ 
 
 /* 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;

 /* 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;

 /* if the scale is empty assume it cannot be used */
 if ((storage_size=H5Dget_storage_size(dsid))<0)
  goto out;

 if (storage_size==0)
  ret = 0;

 return ret;
 
out:
 H5E_BEGIN_TRY {
  H5Sclose(sid);
 } H5E_END_TRY;
 return FAIL;
} 



/*-------------------------------------------------------------------------
 * test error conditions
 *-------------------------------------------------------------------------
 */

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 */
 hsize_t s2_dim[1]  = {DIM2_SIZE};                         /* size of DS 2 dataset */
 hid_t   did;                                              /* dataset ID */
 hid_t   dsid;                                             /* scale ID */
 hid_t   gid;                                              /* group ID */
 hid_t   sid;                                              /* space ID */
 hid_t   sidds;                                            /* 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("test_ds2.h5",H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT))<0)
  goto out;
 /* create a group */
 if ((gid=H5Gcreate(fid,"grp",0))<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=H5Dcreate(fid,"dset_a",H5T_NATIVE_INT,sid,H5P_DEFAULT))<0)
  goto out;
 /* create a dataset for the scale */
 if ((dsid=H5Dcreate(fid,"ds_a",H5T_NATIVE_INT,sidds,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)==SUCCESS)
  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)==SUCCESS)
  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)==SUCCESS)
  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")==SUCCESS)
  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=H5Dcreate(fid,"ds_b",H5T_NATIVE_INT,sidds,H5P_DEFAULT))<0)
  goto out;

 /* open the previous written "ds_a" */
 if ((did = H5Dopen(fid,"ds_a"))<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 = H5Dopen(fid,"dset_a"))<0)
  goto out;

  /* open the previous written "ds_a" */
 if ((dsid = H5Dopen(fid,"ds_a"))<0)
  goto out;

 /* try to attach "ds_a" to "dset_a", not valid */
 if(H5DSattach_scale(did,dsid,0)==SUCCESS)
  goto out;

 /* close */
 if (H5Dclose(dsid)<0)
  goto out;
 if (H5Dclose(did)<0)
  goto out;

 PASSED();

/*-------------------------------------------------------------------------
 * try to attach a dataset that is a scale 
 *-------------------------------------------------------------------------
 */ 

 TESTING2("attach to a dataset that is a scale");

 /* open the previous written "ds_b", that is a scale */
 if ((dsid = H5Dopen(fid,"ds_b"))<0)
  goto out;

 /* open the previous written "ds_a" */
 if ((did = H5Dopen(fid,"ds_a"))<0)
  goto out;

 /* try to attach "ds_a" to "ds_b", not valid */
 if(H5DSattach_scale(dsid,did,0)==SUCCESS)
  goto out;

 /* close */
 if (H5Dclose(dsid)<0)
  goto out;
 if (H5Dclose(did)<0)
  goto out;

 PASSED();

/*-------------------------------------------------------------------------
 * try to attach a scale to an image, pallete or table
 *-------------------------------------------------------------------------
 */ 

 TESTING2("attach to a dataset that is a reserved class dataset");

 /* make an image */
 if (H5IMmake_image_8bit(fid,"image",100,50,NULL)<0)
  goto out;

 /* make a palette */
 if (H5IMmake_palette(fid,"pallete",pal_dims,NULL)<0)
  goto out;

 /* open the previous written "ds_b" */
 if ((dsid = H5Dopen(fid,"ds_b"))<0)
  goto out;

 /* open the image dataset */
 if ((did = H5Dopen(fid,"image"))<0)
  goto out;

 /* try to attach "ds_a" to the image, not valid */
 if(H5DSattach_scale(did,dsid,0)==SUCCESS)
  goto out;

 /* close */
 if (H5Dclose(dsid)<0)
  goto out;
 if (H5Dclose(did)<0)
  goto out;

 PASSED();

 /* close */
 if (H5Fclose(fid)<0)
  goto out;

 return 0;
  
 /* 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;
}




/*-------------------------------------------------------------------------
 * test several rank and types
 *-------------------------------------------------------------------------
 */

static int test_rank(void)
{
 hid_t   fid;                                              /* file ID */
 hid_t   did;                                              /* dataset ID */
 hid_t   dsid;                                             /* scale ID */
 int     rank     = 3;                                     /* rank of data dataset */
 hsize_t dims[3]  = {DIM1_SIZE,DIM2_SIZE,DIM3_SIZE};       /* 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     i;
  
 printf("Testing ranks\n");

/*-------------------------------------------------------------------------
 * create a file, a dataset, scales
 *-------------------------------------------------------------------------
 */  

 /* create a file using default properties */
 if ((fid=H5Fcreate("test_ds3.h5",H5F_ACC_TRUNC,H5P_DEFAULT,H5P_DEFAULT))<0)
  goto out;
 
 /* make a dataset */
 if (H5LTmake_dataset_int(fid,"dset_a",rank,dims,NULL)<0)
  goto out;

 /* make scale datasets  */
 for (i=0; i<rank; i++) 
 {
  sprintf(name,"ds_a_%d",i);
  if (H5LTmake_dataset_int(fid,name,(rank-i),&dims[i],NULL)<0)
   goto out;
 }

 
/*-------------------------------------------------------------------------
 * attach 
 *-------------------------------------------------------------------------
 */ 

 TESTING2("attach");
 
 if ((did = H5Dopen(fid,"dset_a"))<0)
  goto out;
 
 for (i=0; i<rank; i++) 
 {
  sprintf(name,"ds_a_%d",i);
  if((dsid = H5Dopen(fid,name))<0)
   goto out;
  if(H5DSattach_scale(did,dsid,i)<0)
   goto out;
  if (H5DSis_attached(did,dsid,i)<=0)
   goto out;
  if (H5Dclose(dsid)<0)
   goto out;
 }
 
 if (H5Dclose(did)<0)
  goto out;
 
 PASSED();

 
/*-------------------------------------------------------------------------
 * detach 
 *-------------------------------------------------------------------------
 */ 

 TESTING2("detach");
 
 if ((did = H5Dopen(fid,"dset_a"))<0)
  goto out;
 
 for (i=0; i<rank; i++) 
 {
  sprintf(name,"ds_a_%d",i);
  if((dsid = H5Dopen(fid,name))<0)
   goto out;
  if(H5DSdetach_scale(did,dsid,i)<0)
   goto out;
  if (H5DSis_attached(did,dsid,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 = H5Dopen(fid,"dset_a"))<0)
  goto out;
 
 for (i=0; i<rank; i++) 
 {
  sprintf(name,"ds_a_%d",i);
  if((dsid = H5Dopen(fid,name))<0)
   goto out;
  if (H5DSset_scale(dsid,name)<0)
   goto out;
  if(H5DSattach_scale(did,dsid,i)<0)
   goto out;
  if (H5DSis_attached(did,dsid,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,i,name)<0)
   goto out;
  if (H5DSget_label(did,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();


/*-------------------------------------------------------------------------
 * create a dataset and attach only to 1 dimension
 *-------------------------------------------------------------------------
 */  

 TESTING2("attach only to 1 dimension");
 
 /* make a dataset */
 if (H5LTmake_dataset_int(fid,"dset_b",rank,dims,NULL)<0)
  goto out;

 if ((did = H5Dopen(fid,"dset_b"))<0)
  goto out;

 /* attach a DS to dimension 1 */
 sprintf(name,"ds_a_%d",1);
 if((dsid = H5Dopen(fid,name))<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,i)<0)
    goto out;
  }
  else
  {
   if(H5DSdetach_scale(did,dsid,i)!=FAIL)
    goto out;
  }
 }

 if (H5Dclose(dsid)<0)
  goto out;
 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;
}