Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
SMPI: Kill the global list of senders and receivers
[simgrid.git] / src / smpi / smpi_base.c
1 #include "private.h"
2 #include "xbt/time.h"
3
4 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(smpi_base, smpi,
5                                 "Logging specific to SMPI (base)");
6 XBT_LOG_EXTERNAL_CATEGORY(smpi_base);
7 XBT_LOG_EXTERNAL_CATEGORY(smpi_bench);
8 XBT_LOG_EXTERNAL_CATEGORY(smpi_kernel);
9 XBT_LOG_EXTERNAL_CATEGORY(smpi_mpi);
10 XBT_LOG_EXTERNAL_CATEGORY(smpi_receiver);
11 XBT_LOG_EXTERNAL_CATEGORY(smpi_sender);
12 XBT_LOG_EXTERNAL_CATEGORY(smpi_util);
13
14 smpi_mpi_global_t smpi_mpi_global = NULL;
15
16 void smpi_mpi_land_func(void *a, void *b, int *length,
17                         MPI_Datatype * datatype);
18
19 void smpi_mpi_land_func(void *a, void *b, int *length,
20                         MPI_Datatype * datatype)
21 {
22   int i;
23   if (*datatype == smpi_mpi_global->mpi_int) {
24     int *x = a, *y = b;
25     for (i = 0; i < *length; i++) {
26       y[i] = x[i] && y[i];
27     }
28   }
29 }
30
31 void smpi_mpi_sum_func(void *a, void *b, int *length,
32                        MPI_Datatype * datatype);
33
34 void smpi_mpi_sum_func(void *a, void *b, int *length, MPI_Datatype * datatype)
35 {
36   int i;
37   if (*datatype == smpi_mpi_global->mpi_int) {
38     int *x = a, *y = b;
39     for (i = 0; i < *length; i++) {
40       y[i] = x[i] + y[i];
41     }
42   }
43 }
44
45 int smpi_mpi_comm_rank(smpi_mpi_communicator_t comm)
46 {
47   return comm->index_to_rank_map[smpi_host_index()];
48 }
49
50 void smpi_process_init()
51 {
52   smx_host_t host;
53   int i;
54   smpi_host_data_t hdata;
55
56   smpi_global->running_hosts_count++;
57
58   // initialize some local variables
59   host = SIMIX_host_self();
60
61   hdata = xbt_new(s_smpi_host_data_t, 1);
62   SIMIX_host_set_data(host, hdata);
63   SIMIX_process_set_data(SIMIX_process_self(),hdata);
64
65   for (i = 0; i < smpi_global->host_count && host != smpi_global->hosts[i]; i++);
66
67   hdata->index = i;
68   hdata->mutex = SIMIX_mutex_init();
69   hdata->cond = SIMIX_cond_init();
70
71   hdata->pending_recv_request_queue = xbt_fifo_new();
72
73   hdata->main = SIMIX_process_self();
74   hdata->sender = SIMIX_process_create("smpi_sender",
75           smpi_sender, hdata,
76           SIMIX_host_get_name(SIMIX_host_self()), 0, NULL,
77           /*props */ NULL);
78   hdata->receiver = SIMIX_process_create("smpi_receiver",
79           smpi_receiver, hdata,
80           SIMIX_host_get_name(SIMIX_host_self()), 0, NULL,
81           /*props */ NULL);
82
83   smpi_global->main_processes[hdata->index] = SIMIX_process_self();
84   return;
85 }
86
87 void smpi_process_finalize()
88 {
89   int i;
90   smpi_host_data_t hdata =  SIMIX_host_get_data(SIMIX_host_self());
91
92   i = --smpi_global->running_hosts_count;
93
94   SIMIX_mutex_destroy(smpi_host_mutex());
95   SIMIX_cond_destroy(smpi_host_cond());
96   xbt_fifo_free(hdata->pending_recv_request_queue);
97
98   if (0 >= i) {
99
100     // wake up senders/receivers
101           /* MQ: (FIXME) Don't do so: it breaks since some hosts are already gone
102     for (i = 0; i < smpi_global->host_count; i++) {
103       smpi_host_data_t remote_hdata =  SIMIX_process_get_data(smpi_global->main_processes[i]);
104
105       if (SIMIX_process_is_suspended(remote_hdata->sender))
106         SIMIX_process_resume(remote_hdata->sender);
107
108       if (SIMIX_process_is_suspended(remote_hdata->receiver))
109         SIMIX_process_resume(remote_hdata->receiver);
110     }*/
111
112     SIMIX_mutex_destroy(smpi_mpi_global->mpi_comm_world->barrier_mutex);
113     SIMIX_cond_destroy(smpi_mpi_global->mpi_comm_world->barrier_cond);
114     xbt_free(smpi_mpi_global->mpi_comm_world);
115
116     xbt_free(smpi_mpi_global->mpi_byte);
117     xbt_free(smpi_mpi_global->mpi_int);
118     xbt_free(smpi_mpi_global->mpi_double);
119
120     xbt_free(smpi_mpi_global->mpi_land);
121     xbt_free(smpi_mpi_global->mpi_sum);
122
123     xbt_free(smpi_mpi_global);
124   }
125 }
126
127 int smpi_mpi_barrier(smpi_mpi_communicator_t comm)
128 {
129
130   SIMIX_mutex_lock(comm->barrier_mutex);
131   ++comm->barrier_count;
132   if (comm->barrier_count > comm->size) {       // only happens on second barrier...
133     comm->barrier_count = 0;
134   } else if (comm->barrier_count == comm->size) {
135     SIMIX_cond_broadcast(comm->barrier_cond);
136   }
137   while (comm->barrier_count < comm->size) {
138     SIMIX_cond_wait(comm->barrier_cond, comm->barrier_mutex);
139   }
140   SIMIX_mutex_unlock(comm->barrier_mutex);
141
142   return MPI_SUCCESS;
143 }
144
145 int smpi_mpi_isend(smpi_mpi_request_t request)
146 {
147         smpi_host_data_t hdata =  SIMIX_host_get_data(SIMIX_host_self());
148   int retval = MPI_SUCCESS;
149
150   if (NULL == request) {
151     retval = MPI_ERR_INTERN;
152   } else {
153     xbt_fifo_push(smpi_global->pending_send_request_queues[hdata->index], request);
154
155     if (SIMIX_process_is_suspended(hdata->sender)) {
156       SIMIX_process_resume(hdata->sender);
157     }
158   }
159
160   return retval;
161 }
162
163 int smpi_mpi_irecv(smpi_mpi_request_t request)
164 {
165   int retval = MPI_SUCCESS;
166   smpi_host_data_t hdata =  SIMIX_host_get_data(SIMIX_host_self());
167
168   if (NULL == request) {
169     retval = MPI_ERR_INTERN;
170   } else {
171     xbt_fifo_push(hdata->pending_recv_request_queue, request);
172
173     if (SIMIX_process_is_suspended(hdata->receiver)) {
174       SIMIX_process_resume(hdata->receiver);
175     }
176   }
177
178   return retval;
179 }
180
181 int smpi_mpi_wait(smpi_mpi_request_t request, smpi_mpi_status_t * status)
182 {
183   int retval = MPI_SUCCESS;
184
185   if (NULL == request) {
186     retval = MPI_ERR_INTERN;
187   } else {
188     SIMIX_mutex_lock(request->mutex);
189     while (!request->completed) {
190       SIMIX_cond_wait(request->cond, request->mutex);
191     }
192     if (NULL != status) {
193       status->MPI_SOURCE = request->src;
194       status->MPI_TAG = request->tag;
195       status->MPI_ERROR = MPI_SUCCESS;
196     }
197     SIMIX_mutex_unlock(request->mutex);
198   }
199
200   return retval;
201 }