Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Bug fix : we could get stuck into an infinite loop for floating point precision reason.
[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 "surf_private.h"
9 #include "xbt/module.h"
10
11 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_global, surf,
12                                 "Logging specific to the SURF global module");
13
14
15 static double NOW = 0;
16
17 xbt_dynar_t resource_list = NULL;
18 tmgr_history_t history = NULL;
19 lmm_system_t maxmin_system = NULL;
20 xbt_dynar_t surf_path = NULL;
21
22 double generic_maxmin_share_resources(xbt_swag_t running_actions,
23                                        size_t offset)
24 {
25   return  generic_maxmin_share_resources2(running_actions, offset,
26                                           maxmin_system);
27 }
28
29 double generic_maxmin_share_resources2(xbt_swag_t running_actions,
30                                        size_t offset,
31                                        lmm_system_t sys)
32 {
33   surf_action_t action = NULL;
34   double min = -1;
35   double value = -1;
36 #define VARIABLE(action) (*((lmm_variable_t*)(((char *) (action)) + (offset))))
37
38   lmm_solve(sys);
39
40   xbt_swag_foreach(action, running_actions) {
41     value = lmm_variable_getvalue(VARIABLE(action));
42     if ((value > 0) || (action->max_duration >= 0))
43       break;
44   }
45
46   if (!action)
47     return -1.0;
48
49   if (value > 0) {
50     min = value = action->remains / value;
51     if ((action->max_duration >= 0) && (action->max_duration < min))
52       min = action->max_duration;
53   } else
54     min = action->max_duration;
55
56
57   for (action = xbt_swag_getNext(action, running_actions->offset);
58        action;
59        action = xbt_swag_getNext(action, running_actions->offset)) {
60     value = lmm_variable_getvalue(VARIABLE(action));
61     if (value > 0) {
62       value = action->remains / value;
63       if (value < min)
64         min = value;
65     }
66     if ((action->max_duration >= 0) && (action->max_duration < min))
67       min = action->max_duration;
68   }
69 #undef VARIABLE
70   return min;
71 }
72
73 e_surf_action_state_t surf_action_get_state(surf_action_t action)
74 {
75   surf_action_state_t action_state =
76       &(action->resource_type->common_public->states);
77
78   if (action->state_set == action_state->ready_action_set)
79     return SURF_ACTION_READY;
80   if (action->state_set == action_state->running_action_set)
81     return SURF_ACTION_RUNNING;
82   if (action->state_set == action_state->failed_action_set)
83     return SURF_ACTION_FAILED;
84   if (action->state_set == action_state->done_action_set)
85     return SURF_ACTION_DONE;
86   return SURF_ACTION_NOT_IN_THE_SYSTEM;
87 }
88
89 void surf_action_free(surf_action_t * action)
90 {
91   (*action)->resource_type->common_public->action_cancel(*action);
92   xbt_free(*action);
93   *action = NULL;
94 }
95
96 void surf_action_change_state(surf_action_t action,
97                               e_surf_action_state_t state)
98 {
99   surf_action_state_t action_state =
100       &(action->resource_type->common_public->states);
101
102   xbt_swag_remove(action, action->state_set);
103
104   if (state == SURF_ACTION_READY)
105     action->state_set = action_state->ready_action_set;
106   else if (state == SURF_ACTION_RUNNING)
107     action->state_set = action_state->running_action_set;
108   else if (state == SURF_ACTION_FAILED)
109     action->state_set = action_state->failed_action_set;
110   else if (state == SURF_ACTION_DONE)
111     action->state_set = action_state->done_action_set;
112   else
113     action->state_set = NULL;
114
115   if (action->state_set)
116     xbt_swag_insert(action, action->state_set);
117 }
118
119 void surf_action_set_data(surf_action_t action,
120                           void *data)
121 {
122   action->data=data;
123 }
124
125 void surf_init(int *argc, char **argv)
126 {
127   int i,j;
128   char *opt;
129
130   xbt_init(argc, argv);
131   if (!surf_path) {
132     const char *initial_path = "./";
133     surf_path = xbt_dynar_new(sizeof(char*), NULL);
134     xbt_dynar_push(surf_path,&initial_path);
135
136     for (i=1; i<*argc; i++) {
137       if (!strncmp(argv[i],"--surf-path=",strlen("--surf-path="))) {
138         opt=strchr(argv[i],'=');
139         opt++;
140         xbt_dynar_push(surf_path,&opt);
141         /*remove this from argv*/
142         for (j=i+1; j<*argc; j++) {
143           argv[j-1] = argv[j];
144         } 
145         argv[j-1] = NULL;
146         (*argc)--;
147         i--; /* compensate effect of next loop incrementation */
148       }
149     }
150   }
151   if (!resource_list)
152     resource_list = xbt_dynar_new(sizeof(surf_resource_private_t), NULL);
153   if (!history)
154     history = tmgr_history_new();
155   if (!maxmin_system)
156     maxmin_system = lmm_system_new();
157 }
158
159 static char* path_name = NULL;
160 FILE *surf_fopen(const char *name, const char *mode)
161 {
162   int i; 
163   char* path = NULL;
164   FILE *file = NULL;
165
166   xbt_assert0(surf_path,"surf_init has to be called before using surf_fopen");
167   if(!path_name) path_name=xbt_new0(char,strlen(name)+1);
168
169   xbt_dynar_foreach(surf_path,i,path) {
170     if(strlen(path_name)<strlen(path)+strlen(name)+2) 
171       path_name=xbt_realloc(path_name,strlen(path)+strlen(name)+2);
172     strcpy(path_name, path);
173     strcat(path_name,"/");
174     strcat(path_name,name);
175     file = fopen(path_name,mode);
176     if(file) return file;
177   }
178   return file;
179 }
180
181 void surf_finalize(void)
182 {
183   int i;
184   surf_resource_t resource = NULL;
185
186   xbt_dynar_foreach(resource_list, i, resource) {
187     resource->common_private->finalize();
188   }
189
190   if (maxmin_system) {
191     lmm_system_free(maxmin_system);
192     maxmin_system = NULL;
193   }
194   if (history) {
195     tmgr_history_free(history);
196     history = NULL;
197   }
198   if (resource_list)
199     xbt_dynar_free(&resource_list);
200
201   if(surf_path) 
202     xbt_dynar_free(&surf_path);
203
204   tmgr_finalize();
205   surf_parse_lex_destroy();
206   if(path_name) {
207     xbt_free(path_name);
208     path_name = NULL;
209   }
210 }
211
212 double surf_solve(void)
213 {
214   static int first_run = 1;
215
216   double min = -1.0;
217   double next_event_date = -1.0;
218   double resource_next_action_end = -1.0;
219   double value = -1.0;
220   surf_resource_object_t resource_obj = NULL;
221   surf_resource_t resource = NULL;
222   tmgr_trace_event_t event = NULL;
223   int i;
224
225   if (first_run) {
226     while ((next_event_date = tmgr_history_next_date(history)) != -1.0) {
227       if (next_event_date > NOW)
228         break;
229       while ((event =
230               tmgr_history_get_next_event_leq(history, next_event_date,
231                                               &value,
232                                               (void **) &resource_obj))) {
233         resource_obj->resource->common_private->
234             update_resource_state(resource_obj, event, value);
235       }
236     }
237     xbt_dynar_foreach(resource_list, i, resource) {
238       resource->common_private->update_actions_state(NOW, 0.0);
239     }
240     first_run = 0;
241     return 0.0;
242   }
243
244   min = -1.0;
245
246   xbt_dynar_foreach(resource_list, i, resource) {
247     resource_next_action_end =
248         resource->common_private->share_resources(NOW);
249     if (((min < 0.0) || (resource_next_action_end < min))
250         && (resource_next_action_end >= 0.0))
251       min = resource_next_action_end;
252   }
253
254   if (min < 0.0)
255     return -1.0;
256
257   while ((next_event_date = tmgr_history_next_date(history)) != -1.0) {
258     if (next_event_date > NOW + min)
259       break;
260     while ((event =
261             tmgr_history_get_next_event_leq(history, next_event_date,
262                                             &value,
263                                             (void **) &resource_obj))) {
264       if (resource_obj->resource->common_private->
265           resource_used(resource_obj)) {
266         min = next_event_date - NOW;
267       }
268       /* update state of resource_obj according to new value. Does not touch lmm.
269          It will be modified if needed when updating actions */
270       resource_obj->resource->common_private->
271           update_resource_state(resource_obj, event, value);
272     }
273   }
274
275
276   xbt_dynar_foreach(resource_list, i, resource) {
277     resource->common_private->update_actions_state(NOW, min);
278   }
279
280   NOW = NOW + min;
281
282   return min;
283 }
284
285 double surf_get_clock(void)
286 {
287   return NOW;
288 }