| 1 | #include "mpi.h"
|
|---|
| 2 | #include "assert.h"
|
|---|
| 3 |
|
|---|
| 4 | int many_to_one () {
|
|---|
| 5 | int size;
|
|---|
| 6 | int rank;
|
|---|
| 7 | MPI_Comm comm = MPI_COMM_WORLD;
|
|---|
| 8 | MPI_Datatype type = MPI_INT;
|
|---|
| 9 |
|
|---|
| 10 | MPI_Comm_size (comm, &size);
|
|---|
| 11 | MPI_Comm_rank (comm, &rank);
|
|---|
| 12 |
|
|---|
| 13 | if (rank != 0) {
|
|---|
| 14 | int dat[3] = {rank, rank + 1, rank + 2};
|
|---|
| 15 | MPI_Request reqs[3];
|
|---|
| 16 |
|
|---|
| 17 | for (int i = 0; i < 3; i++)
|
|---|
| 18 | MPI_Isend(dat + i, 1, MPI_INT, 0, i, comm, reqs + i);
|
|---|
| 19 | MPI_Waitall(3, reqs, MPI_STATUS_IGNORE);
|
|---|
| 20 | } else {
|
|---|
| 21 | int ub = 3 * (size - 1);
|
|---|
| 22 | int dat[ub];
|
|---|
| 23 | MPI_Request reqs[ub];
|
|---|
| 24 | MPI_Status statuses[ub];
|
|---|
| 25 |
|
|---|
| 26 | for (int i = 0; i < ub; i++) {
|
|---|
| 27 | int nd = $choose_int(2);
|
|---|
| 28 |
|
|---|
| 29 | if (nd)
|
|---|
| 30 | MPI_Irecv(dat + i, 1, MPI_INT, MPI_ANY_SOURCE, MPI_ANY_TAG, comm, reqs + i);
|
|---|
| 31 | else {
|
|---|
| 32 | MPI_Recv(dat + i, 1, MPI_INT, MPI_ANY_SOURCE, MPI_ANY_TAG, comm, statuses + i);
|
|---|
| 33 | reqs[i] = MPI_REQUEST_NULL;
|
|---|
| 34 | }
|
|---|
| 35 | }
|
|---|
| 36 | MPI_Waitall(ub, reqs, statuses);
|
|---|
| 37 | for (int i = 0; i < ub; i++) {
|
|---|
| 38 | int src = statuses[i].MPI_SOURCE;
|
|---|
| 39 | int tag = statuses[i].MPI_TAG;
|
|---|
| 40 |
|
|---|
| 41 | assert (dat[i] == src + tag);
|
|---|
| 42 | }
|
|---|
| 43 | }
|
|---|
| 44 | return 0;
|
|---|
| 45 | }
|
|---|
| 46 |
|
|---|
| 47 | int main() {
|
|---|
| 48 | MPI_Init(NULL, NULL);
|
|---|
| 49 | many_to_one();
|
|---|
| 50 | MPI_Finalize();
|
|---|
| 51 | return 0;
|
|---|
| 52 | }
|
|---|