Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Merge branch 'clean_events' of github.com:Takishipp/simgrid into clean_events
[simgrid.git] / teshsuite / smpi / isp / umpire / intercomm_create-deadlock3.c
1 /* -*- Mode: C; -*- */
2 /* Creator: Bronis R. de Supinski (bronis@llnl.gov) Fri Mar  17 2000 */
3 /* no-error.c -- do some MPI calls without any errors */
4
5 #include <stdio.h>
6 #include "mpi.h"
7
8 #define buf_size 128
9
10 #define INTERCOMM_CREATE_TAG 666
11
12 int
13 main (int argc, char **argv)
14 {
15   int nprocs = -1;
16   int rank = -1;
17   char processor_name[128];
18   int namelen = 128;
19   int buf0[buf_size];
20   int buf1[buf_size];
21   MPI_Status status;
22   MPI_Comm temp, intercomm;
23   int trank, tnprocs;
24   int drank, dnprocs, rleader, rnprocs;
25
26   /* init */
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);
32   fflush (stdout);
33
34   MPI_Barrier (MPI_COMM_WORLD);
35
36   if (nprocs < 4) {
37     printf ("not enough tasks\n");
38   }
39   else {
40     /* need to make split communicator temporarily... */
41     MPI_Comm_split (MPI_COMM_WORLD, rank % 2, nprocs - rank, &temp);
42
43     if (temp != MPI_COMM_NULL) {
44       MPI_Comm_size (temp, &tnprocs);
45       MPI_Comm_rank (temp, &trank);
46
47       /* create an intercommunicator temporarily so can merge it... */
48       rleader = ((rank + nprocs) % 2) ?  nprocs - 2 : nprocs - 1;
49
50       if (rank == 1) {
51         MPI_Recv (buf0, buf_size, MPI_INT, 0, 0, MPI_COMM_WORLD, &status);
52       }
53
54       MPI_Intercomm_create (temp, 0, MPI_COMM_WORLD, rleader,
55                             INTERCOMM_CREATE_TAG, &intercomm);
56
57       if (rank == 0) {
58         memset (buf0, 0, buf_size);
59         
60         MPI_Send (buf0, buf_size, MPI_INT, 1, 0, MPI_COMM_WORLD);
61       }
62
63       MPI_Comm_free (&temp);
64
65       if (intercomm != MPI_COMM_NULL) {
66         MPI_Comm_size (intercomm, &dnprocs);
67         MPI_Comm_rank (intercomm, &drank);
68         MPI_Comm_remote_size (intercomm, &rnprocs);
69
70         if (rnprocs > drank) {
71           if (rank % 2) {
72             memset (buf1, 1, buf_size);
73
74             MPI_Recv (buf0, buf_size, MPI_INT, drank, 0, intercomm, &status);
75
76             MPI_Send (buf1, buf_size, MPI_INT, drank, 0, intercomm);
77           }
78           else {
79             memset (buf0, 0, buf_size);
80         
81             MPI_Send (buf0, buf_size, MPI_INT, drank, 0, intercomm);
82         
83             MPI_Recv (buf1, buf_size, MPI_INT, drank, 0, intercomm, &status);
84           }
85         }
86         else {
87           printf ("(%d) Intercomm too small (lrank = %d; remote size = %d)\n",
88                   rank, drank, rnprocs);
89         }
90
91         MPI_Comm_free (&intercomm);
92       }
93       else {
94         printf ("(%d) Got MPI_COMM_NULL\n", rank);
95       }
96     }
97     else {
98       printf ("(%d) MPI_Comm_split got MPI_COMM_NULL\n", rank);
99     }
100   }
101
102   MPI_Barrier (MPI_COMM_WORLD);
103
104   MPI_Finalize ();
105   printf ("(%d) Finished normally\n", rank);
106 }
107
108 /* EOF */