hdf5/test/th5s.c

1628 lines
62 KiB
C
Raw Normal View History

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* Copyright by The HDF Group. *
* Copyright by the Board of Trustees of the University of Illinois. *
* All rights reserved. *
* *
* This file is part of HDF5. The full HDF5 copyright notice, including *
* terms governing use, modification, and redistribution, is contained in *
* the files COPYING and Copyright.html. COPYING can be found at the root *
* of the source code distribution tree; Copyright.html can be found at the *
* root level of an installed copy of the electronic HDF5 document set and *
* is linked from the top-level documents page. It can also be found at *
* http://hdfgroup.org/HDF5/doc/Copyright.html. If you do not have *
* access to either file, you may request a copy from help@hdfgroup.org. *
* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/***********************************************************
*
* Test program: th5s
*
* Test the dataspace functionality
*
*************************************************************/
#include "testhdf5.h"
#include "H5private.h"
#include "H5Bprivate.h"
#include "H5Sprivate.h"
#include "H5Pprivate.h"
#define TESTFILE "th5s.h5"
#define DATAFILE "th5s1.h5"
#define NULLFILE "th5s2.h5"
#define BASICFILE "th5s3.h5"
#define BASICDATASET "basic_dataset"
#define BASICDATASET2 "basic_dataset2"
#define BASICATTR "basic_attribute"
#define NULLDATASET "null_dataset"
#define NULLATTR "null_attribute"
/* 3-D dataset with fixed dimensions */
#define SPACE1_NAME "Space1"
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
#define SPACE1_RANK 3
#define SPACE1_DIM1 3
#define SPACE1_DIM2 15
#define SPACE1_DIM3 13
/* 4-D dataset with one unlimited dimension */
#define SPACE2_NAME "Space2"
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
#define SPACE2_RANK 4
#define SPACE2_DIM1 0
#define SPACE2_DIM2 15
#define SPACE2_DIM3 13
#define SPACE2_DIM4 23
#define SPACE2_MAX1 H5S_UNLIMITED
#define SPACE2_MAX2 15
#define SPACE2_MAX3 13
#define SPACE2_MAX4 23
/* Scalar dataset with simple datatype */
#define SPACE3_NAME "Scalar1"
#define SPACE3_RANK 0
unsigned space3_data=65;
/* Scalar dataset with compound datatype */
#define SPACE4_NAME "Scalar2"
#define SPACE4_RANK 0
#define SPACE4_FIELDNAME1 "c1"
#define SPACE4_FIELDNAME2 "u"
#define SPACE4_FIELDNAME3 "f"
#define SPACE4_FIELDNAME4 "c2"
size_t space4_field1_off=0;
size_t space4_field2_off=0;
size_t space4_field3_off=0;
size_t space4_field4_off=0;
struct space4_struct {
char c1;
unsigned u;
float f;
char c2;
} space4_data={'v',987123,(float)-3.14,'g'}; /* Test data for 4th dataspace */
/****************************************************************
**
** test_h5s_basic(): Test basic H5S (dataspace) code.
**
****************************************************************/
static void
test_h5s_basic(void)
{
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
hid_t fid1; /* HDF5 File IDs */
hid_t sid1, sid2; /* Dataspace ID */
hid_t dset1; /* Dataset ID */
hid_t aid1; /* Attribute ID */
int rank; /* Logical rank of dataspace */
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
hsize_t dims1[] = {SPACE1_DIM1, SPACE1_DIM2, SPACE1_DIM3};
hsize_t dims2[] = {SPACE2_DIM1, SPACE2_DIM2, SPACE2_DIM3,
SPACE2_DIM4};
hsize_t dims3[H5S_MAX_RANK+1];
hsize_t max2[] = {SPACE2_MAX1, SPACE2_MAX2, SPACE2_MAX3,
SPACE2_MAX4};
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
hsize_t tdims[4]; /* Dimension array to test with */
hsize_t tmax[4];
hssize_t n; /* Number of dataspace elements */
[svn-r338] Changes since 19980407 ---------------------- ./src/H5B.c ./src/H5D.c ./src/H5Dprivate.h ./src/H5Dpublic.h ./src/H5F.c ./src/H5Farray.c ./src/H5Fcore.c ./src/H5Ffamily.c ./src/H5Fistore.c ./src/H5Flow.c ./src/H5Fprivate.h ./src/H5Fpublic.h ./src/H5Fsec2.c ./src/H5Fsplit.c ./src/H5Fstdio.c ./src/H5G.c ./src/H5Gent.c ./src/H5Gnode.c ./src/H5HG.c ./src/H5HL.c ./src/H5MF.c ./src/H5MFprivate.h ./src/H5O.c ./src/H5Ocont.c ./src/H5Odtype.c ./src/H5Oefl.c ./src/H5Olayout.c ./src/H5Oname.c ./src/H5Oprivate.h ./src/H5Osdspace.c ./src/H5Oshared.c ./src/H5Ostab.c ./src/H5P.c ./src/H5Ppublic.h ./src/H5S.c ./src/H5Sprivate.h ./src/H5Spublic.h ./src/H5Ssimp.c ./src/H5Tconv.c ./src/H5Tpkg.h ./src/H5V.c ./src/H5Vprivate.h ./src/H5private.h ./src/H5public.h ./src/h5ls.c ./test/cmpd_dset.c ./test/dsets.c ./test/extend.c ./test/external.c ./test/hyperslab.c ./test/iopipe.c ./test/istore.c ./test/shtype.c ./test/tfile.c ./test/th5s.c Anything having to do with the size of a dataset now uses the types `hsize_t' and `hssize_t' which must be the same size and at least as large as `size_t'. This isn't fully tested yet, so hsize_t and hssize_t are defined as size_t and ssize_t in H5public.h. Setting them to larger values will trip up gcc versions less than 2.8.1 on x86 platforms. Documented unused function formals with `__unused__' before the formal name. This also has the effect of supressing warning messages for gcc since it's defined to be `__attribute__((unused))' in the H5private.h file. ./src/debug.c ./src/h5ls.c If the file name contains a `%' then the file is opened as a file family with H5P_DEFAULT for the file member access property list. ./src/h5ls.c The group name is optional, defaulting to `/'. ./src/hdf5.h Added some missing public header files.
1998-04-09 05:43:02 +08:00
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Dataspace Manipulation\n"));
sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL);
CHECK(sid1, FAIL, "H5Screate_simple");
n = H5Sget_simple_extent_npoints(sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, SPACE1_DIM1 * SPACE1_DIM2 * SPACE1_DIM3,
"H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE1_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid1, tdims, NULL);
CHECK(rank, FAIL, "H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tdims, dims1, SPACE1_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
sid2 = H5Screate_simple(SPACE2_RANK, dims2, max2);
CHECK(sid2, FAIL, "H5Screate_simple");
n = H5Sget_simple_extent_npoints(sid2);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, SPACE2_DIM1 * SPACE2_DIM2 * SPACE2_DIM3 * SPACE2_DIM4,
"H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid2);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE2_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid2, tdims, tmax);
CHECK(rank, FAIL, "H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tdims, dims2, SPACE2_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tmax, max2, SPACE2_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
/* Change max dims from zero to non-zero and back again */
ret = H5Sset_extent_simple(sid1, SPACE1_RANK, dims1, max2);
CHECK(ret, FAIL, "H5Sset_extent_simple");
ret = H5Sset_extent_simple(sid1, SPACE1_RANK, dims1, NULL);
CHECK(ret, FAIL, "H5Sset_extent_simple");
rank = H5Sget_simple_extent_dims(sid1, tdims, tmax);
CHECK(rank, FAIL, "H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tdims, dims1, SPACE1_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tmax, dims1, SPACE1_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
1998-03-13 05:04:53 +08:00
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
1998-03-13 05:04:53 +08:00
ret = H5Sclose(sid2);
CHECK(ret, FAIL, "H5Sclose");
/*
* Check to be sure we can't create a simple data space that has too many
* dimensions.
*/
H5E_BEGIN_TRY {
sid1 = H5Screate_simple(H5S_MAX_RANK+1, dims3, NULL);
} H5E_END_TRY;
VERIFY(sid1, FAIL, "H5Screate_simple");
/*
* Try reading a file that has been prepared that has a dataset with a
* higher dimensionality than what the library can handle.
*
* If this test fails and the H5S_MAX_RANK variable has changed, follow
* the instructions in space_overflow.c for regenerating the th5s.h5 file.
*/
{
char testfile[512]="";
2002-11-20 20:52:49 +08:00
char *srcdir = HDgetenv("srcdir");
if (srcdir && ((HDstrlen(srcdir) + HDstrlen(TESTFILE) + 1) < sizeof(testfile))){
2002-11-20 20:52:49 +08:00
HDstrcpy(testfile, srcdir);
HDstrcat(testfile, "/");
}
2002-11-20 20:52:49 +08:00
HDstrcat(testfile, TESTFILE);
fid1 = H5Fopen(testfile, H5F_ACC_RDONLY, H5P_DEFAULT);
CHECK_I(fid1, "H5Fopen");
if (fid1 >= 0){
dset1 = H5Dopen2(fid1, "dset", H5P_DEFAULT);
VERIFY(dset1, FAIL, "H5Dopen2");
ret = H5Fclose(fid1);
CHECK_I(ret, "H5Fclose");
}
else
printf("***cannot open the pre-created H5S_MAX_RANK test file (%s)\n",
testfile);
}
/* Verify that incorrect dimensions don't work */
dims1[0]=0;
sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL);
VERIFY(sid1, FAIL, "H5Screate_simple");
dims1[0] = H5S_UNLIMITED;
sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL);
VERIFY(sid1, FAIL, "H5Screate_simple");
dims1[0]=0;
sid1 = H5Screate(H5S_SIMPLE);
CHECK(sid1, FAIL, "H5Screate");
ret = H5Sset_extent_simple(sid1,SPACE1_RANK,dims1,NULL);
VERIFY(ret, FAIL, "H5Sset_extent_simple");
ret = H5Sclose(sid1);
CHECK_I(ret, "H5Sclose");
dims1[0] = H5S_UNLIMITED;
sid1 = H5Screate(H5S_SIMPLE);
CHECK(sid1, FAIL, "H5Screate");
ret = H5Sset_extent_simple(sid1,SPACE1_RANK,dims1,NULL);
VERIFY(ret, FAIL, "H5Sset_extent_simple");
ret = H5Sclose(sid1);
CHECK_I(ret, "H5Sclose");
/*
* Try writing simple dataspaces without setting their extents
*/
/* Create the file */
fid1 = H5Fcreate(BASICFILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
CHECK(fid1, FAIL, "H5Fcreate");
dims1[0]=SPACE1_DIM1;
sid1 = H5Screate(H5S_SIMPLE);
CHECK(sid1, FAIL, "H5Screate");
sid2 = H5Screate_simple(1, dims1, dims1);
CHECK(sid2, FAIL, "H5Screate");
/* This dataset's space has no extent; it should not be created */
H5E_BEGIN_TRY {
dset1 = H5Dcreate2(fid1, BASICDATASET, H5T_NATIVE_INT, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
} H5E_END_TRY
VERIFY(dset1, FAIL, "H5Dcreate2");
dset1 = H5Dcreate2(fid1, BASICDATASET2, H5T_NATIVE_INT, sid2, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
CHECK(dset1, FAIL, "H5Dcreate2");
/* Try some writes with the bad dataspace (sid1) */
H5E_BEGIN_TRY {
ret = H5Dwrite(dset1, H5T_NATIVE_INT, sid1, H5S_ALL, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dwrite");
H5E_BEGIN_TRY {
ret = H5Dwrite(dset1, H5T_NATIVE_INT, H5S_ALL, sid1, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dwrite");
H5E_BEGIN_TRY {
ret = H5Dwrite(dset1, H5T_NATIVE_INT, sid1, sid1, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dwrite");
/* Try to iterate using the bad dataspace */
H5E_BEGIN_TRY {
ret = H5Diterate(&n, H5T_NATIVE_INT, sid1, NULL, NULL);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Diterate");
/* Try to fill using the bad dataspace */
H5E_BEGIN_TRY {
ret = H5Dfill(NULL, H5T_NATIVE_INT, &n, H5T_NATIVE_INT, sid1);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dfill");
/* Now use the bad dataspace as the space for an attribute */
H5E_BEGIN_TRY {
aid1 = H5Acreate2(dset1, BASICATTR, H5T_NATIVE_INT, sid1, H5P_DEFAULT, H5P_DEFAULT);
} H5E_END_TRY
VERIFY(aid1, FAIL, "H5Acreate2");
/* Make sure that dataspace reads using the bad dataspace fail */
H5E_BEGIN_TRY {
ret = H5Dread(dset1, H5T_NATIVE_INT, sid1, H5S_ALL, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dread");
H5E_BEGIN_TRY {
ret = H5Dread(dset1, H5T_NATIVE_INT, H5S_ALL, sid1, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dread");
H5E_BEGIN_TRY {
ret = H5Dread(dset1, H5T_NATIVE_INT, sid1, sid1, H5P_DEFAULT, &n);
} H5E_END_TRY
VERIFY(ret, FAIL, "H5Dread");
/* Clean up */
ret = H5Dclose(dset1);
CHECK(ret, FAIL, "H5Dclose");
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(sid2);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Fclose(fid1);
CHECK(ret, FAIL, "H5Fclose");
} /* test_h5s_basic() */
/****************************************************************
**
** test_h5s_null(): Test NULL data space
**
****************************************************************/
static void
test_h5s_null(void)
{
hid_t fid; /* File ID */
hid_t sid; /* Dataspace IDs */
hid_t dset_sid, dset_sid2; /* Dataspace IDs */
hid_t attr_sid; /* Dataspace IDs */
hid_t did; /* Dataset ID */
hid_t attr; /*Attribute ID */
H5S_class_t stype; /* dataspace type */
hssize_t nelem; /* Number of elements */
unsigned uval=2; /* Buffer for writing to dataset */
int val=1; /* Buffer for writing to attribute */
H5S_sel_type sel_type; /* Type of selection currently */
hsize_t dims[1]={10}; /* Dimensions for converting null dataspace to simple */
H5S_class_t space_type; /* Type of dataspace */
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Null Dataspace\n"));
/* Create the file */
fid = H5Fcreate(NULLFILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
CHECK(fid, FAIL, "H5Fcreate");
sid = H5Screate(H5S_NULL);
CHECK(sid, FAIL, "H5Screate");
/* Check that the null dataspace actually has 0 elements */
nelem = H5Sget_simple_extent_npoints(sid);
VERIFY(nelem, 0, "H5Sget_simple_extent_npoints");
/* Check that the dataspace was created with an "all" selection */
sel_type = H5Sget_select_type(sid);
VERIFY(sel_type, H5S_SEL_ALL, "H5Sget_select_type");
/* Check that the null dataspace has 0 elements selected */
nelem = H5Sget_select_npoints(sid);
VERIFY(nelem, 0, "H5Sget_select_npoints");
/* Change to "none" selection */
ret = H5Sselect_none(sid);
CHECK(ret, FAIL, "H5Sselect_none");
/* Check that the null dataspace has 0 elements selected */
nelem = H5Sget_select_npoints(sid);
VERIFY(nelem, 0, "H5Sget_select_npoints");
/* Check to be sure we can't set a hyperslab selection on a null dataspace */
H5E_BEGIN_TRY {
hsize_t start[1]={0};
hsize_t count[1]={0};
ret = H5Sselect_hyperslab(sid, H5S_SELECT_SET, start, NULL, count, NULL);
} H5E_END_TRY;
VERIFY(ret, FAIL, "H5Sselect_hyperslab");
/* Check to be sure we can't set a point selection on a null dataspace */
H5E_BEGIN_TRY {
hsize_t coord[1][1]; /* Coordinates for point selection */
coord[0][0]=0;
ret = H5Sselect_elements(sid, H5S_SELECT_SET, (size_t)1, (const hsize_t *)coord);
} H5E_END_TRY;
VERIFY(ret, FAIL, "H5Sselect_elements");
/* Create first dataset */
did = H5Dcreate2(fid, NULLDATASET, H5T_NATIVE_UINT, sid, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
CHECK(did, FAIL, "H5Dcreate2");
/* Write "nothing" to the dataset */
ret = H5Dwrite(did, H5T_NATIVE_UINT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &uval);
CHECK(ret, FAIL, "H5Dwrite");
/* Write "nothing" to the dataset (with type conversion :-) */
ret = H5Dwrite(did, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &val);
CHECK(ret, FAIL, "H5Dwrite");
/* Try reading from the dataset (make certain our buffer is unmodified) */
ret = H5Dread(did, H5T_NATIVE_UINT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &uval);
CHECK(ret, FAIL, "H5Dread");
VERIFY(uval, 2, "H5Dread");
/* Try reading from the dataset (with type conversion :-) (make certain our buffer is unmodified) */
ret = H5Dread(did, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &val);
CHECK(ret, FAIL, "H5Dread");
VERIFY(val, 1, "H5Dread");
/* Create an attribute for the group */
attr = H5Acreate2(did, NULLATTR, H5T_NATIVE_INT, sid, H5P_DEFAULT, H5P_DEFAULT);
CHECK(attr, FAIL, "H5Acreate2");
/* Write "nothing" to the attribute */
ret = H5Awrite(attr, H5T_NATIVE_INT, &val);
CHECK(ret, FAIL, "H5Awrite");
/* Write "nothing" to the attribute (with type conversion :-) */
ret = H5Awrite(attr, H5T_NATIVE_UINT, &uval);
CHECK(ret, FAIL, "H5Awrite");
/* Try reading from the attribute (make certain our buffer is unmodified) */
ret = H5Aread(attr, H5T_NATIVE_INT, &val);
CHECK(ret, FAIL, "H5Aread");
VERIFY(val, 1, "H5Aread");
/* Try reading from the attribute (with type conversion :-) (make certain our buffer is unmodified) */
ret = H5Aread(attr, H5T_NATIVE_UINT, &uval);
CHECK(ret, FAIL, "H5Aread");
VERIFY(uval, 2, "H5Aread");
/* Close attribute */
ret=H5Aclose(attr);
CHECK(ret, FAIL, "H5Aclose");
/* Close the dataset */
ret = H5Dclose(did);
CHECK(ret, FAIL, "H5Dclose");
/* Verify that we've got the right kind of dataspace */
space_type = H5Sget_simple_extent_type(sid);
VERIFY(space_type, H5S_NULL, "H5Sget_simple_extent_type");
/* Convert the null dataspace to a simple dataspace */
ret = H5Sset_extent_simple(sid, 1, dims, NULL);
CHECK(ret, FAIL, "H5Sset_extent_simple");
/* Verify that we've got the right kind of dataspace now */
space_type = H5Sget_simple_extent_type(sid);
VERIFY(space_type, H5S_SIMPLE, "H5Sget_simple_extent_type");
/* Close the dataspace */
ret = H5Sclose(sid);
CHECK(ret, FAIL, "H5Sclose");
/* Close the file */
ret = H5Fclose(fid);
CHECK(ret, FAIL, "H5Fclose");
/*============================================
* Reopen the file to check the data space
*============================================
*/
fid = H5Fopen(NULLFILE, H5F_ACC_RDONLY, H5P_DEFAULT);
CHECK(fid, FAIL, "H5Fopen");
/* Reopen the dataset */
did = H5Dopen2(fid, NULLDATASET, H5P_DEFAULT);
CHECK(did, FAIL, "H5Dopen2");
/* Get the space of the dataset */
dset_sid = H5Dget_space(did);
CHECK(dset_sid, FAIL, "H5Dget_space");
/* Query the NULL dataspace */
dset_sid2 = H5Scopy(dset_sid);
CHECK(dset_sid2, FAIL, "H5Scopy");
/* Verify the class type of dataspace */
stype = H5Sget_simple_extent_type(dset_sid2);
VERIFY(stype, H5S_NULL, "H5Sget_simple_extent_type");
/* Verify there is zero element in the dataspace */
ret = (herr_t)H5Sget_simple_extent_npoints(dset_sid2);
VERIFY(ret, 0, "H5Sget_simple_extent_npoints");
/* Try reading from the dataset (make certain our buffer is unmodified) */
ret = H5Dread(did, H5T_NATIVE_UINT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &uval);
CHECK(ret, FAIL, "H5Dread");
VERIFY(uval, 2, "H5Dread");
/* Close the dataspace */
ret = H5Sclose(dset_sid);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(dset_sid2);
CHECK(ret, FAIL, "H5Sclose");
/* Open the attribute for the dataset */
attr = H5Aopen(did, NULLATTR, H5P_DEFAULT);
CHECK(attr, FAIL, "H5Aopen");
/* Get the space of the dataset */
attr_sid = H5Aget_space(attr);
CHECK(attr_sid, FAIL, "H5Aget_space");
/* Verify the class type of dataspace */
stype = H5Sget_simple_extent_type(attr_sid);
VERIFY(stype, H5S_NULL, "H5Sget_simple_extent_type");
/* Verify there is zero element in the dataspace */
ret = (herr_t)H5Sget_simple_extent_npoints(attr_sid);
VERIFY(ret, 0, "H5Sget_simple_extent_npoints");
/* Close the dataspace */
ret = H5Sclose(attr_sid);
CHECK(ret, FAIL, "H5Sclose");
/* Try reading from the attribute (make certain our buffer is unmodified) */
ret = H5Aread(attr, H5T_NATIVE_INT, &val);
CHECK(ret, FAIL, "H5Aread");
VERIFY(val, 1, "H5Aread");
/* Close attribute */
ret=H5Aclose(attr);
CHECK(ret, FAIL, "H5Aclose");
/* Close the dataset */
ret = H5Dclose(did);
CHECK(ret, FAIL, "H5Dclose");
/* Close the file */
ret = H5Fclose(fid);
CHECK(ret, FAIL, "H5Fclose");
} /* end test_h5s_null() */
/****************************************************************
**
** test_h5s_encode(): Test H5S (dataspace) encoding and decoding.
**
****************************************************************/
static void
test_h5s_encode(void)
{
hid_t sid1, sid2, sid3; /* Dataspace ID */
hid_t decoded_sid1, decoded_sid2, decoded_sid3;
int rank; /* Logical rank of dataspace */
hsize_t dims1[] = {SPACE1_DIM1, SPACE1_DIM2, SPACE1_DIM3};
size_t sbuf_size=0, null_size=0, scalar_size=0;
unsigned char *sbuf=NULL, *null_sbuf=NULL, *scalar_buf=NULL;
hsize_t tdims[4]; /* Dimension array to test with */
hssize_t n; /* Number of dataspace elements */
hsize_t start[] = {0, 0, 0};
hsize_t stride[] = {2, 5, 3};
hsize_t count[] = {2, 2, 2};
hsize_t block[] = {1, 3, 1};
H5S_sel_type sel_type;
H5S_class_t space_type;
hssize_t nblocks;
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Dataspace Encoding and Decoding\n"));
/*-------------------------------------------------------------------------
* Test encoding and decoding of simple dataspace and hyperslab selection.
*-------------------------------------------------------------------------
*/
sid1 = H5Screate_simple(SPACE1_RANK, dims1, NULL);
CHECK(sid1, FAIL, "H5Screate_simple");
ret = H5Sselect_hyperslab(sid1, H5S_SELECT_SET, start, stride, count, block);
CHECK(ret, FAIL, "H5Sselect_hyperslab");
/* Encode simple data space in a buffer */
ret = H5Sencode(sid1, NULL, &sbuf_size);
CHECK(ret, FAIL, "H5Sencode");
if(sbuf_size>0)
sbuf = (unsigned char*)HDcalloc((size_t)1, sbuf_size);
/* Try decoding bogus buffer */
H5E_BEGIN_TRY {
ret = H5Sdecode(sbuf);
} H5E_END_TRY;
VERIFY(ret, FAIL, "H5Sdecode");
ret = H5Sencode(sid1, sbuf, &sbuf_size);
CHECK(ret, FAIL, "H5Sencode");
/* Decode from the dataspace buffer and return an object handle */
decoded_sid1=H5Sdecode(sbuf);
CHECK(decoded_sid1, FAIL, "H5Sdecode");
/* Verify the decoded dataspace */
n = H5Sget_simple_extent_npoints(decoded_sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, SPACE1_DIM1 * SPACE1_DIM2 * SPACE1_DIM3,
"H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(decoded_sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE1_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(decoded_sid1, tdims, NULL);
CHECK(rank, FAIL, "H5Sget_simple_extent_dims");
VERIFY(HDmemcmp(tdims, dims1, SPACE1_RANK * sizeof(hsize_t)), 0,
"H5Sget_simple_extent_dims");
/* Verify hyperslabe selection */
sel_type = H5Sget_select_type(decoded_sid1);
VERIFY(sel_type, H5S_SEL_HYPERSLABS, "H5Sget_select_type");
nblocks = H5Sget_select_hyper_nblocks(decoded_sid1);
VERIFY(nblocks, 2*2*2, "H5Sget_select_hyper_nblocks");
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(decoded_sid1);
CHECK(ret, FAIL, "H5Sclose");
/*-------------------------------------------------------------------------
* Test encoding and decoding of null dataspace.
*-------------------------------------------------------------------------
*/
sid2 = H5Screate(H5S_NULL);
CHECK(sid2, FAIL, "H5Screate");
/* Encode null data space in a buffer */
ret = H5Sencode(sid2, NULL, &null_size);
CHECK(ret, FAIL, "H5Sencode");
if(null_size>0)
null_sbuf = (unsigned char*)HDcalloc((size_t)1, null_size);
ret = H5Sencode(sid2, null_sbuf, &null_size);
CHECK(ret, FAIL, "H5Sencode");
/* Decode from the dataspace buffer and return an object handle */
decoded_sid2=H5Sdecode(null_sbuf);
CHECK(decoded_sid2, FAIL, "H5Sdecode");
/* Verify decoded dataspace */
space_type = H5Sget_simple_extent_type(decoded_sid2);
VERIFY(space_type, H5S_NULL, "H5Sget_simple_extent_type");
ret = H5Sclose(sid2);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(decoded_sid2);
CHECK(ret, FAIL, "H5Sclose");
/*-------------------------------------------------------------------------
* Test encoding and decoding of scalar dataspace.
*-------------------------------------------------------------------------
*/
/* Create scalar dataspace */
sid3 = H5Screate(H5S_SCALAR);
CHECK(sid3, FAIL, "H5Screate_simple");
/* Encode scalar data space in a buffer */
ret = H5Sencode(sid3, NULL, &scalar_size);
CHECK(ret, FAIL, "H5Sencode");
if(scalar_size>0)
scalar_buf = (unsigned char*)HDcalloc((size_t)1, scalar_size);
ret = H5Sencode(sid3, scalar_buf, &scalar_size);
CHECK(ret, FAIL, "H5Sencode");
/* Decode from the dataspace buffer and return an object handle */
decoded_sid3=H5Sdecode(scalar_buf);
CHECK(decoded_sid3, FAIL, "H5Sdecode");
/* Verify extent type */
space_type = H5Sget_simple_extent_type(decoded_sid3);
VERIFY(space_type, H5S_SCALAR, "H5Sget_simple_extent_type");
/* Verify decoded dataspace */
n = H5Sget_simple_extent_npoints(decoded_sid3);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, 1, "H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(decoded_sid3);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, 0, "H5Sget_simple_extent_ndims");
ret = H5Sclose(sid3);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(decoded_sid3);
CHECK(ret, FAIL, "H5Sclose");
HDfree(sbuf);
HDfree(null_sbuf);
HDfree(scalar_buf);
} /* test_h5s_encode() */
/****************************************************************
**
** test_h5s_scalar_write(): Test scalar H5S (dataspace) writing code.
**
****************************************************************/
static void
test_h5s_scalar_write(void)
{
hid_t fid1; /* HDF5 File IDs */
hid_t dataset; /* Dataset ID */
hid_t sid1; /* Dataspace ID */
int rank; /* Logical rank of dataspace */
hsize_t tdims[4]; /* Dimension array to test with */
hssize_t n; /* Number of dataspace elements */
H5S_class_t ext_type; /* Extent type */
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Scalar Dataspace Manipulation during Writing\n"));
/* Create file */
fid1 = H5Fcreate(DATAFILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
CHECK(fid1, FAIL, "H5Fcreate");
/* Create scalar dataspace */
sid1 = H5Screate_simple(SPACE3_RANK, NULL, NULL);
CHECK(sid1, FAIL, "H5Screate_simple");
n = H5Sget_simple_extent_npoints(sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, 1, "H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE3_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid1, tdims, NULL);
VERIFY(rank, 0, "H5Sget_simple_extent_dims");
/* Verify extent type */
ext_type = H5Sget_simple_extent_type(sid1);
VERIFY(ext_type, H5S_SCALAR, "H5Sget_simple_extent_type");
/* Create a dataset */
dataset = H5Dcreate2(fid1, "Dataset1", H5T_NATIVE_UINT, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
CHECK(dataset, FAIL, "H5Dcreate2");
ret = H5Dwrite(dataset, H5T_NATIVE_UINT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &space3_data);
CHECK(ret, FAIL, "H5Dwrite");
/* Close Dataset */
ret = H5Dclose(dataset);
CHECK(ret, FAIL, "H5Dclose");
/* Close scalar dataspace */
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
/* Close file */
ret = H5Fclose(fid1);
CHECK(ret, FAIL, "H5Fclose");
} /* test_h5s_scalar_write() */
/****************************************************************
**
** test_h5s_scalar_read(): Test scalar H5S (dataspace) reading code.
**
****************************************************************/
static void
test_h5s_scalar_read(void)
{
hid_t fid1; /* HDF5 File IDs */
hid_t dataset; /* Dataset ID */
hid_t sid1; /* Dataspace ID */
int rank; /* Logical rank of dataspace */
hsize_t tdims[4]; /* Dimension array to test with */
hssize_t n; /* Number of dataspace elements */
unsigned rdata; /* Scalar data read in */
herr_t ret; /* Generic return value */
H5S_class_t ext_type; /* Extent type */
/* Output message about test being performed */
MESSAGE(5, ("Testing Scalar Dataspace Manipulation during Reading\n"));
/* Create file */
fid1 = H5Fopen(DATAFILE, H5F_ACC_RDWR, H5P_DEFAULT);
CHECK(fid1, FAIL, "H5Fopen");
/* Create a dataset */
dataset = H5Dopen2(fid1, "Dataset1", H5P_DEFAULT);
CHECK(dataset, FAIL, "H5Dopen2");
sid1 = H5Dget_space(dataset);
CHECK(sid1, FAIL, "H5Dget_space");
n = H5Sget_simple_extent_npoints(sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, 1, "H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE3_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid1, tdims, NULL);
VERIFY(rank, 0, "H5Sget_simple_extent_dims");
/* Verify extent type */
ext_type = H5Sget_simple_extent_type(sid1);
VERIFY(ext_type, H5S_SCALAR, "H5Sget_simple_extent_type");
ret = H5Dread(dataset, H5T_NATIVE_UINT, H5S_ALL, H5S_ALL, H5P_DEFAULT, &rdata);
CHECK(ret, FAIL, "H5Dread");
VERIFY(rdata, space3_data, "H5Dread");
/* Close Dataset */
ret = H5Dclose(dataset);
CHECK(ret, FAIL, "H5Dclose");
/* Close scalar dataspace */
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
/* Close file */
ret = H5Fclose(fid1);
CHECK(ret, FAIL, "H5Fclose");
} /* test_h5s_scalar_read() */
/****************************************************************
**
** test_h5s_compound_scalar_write(): Test scalar H5S (dataspace) writing for
** compound datatypes.
**
****************************************************************/
static void
test_h5s_compound_scalar_write(void)
{
hid_t fid1; /* HDF5 File IDs */
hid_t dataset; /* Dataset ID */
hid_t tid1; /* Attribute datatype ID */
hid_t sid1; /* Dataspace ID */
int rank; /* Logical rank of dataspace */
hsize_t tdims[4]; /* Dimension array to test with */
hssize_t n; /* Number of dataspace elements */
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Scalar Dataspace Manipulation for Writing Compound Datatypes\n"));
/* Create file */
fid1 = H5Fcreate(DATAFILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
CHECK(fid1, FAIL, "H5Fcreate");
/* Create the compound datatype. */
tid1 = H5Tcreate (H5T_COMPOUND, sizeof(struct space4_struct));
CHECK(tid1, FAIL, "H5Tcreate");
space4_field1_off=HOFFSET(struct space4_struct, c1);
ret = H5Tinsert(tid1, SPACE4_FIELDNAME1, space4_field1_off,
H5T_NATIVE_SCHAR);
CHECK(ret, FAIL, "H5Tinsert");
space4_field2_off=HOFFSET(struct space4_struct, u);
ret = H5Tinsert(tid1, SPACE4_FIELDNAME2, space4_field2_off,
H5T_NATIVE_UINT);
CHECK(ret, FAIL, "H5Tinsert");
space4_field3_off=HOFFSET(struct space4_struct, f);
ret = H5Tinsert(tid1, SPACE4_FIELDNAME3, space4_field3_off,
H5T_NATIVE_FLOAT);
CHECK(ret, FAIL, "H5Tinsert");
space4_field4_off=HOFFSET(struct space4_struct, c2);
ret = H5Tinsert(tid1, SPACE4_FIELDNAME4, space4_field4_off,
H5T_NATIVE_SCHAR);
CHECK(ret, FAIL, "H5Tinsert");
/* Create scalar dataspace */
sid1 = H5Screate_simple(SPACE3_RANK, NULL, NULL);
CHECK(sid1, FAIL, "H5Screate_simple");
n = H5Sget_simple_extent_npoints(sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, 1, "H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE3_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid1, tdims, NULL);
VERIFY(rank, 0, "H5Sget_simple_extent_dims");
/* Create a dataset */
dataset = H5Dcreate2(fid1, "Dataset1", tid1, sid1, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT);
CHECK(dataset, FAIL, "H5Dcreate2");
ret = H5Dwrite(dataset, tid1, H5S_ALL, H5S_ALL, H5P_DEFAULT, &space4_data);
CHECK(ret, FAIL, "H5Dwrite");
/* Close Dataset */
ret = H5Dclose(dataset);
CHECK(ret, FAIL, "H5Dclose");
/* Close scalar dataspace */
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
/* Close file */
ret = H5Fclose(fid1);
CHECK(ret, FAIL, "H5Fclose");
} /* test_h5s_compound_scalar_write() */
/****************************************************************
**
** test_h5s_compound_scalar_read(): Test scalar H5S (dataspace) reading for
** compound datatypes.
**
****************************************************************/
static void
test_h5s_compound_scalar_read(void)
{
hid_t fid1; /* HDF5 File IDs */
hid_t dataset; /* Dataset ID */
hid_t sid1; /* Dataspace ID */
hid_t type; /* Datatype */
int rank; /* Logical rank of dataspace */
hsize_t tdims[4]; /* Dimension array to test with */
hssize_t n; /* Number of dataspace elements */
struct space4_struct rdata; /* Scalar data read in */
herr_t ret; /* Generic return value */
/* Output message about test being performed */
MESSAGE(5, ("Testing Scalar Dataspace Manipulation for Reading Compound Datatypes\n"));
/* Create file */
fid1 = H5Fopen(DATAFILE, H5F_ACC_RDWR, H5P_DEFAULT);
CHECK(fid1, FAIL, "H5Fopen");
/* Create a dataset */
dataset = H5Dopen2(fid1, "Dataset1", H5P_DEFAULT);
CHECK(dataset, FAIL, "H5Dopen2");
sid1 = H5Dget_space(dataset);
CHECK(sid1, FAIL, "H5Dget_space");
n = H5Sget_simple_extent_npoints(sid1);
CHECK(n, FAIL, "H5Sget_simple_extent_npoints");
VERIFY(n, 1, "H5Sget_simple_extent_npoints");
rank = H5Sget_simple_extent_ndims(sid1);
CHECK(rank, FAIL, "H5Sget_simple_extent_ndims");
VERIFY(rank, SPACE3_RANK, "H5Sget_simple_extent_ndims");
rank = H5Sget_simple_extent_dims(sid1, tdims, NULL);
VERIFY(rank, 0, "H5Sget_simple_extent_dims");
type=H5Dget_type(dataset);
CHECK(type, FAIL, "H5Dget_type");
ret = H5Dread(dataset, type, H5S_ALL, H5S_ALL, H5P_DEFAULT, &rdata);
CHECK(ret, FAIL, "H5Dread");
if(HDmemcmp(&space4_data,&rdata,sizeof(struct space4_struct))) {
printf("scalar data different: space4_data.c1=%c, read_data4.c1=%c\n",space4_data.c1,rdata.c1);
printf("scalar data different: space4_data.u=%u, read_data4.u=%u\n",space4_data.u,rdata.u);
printf("scalar data different: space4_data.f=%f, read_data4.f=%f\n",space4_data.f,rdata.f);
TestErrPrintf("scalar data different: space4_data.c1=%c, read_data4.c1=%c\n",space4_data.c1,rdata.c2);
} /* end if */
/* Close datatype */
ret = H5Tclose(type);
CHECK(ret, FAIL, "H5Tclose");
/* Close Dataset */
ret = H5Dclose(dataset);
CHECK(ret, FAIL, "H5Dclose");
/* Close scalar dataspace */
ret = H5Sclose(sid1);
CHECK(ret, FAIL, "H5Sclose");
/* Close file */
ret = H5Fclose(fid1);
CHECK(ret, FAIL, "H5Fclose");
} /* test_h5s_compound_scalar_read() */
/* Data arrays for chunk test */
double chunk_data_dbl[50000][3];
float chunk_data_flt[50000][3];
/****************************************************************
**
** test_h5s_chunk(): Exercise chunked I/O, testing when data conversion
** is necessary and the entire chunk read in doesn't fit into the
** conversion buffer
**
****************************************************************/
static void
test_h5s_chunk(void)
{
herr_t status;
hid_t fileID, dsetID;
hid_t plist_id;
hid_t space_id;
hsize_t dims[2];
hsize_t csize[2];
int i,j;
fileID = H5Fcreate(DATAFILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
CHECK(fileID, FAIL, "H5Fcreate");
plist_id = H5Pcreate(H5P_DATASET_CREATE);
CHECK(plist_id, FAIL, "H5Pcreate");
csize[0] = 50000;
csize[1] = 3;
status = H5Pset_chunk(plist_id, 2, csize);
CHECK(status, FAIL, "H5Pset_chunk");
/* Create the data space */
dims[0] = 50000;
dims[1] = 3;
space_id = H5Screate_simple(2, dims, NULL);
CHECK(space_id, FAIL, "H5Screate_simple");
dsetID = H5Dcreate2(fileID, "coords", H5T_NATIVE_FLOAT, space_id, H5P_DEFAULT, plist_id, H5P_DEFAULT);
CHECK(dsetID, FAIL, "H5Dcreate2");
/* Initialize float array */
for(i = 0; i < 50000; i++)
for(j = 0; j < 3; j++)
chunk_data_flt[i][j] = (float)((i + 1) * 2.5 - j * 100.3);
status = H5Dwrite(dsetID, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, chunk_data_flt);
CHECK(status, FAIL, "H5Dwrite");
status = H5Pclose(plist_id);
CHECK(status, FAIL, "H5Pclose");
status = H5Sclose(space_id);
CHECK(status, FAIL, "H5Sclose");
status = H5Dclose(dsetID);
CHECK(status, FAIL, "H5Dclose");
status = H5Fclose(fileID);
CHECK(status, FAIL, "H5Fclose");
/* Reset/initialize the data arrays to read in */
HDmemset(chunk_data_dbl, 0, sizeof(double) * 50000 * 3);
HDmemset(chunk_data_flt, 0, sizeof(float) * 50000 * 3);
fileID = H5Fopen(DATAFILE, H5F_ACC_RDONLY, H5P_DEFAULT);
CHECK(fileID, FAIL, "H5Fopen");
dsetID = H5Dopen2(fileID, "coords", H5P_DEFAULT);
CHECK(dsetID, FAIL, "H5Dopen2");
status= H5Dread (dsetID, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, H5P_DEFAULT, chunk_data_dbl);
CHECK(status, FAIL, "H5Dread");
status= H5Dread (dsetID, H5T_NATIVE_FLOAT, H5S_ALL, H5S_ALL, H5P_DEFAULT, chunk_data_flt);
CHECK(status, FAIL, "H5Dread");
status = H5Dclose(dsetID);
CHECK(status, FAIL, "H5Dclose");
status = H5Fclose(fileID);
CHECK(status, FAIL, "H5Fclose");
for(i=0; i<50000; i++) {
for(j=0; j<3; j++) {
/* Check if the two values are within 0.001% range. */
if(!DBL_REL_EQUAL(chunk_data_dbl[i][j], chunk_data_flt[i][j], 0.00001))
TestErrPrintf("%u: chunk_data_dbl[%d][%d]=%e, chunk_data_flt[%d][%d]=%e\n", (unsigned)__LINE__, i, j, chunk_data_dbl[i][j], i, j, chunk_data_flt[i][j]);
} /* end for */
} /* end for */
} /* test_h5s_chunk() */
/****************************************************************
**
** test_h5s_extent_equal(): Exercise extent comparison code
**
****************************************************************/
static void
test_h5s_extent_equal(void)
{
hid_t null_space; /* Null dataspace */
hid_t scalar_space; /* Scalar dataspace */
hid_t d1_space1, d1_space2, d1_space3, d1_space4; /* 1-D dataspaces */
hid_t d2_space1, d2_space2, d2_space3, d2_space4; /* 2-D dataspaces */
hid_t d3_space1, d3_space2, d3_space3, d3_space4; /* 3-D dataspaces */
hsize_t d1_dims1[1] = {10}, /* 1-D dimensions */
d1_dims2[1] = {20},
d1_dims3[1] = {H5S_UNLIMITED};
hsize_t d2_dims1[2] = {10, 10}, /* 2-D dimensions */
d2_dims2[2] = {20, 20},
d2_dims3[2] = {H5S_UNLIMITED, H5S_UNLIMITED};
hsize_t d3_dims1[3] = {10, 10, 10}, /* 3-D dimensions */
d3_dims2[3] = {20, 20, 20},
d3_dims3[3] = {H5S_UNLIMITED, H5S_UNLIMITED, H5S_UNLIMITED};
htri_t ext_equal; /* Whether two dataspace extents are equal */
herr_t ret; /* Generic error return */
/* Create dataspaces */
null_space = H5Screate(H5S_NULL);
CHECK(null_space, FAIL, "H5Screate");
scalar_space = H5Screate(H5S_SCALAR);
CHECK(scalar_space, FAIL, "H5Screate");
d1_space1 = H5Screate_simple(1, d1_dims1, NULL);
CHECK(d1_space1, FAIL, "H5Screate");
d1_space2 = H5Screate_simple(1, d1_dims2, NULL);
CHECK(d1_space2, FAIL, "H5Screate");
d1_space3 = H5Screate_simple(1, d1_dims1, d1_dims2);
CHECK(d1_space3, FAIL, "H5Screate");
d1_space4 = H5Screate_simple(1, d1_dims1, d1_dims3);
CHECK(d1_space4, FAIL, "H5Screate");
d2_space1 = H5Screate_simple(2, d2_dims1, NULL);
CHECK(d2_space1, FAIL, "H5Screate");
d2_space2 = H5Screate_simple(2, d2_dims2, NULL);
CHECK(d2_space2, FAIL, "H5Screate");
d2_space3 = H5Screate_simple(2, d2_dims1, d2_dims2);
CHECK(d2_space3, FAIL, "H5Screate");
d2_space4 = H5Screate_simple(2, d2_dims1, d2_dims3);
CHECK(d2_space4, FAIL, "H5Screate");
d3_space1 = H5Screate_simple(3, d3_dims1, NULL);
CHECK(d3_space1, FAIL, "H5Screate");
d3_space2 = H5Screate_simple(3, d3_dims2, NULL);
CHECK(d3_space2, FAIL, "H5Screate");
d3_space3 = H5Screate_simple(3, d3_dims1, d3_dims2);
CHECK(d3_space3, FAIL, "H5Screate");
d3_space4 = H5Screate_simple(3, d3_dims1, d3_dims3);
CHECK(d3_space4, FAIL, "H5Screate");
/* Compare all dataspace combinations */
/* Compare null dataspace against all others, including itself */
ext_equal = H5Sextent_equal(null_space, null_space);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(null_space, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare scalar dataspace against all others, including itself */
ext_equal = H5Sextent_equal(scalar_space, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, scalar_space);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(scalar_space, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 1-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d1_space1, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d1_space1);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space1, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare larger 1-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d1_space2, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d1_space2);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space2, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 1-D dataspace w/fixed max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d1_space3, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d1_space3);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space3, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 1-D dataspace w/unlimited max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d1_space4, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d1_space4);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d1_space4, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 2-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d2_space1, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d2_space1);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space1, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare larger 2-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d2_space2, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d2_space2);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space2, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 2-D dataspace w/fixed max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d2_space3, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d2_space3);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space3, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 2-D dataspace w/unlimited max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d2_space4, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d2_space4);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d2_space4, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 3-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d3_space1, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d3_space1);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space1, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare larger 2-D dataspace w/no max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d3_space2, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d3_space2);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space2, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 2-D dataspace w/fixed max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d3_space3, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d3_space3);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space3, d3_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
/* Compare small 2-D dataspace w/unlimited max. dims against all others, including itself */
ext_equal = H5Sextent_equal(d3_space4, null_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, scalar_space);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d1_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d1_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d1_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d1_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d2_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d2_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d2_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d2_space4);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d3_space1);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d3_space2);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d3_space3);
VERIFY(ext_equal, FALSE, "H5Sextent_equal");
ext_equal = H5Sextent_equal(d3_space4, d3_space4);
VERIFY(ext_equal, TRUE, "H5Sextent_equal");
/* Close dataspaces */
ret = H5Sclose(null_space);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(scalar_space);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d1_space1);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d1_space2);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d1_space3);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d1_space4);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d2_space1);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d2_space2);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d2_space3);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d2_space4);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d3_space1);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d3_space2);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d3_space3);
CHECK(ret, FAIL, "H5Sclose");
ret = H5Sclose(d3_space4);
CHECK(ret, FAIL, "H5Sclose");
} /* test_h5s_extent_equal() */
/****************************************************************
**
** test_h5s(): Main H5S (dataspace) testing routine.
**
****************************************************************/
void
test_h5s(void)
{
/* Output message about test being performed */
MESSAGE(5, ("Testing Dataspaces\n"));
test_h5s_basic(); /* Test basic H5S code */
test_h5s_null(); /* Test Null dataspace H5S code */
test_h5s_encode(); /* Test encoding and decoding */
test_h5s_scalar_write(); /* Test scalar H5S writing code */
test_h5s_scalar_read(); /* Test scalar H5S reading code */
test_h5s_compound_scalar_write(); /* Test compound datatype scalar H5S writing code */
test_h5s_compound_scalar_read(); /* Test compound datatype scalar H5S reading code */
/* This test was added later to exercise a bug in chunked I/O */
test_h5s_chunk(); /* Exercise bug fix for chunked I/O */
test_h5s_extent_equal(); /* Test extent comparison code */
} /* test_h5s() */
/*-------------------------------------------------------------------------
* Function: cleanup_h5s
*
* Purpose: Cleanup temporary test files
*
* Return: none
*
* Programmer: Albert Cheng
* July 2, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
void
cleanup_h5s(void)
{
remove(DATAFILE);
remove(NULLFILE);
remove(BASICFILE);
}