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
|
/********************************************************************
*
* Testing for thread safety in H5A (dataset attribute) library
* operations. -- Threaded program --
* ------------------------------------------------------------------
*
* Plan: Attempt to break H5Acreate by making many simultaneous create
* calls.
*
* Claim: N calls to H5Acreate should create N attributes for a dataset
* if threadsafe. If some unprotected shared data exists for the
* dataset (eg, a count of the number of attributes in the
* dataset), there is a small chance that consecutive reads occur
* before a write to that shared variable.
*
* Created: Oct 5 1999
* Programmer: Chee Wai LEE
*
* Modification History
* --------------------
* May 15 2000 - incorporated into library tests (Chee Wai LEE)
*
********************************************************************/
#include "ttsafe.h"
#ifndef H5_HAVE_THREADSAFE
static int dummy; /* just to create a non-empty object file */
#else
#include <stdio.h>
#include <stdlib.h>
#define FILE "ttsafe.h5"
#define DATASETNAME "IntData"
#define NUM_THREADS 16
void *tts_acreate_thread(void *);
typedef struct acreate_data_struct {
hid_t dataset;
hid_t datatype;
hid_t dataspace;
int current_index;
} ttsafe_name_data_t;
void tts_acreate(void) {
/* Pthread definitions
*/
pthread_t threads[NUM_THREADS];
/* HDF5 data definitions
*/
hid_t file, dataset;
hid_t dataspace, datatype;
hid_t attribute;
hsize_t dimsf[1]; /* dataset dimensions */
int data; /* data to write */
int buffer, ret;
int i;
ttsafe_name_data_t *attrib_data;
/* create a hdf5 file using H5F_ACC_TRUNC access,
* default file creation plist and default file
* access plist
*/
file = H5Fcreate(FILE, H5F_ACC_TRUNC, H5P_DEFAULT, H5P_DEFAULT);
/* create a simple dataspace for the dataset
*/
dimsf[0] = 1;
dataspace = H5Screate_simple(1,dimsf,NULL);
/* define datatype for the data using native little endian integers
*/
datatype = H5Tcopy(H5T_NATIVE_INT);
H5Tset_order(datatype, H5T_ORDER_LE);
/* create a new dataset within the file
*/
dataset = H5Dcreate(file, DATASETNAME, datatype, dataspace,
H5P_DEFAULT);
/* initialize data for dataset and write value to dataset
*/
data = NUM_THREADS;
H5Dwrite(dataset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL,
H5P_DEFAULT, &data);
/* simultaneously create a large number of attributes to be
associated with the dataset
*/
for (i=0;i<NUM_THREADS;i++) {
attrib_data = malloc(sizeof(ttsafe_name_data_t));
attrib_data->dataset = dataset;
attrib_data->datatype = datatype;
attrib_data->dataspace = dataspace;
attrib_data->current_index = i;
pthread_create(&threads[i],NULL,tts_acreate_thread,attrib_data);
}
for (i=0;i<NUM_THREADS;i++) {
pthread_join(threads[i],NULL);
}
/* verify the correctness of the test */
for (i=0; i<NUM_THREADS; i++) {
attribute = H5Aopen_name(dataset,gen_name(i));
if (attribute < 0) {
fprintf(stderr,"unable to open appropriate attribute. Test failed!\n");
} else {
ret = H5Aread(attribute, H5T_NATIVE_INT, &buffer);
if ((ret < 0) || (buffer != i)) {
fprintf(stderr,"wrong data values. Test failed!\n");
}
H5Aclose(attribute);
}
}
/* close remaining resources
*/
H5Sclose(dataspace);
H5Tclose(datatype);
H5Dclose(dataset);
H5Fclose(file);
}
void *tts_acreate_thread(void *client_data) {
hid_t attribute;
hsize_t dimsf[1]; /* dataset dimensions */
char *attribute_name;
int *attribute_data; /* data for attributes */
int i;
ttsafe_name_data_t *attrib_data = (ttsafe_name_data_t *)client_data;
/* create attribute
*/
attribute_name = gen_name(attrib_data->current_index);
attribute = H5Acreate(attrib_data->dataset,
attribute_name,
attrib_data->datatype,
attrib_data->dataspace,
H5P_DEFAULT);
/* Write data to the attribute
*/
attribute_data = malloc(sizeof(int));
*attribute_data = attrib_data->current_index;
H5Awrite(attribute,H5T_NATIVE_INT,attribute_data);
H5Aclose(attribute);
return NULL;
}
void cleanup_acreate(void) {
}
#endif /*H5_HAVE_THREADSAFE*/
|
