mirror of
https://github.com/HDFGroup/hdf5.git
synced 2024-11-27 02:10:55 +08:00
cc5b2ff93b
Code cleanup
Description:
Big clean up on high-level library code, including:
- Removing include of HDF5 private header file(!)
- Reduce number of compiler warnings
- Untangle header files, etc.
Platforms tested:
FreeBSD 4.11 (sleipnir) w/C++
Linux 2.4/64 (mir) w/C++ & FORTRAN
2871 lines
76 KiB
C
2871 lines
76 KiB
C
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
|
|
* 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://hdf.ncsa.uiuc.edu/HDF5/doc/Copyright.html. If you do not have *
|
|
* access to either file, you may request a copy from hdfhelp@ncsa.uiuc.edu. *
|
|
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
|
|
|
|
#include <stdlib.h>
|
|
#include <string.h>
|
|
#include "h5hltest.h"
|
|
#include "H5DSpublic.h"
|
|
#include "H5LTpublic.h"
|
|
#include "H5IMpublic.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);
|
|
static herr_t op_bogus(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);
|
|
static int test_types(void);
|
|
static int test_iterators(void);
|
|
static int test_data(void);
|
|
static int read_data( const char* fname, int ndims, hsize_t *dims, float **buf );
|
|
|
|
|
|
#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"
|
|
|
|
#define FILE1 "test_ds1.h5"
|
|
#define FILE2 "test_ds2.h5"
|
|
#define FILE3 "test_ds3.h5"
|
|
#define FILE4 "test_ds4.h5"
|
|
#define FILE5 "test_ds5.h5"
|
|
#define FILE6 "test_ds6.h5"
|
|
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* 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;
|
|
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;
|
|
|
|
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
|
|
* H5DSiterate_scales
|
|
* H5DSget_num_scales
|
|
*
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
|
|
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");
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* 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 = 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)==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 = 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",(size_t)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();
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* 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 = H5Dopen(fid,"dset_e"))<0)
|
|
goto out;
|
|
if ((dsid = H5Dopen(fid,"ds_e_1"))<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 = 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,(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 = 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,(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
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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_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 = H5Dopen(fid,"dset_f"))<0)
|
|
goto out;
|
|
|
|
if ((dsid = H5Dopen(fid,"ds_f_1"))<0)
|
|
goto out;
|
|
if (H5DSattach_scale(did,dsid,DIM0)<0)
|
|
goto out;
|
|
if (H5Dclose(dsid)<0)
|
|
goto out;
|
|
if ((dsid = H5Dopen(fid,"ds_f_11"))<0)
|
|
goto out;
|
|
if (H5DSattach_scale(did,dsid,DIM0)<0)
|
|
goto out;
|
|
if (H5Dclose(dsid)<0)
|
|
goto out;
|
|
if ((dsid = H5Dopen(fid,"ds_f_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_f" */
|
|
if ((did = H5Dopen(fid,"dset_f"))<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 = -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;
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* 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 */
|
|
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 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;
|
|
|
|
/* 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 */
|
|
storage_size=H5Dget_storage_size(dsid);
|
|
|
|
if (storage_size==0)
|
|
ret = 0;
|
|
|
|
return ret;
|
|
|
|
out:
|
|
H5E_BEGIN_TRY {
|
|
H5Sclose(sid);
|
|
} H5E_END_TRY;
|
|
return FAIL;
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* Function: op_bogus
|
|
*
|
|
* Purpose: example operator function used by H5DSiterate_scales, that does nothing
|
|
*
|
|
* 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 op_bogus(hid_t dset, unsigned dim, hid_t scale_id, void *visitor_data)
|
|
{
|
|
/* Stop compiler from whining about "unused parameters" */
|
|
dset = dset;
|
|
dim = dim;
|
|
scale_id = scale_id;
|
|
visitor_data = visitor_data;
|
|
|
|
/* define a default zero value for return. This will cause the iterator to continue */
|
|
return 0;
|
|
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* 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 */
|
|
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=H5Gcreate(fid,"grp",(size_t)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)==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=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)==SUCCEED)
|
|
goto out;
|
|
|
|
/* close */
|
|
if (H5Dclose(dsid)<0)
|
|
goto out;
|
|
if (H5Dclose(did)<0)
|
|
goto out;
|
|
|
|
/* open the previous written "ds_a" */
|
|
if ((did = H5Dopen(fid,"ds_a"))<0)
|
|
goto out;
|
|
|
|
/* open the previous written "ds_b" */
|
|
if ((dsid = H5Dopen(fid,"ds_b"))<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;
|
|
|
|
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)==SUCCEED)
|
|
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",(hsize_t)100,(hsize_t)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)==SUCCEED)
|
|
goto out;
|
|
|
|
/* close */
|
|
if (H5Dclose(dsid)<0)
|
|
goto out;
|
|
if (H5Dclose(did)<0)
|
|
goto out;
|
|
|
|
PASSED();
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* is scale
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
TESTING2("is scale");
|
|
|
|
/* open a non scale dataset */
|
|
if ((did = H5Dopen(fid,"dset_a"))<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 = H5Gopen(fid,"grp"))<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 = H5Dopen(fid,"ds_a"))<0)
|
|
goto out;
|
|
|
|
/* open the previous written "dset_a" */
|
|
if ((did = H5Dopen(fid,"dset_a"))<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 = H5Dopen(fid,"ds_a"))<0)
|
|
goto out;
|
|
|
|
/* open the group. */
|
|
if ((gid = H5Gopen(fid,"grp"))<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
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
TESTING2("detach scale when scale is group");
|
|
|
|
/* open the previous written "dset_a" */
|
|
if ((did = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
/* open the group. */
|
|
if ((gid = H5Gopen(fid,"grp"))<0)
|
|
goto out;
|
|
|
|
/* 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 */
|
|
out:
|
|
H5E_BEGIN_TRY {
|
|
H5Sclose(sid);
|
|
H5Sclose(sidds);
|
|
H5Dclose(did);
|
|
H5Dclose(dsid);
|
|
H5Gclose(gid);
|
|
H5Fclose(fid);
|
|
} H5E_END_TRY;
|
|
H5_FAILED();
|
|
return FAIL;
|
|
}
|
|
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* test iterators
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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=H5Gcreate(fid,"grp",(size_t)0))<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 = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
/* try to iterate trough the 1st dimension of "dset_a", return error */
|
|
if (H5DSiterate_scales(did,0,NULL,verifiy_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 = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
/* try to iterate trough the 3rd dimension of "dset_a", return error */
|
|
if (H5DSiterate_scales(did,3,NULL,verifiy_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 = H5Dopen(fid,"dset_a"))<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 = H5Dopen(fid,dname))<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();
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* iterate on group
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
TESTING2("iterate on group ");
|
|
|
|
/* open */
|
|
if ((gid = H5Gopen(fid,"grp"))<0)
|
|
goto out;
|
|
|
|
/* try to iterate, return error */
|
|
if (H5DSiterate_scales(gid,0,NULL,verifiy_scale,NULL)==SUCCEED)
|
|
goto out;
|
|
|
|
/* close */
|
|
if (H5Gclose(gid)<0)
|
|
goto out;
|
|
|
|
PASSED();
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* iterate in deleted scales
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
TESTING2("iterate in deleted scales ");
|
|
|
|
if (H5Gunlink(fid,"ds_0")<0)
|
|
goto out;
|
|
|
|
/* open the previously written "dset_a" */
|
|
if ((did = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
/* iterate */
|
|
if (H5DSiterate_scales(did,0,NULL,op_bogus,NULL)==SUCCEED)
|
|
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 {
|
|
H5Gclose(gid);
|
|
H5Fclose(fid);
|
|
} H5E_END_TRY;
|
|
H5_FAILED();
|
|
return FAIL;
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* test several ranks
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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 = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
for (i=0; i<3; i++)
|
|
{
|
|
sprintf(name,"ds_a_%d",i);
|
|
if((dsid = H5Dopen(fid,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 (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<3; i++)
|
|
{
|
|
sprintf(name,"ds_a_%d",i);
|
|
if((dsid = H5Dopen(fid,name))<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 = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
for (i=0; i<3; 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,(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=H5Dcreate(fid,"dset_b",H5T_NATIVE_INT,sid,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=H5Dcreate(fid,"ds_b_1",H5T_NATIVE_FLOAT,sidds,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;
|
|
}
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* attach scales with several datatypes
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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 = H5Dopen(fid,"dset_a"))<0)
|
|
goto out;
|
|
|
|
/* 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;
|
|
/* 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 = H5Dopen(fid,DS_2_NAME))<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 = H5Dopen(fid,"dset_b"))<0)
|
|
goto out;
|
|
/* get the DS dataset id */
|
|
if ((dsid = H5Dopen(fid,"ds_b_1"))<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 = H5Dopen(fid,"ds_b_2"))<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;
|
|
}
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* read realistic data and generate an HDF5 file with dimension scales
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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=H5Dcreate(fid,"data",H5T_NATIVE_FLOAT,sid,dcpl))<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 = H5Dopen(fid,"data"))<0)
|
|
goto out;
|
|
|
|
/* get the DS dataset id */
|
|
if ((dsid = H5Dopen(fid,"lat"))<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 = H5Dopen(fid,"lon"))<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;
|
|
}
|
|
|
|
|
|
|
|
/*-------------------------------------------------------------------------
|
|
* read_data
|
|
* utility function to read ASCII data
|
|
* the files have a header of the type
|
|
*
|
|
* dimension i
|
|
* n
|
|
*
|
|
* followed by the data
|
|
*
|
|
*-------------------------------------------------------------------------
|
|
*/
|
|
|
|
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;
|
|
|
|
}
|
|
|
|
|