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
|
/* This is part of the netCDF package. Copyright 2005 University
Corporation for Atmospheric Research/Unidata See COPYRIGHT file for
conditions of use.
This program excersizes HDF5 variable length array code.
@author Ed Hartnett
*/
#include <nc_tests.h>
#include <hdf5.h>
#include <nc_logging.h>
#define FILE_NAME "tst_h_vl.nc"
int
main()
{
printf("\n*** Checking HDF5 VLEN types even more.\n");
printf("*** Reading file created by netcdf-4 tests tst_vl...");
{
hid_t fileid, grpid, fapl_id, hdf_typeid = 0, base_hdf_typeid = 0, attid = 0;
hid_t file_typeid, spaceid, native_typeid;
H5O_info_t obj_info;
char obj_name[NC_MAX_NAME + 1];
hsize_t num_obj, i;
int obj_class;
H5T_class_t class;
size_t size, type_size;
hssize_t att_npoints;
int att_ndims;
hsize_t dims[1];
void *vldata;
/* Open the file and read the vlen data. */
if ((fapl_id = H5Pcreate(H5P_FILE_ACCESS)) < 0) ERR;
if (H5Pset_fclose_degree(fapl_id, H5F_CLOSE_SEMI)) ERR;
if (H5Pset_libver_bounds(fapl_id, H5F_LIBVER_LATEST, H5F_LIBVER_LATEST) < 0) ERR;
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDONLY, fapl_id)) < 0) ERR;
if ((grpid = H5Gopen(fileid, "/")) < 0) ERR;
if (H5Gget_num_objs(grpid, &num_obj) < 0) ERR;
for (i = 0; i < num_obj; i++)
{
if (H5Oget_info_by_idx(grpid, ".", H5_INDEX_CRT_ORDER, H5_ITER_INC,
i, &obj_info, H5P_DEFAULT) < 0) ERR_RET;
obj_class = obj_info.type;
if ((size = H5Lget_name_by_idx(grpid, ".", H5_INDEX_CRT_ORDER, H5_ITER_INC, i,
NULL, 0, H5P_DEFAULT)) < 0) ERR_RET;
if (size > NC_MAX_NAME) ERR_RET;
if (H5Lget_name_by_idx(grpid, ".", H5_INDEX_CRT_ORDER, H5_ITER_INC, i,
obj_name, size+1, H5P_DEFAULT) < 0) ERR_RET;
/*printf("nc4_rec_read_metadata: encountered HDF5 object obj_class %d obj_name %s.\n",
obj_class, obj_name);*/
/* Deal with groups and datasets. */
switch (obj_class)
{
case H5O_TYPE_GROUP:
break;
case H5O_TYPE_DATASET:
break;
case H5O_TYPE_NAMED_DATATYPE:
if ((hdf_typeid = H5Topen2(grpid, obj_name, H5P_DEFAULT)) < 0) ERR_RET;
if ((class = H5Tget_class(hdf_typeid)) < 0) ERR_RET;
switch (class)
{
case H5T_VLEN:
/* Find the base type of this vlen (i.e. what is this a * vlen of?) and its size. */
if (!(base_hdf_typeid = H5Tget_super(hdf_typeid))) ERR_RET;
if (!(type_size = H5Tget_size(base_hdf_typeid))) ERR_RET;
break;
default:
LOG((0, "unknown datatype class"));
return NC_EBADCLASS;
}
break;
case H5G_LINK:
break;
default:
LOG((0, "Unknown object class %d in nc4_rec_read_metadata!",
obj_class));
}
}
num_obj = H5Aget_num_attrs(grpid);
for (i = 0; i < num_obj; i++)
{
if (attid > 0)
H5Aclose(attid);
if ((attid = H5Aopen_idx(grpid, (unsigned int)i)) < 0) ERR_RET;
if (H5Aget_name(attid, NC_MAX_NAME + 1, obj_name) < 0) ERR_RET;
LOG((4, "reading attribute of _netCDF group, named %s", obj_name));
if ((attid && H5Aclose(attid))) ERR_RET;
}
/* Open the HDF5 attribute. */
if ((attid = H5Aopen_name(grpid, obj_name)) < 0) ERR;
/* Get type of attribute in file. */
if ((file_typeid = H5Aget_type(attid)) < 0) ERR;
if ((class = H5Tget_class(file_typeid)) < 0) ERR;
if (class != H5T_VLEN) ERR;
/* Get len. */
if ((spaceid = H5Aget_space(attid)) < 0) ERR;
if ((att_ndims = H5Sget_simple_extent_ndims(spaceid)) < 0) ERR;
if (att_ndims != 1) ERR;
if ((att_npoints = H5Sget_simple_extent_npoints(spaceid)) < 0) ERR;
if (H5Sget_simple_extent_dims(spaceid, dims, NULL) < 0) ERR;
if (dims[0] != 3) ERR;
if ((native_typeid = H5Tget_native_type(file_typeid, H5T_DIR_DEFAULT)) < 0) ERR;
if (!(vldata = malloc((unsigned int)(dims[0] * sizeof(hvl_t))))) ERR;
if (H5Aread(attid, native_typeid, vldata) < 0) ERR;
if (H5Tclose(native_typeid) < 0) ERR;
if (H5Tclose(file_typeid) < 0) ERR;
if (hdf_typeid && H5Tclose(hdf_typeid) < 0) ERR;
if (base_hdf_typeid && H5Tclose(base_hdf_typeid) < 0) ERR;
if (H5Aclose(attid) ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
}
SUMMARIZE_ERR;
FINAL_RESULTS;
}
|
