[svn-r132] Examples directory 'examples' was added. It contains two programs

test_write.c    -   creates an HDF5 file, writes first dataset,
                         then closes the file,
                         reopens it and adds the second data set

     test_read.c     -   finds datasets created by test_write.c program
                         and prints out some information about datasets.
This commit is contained in:
Elena Pourmal 1997-10-22 16:15:54 -05:00
parent cca9fb56f7
commit 4fe5ac0999
2 changed files with 228 additions and 0 deletions

98
examples/test_read.c Normal file
View File

@ -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);
}

130
examples/test_write.c Normal file
View File

@ -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);
}