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
|
#include "h5rados_example.h"
int main(int argc, char *argv[]) {
rados_t cluster;
char *pool = "mypool";
hid_t file = -1, dset = -1, fapl = -1, space = -1, type = -1, dcpl = -1, dapl = -1, def_dcpl = -1, def_dapl = -1;
int ndims;
hsize_t dims[2];
htri_t tri_ret;
(void)MPI_Init(&argc, &argv);
if(argc != 3)
PRINTF_ERROR("argc != 3\n");
if(rados_create(&cluster, NULL) < 0)
ERROR;
if(rados_conf_read_file(cluster, CEPH_CONFIG_FILE) < 0)
ERROR;
/* Initialize VOL */
if(H5VLrados_init(cluster, pool) < 0)
ERROR;
/* Set up FAPL */
if((fapl = H5Pcreate(H5P_FILE_ACCESS)) < 0)
ERROR;
if(H5Pset_fapl_rados(fapl, MPI_COMM_WORLD, MPI_INFO_NULL) < 0)
ERROR;
if(H5Pset_all_coll_metadata_ops(fapl, true) < 0)
ERROR;
/* Open file */
if((file = H5Fopen(argv[1], H5F_ACC_RDONLY, fapl)) < 0)
ERROR;
printf("Opening dataset\n");
/* Open dataset */
if((dset = H5Dopen2(file, argv[2], H5P_DEFAULT)) < 0)
ERROR;
/* Check dataset dataspace */
if((space = H5Dget_space(dset)) < 0)
ERROR;
if((ndims = H5Sget_simple_extent_ndims(space)) < 0)
ERROR;
if(ndims != 2)
PRINTF_ERROR("ndims == %d, expected 2\n", ndims);
if(H5Sget_simple_extent_dims(space, dims, NULL) < 0)
ERROR;
if(dims[0] != 4)
PRINTF_ERROR("dims[0] == %d, expected 4\n", (int)dims[0]);
if(dims[1] != 6)
PRINTF_ERROR("dims[1] == %d, expected 6\n", (int)dims[1]);
/* Check dataset datatype */
if((type = H5Dget_type(dset)) < 0)
ERROR;
if((tri_ret = H5Tequal(type, H5T_NATIVE_INT)) < 0)
ERROR;
if(!tri_ret)
PRINTF_ERROR("datatype does not equal H5T_NATIVE_INT\n");
/* Check DCPL */
if((dcpl = H5Dget_create_plist(dset)) < 0)
ERROR;
if((def_dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0)
ERROR;
if((tri_ret = H5Pequal(dcpl, def_dcpl)) < 0)
ERROR;
if(!tri_ret)
PRINTF_ERROR("DCPL does not equal default\n");
/* Check DAPL */
if((dapl = H5Dget_access_plist(dset)) < 0)
ERROR;
if((def_dapl = H5Pcreate(H5P_DATASET_ACCESS)) < 0)
ERROR;
if((tri_ret = H5Pequal(dapl, def_dapl)) < 0)
ERROR;
if(!tri_ret)
PRINTF_ERROR("DAPL does not equal default\n");
/* Close */
if(H5Dclose(dset) < 0)
ERROR;
if(H5Fclose(file) < 0)
ERROR;
if(H5Pclose(fapl) < 0)
ERROR;
if(H5Sclose(space) < 0)
ERROR;
if(H5Tclose(type) < 0)
ERROR;
if(H5Pclose(dcpl) < 0)
ERROR;
if(H5Pclose(dapl) < 0)
ERROR;
if(H5Pclose(def_dcpl) < 0)
ERROR;
if(H5Pclose(def_dapl) < 0)
ERROR;
printf("Success\n");
(void)MPI_Finalize();
return 0;
error:
H5E_BEGIN_TRY {
H5Dclose(dset);
H5Fclose(file);
H5Pclose(fapl);
H5Sclose(space);
H5Tclose(type);
H5Pclose(dcpl);
H5Pclose(dapl);
H5Pclose(def_dcpl);
H5Pclose(def_dapl);
} H5E_END_TRY;
(void)MPI_Finalize();
return 1;
}
|
