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
|
/*
* Copyright © 1998 NCSA
* All rights reserved.
*
* Programmer: Saurabh Bagchi (bagchi@uiuc.edu)
* Wednesday, August 11, 1999.
*
* Modifications: Saurabh Bagchi (Aug 17, 1999)
* Modified to work with VFL (HDF51.3).
*/
/* Test the following functionality of the GASS driver.
1. Open a remote file for append.
2. Create a new dataset within the file.
3. Create a local memory buffer to hold the data.
4. Write the local data to the remote dataset.
*/
#include <h5test.h>
#ifndef HAVE_GASS
int main(void)
{
printf("Test skipped because GASS driver not available\n");
return 0;
}
#else
#define URL "ftp://acheng:acheng12@eirene/tmp/acheng/junk.w"
#define DATASETNAME "Int1Array"
#define NX 5 /* dataset dimensions */
#define NY 6
#define RANK 2
int main (void)
{
hid_t fapl =-1, file;
hid_t dataspace, datatype, dataset;
hsize_t dimsf[2];
herr_t status = 0;
int data[NX][NY]; /* data to write */
int i, j;
GASS_Info ginf;
/*
* Data and output buffer initialization.
*/
for (j = 0; j < NX; j++) {
for (i = 0; i < NY; i++)
data[j][i] = i*i + j*j;
}
/*
* 0 1 4 9 16 25
* 1 2 5 10 17 26
* 4 5 8 13 20 29
* 9 10 13 18 25 34
* 16 17 20 25 32 41
*/
/* Create access property list and set the driver to GASS */
fapl = H5Pcreate (H5P_FILE_ACCESS);
if (fapl < 0) {
printf (" H5Pcreate failed. \n");
return -1;
}
ginf.block_size = 0;
ginf.max_length =0;
/* ginf = GASS_INFO_NULL; */
status = H5Pset_fapl_gass (fapl, ginf);
if (status < 0) {
printf ("H5Pset_fapl_gass failed. \n");
return -1;
}
/*
* Open an existing file using H5F_ACC_RDWR access,
* and gass file access properties.
*/
file = H5Fopen(URL, H5F_ACC_RDWR, fapl);
if (file < 0) {
printf ("H5Fopen failed. \n");
return -1;
}
/*
* Describe the size of the array and create the data space for fixed
* size dataset.
*/
dimsf[0] = NX;
dimsf[1] = NY;
dataspace = H5Screate_simple(RANK, dimsf, NULL);
if (dataspace < 0) {
printf ("H5Screate failed. \n");
return -1;
}
/*
* Define datatype for the data in the file.
* We will store little endian INT numbers.
*/
datatype = H5Tcopy(H5T_NATIVE_INT);
if (datatype < 0) {
printf ("H5Tcopy failed. \n");
return -1;
}
status = H5Tset_order(datatype, H5T_ORDER_LE);
if (status < 0) {
printf ("H5Tset_order failed. \n");
return -1;
}
/*
* Create a new dataset within the file using defined dataspace and
* datatype and default dataset creation properties.
*/
dataset = H5Dcreate(file, DATASETNAME, datatype, dataspace,
H5P_DEFAULT);
if (dataset < 0) {
printf ("H5Dcreate failed. \n");
return -1;
}
/*
* Write the data to the dataset using default transfer properties.
*/
status = H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL,
H5P_DEFAULT, data);
if (status < 0) {
printf ("H5Dwrite failed. \n");
return -1;
}
/*
* Close/release resources.
*/
H5Sclose(dataspace);
H5Tclose(datatype);
H5Dclose(dataset);
H5Fclose(file);
H5Pclose(fapl);
return 0;
}
#endif
|
