Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Make sure that all the categories are connected with XBT_LOG_CONNECT() too, so that...
[simgrid.git] / src / surf / surf.c
1 /*      $Id$     */
2
3 /* Copyright (c) 2004 Arnaud Legrand. All rights reserved.                  */
4
5 /* This program is free software; you can redistribute it and/or modify it
6  * under the terms of the license (GNU LGPL) which comes with this package. */
7
8 #include <ctype.h>
9
10 #include "surf_private.h"
11 #include "xbt/module.h"
12
13 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_kernel, surf,
14                                 "Logging specific to SURF (kernel)");
15
16 int use_sdp_solver = 0;
17 int use_lagrange_solver = 0;
18
19 /* Additional declarations for Windows potability. */
20
21 #ifndef MAX_DRIVE
22 #define MAX_DRIVE 26
23 #endif
24
25 #ifdef _WIN32
26 #include <windows.h>
27 static const char *disk_drives_letter_table[MAX_DRIVE] = {
28   "A:\\",
29   "B:\\",
30   "C:\\",
31   "D:\\",
32   "E:\\",
33   "F:\\",
34   "G:\\",
35   "H:\\",
36   "I:\\",
37   "J:\\",
38   "K:\\",
39   "L:\\",
40   "M:\\",
41   "N:\\",
42   "O:\\",
43   "P:\\",
44   "Q:\\",
45   "R:\\",
46   "S:\\",
47   "T:\\",
48   "U:\\",
49   "V:\\",
50   "W:\\",
51   "X:\\",
52   "Y:\\",
53   "Z:\\"
54 };
55 #endif                          /* #ifdef _WIN32 */
56
57 /*
58  * Returns the initial path. On Windows the initial path is
59  * the current directory for the current process in the other
60  * case the function returns "./" that represents the current
61  * directory on Unix/Linux platforms.
62  */
63
64 const char *__surf_get_initial_path(void)
65 {
66
67 #ifdef _WIN32
68   unsigned i;
69   char current_directory[MAX_PATH + 1] = { 0 };
70   unsigned int len = GetCurrentDirectory(MAX_PATH + 1, current_directory);
71   char root[4] = { 0 };
72
73   if (!len)
74     return NULL;
75
76   strncpy(root, current_directory, 3);
77
78   for (i = 0; i < MAX_DRIVE; i++) {
79     if (toupper(root[0]) == disk_drives_letter_table[i][0])
80       return disk_drives_letter_table[i];
81   }
82
83   return NULL;
84 #else
85   return "./";
86 #endif
87 }
88
89 /* The __surf_is_absolute_file_path() returns 1 if
90  * file_path is a absolute file path, in the other
91  * case the function returns 0.
92  */
93 int __surf_is_absolute_file_path(const char *file_path)
94 {
95 #ifdef _WIN32
96   WIN32_FIND_DATA wfd = { 0 };
97   HANDLE hFile = FindFirstFile(file_path, &wfd);
98
99   if (INVALID_HANDLE_VALUE == hFile)
100     return 0;
101
102   FindClose(hFile);
103   return 1;
104 #else
105   return (file_path[0] == '/');
106 #endif
107 }
108
109 typedef struct surf_model_object {
110   surf_model_t model;
111 } s_surf_model_object_t, *surf_model_object_t;
112
113 static double NOW = 0;
114
115 xbt_dynar_t model_list = NULL;
116 tmgr_history_t history = NULL;
117 lmm_system_t maxmin_system = NULL;
118 xbt_dynar_t surf_path = NULL;
119 const char *surf_action_state_names[6] = {
120   "SURF_ACTION_READY",
121   "SURF_ACTION_RUNNING",
122   "SURF_ACTION_FAILED",
123   "SURF_ACTION_DONE",
124   "SURF_ACTION_TO_FREE",
125   "SURF_ACTION_NOT_IN_THE_SYSTEM"
126 };
127
128 s_surf_model_description_t surf_network_model_description[surf_network_model_description_size] = {
129   {"CM02", NULL, surf_network_model_init_CM02},
130 #ifdef HAVE_GTNETS
131   {"GTNets", NULL, surf_network_model_init_GTNETS},
132 #endif
133 #ifdef HAVE_SDP
134   {"SDP", NULL, surf_network_model_init_SDP},
135 #endif
136   {"Reno", NULL, surf_network_model_init_Reno},
137   {"Vegas", NULL, surf_network_model_init_Vegas}
138 };
139
140 s_surf_model_description_t surf_cpu_model_description[surf_cpu_model_description_size] = {
141   {"Cas01", NULL, surf_cpu_model_init_Cas01},
142 };
143
144 s_surf_model_description_t surf_workstation_model_description[surf_workstation_model_description_size] = {
145   {"CLM03", NULL, surf_workstation_model_init_CLM03, create_workstations},
146   {"compound", NULL, surf_workstation_model_init_compound, NULL},
147   {"ptask_L07", NULL, surf_workstation_model_init_ptask_L07, NULL}
148 };
149
150 void update_model_description(s_surf_model_description_t * table,
151                                  int table_size,
152                                  const char *name,
153                                  surf_model_t model)
154 {
155   int i = find_model_description(table, table_size, name);
156   table[i].model = model;
157 }
158
159 int find_model_description(s_surf_model_description_t * table,
160                               int table_size, const char *name)
161 {
162   int i;
163   char *name_list = NULL;
164
165   for (i = 0; i < table_size; i++)
166     if (!strcmp(name, table[i].name)) {
167       return i;
168     }
169   name_list = strdup(table[0].name);
170   for (i = 1; i < table_size; i++) {
171     name_list =
172         xbt_realloc(name_list,
173                     strlen(name_list) + strlen(table[i].name) + 2);
174     strcat(name_list, ", ");
175     strcat(name_list, table[i].name);
176   }
177   xbt_assert2(0, "Model '%s' is invalid! Valid models are: %s.", name,
178               name_list);
179 }
180
181 double generic_maxmin_share_resources(xbt_swag_t running_actions,
182                                       size_t offset,
183                                       lmm_system_t sys,
184                                       void (*solve) (lmm_system_t))
185 {
186   surf_action_t action = NULL;
187   double min = -1;
188   double value = -1;
189 #define VARIABLE(action) (*((lmm_variable_t*)(((char *) (action)) + (offset))))
190
191   xbt_assert0(solve, "Give me a real solver function!");
192   solve(sys);
193
194   xbt_swag_foreach(action, running_actions) {
195     value = lmm_variable_getvalue(VARIABLE(action));
196     if ((value > 0) || (action->max_duration >= 0))
197       break;
198   }
199
200   if (!action)
201     return -1.0;
202
203   if (value > 0) {
204     if(action->remains>0) 
205       min = action->remains / value;
206     else 
207       min = 0.0;
208     if ((action->max_duration >= 0) && (action->max_duration < min))
209       min = action->max_duration;
210   } else
211     min = action->max_duration;
212
213   DEBUG5("Found action (%p: duration = %f, remains = %f, value = %f) ! %f",
214          action, action->max_duration, action->remains, value, min);
215
216   for (action = xbt_swag_getNext(action, running_actions->offset);
217        action;
218        action = xbt_swag_getNext(action, running_actions->offset)) {
219     value = lmm_variable_getvalue(VARIABLE(action));
220     if (value > 0) {
221       if(action->remains>0) 
222         value = action->remains / value;
223       else 
224         value = 0.0;
225       if (value < min) {
226         min = value;
227         DEBUG2("Updating min (value) with %p: %f", action, min);
228       }
229     }
230     if ((action->max_duration >= 0) && (action->max_duration < min)) {
231       min = action->max_duration;
232       DEBUG2("Updating min (duration) with %p: %f", action, min);
233     }
234   }
235   DEBUG1("min value : %f", min);
236
237 #undef VARIABLE
238   return min;
239 }
240
241 e_surf_action_state_t surf_action_get_state(surf_action_t action)
242 {
243   surf_action_state_t action_state =
244       &(action->model_type->common_public->states);
245
246   if (action->state_set == action_state->ready_action_set)
247     return SURF_ACTION_READY;
248   if (action->state_set == action_state->running_action_set)
249     return SURF_ACTION_RUNNING;
250   if (action->state_set == action_state->failed_action_set)
251     return SURF_ACTION_FAILED;
252   if (action->state_set == action_state->done_action_set)
253     return SURF_ACTION_DONE;
254   return SURF_ACTION_NOT_IN_THE_SYSTEM;
255 }
256
257 double surf_action_get_start_time(surf_action_t action)
258 {
259   return action->start;
260 }
261
262 double surf_action_get_finish_time(surf_action_t action)
263 {
264   return action->finish;
265 }
266
267 void surf_action_free(surf_action_t * action)
268 {
269   (*action)->model_type->common_public->action_cancel(*action);
270   free(*action);
271   *action = NULL;
272 }
273
274 void surf_action_change_state(surf_action_t action,
275                               e_surf_action_state_t state)
276 {
277   surf_action_state_t action_state =
278       &(action->model_type->common_public->states);
279   XBT_IN2("(%p,%s)", action, surf_action_state_names[state]);
280   xbt_swag_remove(action, action->state_set);
281
282   if (state == SURF_ACTION_READY)
283     action->state_set = action_state->ready_action_set;
284   else if (state == SURF_ACTION_RUNNING)
285     action->state_set = action_state->running_action_set;
286   else if (state == SURF_ACTION_FAILED)
287     action->state_set = action_state->failed_action_set;
288   else if (state == SURF_ACTION_DONE)
289     action->state_set = action_state->done_action_set;
290   else
291     action->state_set = NULL;
292
293   if (action->state_set)
294     xbt_swag_insert(action, action->state_set);
295   XBT_OUT;
296 }
297
298 void surf_action_set_data(surf_action_t action, void *data)
299 {
300   action->data = data;
301 }
302
303 void surf_init(int *argc, char **argv)
304 {
305   int i, j;
306   char *opt;
307
308   const char *initial_path;
309
310   xbt_init(argc, argv);
311   if (!surf_path) {
312
313     /* retrieves the current directory of the current process */
314     initial_path = __surf_get_initial_path();
315
316     xbt_assert0((initial_path),
317                 "__surf_get_initial_path() failed! Can't resolves current Windows directory");
318
319     surf_path = xbt_dynar_new(sizeof(char *), NULL);
320     xbt_dynar_push(surf_path, &initial_path);
321
322     for (i = 1; i < *argc; i++) {
323       if (!strncmp(argv[i], "--surf-path=", strlen("--surf-path="))) {
324         opt = strchr(argv[i], '=');
325         opt++;
326         xbt_dynar_push(surf_path, &opt);
327         /*remove this from argv */
328         for (j = i + 1; j < *argc; j++) {
329           argv[j - 1] = argv[j];
330         }
331         argv[j - 1] = NULL;
332         (*argc)--;
333         i--;                    /* compensate effect of next loop incrementation */
334       }
335     }
336   }
337   if (!model_list)
338     model_list = xbt_dynar_new(sizeof(surf_model_private_t), NULL);
339   if (!history)
340     history = tmgr_history_new();
341 }
342
343 static char *path_name = NULL;
344 FILE *surf_fopen(const char *name, const char *mode)
345 {
346   unsigned int iter;
347   char *path = NULL;
348   FILE *file = NULL;
349   unsigned int path_name_len = 0;       /* don't count '\0' */
350
351   xbt_assert0(name, "Need a non-NULL file name");
352
353   xbt_assert0(surf_path,
354               "surf_init has to be called before using surf_fopen");
355
356   if (__surf_is_absolute_file_path(name)) {     /* don't mess with absolute file names */
357     return fopen(name, mode);
358
359   } else {                      /* search relative files in the path */
360
361     if (!path_name) {
362       path_name_len = strlen(name);
363       path_name = xbt_new0(char, path_name_len + 1);
364     }
365
366     xbt_dynar_foreach(surf_path, iter, path) {
367       if (path_name_len < strlen(path) + strlen(name) + 1) {
368         path_name_len = strlen(path) + strlen(name) + 1;        /* plus '/' */
369         path_name = xbt_realloc(path_name, path_name_len + 1);
370       }
371       #ifdef WIN32
372       sprintf(path_name, "%s\\%s", path, name);
373       #else
374       sprintf(path_name, "%s/%s", path, name);
375       #endif
376       file = fopen(path_name, mode);
377       if (file)
378         return file;
379     }
380   }
381   return file;
382 }
383
384 void surf_exit(void)
385 {
386   unsigned int iter;
387   surf_model_t model = NULL;
388
389   xbt_dynar_foreach(model_list, iter, model) {
390     model->common_private->finalize();
391   }
392
393   if (maxmin_system) {
394     lmm_system_free(maxmin_system);
395     maxmin_system = NULL;
396   }
397   if (history) {
398     tmgr_history_free(history);
399     history = NULL;
400   }
401   if (model_list)
402     xbt_dynar_free(&model_list);
403
404   if (surf_path)
405     xbt_dynar_free(&surf_path);
406
407   tmgr_finalize();
408   surf_parse_lex_destroy();
409   if (path_name) {
410     free(path_name);
411     path_name = NULL;
412   }
413   surf_parse_free_callbacks();
414   xbt_dict_free(&route_table);
415   NOW=0; /* Just in case the user plans to restart the simulation afterward */
416   xbt_exit();
417 }
418
419 double surf_solve(void)
420 {
421   static int first_run = 1;
422
423   double min = -1.0;
424   double next_event_date = -1.0;
425   double model_next_action_end = -1.0;
426   double value = -1.0;
427   surf_model_object_t model_obj = NULL;
428   surf_model_t model = NULL;
429   tmgr_trace_event_t event = NULL;
430   unsigned int iter;
431
432   if (first_run) {
433     DEBUG0
434         ("First Run! Let's \"purge\" events and put models in the right state");
435     while ((next_event_date = tmgr_history_next_date(history)) != -1.0) {
436       if (next_event_date > NOW)
437         break;
438       while ((event =
439               tmgr_history_get_next_event_leq(history, next_event_date,
440                                               &value,
441                                               (void **) &model_obj))) {
442         model_obj->model->common_private->
443             update_resource_state(model_obj, event, value);
444       }
445     }
446     xbt_dynar_foreach(model_list, iter, model) {
447       model->common_private->update_actions_state(NOW, 0.0);
448     }
449     first_run = 0;
450     return 0.0;
451   }
452
453   min = -1.0;
454
455   DEBUG0("Looking for next action end");
456   xbt_dynar_foreach(model_list, iter, model) {
457     DEBUG1("Running for Resource [%s]", model->common_public->name);
458     model_next_action_end =
459         model->common_private->share_resources(NOW);
460     DEBUG2("Resource [%s] : next action end = %f",
461            model->common_public->name, model_next_action_end);
462     if (((min < 0.0) || (model_next_action_end < min))
463         && (model_next_action_end >= 0.0))
464       min = model_next_action_end;
465   }
466   DEBUG1("Next action end : %f", min);
467
468   if (min < 0.0)
469     return -1.0;
470
471   DEBUG0("Looking for next event");
472   while ((next_event_date = tmgr_history_next_date(history)) != -1.0) {
473     DEBUG1("Next event : %f", next_event_date);
474     if (next_event_date > NOW + min)
475       break;
476     DEBUG0("Updating models");
477     while ((event =
478             tmgr_history_get_next_event_leq(history, next_event_date,
479                                             &value,
480                                             (void **) &model_obj))) {
481       if (model_obj->model->common_private->
482           resource_used(model_obj)) {
483         min = next_event_date - NOW;
484         DEBUG1
485             ("This event will modify model state. Next event set to %f",
486              min);
487       }
488       /* update state of model_obj according to new value. Does not touch lmm.
489          It will be modified if needed when updating actions */
490       model_obj->model->common_private->
491           update_resource_state(model_obj, event, value);
492     }
493   }
494
495   DEBUG1("Duration set to %f", min);
496
497   NOW = NOW + min;
498
499   xbt_dynar_foreach(model_list, iter, model) {
500     model->common_private->update_actions_state(NOW, min);
501   }
502
503   return min;
504 }
505
506 double surf_get_clock(void)
507 {
508   return NOW;
509 }