summaryrefslogtreecommitdiffstats
path: root/examples/h5_write.c
blob: 00ac680784edb7f1a323f60694aaf26519d6d092 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
/*  
 *  This example writes data to HDF5 file.
 *  Data conversion is performed during write operation.  
 */
 
#include "hdf5.h"

#define FILE        "SDS.h5"
#define DATASETNAME "IntArray" 
#define NX     5                      /* dataset dimensions */
#define NY     6
#define RANK   2

main ()
{
   hid_t       file, dataset;         /* file and dataset handles */
   hid_t       datatype, dataspace;   /* handles */
   size_t      dimsf[2];              /* dataset dimensions */
   herr_t      status;                             
   int32       data[NX][NY];          /* data to write */
   int         i, j;

/* 
 * Data  and output buffer initialization. 
 */

for (j = 0; j < NX; j++) {
    for (i = 0; i < NY; i++)
        data[j][i] = i + j;
}     
                                       /*  0 1 2 3 4 5 
                                           1 2 3 4 5 6
                                           2 3 4 5 6 7
                                           3 4 5 6 7 8
                                           4 5 6 7 8 9   */

/*
 * Create a new file using H5ACC_OVERWRITE access,
 * default file creation properties, and default file
 * access properties.
 */
file = H5Fcreate(FILE, H5ACC_OVERWRITE, H5C_DEFAULT, H5C_DEFAULT);

/*
 * Describe the size of the array and create the data space for fixed
 * size dataset. 
 */
dimsf[0] = NX;
dimsf[1] = NY;
dataspace = H5Pcreate_simple(RANK, dimsf, NULL); 

/* 
 * Define datatype for the data in the file.
 * We will store little endian INT32 numbers.
 */
datatype = H5Tcopy(H5T_NATIVE_INT32);
status = H5Tset_order(datatype, H5T_ORDER_LE);
/*
 * Create a new dataset within the file using defined dataspace and
 * datatype and default dataset creation properties.
 */
dataset = H5Dcreate(file, DATASETNAME, datatype, dataspace,
                    H5C_DEFAULT);

/*
 * Write the data to the dataset using default transfer properties.
 */
status = H5Dwrite(dataset, H5T_NATIVE_INT32, H5P_ALL, H5P_ALL,
                  H5C_DEFAULT, data);

/*
 * Close/release resources.
 */
H5Pclose(dataspace);
H5Tclose(datatype);
H5Dclose(dataset);
H5Fclose(file);
 
}