teshsuite/smpi/mpich3-test/rma/win_dynamic_acc.c

   1 /* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
   2 /*
   3  *
   4  *  (C) 2012 by Argonne National Laboratory.
   5  *      See COPYRIGHT in top-level directory.
   6  */
   7
   8 #include <stdio.h>
   9 #include <stdlib.h>
  10 #include <assert.h>
  11 #include <mpi.h>
  12 #include "mpitest.h"
  13
  14 #define ITER 100
  15
  16 const int verbose = 0;
  17
  18 int main(int argc, char **argv) {
  19     int       i, j, rank, nproc;
  20     int       errors = 0, all_errors = 0;
  21     int       val = 0, one = 1;
  22     MPI_Aint *val_ptrs;
  23     MPI_Win   dyn_win;
  24
  25     MPI_Init(&argc, &argv);
  26
  27     MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  28     MPI_Comm_size(MPI_COMM_WORLD, &nproc);
  29
  30     val_ptrs = malloc(nproc * sizeof(MPI_Aint));
  31     MPI_Get_address(&val, &val_ptrs[rank]);
  32
  33     MPI_Allgather(MPI_IN_PLACE, 0, MPI_DATATYPE_NULL, val_ptrs, 1, MPI_AINT,
  34                   MPI_COMM_WORLD);
  35
  36     MPI_Win_create_dynamic(MPI_INFO_NULL, MPI_COMM_WORLD, &dyn_win);
  37     MPI_Win_attach(dyn_win, &one, sizeof(int));
  38
  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);
  43     }
  44
  45     MPI_Barrier(MPI_COMM_WORLD);
  46
  47     /* Read and verify my data */
  48     if ( val != ITER ) {
  49         errors++;
  50         printf("%d -- Got %d, expected %d\n", rank, val, ITER);
  51     }
  52
  53     MPI_Win_detach(dyn_win, &one);
  54     MPI_Win_free(&dyn_win);
  55
  56     MPI_Reduce(&errors, &all_errors, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD);
  57
  58     if (rank == 0 && all_errors == 0)
  59         printf(" No Errors\n");
  60
  61     free(val_ptrs);
  62     MPI_Finalize();
  63
  64     return 0;
  65 }