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
|
/*
* Copyright (C) by Argonne National Laboratory
* See COPYRIGHT in top-level directory
*/
#include <stdio.h>
#include <assert.h>
#include <pthread.h>
#include <stdlib.h>
#include "yaksa.h"
#define DIMSIZE (80)
int **inbuf_, **outbuf_;
void *thread_fn(void *arg);
void *thread_fn(void *arg)
{
int rc = YAKSA_SUCCESS;
yaksa_type_t vector, vector_vector;
uintptr_t actual;
uintptr_t tid = (uintptr_t) arg;
int *inbuf = inbuf_[tid];
int *outbuf = outbuf_[tid];
yaksa_init(NULL);
rc = yaksa_type_create_vector(3, 2, 3, YAKSA_TYPE__INT, NULL, &vector);
assert(rc == YAKSA_SUCCESS);
rc = yaksa_type_create_vector(5, 1, 10, vector, NULL, &vector_vector);
assert(rc == YAKSA_SUCCESS);
for (int i = 0; i < DIMSIZE * DIMSIZE; i++) {
inbuf[i] = i;
outbuf[i] = -1;
}
yaksa_request_t request;
rc = yaksa_ipack(inbuf, 1, vector_vector, 0, outbuf, DIMSIZE * DIMSIZE * sizeof(int), &actual,
NULL, YAKSA_OP__REPLACE, &request);
assert(rc == YAKSA_SUCCESS);
assert(actual == 5 * 3 * 2 * sizeof(int));
rc = yaksa_request_wait(request);
assert(rc == YAKSA_SUCCESS);
int idx = 0;
int val = 0;
for (int i = 0; i < 5; i++) {
for (int j = 0; j < 3; j++) {
for (int k = 0; k < 2; k++) {
inbuf[val] = -1;
if (outbuf[idx] != val)
fprintf(stderr, "outbuf[%d] = %d instead of %d\n", idx, outbuf[idx], val);
idx++;
val++;
}
val++;
}
val += 71;
}
uintptr_t actual_unpack_bytes;
rc = yaksa_iunpack(outbuf, actual, inbuf, 1, vector_vector, 0, &actual_unpack_bytes,
NULL, YAKSA_OP__REPLACE, &request);
assert(rc == YAKSA_SUCCESS);
rc = yaksa_request_wait(request);
assert(rc == YAKSA_SUCCESS);
idx = 0;
for (int i = 0; i < DIMSIZE; i++) {
for (int j = 0; j < DIMSIZE; j++) {
if (inbuf[idx] != idx)
fprintf(stderr, "inbuf[%d] = %d instead of %d\n", idx, inbuf[idx], idx);
idx++;
}
}
yaksa_type_free(vector_vector);
yaksa_type_free(vector);
yaksa_finalize();
return NULL;
}
#define MAX_THREADS (16)
#define ITERS 10000
int main()
{
pthread_t thread[MAX_THREADS];
inbuf_ = (int **) malloc(MAX_THREADS * sizeof(int *));
outbuf_ = (int **) malloc(MAX_THREADS * sizeof(int *));
for (int i = 0; i < MAX_THREADS; i++) {
inbuf_[i] = (int *) malloc(DIMSIZE * DIMSIZE * sizeof(int));
outbuf_[i] = (int *) malloc(DIMSIZE * DIMSIZE * sizeof(int));
}
for (int j = 0; j < ITERS; j++) {
for (uintptr_t i = 0; i < MAX_THREADS; i++)
pthread_create(&thread[i], NULL, thread_fn, (void *) i);
for (int i = 0; i < MAX_THREADS; i++)
pthread_join(thread[i], NULL);
}
for (int i = 0; i < MAX_THREADS; i++) {
free(inbuf_[i]);
free(outbuf_[i]);
}
free(inbuf_);
free(outbuf_);
return 0;
}
|
