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
|
/*
* Copyright (C) 1998 NCSA
* All rights reserved.
*
* Programmer: Robb Matzke <matzke@llnl.gov>
* Wednesday, January 21, 1998
*
* Purpose: Simple data space functions.
*/
#include <H5private.h>
#include <H5Eprivate.h>
#include <H5Pprivate.h>
/* Interface initialization */
#define PABLO_MASK H5P_simp_mask
#define INTERFACE_INIT NULL
static intn interface_initialize_g = FALSE;
/*-------------------------------------------------------------------------
* Function: H5P_simp_init
*
* Purpose: Generates element numbering information for the data
* spaces involved in a data space conversion.
*
* Return: Success: Number of elements that can be efficiently
* transferred at a time.
*
* Failure: Zero
*
* Programmer: Robb Matzke
* Wednesday, January 21, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
size_t
H5P_simp_init (const struct H5O_layout_t *layout, const H5P_t *mem_space,
const H5P_t *file_space, H5P_number_t *numbering/*out*/)
{
size_t nelmts;
FUNC_ENTER (H5P_simp_init, 0);
/* Check args */
assert (layout);
assert (mem_space && H5P_SIMPLE==mem_space->type);
assert (file_space && H5P_SIMPLE==file_space->type);
assert (numbering);
/* Numbering is implied by the hyperslab, C order */
HDmemset (numbering, 0, sizeof(H5P_number_t));
/* Data can be efficiently copied at any size */
nelmts = H5P_get_npoints (file_space);
FUNC_LEAVE (nelmts);
}
/*-------------------------------------------------------------------------
* Function: H5P_simp_fgath
*
* Purpose: Gathers data points from file F and accumulates them in the
* type conversion buffer BUF. The LAYOUT argument describes
* how the data is stored on disk. ELMT_SIZE is the size in
* bytes of a datum which this function treats as opaque.
* FILE_SPACE describes the data space of the dataset on disk
* and the elements that have been selected for reading (via
* hyperslab, etc) and NUMBERING describes how those elements
* are numbered (initialized by the H5P_*_init() call). This
* function will copy at most NELMTS elements beginning at the
* element numbered START.
*
* Return: Success: Number of elements copied.
*
* Failure: 0
*
* Programmer: Robb Matzke
* Wednesday, January 21, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
size_t
H5P_simp_fgath (H5F_t *f, const struct H5O_layout_t *layout,
size_t elmt_size, const H5P_t *file_space,
const H5P_number_t *numbering, intn start, intn nelmts,
void *buf/*out*/)
{
size_t offset[H5O_LAYOUT_NDIMS]; /*offset of hyperslab */
size_t size[H5O_LAYOUT_NDIMS]; /*size of hyperslab */
intn i; /*counters */
FUNC_ENTER (H5P_simp_fgath, 0);
/* Check args */
assert (f);
assert (layout);
assert (elmt_size>0);
assert (file_space);
assert (numbering);
assert (nelmts>0);
assert (buf);
/*
* The prototype doesn't support strip mining.
*/
assert (0==start);
assert (nelmts==H5P_get_npoints (file_space));
/*
* Quincey, this is where we look at FILE_SPACE to decide what the
* hyperslab is to read from disk. For now, since the H5P interface
* doesn't support hyperslabs, we'll assume the caller is asking for the
* entire array. --RPM
*/
assert (nelmts == H5P_get_npoints (file_space));
for (i=0; i<layout->ndims; i++) offset[i] = 0;
i = H5P_get_dims (file_space, size);
assert (i+1 == layout->ndims);
size[i] = elmt_size;
/*
* Gather from file.
*/
if (H5F_arr_read (f, layout, size, offset, offset, size, buf/*out*/)<0) {
HRETURN_ERROR (H5E_DATASPACE, H5E_READERROR, 0, "read error");
}
FUNC_LEAVE (nelmts);
}
/*-------------------------------------------------------------------------
* Function: H5P_simp_mscat
*
* Purpose: Scatters data points from the type conversion buffer
* TCONV_BUF to the application buffer BUF. Each element is
* ELMT_SIZE bytes and they are organized in application memory
* according to MEM_SPACE. The NUMBERING information together
* with START and NELMTS describe how the elements stored in
* TCONV_BUF are globally numbered.
*
* Return: Success: SUCCEED
*
* Failure: FAIL
*
* Programmer: Robb Matzke
* Wednesday, January 21, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
herr_t
H5P_simp_mscat (const void *tconv_buf, size_t elmt_size,
const H5P_t *mem_space, const H5P_number_t *numbering,
intn start, intn nelmts, void *buf/*out*/)
{
FUNC_ENTER (H5P_simp_mscat, FAIL);
/* Check args */
assert (tconv_buf);
assert (elmt_size>0);
assert (mem_space && H5P_SIMPLE==mem_space->type);
assert (numbering);
assert (nelmts>0);
assert (buf);
/*
* The prototype doesn't support strip mining.
*/
assert (0==start);
assert (nelmts==H5P_get_npoints (mem_space));
/*
* Quincey, this is where we look at the hyperslab spec of MEM_SPACE to
* figure out how to scatter. For now we just assume that data points
* are copied directly from TCONV_BUF to BUF.
*/
HDmemcpy (buf, tconv_buf, nelmts*elmt_size);
FUNC_LEAVE (SUCCEED);
}
/*-------------------------------------------------------------------------
* Function: H5P_simp_mgath
*
* Purpose: Gathers dataset elements from application memory and copies
* them into the data type conversion buffer. NOT IMPLEMENTED
* YET.
*
* Return: Success: Number of elements copied.
*
* Failure: 0
*
* Programmer: Robb Matzke
* Wednesday, January 21, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
size_t
H5P_simp_mgath (void)
{
FUNC_ENTER (H5P_simp_mgath, 0);
HRETURN_ERROR (H5E_DATASPACE, H5E_UNSUPPORTED, 0,
"not implemented yet");
FUNC_LEAVE (0);
}
/*-------------------------------------------------------------------------
* Function: H5P_simp_fscat
*
* Purpose: Scatters dataset elements from the type conversion buffer to
* the file. NOT IMPLEMENTED YET.
*
* Return: Success: SUCCEED
*
* Failure: FAIL
*
* Programmer: Robb Matzke
* Wednesday, January 21, 1998
*
* Modifications:
*
*-------------------------------------------------------------------------
*/
herr_t
H5P_simp_fscat (void)
{
FUNC_ENTER (H5P_simp_fscat, FAIL);
HRETURN_ERROR (H5E_DATASPACE, H5E_UNSUPPORTED, FAIL,
"not implemented yet");
FUNC_LEAVE (FAIL);
}
|
