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

   1 /* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
   2 /*
   3  *
   4  *  (C) 2003 by Argonne National Laboratory.
   5  *      See COPYRIGHT in top-level directory.
   6  */
   7 #include "mpi.h"
   8 #include <stdio.h>
   9 #include <stdlib.h>
  10 #include "mpitest.h"
  11
  12 #ifndef MAX_INT
  13 #define MAX_INT 0x7fffffff
  14 #endif
  15
  16 /*
  17 static char MTEST_Descrip[] = "Test MPI_Accumulate with fence";
  18 */
  19
  20 int main( int argc, char *argv[] )
  21 {
  22     int errs = 0;
  23     int rank, size, source;
  24     int minsize = 2, count, i;
  25     MPI_Comm      comm;
  26     MPI_Win       win;
  27     int           *winbuf, *sbuf;
  28
  29     MTest_Init( &argc, &argv );
  30
  31     /* The following illustrates the use of the routines to
  32        run through a selection of communicators and datatypes.
  33        Use subsets of these for tests that do not involve combinations
  34        of communicators, datatypes, and counts of datatypes */
  35     while (MTestGetIntracommGeneral( &comm, minsize, 1 )) {
  36         if (comm == MPI_COMM_NULL) continue;
  37         /* Determine the sender and receiver */
  38         MPI_Comm_rank( comm, &rank );
  39         MPI_Comm_size( comm, &size );
  40         source = 0;
  41
  42         for (count = 1; count < 65000; count = count * 2) {
  43             /* We compare with an integer value that can be as large as
  44                size * (count * count + (1/2)*(size-1))
  45                For large machines (size large), this can exceed the
  46                maximum integer for some large values of count.  We check
  47                that in advance and break this loop if the above value
  48                would exceed MAX_INT.  Specifically,
  49
  50                size*count*count + (1/2)*size*(size-1) > MAX_INT
  51                count*count > (MAX_INT/size - (1/2)*(size-1))
  52             */
  53             if (count * count > (MAX_INT/size - (size-1)/2)) break;
  54             winbuf = (int *)malloc( count * sizeof(int) );
  55             sbuf   = (int *)malloc( count * sizeof(int) );
  56
  57             for (i=0; i<count; i++) winbuf[i] = 0;
  58             for (i=0; i<count; i++) sbuf[i] = rank + i * count;
  59             MPI_Win_create( winbuf, count * sizeof(int), sizeof(int),
  60                             MPI_INFO_NULL, comm, &win );
  61             MPI_Win_fence( 0, win );
  62             MPI_Accumulate( sbuf, count, MPI_INT, source, 0, count, MPI_INT,
  63                                 MPI_SUM, win );
  64             MPI_Win_fence( 0, win );
  65             if (rank == source) {
  66                 /* Check the results */
  67                 for (i=0; i<count; i++) {
  68                     int result = i * count * size + (size*(size-1))/2;
  69                     if (winbuf[i] != result) {
  70                         if (errs < 10) {
  71                             fprintf( stderr, "Winbuf[%d] = %d, expected %d (count = %d, size = %d)\n",
  72                                      i, winbuf[i], result, count, size );
  73                         }
  74                         errs++;
  75                     }
  76                 }
  77             }
  78             free( winbuf );
  79             free( sbuf );
  80             MPI_Win_free( &win );
  81         }
  82         MTestFreeComm(&comm);
  83     }
  84
  85     MTest_Finalize( errs );
  86     MPI_Finalize();
  87     return 0;
  88 }