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
|
/*
* Copyright (C) by Argonne National Laboratory
* See COPYRIGHT in top-level directory
*/
#include <stdio.h>
#include <stdlib.h>
#include "mpi.h"
#include "mpithreadtest.h"
#include "mpitest.h"
int rank;
MTEST_THREAD_RETURN_TYPE run_test(void *arg);
MTEST_THREAD_RETURN_TYPE run_test(void *arg)
{
MPI_Status reqstat;
int i, j;
int peer = rank ? 0 : 1;
for (j = 0; j < 10; j++) {
if (rank % 2) {
for (i = 0; i < 16; i++)
MPI_Send(NULL, 0, MPI_CHAR, peer, j, MPI_COMM_WORLD);
for (i = 0; i < 16; i++)
MPI_Recv(NULL, 0, MPI_CHAR, peer, j, MPI_COMM_WORLD, &reqstat);
} else {
for (i = 0; i < 16; i++)
MPI_Recv(NULL, 0, MPI_CHAR, peer, j, MPI_COMM_WORLD, &reqstat);
for (i = 0; i < 16; i++)
MPI_Send(NULL, 0, MPI_CHAR, peer, j, MPI_COMM_WORLD);
}
}
return (MTEST_THREAD_RETURN_TYPE) NULL;
}
int main(int argc, char **argv)
{
int zero = 0, pmode, nprocs;
MTest_Init_thread(&argc, &argv, MPI_THREAD_MULTIPLE, &pmode);
if (pmode != MPI_THREAD_MULTIPLE) {
fprintf(stderr, "Thread Multiple not supported by the MPI implementation\n");
MPI_Abort(MPI_COMM_WORLD, -1);
}
MPI_Comm_size(MPI_COMM_WORLD, &nprocs);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
if (nprocs != 2) {
fprintf(stderr, "Need two processes\n");
MPI_Abort(MPI_COMM_WORLD, -1);
}
MTest_Start_thread(run_test, NULL);
run_test(&zero);
MTest_Join_threads();
MTest_Finalize(0);
return 0;
}
|
