-
-
Notifications
You must be signed in to change notification settings - Fork 265
/
h5ex_d_unlimgzip.c
295 lines (257 loc) · 8.21 KB
/
h5ex_d_unlimgzip.c
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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
/************************************************************
This example shows how to create and extend an unlimited
dataset with gzip compression. The program first writes
integers to a gzip compressed dataset with dataspace
dimensions of DIM0xDIM1, then closes the file. Next, it
reopens the file, reads back the data, outputs it to the
screen, extends the dataset, and writes new data to the
extended portions of the dataset. Finally it reopens the
file again, reads back the data, and outputs it to the
screen.
************************************************************/
#include "hdf5.h"
#include <stdio.h>
#include <stdlib.h>
#define FILE "h5ex_d_unlimgzip.h5"
#define DATASET "DS1"
#define DIM0 4
#define DIM1 7
#define EDIM0 6
#define EDIM1 10
#define CHUNK0 4
#define CHUNK1 4
int
main(void)
{
hid_t file = H5I_INVALID_HID;
hid_t space = H5I_INVALID_HID;
hid_t dset = H5I_INVALID_HID;
hid_t dcpl = H5I_INVALID_HID;
herr_t status;
htri_t avail;
H5Z_filter_t filter_type;
hsize_t dims[2] = {DIM0, DIM1};
hsize_t extdims[2] = {EDIM0, EDIM1};
hsize_t maxdims[2];
hsize_t chunk[2] = {CHUNK0, CHUNK1};
hsize_t start[2];
hsize_t count[2];
size_t nelmts;
unsigned int flags, filter_info;
int wdata[DIM0][DIM1]; /* Write buffer */
int wdata2[EDIM0][EDIM1]; /* Write buffer for extension */
int **rdata = NULL; /* Read buffer */
int ndims;
hsize_t i, j;
/*
* Check if gzip compression is available and can be used for both
* compression and decompression. Normally we do not perform error
* checking in these examples for the sake of clarity, but in this
* case we will make an exception because this filter is an
* optional part of the hdf5 library.
*/
avail = H5Zfilter_avail(H5Z_FILTER_DEFLATE);
if (!avail) {
printf("gzip filter not available.\n");
return 1;
}
status = H5Zget_filter_info(H5Z_FILTER_DEFLATE, &filter_info);
if (!(filter_info & H5Z_FILTER_CONFIG_ENCODE_ENABLED) ||
!(filter_info & H5Z_FILTER_CONFIG_DECODE_ENABLED)) {
printf("gzip filter not available for encoding and decoding.\n");
return 1;
}
/*
* Initialize data.
*/
for (i = 0; i < DIM0; i++)
for (j = 0; j < DIM1; j++)
wdata[i][j] = i * j - j;
/*
* Create a new file using the default properties.
*/
file = H5Fcreate(FILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
/*
* Create dataspace with unlimited dimensions.
*/
maxdims[0] = H5S_UNLIMITED;
maxdims[1] = H5S_UNLIMITED;
space = H5Screate_simple(2, dims, maxdims);
/*
* Create the dataset creation property list, add the gzip
* compression filter and set the chunk size.
*/
dcpl = H5Pcreate(H5P_DATASET_CREATE);
status = H5Pset_deflate(dcpl, 9);
status = H5Pset_chunk(dcpl, 2, chunk);
/*
* Create the compressed unlimited dataset.
*/
dset = H5Dcreate(file, DATASET, H5T_STD_I32LE, space, H5P_DEFAULT, dcpl, H5P_DEFAULT);
/*
* Write the data to the dataset.
*/
status = H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, wdata[0]);
/*
* Close and release resources.
*/
status = H5Pclose(dcpl);
status = H5Dclose(dset);
status = H5Sclose(space);
status = H5Fclose(file);
/*
* In this next section we read back the data, extend the dataset,
* and write new data to the extended portions.
*/
/*
* Open file and dataset using the default properties.
*/
file = H5Fopen(FILE, H5F_ACC_RDWR, H5P_DEFAULT);
dset = H5Dopen(file, DATASET, H5P_DEFAULT);
/*
* Get dataspace and allocate memory for read buffer. This is a
* two dimensional dataset so the dynamic allocation must be done
* in steps.
*/
space = H5Dget_space(dset);
ndims = H5Sget_simple_extent_dims(space, dims, NULL);
/*
* Allocate array of pointers to rows.
*/
rdata = (int **)malloc(dims[0] * sizeof(int *));
/*
* Allocate space for integer data.
*/
rdata[0] = (int *)malloc(dims[0] * dims[1] * sizeof(int));
/*
* Set the rest of the pointers to rows to the correct addresses.
*/
for (i = 1; i < dims[0]; i++)
rdata[i] = rdata[0] + i * dims[1];
/*
* Read the data using the default properties.
*/
status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata[0]);
/*
* Output the data to the screen.
*/
printf("Dataset before extension:\n");
for (i = 0; i < dims[0]; i++) {
printf(" [");
for (j = 0; j < dims[1]; j++)
printf(" %3d", rdata[i][j]);
printf("]\n");
}
status = H5Sclose(space);
/*
* Extend the dataset.
*/
status = H5Dset_extent(dset, extdims);
/*
* Retrieve the dataspace for the newly extended dataset.
*/
space = H5Dget_space(dset);
/*
* Initialize data for writing to the extended dataset.
*/
for (i = 0; i < EDIM0; i++)
for (j = 0; j < EDIM1; j++)
wdata2[i][j] = j;
/*
* Select the entire dataspace.
*/
status = H5Sselect_all(space);
/*
* Subtract a hyperslab reflecting the original dimensions from the
* selection. The selection now contains only the newly extended
* portions of the dataset.
*/
start[0] = 0;
start[1] = 0;
count[0] = dims[0];
count[1] = dims[1];
status = H5Sselect_hyperslab(space, H5S_SELECT_NOTB, start, NULL, count, NULL);
/*
* Write the data to the selected portion of the dataset.
*/
status = H5Dwrite(dset, H5T_NATIVE_INT, H5S_ALL, space, H5P_DEFAULT, wdata2[0]);
/*
* Close and release resources.
*/
free(rdata[0]);
free(rdata);
status = H5Dclose(dset);
status = H5Sclose(space);
status = H5Fclose(file);
/*
* Now we simply read back the data and output it to the screen.
*/
/*
* Open file and dataset using the default properties.
*/
file = H5Fopen(FILE, H5F_ACC_RDONLY, H5P_DEFAULT);
dset = H5Dopen(file, DATASET, H5P_DEFAULT);
/*
* Retrieve dataset creation property list.
*/
dcpl = H5Dget_create_plist(dset);
/*
* Retrieve and print the filter type. Here we only retrieve the
* first filter because we know that we only added one filter.
*/
nelmts = 0;
filter_type = H5Pget_filter(dcpl, 0, &flags, &nelmts, NULL, 0, NULL, &filter_info);
printf("\nFilter type is: ");
switch (filter_type) {
case H5Z_FILTER_DEFLATE:
printf("H5Z_FILTER_DEFLATE\n");
break;
case H5Z_FILTER_SHUFFLE:
printf("H5Z_FILTER_SHUFFLE\n");
break;
case H5Z_FILTER_FLETCHER32:
printf("H5Z_FILTER_FLETCHER32\n");
break;
case H5Z_FILTER_SZIP:
printf("H5Z_FILTER_SZIP\n");
break;
case H5Z_FILTER_NBIT:
printf("H5Z_FILTER_NBIT\n");
break;
case H5Z_FILTER_SCALEOFFSET:
printf("H5Z_FILTER_SCALEOFFSET\n");
}
/*
* Get dataspace and allocate memory for the read buffer as before.
*/
space = H5Dget_space(dset);
ndims = H5Sget_simple_extent_dims(space, dims, NULL);
rdata = (int **)malloc(dims[0] * sizeof(int *));
rdata[0] = (int *)malloc(dims[0] * dims[1] * sizeof(int));
for (i = 1; i < dims[0]; i++)
rdata[i] = rdata[0] + i * dims[1];
/*
* Read the data using the default properties.
*/
status = H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, rdata[0]);
/*
* Output the data to the screen.
*/
printf("Dataset after extension:\n");
for (i = 0; i < dims[0]; i++) {
printf(" [");
for (j = 0; j < dims[1]; j++)
printf(" %3d", rdata[i][j]);
printf("]\n");
}
/*
* Close and release resources.
*/
free(rdata[0]);
free(rdata);
status = H5Pclose(dcpl);
status = H5Dclose(dset);
status = H5Sclose(space);
status = H5Fclose(file);
return 0;
}