Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
complete reorganisation of examples/smpi/NAS
[simgrid.git] / src / simdag / sd_dotloader.cpp
1 /* Copyright (c) 2009-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/simdag/simdag_private.h"
8 #include "simgrid/simdag.h"
9 #include "xbt/file.h"
10 #include <string.h>
11
12 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(sd_dotparse, sd, "Parsing DOT files");
13
14 #ifdef HAVE_GRAPHVIZ
15 #include <graphviz/cgraph.h>
16 #endif
17
18 typedef enum {
19   sequential =0,
20   parallel
21 } seq_par_t;
22
23 xbt_dynar_t SD_dotload_generic(const char * filename, seq_par_t seq_or_par, bool schedule);
24
25 static void dot_task_p_free(void *task) {
26   SD_task_destroy(*(SD_task_t *)task);
27 }
28
29 /** @brief loads a DOT file describing a DAG
30  * 
31  * See http://www.graphviz.org/doc/info/lang.html  for more details.
32  * The size attribute of a node describes:
33  *   - for a compute task: the amount of flops to execute
34  *   - for a communication task : the amount of bytes to transfer
35  * If this attribute is ommited, the default value is zero.
36  */
37 xbt_dynar_t SD_dotload(const char *filename) {
38   return SD_dotload_generic(filename, sequential, false);
39 }
40
41 xbt_dynar_t SD_PTG_dotload(const char * filename) {
42   return SD_dotload_generic(filename, parallel, false);
43 }
44
45 xbt_dynar_t SD_dotload_with_sched(const char *filename) {
46   return SD_dotload_generic(filename, sequential, true);
47 }
48
49 static int edge_compare(const void *a, const void *b)
50 {
51   unsigned va = AGSEQ(*(Agedge_t **)a);
52   unsigned vb = AGSEQ(*(Agedge_t **)b);
53   return va == vb ? 0 : (va < vb ? -1 : 1);
54 }
55
56 xbt_dynar_t SD_dotload_generic(const char * filename, seq_par_t seq_or_par, bool schedule){
57   xbt_assert(filename, "Unable to use a null file descriptor\n");
58   FILE *in_file = fopen(filename, "r");
59   xbt_assert(in_file != NULL, "Failed to open file: %s", filename);
60
61   unsigned int i;
62   SD_task_t root, end, task;
63   xbt_dict_t computers;
64   xbt_dynar_t computer = NULL;
65   xbt_dict_cursor_t dict_cursor;
66   bool schedule_success = true;
67
68   xbt_dict_t jobs = xbt_dict_new_homogeneous(NULL);
69   xbt_dynar_t result = xbt_dynar_new(sizeof(SD_task_t), dot_task_p_free);
70
71   Agraph_t * dag_dot = agread(in_file, NIL(Agdisc_t *));
72
73   if (schedule)
74     computers = xbt_dict_new_homogeneous(NULL);
75
76   /* Create all the nodes */
77   Agnode_t *node = NULL;
78   for (node = agfstnode(dag_dot); node; node = agnxtnode(dag_dot, node)) {
79     char *name = agnameof(node);
80     double amount = atof(agget(node, (char *) "size"));
81
82     if (!(task = (SD_task_t)xbt_dict_get_or_null(jobs, name))) {
83       if (seq_or_par == sequential){
84         XBT_DEBUG("See <job id=%s amount =%.0f>", name, amount);
85         task = SD_task_create_comp_seq(name, NULL , amount);
86       } else {
87         double alpha = atof(agget(node, (char *) "alpha"));
88         XBT_DEBUG("See <job id=%s amount =%.0f alpha = %.3f>", name, amount, alpha);
89         task = SD_task_create_comp_par_amdahl(name, NULL , amount, alpha);
90       }
91
92       xbt_dict_set(jobs, name, task, NULL);
93
94       if (strcmp(name,"root") && strcmp(name,"end"))
95         xbt_dynar_push(result, &task);
96
97       if((seq_or_par == sequential) &&
98         ((schedule && schedule_success) || XBT_LOG_ISENABLED(sd_dotparse, xbt_log_priority_verbose))){
99         /* try to take the information to schedule the task only if all is right*/
100         char *char_performer = agget(node, (char *) "performer");
101         char *char_order = agget(node, (char *) "order");
102         /* Tasks will execute on in a given "order" on a given set of "performer" hosts */
103         int performer = ((!char_performer || !strcmp(char_performer,"")) ? -1:atoi(char_performer));
104         int order = ((!char_order || !strcmp(char_order, ""))? -1:atoi(char_order));
105
106         if((performer != -1 && order != -1) && performer < (int) sg_host_count()){
107           /* required parameters are given and less performers than hosts are required */
108           XBT_DEBUG ("Task '%s' is scheduled on workstation '%d' in position '%d'", task->name, performer, order);
109           if(!(computer = (xbt_dynar_t) xbt_dict_get_or_null(computers, char_performer))){
110             computer = xbt_dynar_new(sizeof(SD_task_t), NULL);
111             xbt_dict_set(computers, char_performer, computer, NULL);
112           }
113
114           if((unsigned int)order < xbt_dynar_length(computer)){
115             SD_task_t *task_test = (SD_task_t *)xbt_dynar_get_ptr(computer,order);
116             if(*task_test && *task_test != task){
117               /* the user gave the same order to several tasks */
118               schedule_success = false;
119               XBT_VERB("Task '%s' wants to start on performer '%s' at the same position '%s' as task '%s'",
120                        (*task_test)->name, char_performer, char_order, task->name);
121               continue;
122             }
123           }
124           /* the parameter seems to be ok */
125           xbt_dynar_set_as(computer, order, SD_task_t, task);
126         } else {
127           /* one of required parameters is not given */
128           schedule_success = false;
129           XBT_VERB("The schedule is ignored, task '%s' can not be scheduled on %d hosts", task->name, performer);
130         }
131       }
132     } else {
133       XBT_WARN("Task '%s' is defined more than once", name);
134     }
135   }
136
137   /*Check if 'root' and 'end' nodes have been explicitly declared.  If not, create them. */
138   if (!(root = (SD_task_t)xbt_dict_get_or_null(jobs, "root")))
139     root = (seq_or_par == sequential?SD_task_create_comp_seq("root", NULL, 0):
140                                      SD_task_create_comp_par_amdahl("root", NULL, 0, 0));
141
142   SD_task_set_state(root, SD_SCHEDULABLE);   /* by design the root task is always SCHEDULABLE */
143   xbt_dynar_insert_at(result, 0, &root);     /* Put it at the beginning of the dynar */
144
145   if (!(end = (SD_task_t)xbt_dict_get_or_null(jobs, "end")))
146     end = (seq_or_par == sequential?SD_task_create_comp_seq("end", NULL, 0):
147                                     SD_task_create_comp_par_amdahl("end", NULL, 0, 0));
148
149   /* Create edges */
150   xbt_dynar_t edges = xbt_dynar_new(sizeof(Agedge_t*), NULL);
151   for (node = agfstnode(dag_dot); node; node = agnxtnode(dag_dot, node)) {
152     Agedge_t * edge;
153     xbt_dynar_reset(edges);
154     for (edge = agfstout(dag_dot, node); edge; edge = agnxtout(dag_dot, edge))
155       xbt_dynar_push_as(edges, Agedge_t *, edge);
156
157     /* Be sure edges are sorted */
158     xbt_dynar_sort(edges, edge_compare);
159
160     xbt_dynar_foreach(edges, i, edge) {
161       SD_task_t src, dst;
162       char *src_name=agnameof(agtail(edge)), *dst_name=agnameof(aghead(edge));
163       double size = atof(agget(edge, (char *) "size"));
164
165       src = (SD_task_t)xbt_dict_get_or_null(jobs, src_name);
166       dst = (SD_task_t)xbt_dict_get_or_null(jobs, dst_name);
167
168       if (size > 0) {
169         char *name = (char*)xbt_malloc((strlen(src_name)+strlen(dst_name)+6)*sizeof(char));
170         sprintf(name, "%s->%s", src_name, dst_name);
171         XBT_DEBUG("See <transfer id=%s amount = %.0f>", name, size);
172         if (!(task = (SD_task_t)xbt_dict_get_or_null(jobs, name))) {
173           if (seq_or_par == sequential)
174             task = SD_task_create_comm_e2e(name, NULL , size);
175           else
176             task = SD_task_create_comm_par_mxn_1d_block(name, NULL , size);
177           SD_task_dependency_add(NULL, NULL, src, task);
178           SD_task_dependency_add(NULL, NULL, task, dst);
179           xbt_dict_set(jobs, name, task, NULL);
180           xbt_dynar_push(result, &task);
181         } else {
182           XBT_WARN("Task '%s' is defined more than once", name);
183         }
184         xbt_free(name);
185       } else {
186         SD_task_dependency_add(NULL, NULL, src, dst);
187       }
188     }
189   }
190   xbt_dynar_free(&edges);
191
192   XBT_DEBUG("All tasks have been created, put %s at the end of the dynar", end->name);
193   xbt_dynar_push(result, &end);
194
195   /* Connect entry tasks to 'root', and exit tasks to 'end'*/
196   xbt_dynar_foreach (result, i, task){
197     if (xbt_dynar_is_empty(task->tasks_before) && task != root) {
198       XBT_DEBUG("Task '%s' has no source. Add dependency from 'root'", task->name);
199       SD_task_dependency_add(NULL, NULL, root, task);
200     }
201
202     if (xbt_dynar_is_empty(task->tasks_after) && task != end) {
203       XBT_DEBUG("Task '%s' has no destination. Add dependency to 'end'", task->name);
204       SD_task_dependency_add(NULL, NULL, task, end);
205     }
206   }
207
208   agclose(dag_dot);
209   xbt_dict_free(&jobs);
210   fclose(in_file);
211
212   if(schedule){
213     char *computer_name;
214     if (schedule_success) {
215       const sg_host_t *workstations = sg_host_list ();
216       xbt_dict_foreach(computers,dict_cursor,computer_name,computer){
217         SD_task_t previous_task = NULL;
218         xbt_dynar_foreach(computer, i, task){
219           /* add dependency between the previous and the task to avoid parallel execution */
220           if(task){
221             if(previous_task && !SD_task_dependency_exists(previous_task, task))
222               SD_task_dependency_add(NULL, NULL, previous_task, task);
223
224             SD_task_schedulel(task, 1, workstations[atoi(computer_name)]);
225             previous_task = task;
226           }
227         }
228         xbt_dynar_free(&computer);
229       }
230       xbt_dict_free(&computers);
231     } else {
232       XBT_WARN("The scheduling is ignored");
233       xbt_dict_foreach(computers,dict_cursor,computer_name,computer)
234         xbt_dynar_free(&computer);
235       xbt_dict_free(&computers);
236       xbt_dynar_free(&result);
237       result = NULL;
238     }
239   }
240
241   if (result && !acyclic_graph_detail(result)) {
242     XBT_ERROR("The DOT described in %s is not a DAG. It contains a cycle.", xbt_basename((char*)filename));
243     xbt_dynar_free(&result);
244     result = NULL;
245   }
246   return result;
247 }