2 /* Creator: Bronis R. de Supinski (bronis@llnl.gov) Fri Mar 17 2000 */
3 /* no-error.c -- do some MPI calls without any errors */
10 #define INTERCOMM_CREATE_TAG 666
13 main (int argc, char **argv)
17 char processor_name[128];
22 MPI_Comm temp, intercomm;
24 int drank, dnprocs, rleader, rnprocs;
27 MPI_Init (&argc, &argv);
28 MPI_Comm_size (MPI_COMM_WORLD, &nprocs);
29 MPI_Comm_rank (MPI_COMM_WORLD, &rank);
30 MPI_Get_processor_name (processor_name, &namelen);
31 printf ("(%d) is alive on %s\n", rank, processor_name);
34 MPI_Barrier (MPI_COMM_WORLD);
37 printf ("not enough tasks\n");
40 /* need to make split communicator temporarily... */
41 MPI_Comm_split (MPI_COMM_WORLD, rank % 2, nprocs - rank, &temp);
43 if (temp != MPI_COMM_NULL) {
44 MPI_Comm_size (temp, &tnprocs);
45 MPI_Comm_rank (temp, &trank);
47 /* create an intercommunicator temporarily so can merge it... */
48 rleader = ((rank + nprocs) % 2) ? nprocs - 2 : nprocs - 1;
50 if ((rank % 2) && (trank == 0)) {
51 MPI_Recv (buf0, buf_size, MPI_INT, 2, 0, MPI_COMM_WORLD, &status);
54 MPI_Recv (buf0, buf_size, MPI_INT, 1, 0, MPI_COMM_WORLD, &status);
57 MPI_Intercomm_create (temp, 0, MPI_COMM_WORLD, rleader,
58 INTERCOMM_CREATE_TAG, &intercomm);
61 memset (buf0, 0, buf_size);
63 MPI_Send (buf0, buf_size, MPI_INT, rleader, 0, MPI_COMM_WORLD);
66 memset (buf0, 0, buf_size);
68 MPI_Send (buf0, buf_size, MPI_INT, 0, 0, MPI_COMM_WORLD);
71 MPI_Comm_free (&temp);
73 if (intercomm != MPI_COMM_NULL) {
74 MPI_Comm_size (intercomm, &dnprocs);
75 MPI_Comm_rank (intercomm, &drank);
76 MPI_Comm_remote_size (intercomm, &rnprocs);
78 if (rnprocs > drank) {
80 memset (buf1, 1, buf_size);
82 MPI_Recv (buf0, buf_size, MPI_INT, drank, 0, intercomm, &status);
84 MPI_Send (buf1, buf_size, MPI_INT, drank, 0, intercomm);
87 memset (buf0, 0, buf_size);
89 MPI_Send (buf0, buf_size, MPI_INT, drank, 0, intercomm);
91 MPI_Recv (buf1, buf_size, MPI_INT, drank, 0, intercomm, &status);
95 printf ("(%d) Intercomm too small (lrank = %d; remote size = %d)\n",
96 rank, drank, rnprocs);
99 MPI_Comm_free (&intercomm);
102 printf ("(%d) Got MPI_COMM_NULL\n", rank);
106 printf ("(%d) MPI_Comm_split got MPI_COMM_NULL\n", rank);
110 MPI_Barrier (MPI_COMM_WORLD);
113 printf ("(%d) Finished normally\n", rank);