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.
16 const int verbose = 0;
18 int main(int argc, char **argv) {
19 int i, j, rank, nproc;
20 int errors = 0, all_errors = 0;
25 MPI_Init(&argc, &argv);
27 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
28 MPI_Comm_size(MPI_COMM_WORLD, &nproc);
30 val_ptrs = malloc(nproc * sizeof(MPI_Aint));
31 MPI_Get_address(&val, &val_ptrs[rank]);
33 MPI_Allgather(MPI_IN_PLACE, 0, MPI_DATATYPE_NULL, val_ptrs, 1, MPI_AINT,
36 MPI_Win_create_dynamic(MPI_INFO_NULL, MPI_COMM_WORLD, &dyn_win);
37 MPI_Win_attach(dyn_win, &one, sizeof(int));
39 for (i = 0; i < ITER; i++) {
40 MPI_Win_fence(MPI_MODE_NOPRECEDE, dyn_win);
41 MPI_Accumulate(&one, 1, MPI_INT, i%nproc, val_ptrs[i%nproc], 1, MPI_INT, MPI_SUM, dyn_win);
42 MPI_Win_fence(MPI_MODE_NOSUCCEED, dyn_win);
45 MPI_Barrier(MPI_COMM_WORLD);
47 /* Read and verify my data */
50 printf("%d -- Got %d, expected %d\n", rank, val, ITER);
53 MPI_Win_detach(dyn_win, &one);
54 MPI_Win_free(&dyn_win);
56 MPI_Reduce(&errors, &all_errors, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD);
58 if (rank == 0 && all_errors == 0)
59 printf(" No Errors\n");