1 /* -*- Mode: C; c-basic-offset:4 ; indent-tabs-mode:nil ; -*- */
4 * (C) 2003 by Argonne National Laboratory.
5 * See COPYRIGHT in top-level directory.
11 #include "mpicolltest.h"
14 static char MTEST_Descrip[] = "Simple intercomm gather test";
17 int main(int argc, char *argv[])
21 int leftGroup, i, count, rank, rsize, size;
23 MPI_Datatype datatype;
25 MTest_Init(&argc, &argv);
28 /* Get an intercommunicator */
29 while (MTestGetIntercomm(&comm, &leftGroup, 4)) {
30 if (comm == MPI_COMM_NULL)
32 MPI_Comm_rank(comm, &rank);
33 MPI_Comm_remote_size(comm, &rsize);
34 MPI_Comm_size(comm, &size);
35 /* To improve reporting of problems about operations, we
36 * change the error handler to errors return */
37 MPI_Comm_set_errhandler(comm, MPI_ERRORS_RETURN);
38 for (count = 1; count < 65000; count = 2 * count) {
40 buf = (int *) malloc(count * rsize * sizeof(int));
41 for (i = 0; i < count * rsize; i++)
44 err = MTest_Gather(NULL, 0, datatype,
46 (rank == 0) ? MPI_ROOT : MPI_PROC_NULL, comm);
51 /* Test that no other process in this group received the
54 for (i = 0; i < count; i++) {
61 /* Check for the correct data */
62 for (i = 0; i < count * rsize; i++) {
70 /* In the right group */
71 buf = (int *) malloc(count * sizeof(int));
72 for (i = 0; i < count; i++)
73 buf[i] = rank * count + i;
74 err = MTest_Gather(buf, count, datatype, NULL, 0, datatype, 0, comm);