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
|
#include <dx/dx.h>
/* necessary for hdf library routines */
/* IMPORTANT NOTE: Data Explorer does not provide HDF libraries. This
example program is for illustrative purposes only. Because we do
not provide hdf libraries, it will not be possible for you to compile
and link this module unless you have the HDF libraries yourself
*/
#include <df.h>
#define MAXRANK 3
Error m_SimpleImport(Object *in, Object *out)
{
Array a=NULL;
Field f=NULL;
char *filename;
int dims, counts[MAXRANK], numelements, i, j;
float deltas[MAXRANK*MAXRANK], origins[MAXRANK], *data;
/* extract the file name from in[0] */
if (!in[0]) {
DXSetError(ERROR_BAD_PARAMETER,"missing filename");
goto error;
}
else if (!DXExtractString(in[0], &filename)) {
DXSetError(ERROR_BAD_PARAMETER, "filename must be a string");
goto error;
}
/* check to see that the file is accessible, and is an hdf file */
/* DFishdf is an hdf library routine */
if (DFishdf(filename) != 0) {
DXSetError(ERROR_BAD_PARAMETER,
"file \"%s\" is either not accessible, or is not an hdf file",
filename);
goto error;
}
/* initialize */
/* DFSDrestart is an HDF library routine */
DFSDrestart();
/* this returns the dimensionality of the grid (1D, 2D, etc) in dims. The
number of positions in each dimension is returned in the array counts */
/* DFSDgetdims is an HDF library routine */
DFSDgetdims(filename, &dims, counts, MAXRANK);
/* make a new data array (scalar) */
a = DXNewArray(TYPE_FLOAT, CATEGORY_REAL, 0);
if (!a)
goto error;
/* figure out how many elements there are in the data array */
for (i=0, numelements=1; i<dims; numelements *= counts[i], i++)
;
/* allocate space in the data array */
if (!DXAddArrayData(a, 0, numelements, NULL))
goto error;
/* get a pointer to the data array */
data = (float *)DXGetArrayData(a);
if (!data)
goto error;
/* read the data from the HDF file to the data array */
/* DFSDgetdata is an HDF library routine */
DFSDgetdata(filename, dims, counts, data);
/* create a new field */
f = DXNewField();
if (!f)
goto error;
/* set the dependency of the data to be on positions */
if (!DXSetStringAttribute((Object)a, "dep", "positions"))
goto error;
/* set the data array as the data component of f */
if (!DXSetComponentValue(f, "data", (Object)a))
goto error;
a=NULL;
/* create the connections array */
/* DXMakeGridConnections will set the element type */
a = DXMakeGridConnectionsV(dims, counts);
if (!a)
goto error;
if (!DXSetComponentValue(f, "connections", (Object)a))
goto error;
a=NULL;
/* now create the positions array; origin 0 and deltas 1 in each dimension */
for (i=0; i<dims; i++) {
origins[i] = 0.0;
for (j=0; j<dims; j++) {
if (i==j)
deltas[i*dims + j] = 1.0;
else
deltas[i*dims + j] = 0.0;
}
}
a = DXMakeGridPositionsV(dims, counts, origins, deltas);
if (!a)
goto error;
if (!DXSetComponentValue(f, "positions", (Object)a))
goto error;
a=NULL;
/* EndField sets default attributes (such as setting the connections
attribute to be "ref" positions), and creates the bounding box. */
if (!DXEndField(f))
goto error;
out[0]=f;
return OK;
error:
DXDelete((Object)f);
DXDelete((Object)a);
return ERROR;
}
|
