+++ /dev/null
-#include "mpi.h"
-#include <stdio.h>
-#include "test.h"
-
-MPI_Comm GetNextComm( void );
-void addem ( int *, int *, int *, MPI_Datatype * );
-void assoc ( int *, int *, int *, MPI_Datatype * );
-
-void addem( int *invec, int *inoutvec, int *len, MPI_Datatype *dtype)
-{
- int i;
- for ( i=0; i<*len; i++ )
- inoutvec[i] += invec[i];
-}
-
-#define BAD_ANSWER 100000
-
-/*
- The operation is inoutvec[i] = invec[i] op inoutvec[i]
- (see 4.9.4). The order is important.
-
- Note that the computation is in process rank (in the communicator)
- order, independant of the root.
- */
-void assoc( int *invec, int *inoutvec, int *len, MPI_Datatype *dtype)
-{
- int i;
- for ( i=0; i<*len; i++ ) {
- if (inoutvec[i] <= invec[i] ) {
- int rank;
- MPI_Comm_rank( MPI_COMM_WORLD, &rank );
- fprintf( stderr, "[%d] inout[0] = %d, in[0] = %d\n",
- rank, inoutvec[0], invec[0] );
- inoutvec[i] = BAD_ANSWER;
- }
- else
- inoutvec[i] = invec[i];
- }
-}
-
-MPI_Comm GetNextComm( void )
-{
- MPI_Comm comm = MPI_COMM_NULL;
- static int idx = 0;
- int size, rank;
- MPI_Comm_size( MPI_COMM_WORLD, &size );
- MPI_Comm_rank( MPI_COMM_WORLD, &rank );
-
- switch (idx) {
- case 0:
- MPI_Comm_dup( MPI_COMM_WORLD, &comm );
- break;
- case 1:
- /* invert the rank order */
- MPI_Comm_split( MPI_COMM_WORLD, 0, size - rank, &comm );
- break;
- case 2:
- /* Divide into subsets */
- MPI_Comm_split( MPI_COMM_WORLD, rank < (size/2), rank, &comm );
- break;
- case 3:
- /* Another division */
- MPI_Comm_split( MPI_COMM_WORLD, rank < (size/3), size-rank, &comm );
- break;
- case 4:
- /* odd and even */
- MPI_Comm_split( MPI_COMM_WORLD, (rank % 2) == 0, rank, &comm );
- break;
- case 5:
- /* Last case: startover */
- idx = -1;
- break;
- }
- idx++;
- return comm;
-}
-
-int main( int argc, char **argv )
-{
- int rank, size, i;
- int data;
- int errors=0;
- int result = -100;
- int correct_result;
- MPI_Op op_assoc, op_addem;
- MPI_Comm comm;
-
- MPI_Init( &argc, &argv );
- MPI_Op_create( (MPI_User_function *)assoc, 0, &op_assoc );
- MPI_Op_create( (MPI_User_function *)addem, 1, &op_addem );
-
- /* Run this for a variety of communicator sizes */
- while ((comm = GetNextComm()) != MPI_COMM_NULL) {
- MPI_Comm_rank( comm, &rank );
- MPI_Comm_size( comm, &size );
-
- data = rank;
-
- correct_result = 0;
- for (i=0;i<=rank;i++)
- correct_result += i;
-
- MPI_Scan ( &data, &result, 1, MPI_INT, MPI_SUM, comm );
- if (result != correct_result) {
- fprintf( stderr, "[%d] Error suming ints with scan\n", rank );
- errors++;
- }
-
- MPI_Scan ( &data, &result, 1, MPI_INT, MPI_SUM, comm );
- if (result != correct_result) {
- fprintf( stderr, "[%d] Error summing ints with scan (2)\n", rank );
- errors++;
- }
-
- data = rank;
- result = -100;
- MPI_Scan ( &data, &result, 1, MPI_INT, op_addem, comm );
- if (result != correct_result) {
- fprintf( stderr, "[%d] Error summing ints with scan (userop)\n",
- rank );
- errors++;
- }
-
- MPI_Scan ( &data, &result, 1, MPI_INT, op_addem, comm );
- if (result != correct_result) {
- fprintf( stderr, "[%d] Error summing ints with scan (userop2)\n",
- rank );
- errors++;
- }
-/* result = -100;*/
-/* data = rank;*/
-/* MPI_Scan ( &data, &result, 1, MPI_INT, op_assoc, comm );*/
-/* if (result == BAD_ANSWER) {*/
-/* fprintf( stderr, "[%d] Error scanning with non-commutative op\n",*/
-/* rank );*/
-/* errors++;*/
-/* }*/
- MPI_Comm_free( &comm );
- }
-
- MPI_Op_free( &op_assoc );
- MPI_Op_free( &op_addem );
-
- if (errors) {
- MPI_Comm_rank( MPI_COMM_WORLD, &rank );
- printf( "[%d] done with ERRORS(%d)!\n", rank, errors );
- }
-
- Test_Waitforall( );
- MPI_Finalize();
- return errors;
-}