diff options
| -rw-r--r-- | examples/test_read.c | 98 | ||||
| -rw-r--r-- | examples/test_write.c | 130 |
2 files changed, 228 insertions, 0 deletions
diff --git a/examples/test_read.c b/examples/test_read.c new file mode 100644 index 0000000..7f70519 --- /dev/null +++ b/examples/test_read.c @@ -0,0 +1,98 @@ + +/* This program finds two datasets in HDF5 file created by test_write.c + program */ + +#include "hdf5.h" + +/* Define datasets names */ + + +#define DATA1_NAME "ShortIntegers" +#define DATA2_NAME "DoubleFloats" + +main(void) + +{ + + hid_t file_id; /* HDF5 file ID */ + + hid_t t1_id; /* Type ID for the first dataset */ + hid_t d1_id; /* Dimensionality ID for the first dataset */ + hid_t data1_id; /* Data ID for the first dataset */ + hid_t data1_type; /* Datatype of the first dataset */ + uint8 data1_len; /* Length in bytes for data1_type */ + uint8 data1_arch; /* Architecture format data1_type is in */ + uint32 data1_rank; /* Rank of the first dataset */ + + hid_t t2_id; /* Type ID for the second dataset */ + hid_t d2_id; /* Dimensionality ID for the second dataset */ + hid_t data2_id; /* Data ID for the second dataset */ + hid_t data2_type; /* Datatype of the first dataset */ + uint8 data2_len; /* Length in bytes for data2_type */ + uint8 data2_arch; /* Architecture format data2_type is in */ + uint32 data2_rank; /* Rank of the second dataset */ + + + herr_t ret; + hbool_t ret_file; + +/* Check if we have HDF5 file */ + + ret_file = H5Fis_hdf5("HDF5_file"); + if(ret_file == BTRUE) + printf("Is HDF5? yes\n"); + +/* Open HDF5 file */ + + file_id = H5Fopen("HDF5_file", 0, 0); + + +/* Find first dataset and return its dataset ID*/ + + data1_id = H5Mfind_name(file_id,H5_DATASET,DATA1_NAME); + +/* Return datatype and dimensionality IDs for the found dataset */ + + ret = H5Dget_info(data1_id, &t1_id, &d1_id); + +/* Get type, base length and architecture for the first data set */ + + ret = H5Tget_type(t1_id, &data1_type, &data1_len, &data1_arch); + printf("Base length and architecture for the first data set\n"); + printf("%u %u \n", data1_len, data1_arch); + +/* Get rank for the first data set */ + + data1_rank = H5Pget_lrank(d1_id); + printf("data1_rank %u\n", data1_rank); + +/* Find second dataset and return its dataset ID*/ + + data2_id = H5Mfind_name(file_id,H5_DATASET,DATA2_NAME); + +/* Return datatype and dimensionality IDs for the found dataset */ + + ret = H5Dget_info(data2_id, &t2_id, &d2_id); + +/* Get type, base length and architecture for the second data set */ + + ret = H5Tget_type(t2_id, &data2_type, &data2_len, &data2_arch); + printf(" Base length and architecture for the second data set\n"); + printf("%i %i \n", data2_len, data2_arch); + +/* Get rank for the second data set */ + + data2_rank = H5Pget_lrank(d2_id); + printf("data2_rank %u\n", data2_rank); + ret = H5Mrelease(t1_id); + ret = H5Mrelease(d1_id); + ret = H5Mrelease(data1_id); + + ret = H5Mrelease(t2_id); + ret = H5Mrelease(d2_id); + ret = H5Mrelease(data2_id); + + H5Fclose(file_id); + + +} diff --git a/examples/test_write.c b/examples/test_write.c new file mode 100644 index 0000000..15421b7 --- /dev/null +++ b/examples/test_write.c @@ -0,0 +1,130 @@ + +/* This program creates two datasets and stores them in HDF5 file */ + +#include "hdf5.h" + +/* Define data type, dimensionlity, architecture and name for the first + dataset */ + + +#define DATA1_BASE H5T_INT +#define DATA1_LEN 2 +#define DATA1_RANK 2 +#define DATA1_DIM1 10 +#define DATA1_DIM2 5 +#define DATA1_ARCH H5T_BIGENDIAN +#define DATA1_NAME "ShortIntegers" + +/* Define data type, dimensionality, architecture and name for the second + dataset */ + + +#define DATA2_BASE H5T_FLOAT +#define DATA2_LEN 8 +#define DATA2_RANK 1 +#define DATA2_DIM1 8 +#define DATA2_ARCH H5T_LITTLEENDIAN +#define DATA2_NAME "DoubleFloats" + +main(void) + +{ + + hid_t file1_id, file2_id; /* HDF5 file IDs */ + + hid_t t1_id; /* Type ID for the first dataset */ + hid_t d1_id; /* Dimensionality ID for the first dataset */ + hid_t data1_id; /* Data ID for the first dataset */ + + uint32 dims1[DATA1_RANK] = {DATA1_DIM1, DATA1_DIM2}; + int16 data1[DATA1_DIM1][DATA1_DIM2]; + + + hid_t t2_id; /* Type ID for the second dataset */ + hid_t d2_id; /* Dimensionality ID for the second dataset */ + hid_t data2_id; /* Data ID for the second dataset */ + + uint32 dims2[DATA2_RANK] = {DATA2_DIM1}; + float64 data2[DATA2_DIM1]; + + herr_t ret; + + intn k,l; + +/* Initialize datasets */ + + for (k=0; k<DATA1_DIM1; k++) { + for (l=0; l<DATA1_DIM2; l++) { + data1[k][l] = (k+l); + } + } + + + for (k=0; k<DATA2_DIM1; k++) { + data2[k] = (double) (k); + } + +/* Create HDF5 file */ + + file1_id = H5Fcreate("HDF5_file", H5ACC_OVERWRITE,0,0); + +/* Create and initialize datatype object for the first dataset */ + + t1_id = H5Mcreate(file1_id,H5_DATATYPE,"Short Integer data type"); + ret = H5Tset_type(t1_id, DATA1_BASE, DATA1_LEN, DATA1_ARCH); + +/* Create and initialize dimensionality object for the first dataset */ + + d1_id = H5Mcreate(file1_id, H5_DATASPACE, "2D"); + ret = H5Pset_space(d1_id,DATA1_RANK, dims1); + +/* Create and initialize first dataset object */ + + data1_id = H5Mcreate(file1_id, H5_DATASET, DATA1_NAME); + ret = H5Dset_info (data1_id, t1_id, d1_id); + +/* Write first dataset into HDF5_file */ + + ret = H5Dwrite(data1_id, H5P_ALL, data1); + +/* Release all ID's and close the file */ + + ret = H5Mrelease(t1_id); + ret = H5Mrelease(d1_id); + ret = H5Mrelease(data1_id); + + H5Fclose(file1_id); + +/* Reopen the file to add second dataset */ + + file2_id = H5Fopen("HDF5_file", H5ACC_WRITE, 0); + + +/* Create and initialize datatype object for the second dataset */ + + t2_id = H5Mcreate(file2_id,H5_DATATYPE,"Double float data type"); + ret = H5Tset_type(t2_id, DATA2_BASE, DATA2_LEN, DATA2_ARCH); + +/* Create and initialize dimensionality object for the second dataset */ + + d2_id = H5Mcreate(file2_id, H5_DATASPACE, "1D"); + ret = H5Pset_space(d2_id,DATA2_RANK, dims2); + +/* Create and initialize second dataset object */ + + data2_id = H5Mcreate(file2_id, H5_DATASET, DATA2_NAME); + ret = H5Dset_info (data2_id, t2_id, d2_id); + +/* Write second dataset into HDF5_file */ + + ret = H5Dwrite(data2_id, H5P_ALL, data2); + +/* Release all ID's and close the file */ + + ret = H5Mrelease(t2_id); + ret = H5Mrelease(d2_id); + ret = H5Mrelease(data2_id); + + ret = H5Fclose(file2_id); + +} |