Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
fc9d793303836bbe41fab59313f117b7199a0d16
[simgrid.git] / src / smpi / smpi_global.c
1 #include <stdio.h>
2
3 #include "private.h"
4
5 XBT_LOG_NEW_CATEGORY(smpi, "All SMPI categories");
6
7 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(smpi_kernel, smpi,
8                                 "Logging specific to SMPI (kernel)");
9
10 smpi_global_t smpi_global = NULL;
11
12 void *smpi_request_new(void);
13
14 void *smpi_request_new()
15 {
16   smpi_mpi_request_t request = xbt_new(s_smpi_mpi_request_t, 1);
17
18   request->buf = NULL;
19   request->completed = 0;
20   request->mutex = SIMIX_mutex_init();
21   request->cond = SIMIX_cond_init();
22   request->data = NULL;
23   request->forward = 0;
24
25   return request;
26 }
27
28 void smpi_request_free(void *pointer);
29
30 void smpi_request_free(void *pointer)
31 {
32
33   smpi_mpi_request_t request = pointer;
34
35   SIMIX_cond_destroy(request->cond);
36   SIMIX_mutex_destroy(request->mutex);
37   xbt_free(request);
38
39   return;
40 }
41
42 void smpi_request_reset(void *pointer);
43
44 void smpi_request_reset(void *pointer)
45 {
46   smpi_mpi_request_t request = pointer;
47
48   request->buf = NULL;
49   request->completed = 0;
50   request->data = NULL;
51   request->forward = 0;
52
53   return;
54 }
55
56
57 void *smpi_message_new(void);
58
59 void *smpi_message_new()
60 {
61   smpi_received_message_t message = xbt_new(s_smpi_received_message_t, 1);
62   message->buf = NULL;
63   return message;
64 }
65
66 void smpi_message_free(void *pointer);
67
68 void smpi_message_free(void *pointer)
69 {
70   xbt_free(pointer);
71   return;
72 }
73
74 void smpi_message_reset(void *pointer);
75
76 void smpi_message_reset(void *pointer)
77 {
78   smpi_received_message_t message = pointer;
79   message->buf = NULL;
80   return;
81 }
82
83 int smpi_create_request(void *buf, int count, smpi_mpi_datatype_t datatype,
84                         int src, int dst, int tag,
85                         smpi_mpi_communicator_t comm,
86                         smpi_mpi_request_t * requestptr)
87 {
88   int retval = MPI_SUCCESS;
89
90   smpi_mpi_request_t request = NULL;
91
92   // parameter checking prob belongs in smpi_mpi, but this is less repeat code
93   if (NULL == buf) {
94     retval = MPI_ERR_INTERN;
95   } else if (0 > count) {
96     retval = MPI_ERR_COUNT;
97   } else if (NULL == datatype) {
98     retval = MPI_ERR_TYPE;
99   } else if (MPI_ANY_SOURCE != src && (0 > src || comm->size <= src)) {
100     retval = MPI_ERR_RANK;
101   } else if (0 > dst || comm->size <= dst) {
102     retval = MPI_ERR_RANK;
103   } else if (MPI_ANY_TAG != tag && 0 > tag) {
104     retval = MPI_ERR_TAG;
105   } else if (NULL == comm) {
106     retval = MPI_ERR_COMM;
107   } else if (NULL == requestptr) {
108     retval = MPI_ERR_ARG;
109   } else {
110     request = xbt_mallocator_get(smpi_global->request_mallocator);
111     request->comm = comm;
112     request->src = src;
113     request->dst = dst;
114     request->tag = tag;
115     request->buf = buf;
116     request->datatype = datatype;
117     request->count = count;
118
119     *requestptr = request;
120   }
121   return retval;
122 }
123 /* FIXME: understand what they do and put the prototypes in a header file (live in smpi_base.c) */
124 void smpi_mpi_land_func(void *a, void *b, int *length,
125                         MPI_Datatype * datatype);
126 void smpi_mpi_sum_func(void *a, void *b, int *length,
127                        MPI_Datatype * datatype);
128
129 void smpi_global_init()
130 {
131   int i;
132
133   /* Connect our log channels: that must be done manually under windows */
134 #ifdef XBT_LOG_CONNECT
135   XBT_LOG_CONNECT(smpi_base, smpi);
136   XBT_LOG_CONNECT(smpi_bench, smpi);
137   XBT_LOG_CONNECT(smpi_kernel, smpi);
138   XBT_LOG_CONNECT(smpi_mpi, smpi);
139   XBT_LOG_CONNECT(smpi_receiver, smpi);
140   XBT_LOG_CONNECT(smpi_sender, smpi);
141   XBT_LOG_CONNECT(smpi_util, smpi);
142 #endif
143
144   smpi_global = xbt_new(s_smpi_global_t, 1);
145   // config variable
146   smpi_global->reference_speed = SMPI_DEFAULT_SPEED;
147
148   // mallocators
149   smpi_global->request_mallocator =
150     xbt_mallocator_new(SMPI_REQUEST_MALLOCATOR_SIZE, smpi_request_new,
151                        smpi_request_free, smpi_request_reset);
152   smpi_global->message_mallocator =
153     xbt_mallocator_new(SMPI_MESSAGE_MALLOCATOR_SIZE, smpi_message_new,
154                        smpi_message_free, smpi_message_reset);
155
156   smpi_global->process_count = SIMIX_process_count();
157   fprintf(stderr,"There is %d processes\n",smpi_global->process_count);
158
159   // sender/receiver processes
160   smpi_global->main_processes = xbt_new(smx_process_t, smpi_global->process_count);
161
162   // timers
163   smpi_global->timer = xbt_os_timer_new();
164   smpi_global->timer_mutex = SIMIX_mutex_init();
165   smpi_global->timer_cond = SIMIX_cond_init();
166
167   smpi_global->do_once_duration_nodes = NULL;
168   smpi_global->do_once_duration = NULL;
169   smpi_global->do_once_mutex = SIMIX_mutex_init();
170
171
172   smpi_mpi_global = xbt_new(s_smpi_mpi_global_t, 1);
173
174   // global communicator
175   smpi_mpi_global->mpi_comm_world = xbt_new(s_smpi_mpi_communicator_t, 1);
176   smpi_mpi_global->mpi_comm_world->size = smpi_global->process_count;
177   smpi_mpi_global->mpi_comm_world->barrier_count = 0;
178   smpi_mpi_global->mpi_comm_world->barrier_mutex = SIMIX_mutex_init();
179   smpi_mpi_global->mpi_comm_world->barrier_cond = SIMIX_cond_init();
180   smpi_mpi_global->mpi_comm_world->rank_to_index_map =
181     xbt_new(int, smpi_global->process_count);
182   smpi_mpi_global->mpi_comm_world->index_to_rank_map =
183     xbt_new(int, smpi_global->process_count);
184   for (i = 0; i < smpi_global->process_count; i++) {
185     smpi_mpi_global->mpi_comm_world->rank_to_index_map[i] = i;
186     smpi_mpi_global->mpi_comm_world->index_to_rank_map[i] = i;
187   }
188
189   // mpi datatypes
190   smpi_mpi_global->mpi_byte = xbt_new(s_smpi_mpi_datatype_t, 1);
191   smpi_mpi_global->mpi_byte->size = (size_t) 1;
192   smpi_mpi_global->mpi_int = xbt_new(s_smpi_mpi_datatype_t, 1);
193   smpi_mpi_global->mpi_int->size = sizeof(int);
194   smpi_mpi_global->mpi_double = xbt_new(s_smpi_mpi_datatype_t, 1);
195   smpi_mpi_global->mpi_double->size = sizeof(double);
196
197   // mpi operations
198   smpi_mpi_global->mpi_land = xbt_new(s_smpi_mpi_op_t, 1);
199   smpi_mpi_global->mpi_land->func = smpi_mpi_land_func;
200   smpi_mpi_global->mpi_sum = xbt_new(s_smpi_mpi_op_t, 1);
201   smpi_mpi_global->mpi_sum->func = smpi_mpi_sum_func;
202
203 }
204
205 void smpi_global_destroy()
206 {
207   smpi_do_once_duration_node_t curr, next;
208
209   // processes
210   xbt_free(smpi_global->main_processes);
211
212   // mallocators
213   xbt_mallocator_free(smpi_global->request_mallocator);
214   xbt_mallocator_free(smpi_global->message_mallocator);
215
216   xbt_os_timer_free(smpi_global->timer);
217   SIMIX_mutex_destroy(smpi_global->timer_mutex);
218   SIMIX_cond_destroy(smpi_global->timer_cond);
219
220   for (curr = smpi_global->do_once_duration_nodes; NULL != curr; curr = next) {
221     next = curr->next;
222     xbt_free(curr->file);
223     xbt_free(curr);
224   }
225
226   SIMIX_mutex_destroy(smpi_global->do_once_mutex);
227
228   xbt_free(smpi_global);
229   smpi_global = NULL;
230
231   /* free smpi_mpi_global */
232   SIMIX_mutex_destroy(smpi_mpi_global->mpi_comm_world->barrier_mutex);
233   SIMIX_cond_destroy(smpi_mpi_global->mpi_comm_world->barrier_cond);
234   xbt_free(smpi_mpi_global->mpi_comm_world);
235
236   xbt_free(smpi_mpi_global->mpi_byte);
237   xbt_free(smpi_mpi_global->mpi_int);
238   xbt_free(smpi_mpi_global->mpi_double);
239
240   xbt_free(smpi_mpi_global->mpi_land);
241   xbt_free(smpi_mpi_global->mpi_sum);
242
243   xbt_free(smpi_mpi_global);
244
245 }
246
247 int smpi_process_index()
248 {
249   smpi_process_data_t pdata = (smpi_process_data_t) SIMIX_process_get_data(SIMIX_process_self());
250   return pdata->index;
251 }
252
253 smx_mutex_t smpi_process_mutex()
254 {
255   smpi_process_data_t pdata = (smpi_process_data_t) SIMIX_process_get_data(SIMIX_process_self());
256   return pdata->mutex;
257 }
258
259 smx_cond_t smpi_process_cond()
260 {
261   smpi_process_data_t pdata = (smpi_process_data_t) SIMIX_process_get_data(SIMIX_process_self());
262   return pdata->cond;
263 }
264
265 int smpi_run_simulation(int *argc, char **argv)
266 {
267   smx_cond_t cond = NULL;
268   smx_action_t action = NULL;
269
270   xbt_fifo_t actions_failed = xbt_fifo_new();
271   xbt_fifo_t actions_done = xbt_fifo_new();
272
273   srand(SMPI_RAND_SEED);
274
275   SIMIX_global_init(argc, argv);
276
277   // parse the platform file: get the host list
278   SIMIX_create_environment(argv[1]);
279
280   SIMIX_function_register("smpi_simulated_main", smpi_simulated_main);
281   SIMIX_launch_application(argv[2]);
282
283   // must initialize globals between creating environment and launching app....
284   smpi_global_init();
285
286   /* Prepare to display some more info when dying on Ctrl-C pressing */
287   // FIXME: doesn't work
288   //signal(SIGINT, inthandler);
289
290   /* Clean IO before the run */
291   fflush(stdout);
292   fflush(stderr);
293   SIMIX_init();
294
295   while (SIMIX_solve(actions_done, actions_failed) != -1.0) {
296     while ((action = xbt_fifo_pop(actions_failed))) {
297       DEBUG1("** %s failed **", action->name);
298       while ((cond = xbt_fifo_pop(action->cond_list))) {
299         SIMIX_cond_broadcast(cond);
300       }
301     }
302     while ((action = xbt_fifo_pop(actions_done))) {
303       DEBUG1("** %s done **", action->name);
304       while ((cond = xbt_fifo_pop(action->cond_list))) {
305         SIMIX_cond_broadcast(cond);
306       }
307     }
308   }
309
310   // FIXME: cleanup incomplete
311   xbt_fifo_free(actions_failed);
312   xbt_fifo_free(actions_done);
313
314   INFO1("simulation time %g", SIMIX_get_clock());
315
316   smpi_global_destroy();
317
318   SIMIX_clean();
319
320   return 0;
321 }