1 /* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
4 * (C) 2012 by Argonne National Laboratory.
5 * See COPYRIGHT in top-level directory.
14 #define ITER_PER_RANK 25
16 const int verbose = 0;
18 int main(int argc, char **argv)
21 int errors = 0, all_errors = 0;
27 MPI_Init(&argc, &argv);
29 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
30 MPI_Comm_size(MPI_COMM_WORLD, &nproc);
32 iter = ITER_PER_RANK * nproc;
34 val_ptrs = malloc(nproc * sizeof(MPI_Aint));
35 MPI_Get_address(&val, &val_ptrs[rank]);
37 MPI_Allgather(MPI_IN_PLACE, 0, MPI_DATATYPE_NULL, val_ptrs, 1, MPI_AINT, MPI_COMM_WORLD);
39 MPI_Win_create_dynamic(MPI_INFO_NULL, MPI_COMM_WORLD, &dyn_win);
40 MPI_Win_attach(dyn_win, &val, sizeof(int));
42 for (i = 0; i < iter; i++) {
43 MPI_Win_fence(MPI_MODE_NOPRECEDE, dyn_win);
44 MPI_Accumulate(&one, 1, MPI_INT, i % nproc, val_ptrs[i % nproc], 1, MPI_INT, MPI_SUM,
46 MPI_Win_fence(MPI_MODE_NOSUCCEED, dyn_win);
49 MPI_Barrier(MPI_COMM_WORLD);
51 /* Read and verify my data */
54 printf("%d -- Got %d, expected %d\n", rank, val, iter);
57 MPI_Win_detach(dyn_win, &val);
58 MPI_Win_free(&dyn_win);
60 MPI_Reduce(&errors, &all_errors, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD);
62 if (rank == 0 && all_errors == 0)
63 printf(" No Errors\n");