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
|
/* -*- c-file-style: "stroustrup" -*- */
#include "hdf5.h"
#include <stdio.h>
#include <stdlib.h>
int
main(void)
{
int ret_val = EXIT_SUCCESS;
//! <!-- [create] -->
{
__label__ fail_insert, fail_dtype, fail_file;
hid_t file, dtype;
if ((file = H5Fcreate("t1.h5", H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_file;
}
// create a compound datatype with room for real and imaginary parts
if ((dtype = H5Tcreate(H5T_COMPOUND, 2 * sizeof(double))) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_dtype;
}
// add the real part now and the imaginary part later
if (H5Tinsert(dtype, "re", 0, H5T_IEEE_F64LE) < 0) {
ret_val = EXIT_FAILURE;
goto fail_insert;
}
// commit the datatype definition to the file
if (H5Tcommit(file, "pre-complex", dtype, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 0) {
ret_val = EXIT_FAILURE;
}
fail_insert:
H5Tclose(dtype);
fail_dtype:
H5Fclose(file);
fail_file:;
}
//! <!-- [create] -->
//! <!-- [read] -->
{
__label__ fail_dtype, fail_file;
hid_t file, dtype;
if ((file = H5Fopen("t1.h5", H5F_ACC_RDONLY, H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_file;
}
// open the datatype object stored in the file
if ((dtype = H5Topen(file, "pre-complex", H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_dtype;
}
switch (H5Tget_class(dtype)) { // this time we are only interested in compounds
case H5T_COMPOUND:
printf("Record size: %lu bytes\n", H5Tget_size(dtype));
printf("Record has %d field(s).\n", H5Tget_nmembers(dtype));
break;
default:
break;
}
H5Tclose(dtype);
fail_dtype:
H5Fclose(file);
fail_file:;
}
//! <!-- [read] -->
//! <!-- [update] -->
{
__label__ fail_insert, fail_clone, fail_dtype, fail_file;
hid_t file, dtype, clone;
if ((file = H5Fopen("t1.h5", H5F_ACC_RDWR, H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_file;
}
// open the datatype object stored in the file
if ((dtype = H5Topen(file, "pre-complex", H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_dtype;
}
// the original datatype object is immutable and we need to clone it
if ((clone = H5Tcopy(dtype)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_clone;
}
// remember that the original has enough room for real and imaginary parts
// add the imaginary part
if (H5Tinsert(clone, "im", sizeof(double), H5T_IEEE_F64LE) < 0) {
ret_val = EXIT_FAILURE;
goto fail_insert;
}
// commit the "updated" datatype definition to the file
if (H5Tcommit(file, "complex", clone, H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT) < 0) {
ret_val = EXIT_FAILURE;
}
fail_insert:
H5Tclose(clone);
fail_clone:
H5Tclose(dtype);
fail_dtype:
H5Fclose(file);
fail_file:;
}
//! <!-- [update] -->
//! <!-- [delete] -->
{
__label__ fail_file;
hid_t file;
if ((file = H5Fopen("t1.h5", H5F_ACC_RDWR, H5P_DEFAULT)) == H5I_INVALID_HID) {
ret_val = EXIT_FAILURE;
goto fail_file;
}
// delete the "pre-complex" datatype by unlinking
if (H5Ldelete(file, "pre-complex", H5P_DEFAULT) < 0) {
ret_val = EXIT_FAILURE;
}
H5Fclose(file);
fail_file:;
}
//! <!-- [delete] -->
return ret_val;
}
|
