Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
attempt to replace internal dynar by std::vector
[simgrid.git] / src / simdag / sd_task.cpp
1 /* Copyright (c) 2006-2016. The SimGrid Team.
2  * All rights reserved.                                                     */
3
4 /* This program is free software; you can redistribute it and/or modify it
5  * under the terms of the license (GNU LGPL) which comes with this package. */
6
7 #include "src/surf/HostImpl.hpp"
8 #include "src/surf/surf_interface.hpp"
9 #include "src/simdag/simdag_private.h"
10 #include "simgrid/simdag.h"
11 #include "src/instr/instr_private.h"
12
13 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(sd_task, sd, "Logging specific to SimDag (task)");
14
15 /* Destroys a dependency between two tasks. */
16 static void __SD_task_dependency_destroy(void *dependency)
17 {
18   xbt_free(((SD_dependency_t)dependency)->name);
19   xbt_free(dependency);
20 }
21
22 /* Remove all dependencies associated with a task. This function is called when the task is destroyed. */
23 static void __SD_task_remove_dependencies(SD_task_t task)
24 {
25   /* we must destroy the dependencies carefully (with SD_dependency_remove) because each one is stored twice */
26   SD_dependency_t dependency;
27   while (xbt_dynar_is_empty(task->tasks_before) == 0) {
28     xbt_dynar_get_cpy(task->tasks_before, 0, &dependency);
29     SD_task_dependency_remove(dependency->src, dependency->dst);
30   }
31
32   while (xbt_dynar_is_empty(task->tasks_after) == 0) {
33     xbt_dynar_get_cpy(task->tasks_after, 0, &dependency);
34     SD_task_dependency_remove(dependency->src, dependency->dst);
35   }
36 }
37
38 /* Destroys the data memorized by SD_task_schedule. Task state must be SD_SCHEDULED or SD_RUNNABLE. */
39 static void __SD_task_destroy_scheduling_data(SD_task_t task)
40 {
41   if (task->state != SD_SCHEDULED && task->state != SD_RUNNABLE)
42     THROWF(arg_error, 0, "Task '%s' must be SD_SCHEDULED or SD_RUNNABLE", SD_task_get_name(task));
43
44   xbt_free(task->flops_amount);
45   xbt_free(task->bytes_amount);
46   task->flops_amount = nullptr;
47   task->bytes_amount = nullptr;
48 }
49
50 void* SD_task_new_f()
51 {
52   SD_task_t task = xbt_new0(s_SD_task_t, 1);
53   task->tasks_before = xbt_dynar_new(sizeof(SD_dependency_t), nullptr);
54   task->tasks_after = xbt_dynar_new(sizeof(SD_dependency_t), nullptr);
55
56   return task;
57 }
58
59 void SD_task_recycle_f(void *t)
60 {
61   SD_task_t task = static_cast<SD_task_t>(t);
62
63   /* Reset the content */
64   task->kind = SD_TASK_NOT_TYPED;
65   task->state= SD_NOT_SCHEDULED;
66   sd_global->initial_task_set->push_back(task);
67
68   task->marked = 0;
69
70   task->start_time = -1.0;
71   task->finish_time = -1.0;
72   task->surf_action = nullptr;
73   task->watch_points = 0;
74
75   /* dependencies */
76   xbt_dynar_reset(task->tasks_before);
77   xbt_dynar_reset(task->tasks_after);
78   task->unsatisfied_dependencies = 0;
79   task->is_not_ready = 0;
80
81   /* scheduling parameters */
82   task->host_count = 0;
83   task->host_list = nullptr;
84   task->flops_amount = nullptr;
85   task->bytes_amount = nullptr;
86   task->rate = -1;
87 }
88
89 void SD_task_free_f(void *t)
90 {
91   SD_task_t task = static_cast<SD_task_t>(t);
92
93   xbt_dynar_free(&task->tasks_before);
94   xbt_dynar_free(&task->tasks_after);
95   xbt_free(task);
96 }
97
98 /**
99  * \brief Creates a new task.
100  *
101  * \param name the name of the task (can be \c nullptr)
102  * \param data the user data you want to associate with the task (can be \c nullptr)
103  * \param amount amount of the task
104  * \return the new task
105  * \see SD_task_destroy()
106  */
107 SD_task_t SD_task_create(const char *name, void *data, double amount)
108 {
109   SD_task_t task = static_cast<SD_task_t>(xbt_mallocator_get(sd_global->task_mallocator));
110
111   /* general information */
112   task->data = data;            /* user data */
113   task->name = xbt_strdup(name);
114   task->amount = amount;
115   task->remains = amount;
116
117   return task;
118 }
119
120 static inline SD_task_t SD_task_create_sized(const char *name, void *data, double amount, int ws_count)
121 {
122   SD_task_t task = SD_task_create(name, data, amount);
123   task->bytes_amount = xbt_new0(double, ws_count * ws_count);
124   task->flops_amount = xbt_new0(double, ws_count);
125   task->host_count = ws_count;
126   task->host_list = xbt_new0(sg_host_t, ws_count);
127   return task;
128 }
129
130 /** @brief create a end-to-end communication task that can then be auto-scheduled
131  *
132  * Auto-scheduling mean that the task can be used with SD_task_schedulev(). This allows to specify the task costs at
133  * creation, and decouple them from the scheduling process where you just specify which resource should deliver the
134  * mandatory power.
135  *
136  * A end-to-end communication must be scheduled on 2 hosts, and the amount specified at creation is sent from hosts[0]
137  * to hosts[1].
138  */
139 SD_task_t SD_task_create_comm_e2e(const char *name, void *data, double amount)
140 {
141   SD_task_t res = SD_task_create_sized(name, data, amount, 2);
142   res->bytes_amount[2] = amount;
143   res->kind = SD_TASK_COMM_E2E;
144
145   return res;
146 }
147
148 /** @brief create a sequential computation task that can then be auto-scheduled
149  *
150  * Auto-scheduling mean that the task can be used with SD_task_schedulev(). This allows to specify the task costs at
151  * creation, and decouple them from the scheduling process where you just specify which resource should deliver the
152  * mandatory power.
153  *
154  * A sequential computation must be scheduled on 1 host, and the amount specified at creation to be run on hosts[0].
155  *
156  * \param name the name of the task (can be \c nullptr)
157  * \param data the user data you want to associate with the task (can be \c nullptr)
158  * \param flops_amount amount of compute work to be done by the task
159  * \return the new SD_TASK_COMP_SEQ typed task
160  */
161 SD_task_t SD_task_create_comp_seq(const char *name, void *data, double flops_amount)
162 {
163   SD_task_t res = SD_task_create_sized(name, data, flops_amount, 1);
164   res->flops_amount[0] = flops_amount;
165   res->kind = SD_TASK_COMP_SEQ;
166
167   return res;
168 }
169
170 /** @brief create a parallel computation task that can then be auto-scheduled
171  *
172  * Auto-scheduling mean that the task can be used with SD_task_schedulev(). This allows to specify the task costs at
173  * creation, and decouple them from the scheduling process where you just specify which resource should deliver the
174  * mandatory power.
175  *
176  * A parallel computation can be scheduled on any number of host.
177  * The underlying speedup model is Amdahl's law.
178  * To be auto-scheduled, \see SD_task_distribute_comp_amdahl has to be called first.
179  * \param name the name of the task (can be \c nullptr)
180  * \param data the user data you want to associate with the task (can be \c nullptr)
181  * \param flops_amount amount of compute work to be done by the task
182  * \param alpha purely serial fraction of the work to be done (in [0.;1.[)
183  * \return the new task
184  */
185 SD_task_t SD_task_create_comp_par_amdahl(const char *name, void *data, double flops_amount, double alpha)
186 {
187   xbt_assert(alpha < 1. && alpha >= 0., "Invalid parameter: alpha must be in [0.;1.[");
188
189   SD_task_t res = SD_task_create(name, data, flops_amount);
190   res->alpha = alpha;
191   res->kind = SD_TASK_COMP_PAR_AMDAHL;
192
193   return res;
194 }
195
196 /** @brief create a complex data redistribution task that can then be  auto-scheduled
197  *
198  * Auto-scheduling mean that the task can be used with SD_task_schedulev(). 
199  * This allows to specify the task costs at creation, and decouple them from the scheduling process where you just
200  * specify which resource should communicate.
201  *
202  * A data redistribution can be scheduled on any number of host.
203  * The assumed distribution is a 1D block distribution. Each host owns the same share of the \see amount.
204  * To be auto-scheduled, \see SD_task_distribute_comm_mxn_1d_block has to be  called first.
205  * \param name the name of the task (can be \c nullptr)
206  * \param data the user data you want to associate with the task (can be \c nullptr)
207  * \param amount amount of data to redistribute by the task
208  * \return the new task
209  */
210 SD_task_t SD_task_create_comm_par_mxn_1d_block(const char *name, void *data, double amount)
211 {
212   SD_task_t res = SD_task_create(name, data, amount);
213   res->host_list=nullptr;
214   res->kind = SD_TASK_COMM_PAR_MXN_1D_BLOCK;
215
216   return res;
217 }
218
219 /**
220  * \brief Destroys a task.
221  *
222  * The user data (if any) should have been destroyed first.
223  *
224  * \param task the task you want to destroy
225  * \see SD_task_create()
226  */
227 void SD_task_destroy(SD_task_t task)
228 {
229   XBT_DEBUG("Destroying task %s...", SD_task_get_name(task));
230
231   __SD_task_remove_dependencies(task);
232
233   if (task->state == SD_SCHEDULED || task->state == SD_RUNNABLE)
234     __SD_task_destroy_scheduling_data(task);
235
236   int idx = xbt_dynar_search_or_negative(sd_global->return_set, &task);
237   if (idx >=0) {
238     xbt_dynar_remove_at(sd_global->return_set, idx, nullptr);
239   }
240
241   xbt_free(task->name);
242
243   if (task->surf_action != nullptr)
244     task->surf_action->unref();
245
246   xbt_free(task->host_list);
247   xbt_free(task->bytes_amount);
248   xbt_free(task->flops_amount);
249
250   xbt_mallocator_release(sd_global->task_mallocator,task);
251
252   XBT_DEBUG("Task destroyed.");
253 }
254
255 /**
256  * \brief Returns the user data of a task
257  *
258  * \param task a task
259  * \return the user data associated with this task (can be \c nullptr)
260  * \see SD_task_set_data()
261  */
262 void *SD_task_get_data(SD_task_t task)
263 {
264   return task->data;
265 }
266
267 /**
268  * \brief Sets the user data of a task
269  *
270  * The new data can be \c nullptr. The old data should have been freed first
271  * if it was not \c nullptr.
272  *
273  * \param task a task
274  * \param data the new data you want to associate with this task
275  * \see SD_task_get_data()
276  */
277 void SD_task_set_data(SD_task_t task, void *data)
278 {
279   task->data = data;
280 }
281
282 /**
283  * \brief Sets the rate of a task
284  *
285  * This will change the network bandwidth a task can use. This rate  cannot be dynamically changed. Once the task has
286  * started, this call is ineffective. This rate depends on both the nominal bandwidth on the route onto which the task
287  * is scheduled (\see SD_task_get_current_bandwidth) and the amount of data to transfer.
288  *
289  * To divide the nominal bandwidth by 2, the rate then has to be :
290  *    rate = bandwidth/(2*amount)
291  *
292  * \param task a \see SD_TASK_COMM_E2E task (end-to-end communication)
293  * \param rate the new rate you want to associate with this task.
294  */
295 void SD_task_set_rate(SD_task_t task, double rate)
296 {
297   xbt_assert(task->kind == SD_TASK_COMM_E2E, "The rate can be modified for end-to-end communications only.");
298   if(task->start_time<0) {
299     task->rate = rate;
300   } else {
301     XBT_WARN("Task %p has started. Changing rate is ineffective.", task);
302   }
303 }
304
305 /**
306  * \brief Returns the state of a task
307  *
308  * \param task a task
309  * \return the current \ref e_SD_task_state_t "state" of this task:
310  * #SD_NOT_SCHEDULED, #SD_SCHEDULED, #SD_RUNNABLE, #SD_RUNNING, #SD_DONE or #SD_FAILED
311  * \see e_SD_task_state_t
312  */
313 e_SD_task_state_t SD_task_get_state(SD_task_t task)
314 {
315   return task->state;
316 }
317
318 /* Changes the state of a task. Updates the sd_global->watch_point_reached flag.
319  */
320 void SD_task_set_state(SD_task_t task, e_SD_task_state_t new_state)
321 {
322   std::vector<SD_task_t>::iterator idx;
323   switch (new_state) {
324   case SD_NOT_SCHEDULED:
325   case SD_SCHEDULABLE:
326     if (SD_task_get_state(task) == SD_FAILED){
327       sd_global->initial_task_set->push_back(task);
328       sd_global->completed_task_set->erase(std::remove(sd_global->completed_task_set->begin(),
329                                                        sd_global->completed_task_set->end(), task),
330                                                        sd_global->completed_task_set->end());
331     }
332     break;
333   case SD_SCHEDULED:
334     if (SD_task_get_state(task) == SD_RUNNABLE){
335       sd_global->initial_task_set->push_back(task);
336       sd_global->executable_task_set->erase(std::remove(sd_global->executable_task_set->begin(),
337                                                         sd_global->executable_task_set->end(), task),
338                                                         sd_global->executable_task_set->end());
339     }
340     break;
341   case SD_RUNNABLE:
342     idx = std::find(sd_global->initial_task_set->begin(), sd_global->initial_task_set->end(), task);
343     if (idx != sd_global->initial_task_set->end()) {
344       sd_global->executable_task_set->push_back(*idx);
345       sd_global->initial_task_set->erase(idx);
346     }
347     break;
348   case SD_RUNNING:
349     sd_global->executable_task_set->erase(std::remove(sd_global->executable_task_set->begin(),
350                                                       sd_global->executable_task_set->end(), task),
351                                                       sd_global->executable_task_set->end());
352     break;
353   case SD_DONE:
354     sd_global->completed_task_set->push_back(task);
355     task->finish_time = task->surf_action->getFinishTime();
356     task->remains = 0;
357 #if HAVE_JEDULE
358     jedule_log_sd_event(task);
359 #endif
360     break;
361   case SD_FAILED:
362     sd_global->completed_task_set->push_back(task);
363     break;
364   default:
365     xbt_die( "Invalid state");
366   }
367
368   task->state = new_state;
369
370   if (task->watch_points & new_state) {
371     XBT_VERB("Watch point reached with task '%s'!", SD_task_get_name(task));
372     sd_global->watch_point_reached = 1;
373     SD_task_unwatch(task, new_state);   /* remove the watch point */
374   }
375 }
376
377 /**
378  * \brief Returns the name of a task
379  *
380  * \param task a task
381  * \return the name of this task (can be \c nullptr)
382  */
383 const char *SD_task_get_name(SD_task_t task)
384 {
385   return task->name;
386 }
387
388 /** @brief Allows to change the name of a task */
389 void SD_task_set_name(SD_task_t task, const char *name)
390 {
391   xbt_free(task->name);
392   task->name = xbt_strdup(name);
393 }
394
395 /** @brief Returns the dynar of the parents of a task
396  *
397  * \param task a task
398  * \return a newly allocated dynar comprising the parents of this task
399  */
400
401 xbt_dynar_t SD_task_get_parents(SD_task_t task)
402 {
403   unsigned int i;
404   SD_dependency_t dep;
405
406   xbt_dynar_t parents = xbt_dynar_new(sizeof(SD_task_t), nullptr);
407   xbt_dynar_foreach(task->tasks_before, i, dep) {
408     xbt_dynar_push(parents, &(dep->src));
409   }
410   return parents;
411 }
412
413 /** @brief Returns the dynar of the parents of a task
414  *
415  * \param task a task
416  * \return a newly allocated dynar comprising the parents of this task
417  */
418 xbt_dynar_t SD_task_get_children(SD_task_t task)
419 {
420   unsigned int i;
421   SD_dependency_t dep;
422
423   xbt_dynar_t children = xbt_dynar_new(sizeof(SD_task_t), nullptr);
424   xbt_dynar_foreach(task->tasks_after, i, dep) {
425     xbt_dynar_push(children, &(dep->dst));
426   }
427   return children;
428 }
429
430 /**
431  * \brief Returns the number of workstations involved in a task
432  *
433  * Only call this on already scheduled tasks!
434  * \param task a task
435  */
436 int SD_task_get_workstation_count(SD_task_t task)
437 {
438   return task->host_count;
439 }
440
441 /**
442  * \brief Returns the list of workstations involved in a task
443  *
444  * Only call this on already scheduled tasks!
445  * \param task a task
446  */
447 sg_host_t *SD_task_get_workstation_list(SD_task_t task)
448 {
449   return task->host_list;
450 }
451
452 /**
453  * \brief Returns the total amount of work contained in a task
454  *
455  * \param task a task
456  * \return the total amount of work (computation or data transfer) for this task
457  * \see SD_task_get_remaining_amount()
458  */
459 double SD_task_get_amount(SD_task_t task)
460 {
461   return task->amount;
462 }
463
464 /** @brief Sets the total amount of work of a task
465  * For sequential typed tasks (COMP_SEQ and COMM_E2E), it also sets the appropriate values in the flops_amount and
466  * bytes_amount arrays respectively. Nothing more than modifying task->amount is done for parallel  typed tasks
467  * (COMP_PAR_AMDAHL and COMM_PAR_MXN_1D_BLOCK) as the distribution of the amount of work is done at scheduling time.
468  *
469  * \param task a task
470  * \param amount the new amount of work to execute
471  */
472 void SD_task_set_amount(SD_task_t task, double amount)
473 {
474   task->amount = amount;
475   if (task->kind == SD_TASK_COMP_SEQ)
476     task->flops_amount[0] = amount;
477   if (task->kind == SD_TASK_COMM_E2E)
478     task->bytes_amount[2] = amount;
479 }
480
481 /**
482  * \brief Returns the alpha parameter of a SD_TASK_COMP_PAR_AMDAHL task
483  *
484  * \param task a parallel task assuming Amdahl's law as speedup model
485  * \return the alpha parameter (serial part of a task in percent) for this task
486  */
487 double SD_task_get_alpha(SD_task_t task)
488 {
489   xbt_assert(SD_task_get_kind(task) == SD_TASK_COMP_PAR_AMDAHL, "Alpha parameter is not defined for this kind of task");
490   return task->alpha;
491 }
492
493 /**
494  * \brief Returns the remaining amount work to do till the completion of a task
495  *
496  * \param task a task
497  * \return the remaining amount of work (computation or data transfer) of this task
498  * \see SD_task_get_amount()
499  */
500 double SD_task_get_remaining_amount(SD_task_t task)
501 {
502   if (task->surf_action)
503     return task->surf_action->getRemains();
504   else
505     return task->remains;
506 }
507
508 e_SD_task_kind_t SD_task_get_kind(SD_task_t task)
509 {
510   return task->kind;
511 }
512
513 /** @brief Displays debugging information about a task */
514 void SD_task_dump(SD_task_t task)
515 {
516   unsigned int counter;
517   SD_dependency_t dependency;
518
519   XBT_INFO("Displaying task %s", SD_task_get_name(task));
520   char *statename = bprintf("%s%s%s%s%s%s%s",
521                       (task->state == SD_NOT_SCHEDULED ? " not scheduled" : ""),
522                       (task->state == SD_SCHEDULABLE ? " schedulable" : ""),
523                       (task->state == SD_SCHEDULED ? " scheduled" : ""),
524                       (task->state == SD_RUNNABLE ? " runnable" : " not runnable"),
525                       (task->state == SD_RUNNING ? " running" : ""),
526                       (task->state == SD_DONE ? " done" : ""),
527                       (task->state == SD_FAILED ? " failed" : ""));
528   XBT_INFO("  - state:%s", statename);
529   free(statename);
530
531   if (task->kind != 0) {
532     switch (task->kind) {
533     case SD_TASK_COMM_E2E:
534       XBT_INFO("  - kind: end-to-end communication");
535       break;
536     case SD_TASK_COMP_SEQ:
537       XBT_INFO("  - kind: sequential computation");
538       break;
539     case SD_TASK_COMP_PAR_AMDAHL:
540       XBT_INFO("  - kind: parallel computation following Amdahl's law");
541       break;
542     case SD_TASK_COMM_PAR_MXN_1D_BLOCK:
543       XBT_INFO("  - kind: MxN data redistribution assuming 1D block distribution");
544       break;
545     default:
546       XBT_INFO("  - (unknown kind %d)", task->kind);
547     }
548   }
549
550   if (task->category)
551     XBT_INFO("  - tracing category: %s", task->category);
552
553   XBT_INFO("  - amount: %.0f", SD_task_get_amount(task));
554   if (task->kind == SD_TASK_COMP_PAR_AMDAHL)
555     XBT_INFO("  - alpha: %.2f", task->alpha);
556   XBT_INFO("  - Dependencies to satisfy: %d", task->unsatisfied_dependencies);
557   if (xbt_dynar_is_empty(task->tasks_before) == 0) {
558     XBT_INFO("  - pre-dependencies:");
559     xbt_dynar_foreach(task->tasks_before, counter, dependency) {
560       XBT_INFO("    %s", SD_task_get_name(dependency->src));
561     }
562   }
563   if (xbt_dynar_is_empty(task->tasks_after)== 0) {
564     XBT_INFO("  - post-dependencies:");
565     xbt_dynar_foreach(task->tasks_after, counter, dependency) {
566       XBT_INFO("    %s", SD_task_get_name(dependency->dst));
567     }
568   }
569 }
570
571 /** @brief Dumps the task in dotty formalism into the FILE* passed as second argument */
572 void SD_task_dotty(SD_task_t task, void *out)
573 {
574   unsigned int counter;
575   SD_dependency_t dependency;
576   FILE *fout = static_cast<FILE*>(out);
577   fprintf(fout, "  T%p [label=\"%.20s\"", task, task->name);
578   switch (task->kind) {
579   case SD_TASK_COMM_E2E:
580   case SD_TASK_COMM_PAR_MXN_1D_BLOCK:
581     fprintf(fout, ", shape=box");
582     break;
583   case SD_TASK_COMP_SEQ:
584   case SD_TASK_COMP_PAR_AMDAHL:
585     fprintf(fout, ", shape=circle");
586     break;
587   default:
588     xbt_die("Unknown task type!");
589   }
590   fprintf(fout, "];\n");
591   xbt_dynar_foreach(task->tasks_before, counter, dependency) {
592     fprintf(fout, " T%p -> T%p;\n", dependency->src, dependency->dst);
593   }
594 }
595
596 /**
597  * \brief Adds a dependency between two tasks
598  *
599  * \a dst will depend on \a src, ie \a dst will not start before \a src is finished.
600  * Their \ref e_SD_task_state_t "state" must be #SD_NOT_SCHEDULED, #SD_SCHEDULED or #SD_RUNNABLE.
601  *
602  * \param name the name of the new dependency (can be \c nullptr)
603  * \param data the user data you want to associate with this dependency (can be \c nullptr)
604  * \param src the task which must be executed first
605  * \param dst the task you want to make depend on \a src
606  * \see SD_task_dependency_remove()
607  */
608 void SD_task_dependency_add(const char *name, void *data, SD_task_t src, SD_task_t dst)
609 {
610   bool found = false;
611   SD_dependency_t dependency;
612
613   unsigned long length = xbt_dynar_length(src->tasks_after);
614
615   if (src == dst)
616     THROWF(arg_error, 0, "Cannot add a dependency between task '%s' and itself", SD_task_get_name(src));
617
618   e_SD_task_state_t state = SD_task_get_state(src);
619   if (state != SD_NOT_SCHEDULED && state != SD_SCHEDULABLE && state != SD_RUNNING && state != SD_SCHEDULED &&
620        state != SD_RUNNABLE)
621     THROWF(arg_error, 0, "Task '%s' must be SD_NOT_SCHEDULED, SD_SCHEDULABLE, SD_SCHEDULED, SD_RUNNABLE, or SD_RUNNING",
622            SD_task_get_name(src));
623
624   state = SD_task_get_state(dst);
625   if (state != SD_NOT_SCHEDULED && state != SD_SCHEDULABLE && state != SD_SCHEDULED && state != SD_RUNNABLE)
626     THROWF(arg_error, 0, "Task '%s' must be SD_NOT_SCHEDULED, SD_SCHEDULABLE, SD_SCHEDULED, or SD_RUNNABLE",
627            SD_task_get_name(dst));
628
629   XBT_DEBUG("SD_task_dependency_add: src = %s, dst = %s", SD_task_get_name(src), SD_task_get_name(dst));
630   for (unsigned long i = 0; i < length && !found; i++) {
631     xbt_dynar_get_cpy(src->tasks_after, i, &dependency);
632     found = (dependency->dst == dst);
633     XBT_DEBUG("Dependency %lu: dependency->dst = %s", i, SD_task_get_name(dependency->dst));
634   }
635
636   if (found)
637     THROWF(arg_error, 0, "A dependency already exists between task '%s' and task '%s'",
638            SD_task_get_name(src), SD_task_get_name(dst));
639
640   dependency = xbt_new(s_SD_dependency_t, 1);
641
642   dependency->name = xbt_strdup(name);  /* xbt_strdup is cleaver enough to deal with nullptr args itself */
643   dependency->data = data;
644   dependency->src = src;
645   dependency->dst = dst;
646
647   /* src must be executed before dst */
648   xbt_dynar_push(src->tasks_after, &dependency);
649   xbt_dynar_push(dst->tasks_before, &dependency);
650
651   dst->unsatisfied_dependencies++;
652   dst->is_not_ready++;
653
654   /* if the task was runnable, then dst->tasks_before is not empty anymore, so we must go back to state SD_SCHEDULED */
655   if (SD_task_get_state(dst) == SD_RUNNABLE) {
656     XBT_DEBUG("SD_task_dependency_add: %s was runnable and becomes scheduled!", SD_task_get_name(dst));
657     SD_task_set_state(dst, SD_SCHEDULED);
658   }
659 }
660
661 /**
662  * \brief Returns the name given as input when dependency has been created.
663  *
664  * \param src a task
665  * \param dst a task depending on \a src
666  */
667 const char *SD_task_dependency_get_name(SD_task_t src, SD_task_t dst)
668 {
669   unsigned int i;
670   SD_dependency_t dependency;
671
672   xbt_dynar_foreach(src->tasks_after, i, dependency){
673     if (dependency->dst == dst)
674       return dependency->name;
675   }
676   return nullptr;
677 }
678
679 /**
680  * \brief Indicates whether there is a dependency between two tasks.
681  *
682  * \param src a task
683  * \param dst a task depending on \a src
684  *
685  * If src is nullptr, checks whether dst has any pre-dependency.
686  * If dst is nullptr, checks whether src has any post-dependency.
687  */
688 int SD_task_dependency_exists(SD_task_t src, SD_task_t dst)
689 {
690   xbt_assert(src != nullptr || dst != nullptr, "Invalid parameter: both src and dst are nullptr");
691
692   if (src) {
693     if (dst) {
694       unsigned int counter;
695       SD_dependency_t dependency;
696       xbt_dynar_foreach(src->tasks_after, counter, dependency) {
697         if (dependency->dst == dst)
698           return 1;
699       }
700     } else {
701       return xbt_dynar_length(src->tasks_after);
702     }
703   } else {
704     return xbt_dynar_length(dst->tasks_before);
705   }
706   return 0;
707 }
708
709 /**
710  * \brief Remove a dependency between two tasks
711  *
712  * \param src a task
713  * \param dst a task depending on \a src
714  * \see SD_task_dependency_add()
715  */
716 void SD_task_dependency_remove(SD_task_t src, SD_task_t dst)
717 {
718   unsigned long length;
719   bool found = false;
720   SD_dependency_t dependency;
721
722   /* remove the dependency from src->tasks_after */
723   length = xbt_dynar_length(src->tasks_after);
724
725   for (unsigned long i = 0; i < length && !found; i++) {
726     xbt_dynar_get_cpy(src->tasks_after, i, &dependency);
727     if (dependency->dst == dst) {
728       xbt_dynar_remove_at(src->tasks_after, i, nullptr);
729       found = true;
730     }
731   }
732   if (!found)
733     THROWF(arg_error, 0, "No dependency found between task '%s' and '%s': task '%s' is not a successor of task '%s'",
734            SD_task_get_name(src), SD_task_get_name(dst), SD_task_get_name(dst), SD_task_get_name(src));
735
736   /* remove the dependency from dst->tasks_before */
737   length = xbt_dynar_length(dst->tasks_before);
738   found = false;
739
740   for (unsigned long i = 0; i < length && !found; i++) {
741     xbt_dynar_get_cpy(dst->tasks_before, i, &dependency);
742     if (dependency->src == src) {
743       xbt_dynar_remove_at(dst->tasks_before, i, nullptr);
744       __SD_task_dependency_destroy(dependency);
745       dst->unsatisfied_dependencies--;
746       dst->is_not_ready--;
747       found = true;
748     }
749   }
750   /* should never happen... */
751   xbt_assert(found, "SimDag error: task '%s' is a successor of '%s' but task '%s' is not a predecessor of task '%s'",
752               SD_task_get_name(dst), SD_task_get_name(src), SD_task_get_name(src), SD_task_get_name(dst));
753
754   /* if the task was scheduled and dst->tasks_before is empty now, we can make it runnable */
755   if (dst->unsatisfied_dependencies == 0) {
756     if (SD_task_get_state(dst) == SD_SCHEDULED)
757       SD_task_set_state(dst, SD_RUNNABLE);
758     else
759       SD_task_set_state(dst, SD_SCHEDULABLE);
760   }
761
762   if (dst->is_not_ready == 0)
763     SD_task_set_state(dst, SD_SCHEDULABLE);
764 }
765
766 /**
767  * \brief Returns the user data associated with a dependency between two tasks
768  *
769  * \param src a task
770  * \param dst a task depending on \a src
771  * \return the user data associated with this dependency (can be \c nullptr)
772  * \see SD_task_dependency_add()
773  */
774 void *SD_task_dependency_get_data(SD_task_t src, SD_task_t dst)
775 {
776   bool found = false;
777   SD_dependency_t dependency;
778
779   unsigned long length = xbt_dynar_length(src->tasks_after);
780
781   for (unsigned long i = 0; i < length && !found; i++) {
782     xbt_dynar_get_cpy(src->tasks_after, i, &dependency);
783     found = (dependency->dst == dst);
784   }
785   if (!found)
786     THROWF(arg_error, 0, "No dependency found between task '%s' and '%s'",
787            SD_task_get_name(src), SD_task_get_name(dst));
788   return dependency->data;
789 }
790
791 /**
792  * \brief Adds a watch point to a task
793  *
794  * SD_simulate() will stop as soon as the \ref e_SD_task_state_t "state" of this task becomes the one given in argument.
795  * The watch point is then automatically removed.
796  *
797  * \param task a task
798  * \param state the \ref e_SD_task_state_t "state" you want to watch (cannot be #SD_NOT_SCHEDULED)
799  * \see SD_task_unwatch()
800  */
801 void SD_task_watch(SD_task_t task, e_SD_task_state_t state)
802 {
803   if (state & SD_NOT_SCHEDULED)
804     THROWF(arg_error, 0, "Cannot add a watch point for state SD_NOT_SCHEDULED");
805
806   task->watch_points = task->watch_points | state;
807 }
808
809 /**
810  * \brief Removes a watch point from a task
811  *
812  * \param task a task
813  * \param state the \ref e_SD_task_state_t "state" you no longer want to watch
814  * \see SD_task_watch()
815  */
816 void SD_task_unwatch(SD_task_t task, e_SD_task_state_t state)
817 {
818   xbt_assert(state != SD_NOT_SCHEDULED, "SimDag error: Cannot have a watch point for state SD_NOT_SCHEDULED");
819   task->watch_points = task->watch_points & ~state;
820 }
821
822 /**
823  * \brief Returns an approximative estimation of the execution time of a task.
824  *
825  * The estimation is very approximative because the value returned is the time the task would take if it was executed
826  * now and if it was the only task.
827  *
828  * \param task the task to evaluate
829  * \param workstation_nb number of workstations on which the task would be executed
830  * \param workstation_list the workstations on which the task would be executed
831  * \param flops_amount computation amount for each workstation (i.e., an array of workstation_nb doubles)
832  * \param bytes_amount communication amount between each pair of workstations (i.e., a matrix of
833  *        workstation_nb*workstation_nb doubles)
834  * \see SD_schedule()
835  */
836 double SD_task_get_execution_time(SD_task_t task, int workstation_nb, const sg_host_t *workstation_list,
837                                   const double *flops_amount, const double *bytes_amount)
838 {
839   xbt_assert(workstation_nb > 0, "Invalid parameter");
840   double max_time = 0.0;
841
842   /* the task execution time is the maximum execution time of the parallel tasks */
843   for (int i = 0; i < workstation_nb; i++) {
844     double time = 0.0;
845     if (flops_amount != nullptr)
846       time = flops_amount[i] / workstation_list[i]->speed();
847
848     if (bytes_amount != nullptr)
849       for (int j = 0; j < workstation_nb; j++) {
850         if (bytes_amount[i * workstation_nb + j] !=0 ) {
851           time += (SD_route_get_latency(workstation_list[i], workstation_list[j]) +
852                    bytes_amount[i * workstation_nb + j] /
853                    SD_route_get_bandwidth(workstation_list[i], workstation_list[j]));
854         }
855       }
856
857     if (time > max_time) {
858       max_time = time;
859     }
860   }
861   return max_time;
862 }
863
864 static inline void SD_task_do_schedule(SD_task_t task)
865 {
866   if (SD_task_get_state(task) > SD_SCHEDULABLE)
867     THROWF(arg_error, 0, "Task '%s' has already been scheduled", SD_task_get_name(task));
868
869   if (task->unsatisfied_dependencies == 0)
870     SD_task_set_state(task, SD_RUNNABLE);
871   else
872     SD_task_set_state(task, SD_SCHEDULED);
873 }
874
875 /**
876  * \brief Schedules a task
877  *
878  * The task state must be #SD_NOT_SCHEDULED.
879  * Once scheduled, a task is executed as soon as possible in \see SD_simulate, i.e. when its dependencies are satisfied.
880  *
881  * \param task the task you want to schedule
882  * \param host_count number of hosts on which the task will be executed
883  * \param workstation_list the hosts on which the task will be executed
884  * \param flops_amount computation amount for each hosts (i.e., an array of host_count doubles)
885  * \param bytes_amount communication amount between each pair of hosts (i.e., a matrix of host_count*host_count doubles)
886  * \param rate task execution speed rate
887  * \see SD_task_unschedule()
888  */
889 void SD_task_schedule(SD_task_t task, int host_count, const sg_host_t * workstation_list,
890                       const double *flops_amount, const double *bytes_amount, double rate)
891 {
892   xbt_assert(host_count > 0, "workstation_nb must be positive");
893
894   task->host_count = host_count;
895   task->rate = rate;
896
897   if (flops_amount) {
898     task->flops_amount = static_cast<double*>(xbt_realloc(task->flops_amount, sizeof(double) * host_count));
899     memcpy(task->flops_amount, flops_amount, sizeof(double) * host_count);
900   } else {
901     xbt_free(task->flops_amount);
902     task->flops_amount = nullptr;
903   }
904
905   int communication_nb = host_count * host_count;
906   if (bytes_amount) {
907     task->bytes_amount = static_cast<double*>(xbt_realloc(task->bytes_amount, sizeof(double) * communication_nb));
908     memcpy(task->bytes_amount, bytes_amount, sizeof(double) * communication_nb);
909   } else {
910     xbt_free(task->bytes_amount);
911     task->bytes_amount = nullptr;
912   }
913
914   task->host_list =  static_cast<sg_host_t*>(xbt_realloc(task->host_list, sizeof(sg_host_t) * host_count));
915   memcpy(task->host_list, workstation_list, sizeof(sg_host_t) * host_count);
916
917   SD_task_do_schedule(task);
918 }
919
920 /**
921  * \brief Unschedules a task
922  *
923  * The task state must be #SD_SCHEDULED, #SD_RUNNABLE, #SD_RUNNING or #SD_FAILED.
924  * If you call this function, the task state becomes #SD_NOT_SCHEDULED.
925  * Call SD_task_schedule() to schedule it again.
926  *
927  * \param task the task you want to unschedule
928  * \see SD_task_schedule()
929  */
930 void SD_task_unschedule(SD_task_t task)
931 {
932   if (task->state != SD_SCHEDULED && task->state != SD_RUNNABLE && task->state != SD_RUNNING &&
933       task->state != SD_FAILED)
934     THROWF(arg_error, 0, "Task %s: the state must be SD_SCHEDULED, SD_RUNNABLE, SD_RUNNING or SD_FAILED",
935            SD_task_get_name(task));
936
937   if ((task->state == SD_SCHEDULED || task->state == SD_RUNNABLE)
938       /* if the task is scheduled or runnable */
939       && ((task->kind == SD_TASK_COMP_PAR_AMDAHL) || (task->kind == SD_TASK_COMM_PAR_MXN_1D_BLOCK))) {
940           /* Don't free scheduling data for typed tasks */
941     __SD_task_destroy_scheduling_data(task);
942     xbt_free(task->host_list);
943     task->host_list=nullptr;
944     task->host_count = 0;
945   }
946
947   if (SD_task_get_state(task) == SD_RUNNING)
948     /* the task should become SD_FAILED */
949     task->surf_action->cancel();
950   else {
951     if (task->unsatisfied_dependencies == 0)
952       SD_task_set_state(task, SD_SCHEDULABLE);
953     else
954       SD_task_set_state(task, SD_NOT_SCHEDULED);
955   }
956   task->remains = task->amount;
957   task->start_time = -1.0;
958 }
959
960 /* Runs a task. */
961 void SD_task_run(SD_task_t task)
962 {
963   xbt_assert(SD_task_get_state(task) == SD_RUNNABLE, "Task '%s' is not runnable! Task state: %d",
964              SD_task_get_name(task), (int)SD_task_get_state(task));
965   xbt_assert(task->host_list != nullptr, "Task '%s': workstation_list is nullptr!", SD_task_get_name(task));
966
967   XBT_DEBUG("Running task '%s'", SD_task_get_name(task));
968
969   /* Copy the elements of the task into the action */
970   int host_nb = task->host_count;
971   sg_host_t *hosts = xbt_new(sg_host_t, host_nb);
972
973   for (int i = 0; i < host_nb; i++)
974     hosts[i] =  task->host_list[i];
975
976   double *flops_amount = xbt_new0(double, host_nb);
977   double *bytes_amount = xbt_new0(double, host_nb * host_nb);
978
979   if(task->flops_amount)
980     memcpy(flops_amount, task->flops_amount, sizeof(double) * host_nb);
981   if(task->bytes_amount)
982     memcpy(bytes_amount, task->bytes_amount, sizeof(double) * host_nb * host_nb);
983
984   task->surf_action = surf_host_model->executeParallelTask(host_nb, hosts, flops_amount, bytes_amount, task->rate);
985
986   task->surf_action->setData(task);
987
988   XBT_DEBUG("surf_action = %p", task->surf_action);
989
990   if (task->category)
991     TRACE_surf_action(task->surf_action, task->category);
992
993   __SD_task_destroy_scheduling_data(task);      /* now the scheduling data are not useful anymore */
994   SD_task_set_state(task, SD_RUNNING);
995 }
996
997 /**
998  * \brief Returns the start time of a task
999  *
1000  * The task state must be SD_RUNNING, SD_DONE or SD_FAILED.
1001  *
1002  * \param task: a task
1003  * \return the start time of this task
1004  */
1005 double SD_task_get_start_time(SD_task_t task)
1006 {
1007   if (task->surf_action)
1008     return task->surf_action->getStartTime();
1009   else
1010     return task->start_time;
1011 }
1012
1013 /**
1014  * \brief Returns the finish time of a task
1015  *
1016  * The task state must be SD_RUNNING, SD_DONE or SD_FAILED.
1017  * If the state is not completed yet, the returned value is an estimation of the task finish time. This value can
1018  * vary until the task is completed.
1019  *
1020  * \param task: a task
1021  * \return the start time of this task
1022  */
1023 double SD_task_get_finish_time(SD_task_t task)
1024 {
1025   if (task->surf_action)        /* should never happen as actions are destroyed right after their completion */
1026     return task->surf_action->getFinishTime();
1027   else
1028     return task->finish_time;
1029 }
1030
1031 void SD_task_distribute_comp_amdahl(SD_task_t task, int ws_count)
1032 {
1033   xbt_assert(task->kind == SD_TASK_COMP_PAR_AMDAHL, "Task %s is not a SD_TASK_COMP_PAR_AMDAHL typed task."
1034               "Cannot use this function.", SD_task_get_name(task));
1035   task->flops_amount = xbt_new0(double, ws_count);
1036   task->bytes_amount = xbt_new0(double, ws_count * ws_count);
1037   xbt_free(task->host_list);
1038   task->host_count = ws_count;
1039   task->host_list = xbt_new0(sg_host_t, ws_count);
1040   
1041   for(int i=0;i<ws_count;i++){
1042     task->flops_amount[i] = (task->alpha + (1 - task->alpha)/ws_count) * task->amount;
1043   }
1044
1045
1046
1047 /** @brief Auto-schedules a task.
1048  *
1049  * Auto-scheduling mean that the task can be used with SD_task_schedulev(). This allows to specify the task costs at
1050  * creation, and decouple them from the scheduling process where you just specify which resource should deliver the
1051  * mandatory power.
1052  *
1053  * To be auto-schedulable, a task must be type and created with one of the specialized creation functions.
1054  *
1055  * @todo
1056  * We should create tasks kind for the following categories:
1057  *  - Point to point communication (done)
1058  *  - Sequential computation       (done)
1059  *  - group communication (redistribution, several kinds)
1060  *  - parallel tasks with no internal communication (one kind per speedup    model such as Amdahl)
1061  *  - idem+ internal communication. Task type not enough since we cannot store comm cost alongside to comp one)
1062  */
1063 void SD_task_schedulev(SD_task_t task, int count, const sg_host_t * list)
1064 {
1065   int i, j;
1066   SD_dependency_t dep;
1067   unsigned int cpt;
1068   xbt_assert(task->kind != 0, "Task %s is not typed. Cannot automatically schedule it.", SD_task_get_name(task));
1069   switch (task->kind) {
1070   case SD_TASK_COMP_PAR_AMDAHL:
1071     SD_task_distribute_comp_amdahl(task, count);
1072     /* no break */
1073   case SD_TASK_COMM_E2E:
1074   case SD_TASK_COMP_SEQ:
1075     xbt_assert(task->host_count == count, "Got %d locations, but were expecting %d locations", count,task->host_count);
1076     for (i = 0; i < count; i++)
1077       task->host_list[i] = list[i];
1078     if (SD_task_get_kind(task)== SD_TASK_COMP_SEQ && !task->flops_amount){
1079       /*This task has failed and is rescheduled. Reset the flops_amount*/
1080       task->flops_amount = xbt_new0(double, 1);
1081       task->flops_amount[0] = task->remains;
1082     }
1083     SD_task_do_schedule(task);
1084     break;
1085   default:
1086     xbt_die("Kind of task %s not supported by SD_task_schedulev()", SD_task_get_name(task));
1087   }
1088
1089   if (task->kind == SD_TASK_COMM_E2E) {
1090     XBT_VERB("Schedule comm task %s between %s -> %s. It costs %.f bytes", SD_task_get_name(task),
1091           sg_host_get_name(task->host_list[0]), sg_host_get_name(task->host_list[1]), task->bytes_amount[2]);
1092   }
1093
1094   /* Iterate over all children and parents being COMM_E2E to say where I am located (and start them if runnable) */
1095   if (task->kind == SD_TASK_COMP_SEQ) {
1096     XBT_VERB("Schedule computation task %s on %s. It costs %.f flops", SD_task_get_name(task),
1097           sg_host_get_name(task->host_list[0]), task->flops_amount[0]);
1098
1099     xbt_dynar_foreach(task->tasks_before, cpt, dep) {
1100       SD_task_t before = dep->src;
1101       if (before->kind == SD_TASK_COMM_E2E) {
1102         before->host_list[1] = task->host_list[0];
1103
1104         if (before->host_list[0] && (SD_task_get_state(before) < SD_SCHEDULED)) {
1105           SD_task_do_schedule(before);
1106           XBT_VERB ("Auto-Schedule comm task %s between %s -> %s. It costs %.f bytes", SD_task_get_name(before),
1107                sg_host_get_name(before->host_list[0]), sg_host_get_name(before->host_list[1]), before->bytes_amount[2]);
1108         }
1109       }
1110     }
1111     xbt_dynar_foreach(task->tasks_after, cpt, dep) {
1112       SD_task_t after = dep->dst;
1113       if (after->kind == SD_TASK_COMM_E2E) {
1114         after->host_list[0] = task->host_list[0];
1115         if (after->host_list[1] && (SD_task_get_state(after) < SD_SCHEDULED)) {
1116           SD_task_do_schedule(after);
1117           XBT_VERB ("Auto-Schedule comm task %s between %s -> %s. It costs %.f bytes", SD_task_get_name(after),
1118                sg_host_get_name(after->host_list[0]), sg_host_get_name(after->host_list[1]), after->bytes_amount[2]);
1119         }
1120       }
1121     }
1122   }
1123   /* Iterate over all children and parents being MXN_1D_BLOCK to say where I am located (and start them if runnable) */
1124   if (task->kind == SD_TASK_COMP_PAR_AMDAHL) {
1125     XBT_VERB("Schedule computation task %s on %d workstations. %.f flops will be distributed following Amdahl's Law",
1126           SD_task_get_name(task), task->host_count, task->flops_amount[0]);
1127     xbt_dynar_foreach(task->tasks_before, cpt, dep) {
1128       SD_task_t before = dep->src;
1129       if (before->kind == SD_TASK_COMM_PAR_MXN_1D_BLOCK){
1130         if (!before->host_list){
1131           XBT_VERB("Sender side of Task %s is not scheduled yet", SD_task_get_name(before));
1132           before->host_list = xbt_new0(sg_host_t, count);
1133           before->host_count = count;
1134           XBT_VERB("Fill the workstation list with list of Task '%s'", SD_task_get_name(task));
1135           for (i=0;i<count;i++)
1136             before->host_list[i] = task->host_list[i];
1137         } else {
1138           XBT_VERB("Build communication matrix for task '%s'", SD_task_get_name(before));
1139           int src_nb, dst_nb;
1140           double src_start, src_end, dst_start, dst_end;
1141           src_nb = before->host_count;
1142           dst_nb = count;
1143           before->host_list = static_cast<sg_host_t*>(xbt_realloc(before->host_list, (before->host_count+count)*sizeof(sg_host_t)));
1144           for(i=0; i<count; i++)
1145             before->host_list[before->host_count+i] = task->host_list[i];
1146
1147           before->host_count += count;
1148           xbt_free(before->flops_amount);
1149           xbt_free(before->bytes_amount);
1150           before->flops_amount = xbt_new0(double, before->host_count);
1151           before->bytes_amount = xbt_new0(double, before->host_count* before->host_count);
1152
1153           for(i=0;i<src_nb;i++){
1154             src_start = i*before->amount/src_nb;
1155             src_end = src_start + before->amount/src_nb;
1156             for(j=0; j<dst_nb; j++){
1157               dst_start = j*before->amount/dst_nb;
1158               dst_end = dst_start + before->amount/dst_nb;
1159               XBT_VERB("(%s->%s): (%.2f, %.2f)-> (%.2f, %.2f)", sg_host_get_name(before->host_list[i]),
1160                   sg_host_get_name(before->host_list[src_nb+j]), src_start, src_end, dst_start, dst_end);
1161               if ((src_end <= dst_start) || (dst_end <= src_start)) {
1162                 before->bytes_amount[i*(src_nb+dst_nb)+src_nb+j]=0.0;
1163               } else {
1164                 before->bytes_amount[i*(src_nb+dst_nb)+src_nb+j] = MIN(src_end, dst_end) - MAX(src_start, dst_start);
1165               }
1166               XBT_VERB("==> %.2f", before->bytes_amount[i*(src_nb+dst_nb)+src_nb+j]);
1167             }
1168           }
1169
1170           if (SD_task_get_state(before)< SD_SCHEDULED) {
1171             SD_task_do_schedule(before);
1172             XBT_VERB ("Auto-Schedule redistribution task %s. Send %.f bytes from %d hosts to %d hosts.",
1173                   SD_task_get_name(before),before->amount, src_nb, dst_nb);
1174             }
1175         }
1176       }
1177     }
1178     xbt_dynar_foreach(task->tasks_after, cpt, dep) {
1179       SD_task_t after = dep->dst;
1180       if (after->kind == SD_TASK_COMM_PAR_MXN_1D_BLOCK){
1181         if (!after->host_list){
1182           XBT_VERB("Receiver side of Task '%s' is not scheduled yet", SD_task_get_name(after));
1183           after->host_list = xbt_new0(sg_host_t, count);
1184           after->host_count = count;
1185           XBT_VERB("Fill the workstation list with list of Task '%s'", SD_task_get_name(task));
1186           for (i=0;i<count;i++)
1187             after->host_list[i] = task->host_list[i];
1188         } else {
1189           int src_nb, dst_nb;
1190           double src_start, src_end, dst_start, dst_end;
1191           src_nb = count;
1192           dst_nb = after->host_count;
1193           after->host_list = static_cast<sg_host_t*>(xbt_realloc(after->host_list, (after->host_count+count)*sizeof(sg_host_t)));
1194           for(i=after->host_count - 1; i>=0; i--)
1195             after->host_list[count+i] = after->host_list[i];
1196           for(i=0; i<count; i++)
1197             after->host_list[i] = task->host_list[i];
1198
1199           after->host_count += count;
1200
1201           xbt_free(after->flops_amount);
1202           xbt_free(after->bytes_amount);
1203
1204           after->flops_amount = xbt_new0(double, after->host_count);
1205           after->bytes_amount = xbt_new0(double, after->host_count* after->host_count);
1206
1207           for(i=0;i<src_nb;i++){
1208             src_start = i*after->amount/src_nb;
1209             src_end = src_start + after->amount/src_nb;
1210             for(j=0; j<dst_nb; j++){
1211               dst_start = j*after->amount/dst_nb;
1212               dst_end = dst_start + after->amount/dst_nb;
1213               XBT_VERB("(%d->%d): (%.2f, %.2f)-> (%.2f, %.2f)", i, j, src_start, src_end, dst_start, dst_end);
1214               if ((src_end <= dst_start) || (dst_end <= src_start)) {
1215                 after->bytes_amount[i*(src_nb+dst_nb)+src_nb+j]=0.0;
1216               } else {
1217                 after->bytes_amount[i*(src_nb+dst_nb)+src_nb+j] = MIN(src_end, dst_end)- MAX(src_start, dst_start);
1218               }
1219               XBT_VERB("==> %.2f", after->bytes_amount[i*(src_nb+dst_nb)+src_nb+j]);
1220             }
1221           }
1222
1223           if (SD_task_get_state(after)< SD_SCHEDULED) {
1224             SD_task_do_schedule(after);
1225             XBT_VERB ("Auto-Schedule redistribution task %s. Send %.f bytes from %d hosts to %d hosts.",
1226               SD_task_get_name(after),after->amount, src_nb, dst_nb);
1227           }
1228         }
1229       }
1230     }
1231   }
1232 }
1233
1234 /** @brief autoschedule a task on a list of workstations
1235  *
1236  * This function is very similar to SD_task_schedulev(), but takes the list of workstations to schedule onto as
1237  * separate parameters.
1238  * It builds a proper vector of workstations and then call SD_task_schedulev()
1239  */
1240 void SD_task_schedulel(SD_task_t task, int count, ...)
1241 {
1242   va_list ap;
1243   sg_host_t *list = xbt_new(sg_host_t, count);
1244   va_start(ap, count);
1245   for (int i = 0; i < count; i++) {
1246     list[i] = va_arg(ap, sg_host_t);
1247   }
1248   va_end(ap);
1249   SD_task_schedulev(task, count, list);
1250   free(list);
1251 }