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
|
#include "swmr_common.h"
#define CHUNK_SIZE 50
static int
gen_skeleton(const char *filename, unsigned verbose, int comp_level)
{
hid_t fid; /* File ID for new HDF5 file */
hid_t fcpl; /* File creation property list */
hid_t fapl; /* File access property list */
hid_t dcpl; /* Dataset creation property list */
hid_t tid; /* Datatype for dataset elements */
hid_t sid; /* Dataspace ID */
hsize_t dims = 0; /* Dataset starting dimensions */
hsize_t max_dims = H5S_UNLIMITED; /* Dataset maximum dimensions */
hsize_t chunk_dims = CHUNK_SIZE; /* Chunk dimensions */
unsigned u, v; /* Local index variable */
/* Create file access property list */
if((fapl = H5Pcreate(H5P_FILE_ACCESS)) < 0)
return(-1);
if(H5Pset_libver_bounds(fapl, H5F_LIBVER_LATEST, H5F_LIBVER_LATEST) < 0)
return(-1);
#ifdef QAK
/* Increase the initial size of the metadata cache */
{
H5AC_cache_config_t mdc_config;
mdc_config.version = H5AC__CURR_CACHE_CONFIG_VERSION;
H5Pget_mdc_config(fapl, &mdc_config);
printf("mdc_config.initial_size = %lu\n", (unsigned long)mdc_config.initial_size);
printf("mdc_config.epoch_length = %lu\n", (unsigned long)mdc_config.epoch_length);
mdc_config.set_initial_size = 1;
mdc_config.initial_size = 16 * 1024 * 1024;
/* mdc_config.epoch_length = 5000; */
H5Pset_mdc_config(fapl, &mdc_config);
}
#endif /* QAK */
#ifdef QAK
H5Pset_small_data_block_size(fapl, (hsize_t)(50 * CHUNK_SIZE * DTYPE_SIZE));
#endif /* QAK */
#ifdef QAK
H5Pset_fapl_log(fapl, "append.log", H5FD_LOG_ALL, (size_t)(512 * 1024 * 1024));
#endif /* QAK */
/* Create file creation property list */
if((fcpl = H5Pcreate(H5P_FILE_CREATE)) < 0)
return(-1);
#ifdef QAK
H5Pset_link_phase_change(fcpl, 0, 0);
#endif /* QAK */
/* Emit informational message */
if(verbose)
printf("Creating file\n");
/* Create the file */
if((fid = H5Fcreate(filename, H5F_ACC_TRUNC, fcpl, fapl)) < 0)
return(-1);
/* Close file creation property list */
if(H5Pclose(fcpl) < 0)
return(-1);
/* Close file access property list */
if(H5Pclose(fapl) < 0)
return(-1);
/* Create datatype for creating datasets */
if((tid = create_symbol_datatype()) < 0)
return(-1);
/* Create dataspace for creating datasets */
if((sid = H5Screate_simple(1, &dims, &max_dims)) < 0)
return(-1);
/* Create dataset creation property list */
if((dcpl = H5Pcreate(H5P_DATASET_CREATE)) < 0)
return(-1);
if(H5Pset_chunk(dcpl, 1, &chunk_dims) < 0)
return(-1);
if(comp_level >= 0) {
if(H5Pset_deflate(dcpl, (unsigned)comp_level) < 0)
return(-1);
} /* end if */
/* Emit informational message */
if(verbose)
printf("Creating datasets\n");
/* Create the datasets */
for(u = 0; u < NLEVELS; u++)
for(v = 0; v < symbol_count[u]; v++) {
hid_t dsid; /* Dataset ID */
char name_buf[64];
generate_name(name_buf, u, v);
if((dsid = H5Dcreate2(fid, name_buf, tid, sid, H5P_DEFAULT, dcpl, H5P_DEFAULT)) < 0)
return(-1);
if(H5Dclose(dsid) < 0)
return(-1);
} /* end for */
/* Emit informational message */
if(verbose)
printf("Closing objects\n");
/* Close everythign */
if(H5Pclose(dcpl) < 0)
return(-1);
if(H5Sclose(sid) < 0)
return(-1);
if(H5Tclose(tid) < 0)
return(-1);
if(H5Fclose(fid) < 0)
return(-1);
return(0);
} /* end gen_skeleton() */
static void
usage(void)
{
printf("Usage error!\n");
printf("Usage: swmr_generator [-q] [-c <deflate compression level>]\n");
printf("<deflate compression level> should be -1 (for no compression) or 0-9\n");
printf("Defaults to verbose (no '-q' given) and no compression ('-c -1')\n");
exit(1);
} /* end usage() */
int main(int argc, const char *argv[])
{
int comp_level = (-1); /* Compression level (-1 is no compression) */
unsigned verbose = 1; /* Whether to emit some informational messages */
unsigned u; /* Local index variables */
/* Parse command line options */
if(argc > 1) {
u = 1;
while(u < (unsigned)argc) {
if(argv[u][0] == '-') {
switch(argv[u][1]) {
/* Compress dataset chunks */
case 'c':
comp_level = atoi(argv[u + 1]);
if(comp_level < -1 || comp_level > 9)
usage();
u += 2;
break;
/* Be quiet */
case 'q':
verbose = 0;
u++;
break;
default:
usage();
break;
} /* end switch */
} /* end if */
} /* end while */
} /* end if */
/* Emit informational message */
if(verbose) {
printf("Parameters:\n");
printf("\tcompression level = %d\n", comp_level);
} /* end if */
/* Emit informational message */
if(verbose)
printf("Generating skeleton file: %s\n", FILENAME);
/* Generate file skeleton */
if(gen_skeleton(FILENAME, verbose, comp_level) < 0) {
printf("Error generating skeleton file!\n");
exit(1);
} /* end if */
return(0);
}
|
