8381899cb7b2ceced7d1396784821abb121d2f03
[libfirm] / ir / be / beilpsched.c
1 /*
2  * Copyright (C) 1995-2007 University of Karlsruhe.  All right reserved.
3  *
4  * This file is part of libFirm.
5  *
6  * This file may be distributed and/or modified under the terms of the
7  * GNU General Public License version 2 as published by the Free Software
8  * Foundation and appearing in the file LICENSE.GPL included in the
9  * packaging of this file.
10  *
11  * Licensees holding valid libFirm Professional Edition licenses may use
12  * this file in accordance with the libFirm Commercial License.
13  * Agreement provided with the Software.
14  *
15  * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
16  * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR
17  * PURPOSE.
18  */
19
20 /**
21  * @file
22  * @brief       ILP based instruction scheduling.
23  * @author      Christian Wuerdig
24  * @date        22.10.2006
25  * @version     $Id$
26  *
27  * An ILP scheduler based on
28  * "ILP-based Instruction Scheduling for IA-64"
29  * by Daniel Kaestner and Sebastian Winkel
30  * extended with register pressure constraints by Christian Wuerdig
31  */
32 #ifdef HAVE_CONFIG_H
33 #include "config.h"
34 #endif
35
36 #include "firm_config.h"
37
38 #ifdef WITH_ILP
39
40 #include <math.h>
41
42 #ifndef _WIN32
43 #include <strings.h>
44 #endif /* _WIN32 */
45
46 #include "irnode_t.h"
47 #include "irgwalk.h"
48 #include "irbitset.h"
49 #include "irphase_t.h"
50 #include "height.h"
51 #include "iredges.h"
52 #include "pdeq.h"
53 #include "debug.h"
54 #include "irtools.h"
55 #include "irdump.h"
56 #include "irprintf.h"
57 #include "plist.h"
58 #include "irprintf.h"
59
60 #include <lpp/lpp.h>
61 #include <lpp/lpp_net.h>
62
63 #include <libcore/lc_opts.h>
64 #include <libcore/lc_opts_enum.h>
65 #include <libcore/lc_timing.h>
66
67 #include "be.h"
68 #include "benode_t.h"
69 #include "besched_t.h"
70 #include "beilpsched.h"
71 #include "beutil.h"
72 #include "bestat.h"
73 #include "beirg_t.h"
74
75 typedef struct _ilpsched_options_t {
76         unsigned regpress;
77         unsigned time_limit;
78         char     log_file[1024];
79 } ilpsched_options_t;
80
81 typedef struct _unit_type_info_t {
82         int                            n_units;
83         const be_execution_unit_type_t *tp;
84 } unit_type_info_t;
85
86 /**
87  * holding the ILP variables of the different types
88  */
89 typedef struct _ilp_var_types_t {
90         int *x;   /* x_{nt}^k variables */
91         int *a;   /* a_{nt}^k variables */
92         int *y;   /* y_{nt}^k variables */
93 } ilp_var_types_t;
94
95 /**
96  * Holds alive variables for a node live-in to a block.
97  */
98 typedef struct _ilp_livein_node_t {
99         ir_node  *irn;
100         unsigned max_alive_steps;
101         int      *a;
102 } ilp_livein_node_t;
103
104 /* attributes for a node */
105 typedef struct _ilpsched_node_attr_t {
106         unsigned asap;                     /**< The ASAP scheduling control step */
107         unsigned alap;                     /**< The ALAP scheduling control step */
108         unsigned latency;                  /**< Latency of this node (needed for sorting) */
109         unsigned sched_point;              /**< the step in which the node is finally scheduled */
110         unsigned visit_idx;                /**< Index of the node having visited this node last */
111         unsigned consumer_idx;             /**< Index of the node having counted this node as consumer last */
112         unsigned n_consumer;               /**< Number of consumers */
113         ir_node  **block_consumer;         /**< List of consumer being in the same block */
114         waitq    *projkeeps;               /**< A List of Projs and Keeps belonging to this node */
115         unsigned block_idx     : 30;       /**< A unique per block index */
116         unsigned alap_changed  : 1;        /**< the current ALAP has changed, revisit preds */
117         unsigned is_dummy_node : 1;        /**< this node is assigned to DUMMY unit */
118         bitset_t *transitive_block_nodes;  /**< Set of transitive block nodes (predecessors
119                                                                                         for ASAP, successors for ALAP */
120         unsigned         n_unit_types;     /**< number of allowed execution unit types */
121         unit_type_info_t *type_info;       /**< list of allowed execution unit types */
122         ilp_var_types_t  ilp_vars;         /**< the different ILP variables */
123 } ilpsched_node_attr_t;
124
125 /* attributes for a block */
126 typedef struct _ilpsched_block_attr_t {
127         unsigned block_last_idx;        /**< The highest node index in block so far */
128         unsigned n_interesting_nodes;   /**< The number of nodes interesting for scheduling */
129         unsigned max_steps;             /**< Upper bound for block execution */
130         plist_t  *root_nodes;           /**< A list of nodes having no user in current block */
131         ir_node  *head_ilp_nodes;       /**< A linked list of nodes which will contribute to ILP */
132         pset     *livein_nodes;         /**< A set of nodes which are live-in to this block */
133 } ilpsched_block_attr_t;
134
135 typedef union _ilpsched_attr_ {
136         ilpsched_node_attr_t  node_attr;
137         ilpsched_block_attr_t block_attr;
138 } ilpsched_attr_t;
139
140 /* A irn for the phase and it's attributes (either node or block) */
141 typedef struct {
142         ir_node         *irn;
143         ilpsched_attr_t attr;
144 } be_ilpsched_irn_t;
145
146 /* The ILP scheduling environment */
147 typedef struct {
148         ir_phase             ph;            /**< The phase */
149         ir_graph             *irg;          /**< The current irg */
150         heights_t            *height;       /**< The heights object of the irg */
151         void                 *irg_env;      /**< An environment for the irg scheduling, provided by the backend */
152         void                 *block_env;    /**< An environment for scheduling a block, provided by the backend */
153         const arch_env_t     *arch_env;
154         const arch_isa_t     *isa;          /**< The ISA */
155         const be_main_env_t  *main_env;
156         const be_machine_t   *cpu;          /**< the current abstract machine */
157         ilpsched_options_t   *opts;         /**< the ilp options for current irg */
158         const be_irg_t       *birg;         /**< The birg object */
159         be_options_t         *be_opts;      /**< backend options */
160         const ilp_sched_selector_t *sel;    /**< The ILP sched selector provided by the backend */
161         DEBUG_ONLY(firm_dbg_module_t *dbg);
162 } be_ilpsched_env_t;
163
164 /* convenience macros to handle phase irn data */
165 #define get_ilpsched_irn(ilpsched_env, irn) (phase_get_or_set_irn_data(&(ilpsched_env)->ph, (irn)))
166 #define is_ilpsched_block(node)             (is_Block((node)->irn))
167 #define get_ilpsched_block_attr(block)      (&(block)->attr.block_attr)
168 #define get_ilpsched_node_attr(node)        (&(node)->attr.node_attr)
169
170 /* check if node is considered for ILP scheduling */
171 #define consider_for_sched(isa, irn) \
172         (! (is_Block(irn)            ||  \
173                 is_normal_Proj(isa, irn) ||  \
174                 is_Phi(irn)              ||  \
175                 is_NoMem(irn)            ||  \
176                 is_Unknown(irn)          ||  \
177                 is_End(irn)                  \
178                 ))
179
180 /* gives the valid scheduling time step interval for a node */
181 #define VALID_SCHED_INTERVAL(na) ((na)->alap - (na)->asap + 1)
182
183 /* gives the valid interval where a node can die */
184 #define VALID_KILL_INTERVAL(ba, na) ((ba)->max_steps - (na)->asap + 1)
185
186 /* gives the corresponding ILP variable for given node, unit and time step */
187 #define ILPVAR_IDX(na, unit, control_step) \
188         ((unit) * VALID_SCHED_INTERVAL((na)) + (control_step) - (na)->asap + 1)
189
190 /* gives the corresponding dead nodes ILP variable for given node, unit and time step */
191 #define ILPVAR_IDX_DEAD(ba, na, unit, control_step) \
192         ((unit) * VALID_KILL_INTERVAL((ba), (na)) + (control_step) - (na)->asap + 1)
193
194 /* check if a double value is within an epsilon environment of 0 */
195 #define LPP_VALUE_IS_0(dbl) (fabs((dbl)) <= 1e-10)
196
197 #define ilp_timer_push(t)         lc_timer_push((t))
198 #define ilp_timer_pop()           lc_timer_pop()
199 #define ilp_timer_elapsed_usec(t) lc_timer_elapsed_usec((t))
200
201 /* option variable */
202 static ilpsched_options_t ilp_opts = {
203         1,     /* default is with register pressure constraints */
204         300,   /* 300 sec per block time limit */
205         ""     /* no log file */
206 };
207
208 /* ILP options */
209 static const lc_opt_table_entry_t ilpsched_option_table[] = {
210         LC_OPT_ENT_BOOL("regpress",  "Use register pressure constraints", &ilp_opts.regpress),
211         LC_OPT_ENT_INT("time_limit", "ILP time limit per block", &ilp_opts.time_limit),
212         LC_OPT_ENT_STR("lpp_log",    "LPP logfile (stderr and stdout are supported)", ilp_opts.log_file, sizeof(ilp_opts.log_file)),
213         LC_OPT_LAST
214 };
215
216 /*
217         We need this global variable as we compare nodes dependent on heights,
218         but we cannot pass any information to the qsort compare function.
219 */
220 static heights_t *glob_heights;
221
222 /**
223  * Check if irn is a Proj, which has no execution units assigned.
224  * @return 1 if irn is a Proj having no execution units assigned, 0 otherwise
225  */
226 static INLINE int is_normal_Proj(const arch_isa_t *isa, const ir_node *irn) {
227         return is_Proj(irn) && (arch_isa_get_allowed_execution_units(isa, irn) == NULL);
228 }
229
230 /**
231  * Skips normal Projs.
232  * @return predecessor if irn is a normal Proj, otherwise irn.
233  */
234 static INLINE ir_node *skip_normal_Proj(const arch_isa_t *isa, ir_node *irn) {
235         if (is_normal_Proj(isa, irn))
236                 return get_Proj_pred(irn);
237         return irn;
238 }
239
240 static INLINE int fixed_latency(const ilp_sched_selector_t *sel, ir_node *irn, void *env) {
241         unsigned lat = be_ilp_sched_latency(sel, irn, env);
242         if (lat == 0 && ! is_Proj(irn) && ! be_is_Keep(irn))
243                 lat = 1;
244         return lat;
245 }
246
247 static int cmp_live_in_nodes(const void *a, const void *b) {
248         const ilp_livein_node_t *n1 = a;
249         const ilp_livein_node_t *n2 = b;
250
251         return n1->irn != n2->irn;
252 }
253
254 /**
255  * Compare scheduling time steps of two be_ilpsched_irn's.
256  */
257 static int cmp_ilpsched_irn(const void *a, const void *b) {
258         be_ilpsched_irn_t    *n1   = *(be_ilpsched_irn_t **)a;
259         be_ilpsched_irn_t    *n2   = *(be_ilpsched_irn_t **)b;
260         ilpsched_node_attr_t *n1_a = get_ilpsched_node_attr(n1);
261         ilpsched_node_attr_t *n2_a = get_ilpsched_node_attr(n2);
262
263         if (n1_a->sched_point == n2_a->sched_point) {
264                 ir_node *irn_a = n1->irn;
265                 ir_node *irn_b = n2->irn;
266
267                 if (heights_reachable_in_block(glob_heights, irn_a, irn_b))
268                         return 1;
269                 if (heights_reachable_in_block(glob_heights, irn_b, irn_a))
270                         return -1;
271
272                 /*
273                         Ok, timestep is equal and the nodes are parallel,
274                         so check latency and schedule high latency first.
275                 */
276                 return QSORT_CMP(n2_a->latency, n1_a->latency);
277         }
278         else
279                 return QSORT_CMP(n1_a->sched_point, n2_a->sched_point);
280 }
281
282 /**
283  * In case there is no phase information for irn, initialize it.
284  */
285 static void *init_ilpsched_irn(ir_phase *ph, ir_node *irn, void *old) {
286         be_ilpsched_irn_t *res = old ? old : phase_alloc(ph, sizeof(res[0]));
287
288         if (res == old) {
289                 /* if we have already some data: check for reinitialization */
290
291                 if (! is_Block(irn)) {
292                         ilpsched_node_attr_t *na = get_ilpsched_node_attr(res);
293
294                         if (! na->transitive_block_nodes) {
295                                 ir_node               *block      = get_nodes_block(irn);
296                                 be_ilpsched_irn_t     *block_node = phase_get_or_set_irn_data(ph, block);
297                                 ilpsched_block_attr_t *ba         = get_ilpsched_block_attr(block_node);
298
299                                 /* we are called after the block indices have been build: create bitset */
300                                 na->transitive_block_nodes = bitset_obstack_alloc(phase_obst(ph), ba->block_last_idx);
301                         }
302                         else {
303                                 /* we are called from reinit block data: clear the bitset */
304                                 bitset_clear_all(na->transitive_block_nodes);
305                                 na->visit_idx    = 0;
306                                 na->alap_changed = 1;
307                         }
308                 }
309                 return old;
310         }
311
312         res->irn = irn;
313
314         /* set ilpsched irn attributes (either block or irn) */
315         if (is_Block(irn)) {
316                 ilpsched_block_attr_t *ba = get_ilpsched_block_attr(res);
317
318                 ba->n_interesting_nodes = 0;
319                 ba->block_last_idx      = 0;
320                 ba->root_nodes          = plist_new();
321                 ba->head_ilp_nodes      = NULL;
322                 ba->livein_nodes        = new_pset(cmp_live_in_nodes, 16);
323                 ba->max_steps           = 0;
324         }
325         else {
326                 ilpsched_node_attr_t *na = get_ilpsched_node_attr(res);
327                 memset(na, 0, sizeof(*na));
328         }
329
330         return res;
331 }
332
333 /**
334  * Assign a per block unique number to each node.
335  */
336 static void build_block_idx(ir_node *irn, void *walk_env) {
337         be_ilpsched_env_t     *env = walk_env;
338         be_ilpsched_irn_t     *node, *block_node;
339         ilpsched_node_attr_t  *na;
340         ilpsched_block_attr_t *ba;
341
342         set_irn_link(irn, NULL);
343         if (! consider_for_sched(env->arch_env->isa, irn))
344                 return;
345
346         node       = get_ilpsched_irn(env, irn);
347         na         = get_ilpsched_node_attr(node);
348         block_node = get_ilpsched_irn(env, get_nodes_block(irn));
349         ba         = get_ilpsched_block_attr(block_node);
350
351         na->block_idx = ba->block_last_idx++;
352 }
353
354 /********************************************************
355  *                              __        _
356  *                             / /       | |
357  *   __ _ ___  __ _ _ __      / /    __ _| | __ _ _ __
358  *  / _` / __|/ _` | '_ \    / /    / _` | |/ _` | '_ \
359  * | (_| \__ \ (_| | |_) |  / /    | (_| | | (_| | |_) |
360  *  \__,_|___/\__,_| .__/  /_/      \__,_|_|\__,_| .__/
361  *                 | |                           | |
362  *                 |_|                           |_|
363  ********************************************************/
364
365 /**
366  * Add all nodes having no user in current block to last_nodes list.
367  */
368 static void collect_alap_root_nodes(ir_node *irn, void *walk_env) {
369         ir_node               *block;
370         const ir_edge_t       *edge;
371         be_ilpsched_irn_t     *block_node, *node;
372         ilpsched_block_attr_t *ba;
373         ilpsched_node_attr_t  *na;
374         int                   i, j;
375         be_ilpsched_env_t     *env           = walk_env;
376         int                   has_block_user = 0;
377         unsigned              n_consumer     = 0;
378         ir_edge_kind_t        ekind[2]       = { EDGE_KIND_NORMAL, EDGE_KIND_DEP };
379         ir_node               **consumer;
380         unsigned              idx;
381
382         if (! consider_for_sched(env->arch_env->isa, irn))
383                 return;
384
385         block    = get_nodes_block(irn);
386     idx      = get_irn_idx(irn);
387         consumer = NEW_ARR_F(ir_node *, 0);
388
389         DBG((env->dbg, LEVEL_3, "%+F (%+F) is interesting, examining ... ", irn, block));
390
391         /* check data and dependency out edges */
392         for (i = 0; i < 2 && ! has_block_user; ++i) {
393                 foreach_out_edge_kind(irn, edge, ekind[i]) {
394                         ir_node *user = get_edge_src_irn(edge);
395
396                         if (is_normal_Proj(env->arch_env->isa, user)) {
397                                 const ir_edge_t *user_edge;
398
399                                 if (get_irn_mode(user) == mode_X)
400                                         continue;
401
402                                 /* The ABI ensures, that there will be no ProjT nodes in the graph. */
403                                 for (j = 0; j < 2; ++j) {
404                                         foreach_out_edge_kind(user, user_edge, ekind[j]) {
405                                                 ir_node *real_user = get_edge_src_irn(user_edge);
406
407                                                 if (! is_Phi(real_user) && ! is_Block(real_user)) {
408                                                         be_ilpsched_irn_t    *node = get_ilpsched_irn(env, real_user);
409                                                         ilpsched_node_attr_t *ua   = get_ilpsched_node_attr(node);
410
411                                                         /* skip already visited nodes */
412                                                         if (ua->consumer_idx == idx)
413                                                                 continue;
414
415                                                         /* check if node has user in this block and collect the user if it's a data user */
416                                                         if (get_nodes_block(real_user) == block) {
417                                                                 if (i == 0 && j == 0)
418                                                                         ARR_APP1(ir_node *, consumer, real_user);
419                                                                 has_block_user = 1;
420                                                         }
421
422                                                         /* only count data consumer */
423                                                         if (i == 0)
424                                                                 n_consumer++;
425
426                                                         /* mark user as visited by this node */
427                                                         ua->consumer_idx = idx;
428                                                 }
429                                         }
430                                 }
431                         }
432                         else if (is_Block(user)) {
433                                 continue;
434                         }
435                         else if (! is_Phi(user)) {
436                                 be_ilpsched_irn_t    *node = get_ilpsched_irn(env, user);
437                                 ilpsched_node_attr_t *ua   = get_ilpsched_node_attr(node);
438
439                                 /* skip already visited nodes */
440                                 if (ua->consumer_idx == idx)
441                                         continue;
442
443                                 /* check if node has user in this block and collect the user if it's a data user */
444                                 if (get_nodes_block(user) == block) {
445                                         if (i == 0)
446                                                 ARR_APP1(ir_node *, consumer, user);
447                                         has_block_user = 1;
448                                 }
449
450                                 /* only count data consumer */
451                                 if (i == 0)
452                                         n_consumer++;
453
454                                 /* mark user visited by this node */
455                                 ua->consumer_idx = idx;
456                         }
457                         else if (get_nodes_block(user) != block) {
458                                 n_consumer++;
459                         }
460                 }
461         }
462
463         block_node = get_ilpsched_irn(env, block);
464         ba         = get_ilpsched_block_attr(block_node);
465
466         ba->n_interesting_nodes++;
467
468         /* current irn has no user inside this block, add to queue */
469         if (! has_block_user) {
470                 DB((env->dbg, LEVEL_3, "root node\n"));
471                 plist_insert_back(ba->root_nodes, irn);
472         }
473         else {
474                 DB((env->dbg, LEVEL_3, "normal node\n"));
475         }
476
477         /* record number of all consumer and the consumer within the same block */
478         node = get_ilpsched_irn(env, irn);
479         na   = get_ilpsched_node_attr(node);
480         na->n_consumer     = n_consumer;
481         na->block_consumer = NEW_ARR_D(ir_node *, phase_obst(&env->ph), ARR_LEN(consumer));
482         memcpy(na->block_consumer, consumer, ARR_LEN(consumer) * sizeof(na->block_consumer[0]));
483         DEL_ARR_F(consumer);
484 }
485
486 /**
487  * Calculate the ASAP scheduling step for current irn.
488  */
489 static void calculate_irn_asap(ir_node *irn, void *walk_env) {
490         be_ilpsched_env_t     *env = walk_env;
491         int                   i;
492         ir_node               *block;
493         be_ilpsched_irn_t     *node, *block_node;
494         ilpsched_node_attr_t  *na;
495         ilpsched_block_attr_t *ba;
496
497         /* These nodes are handled separate */
498         if (! consider_for_sched(env->arch_env->isa, irn))
499                 return;
500
501         DBG((env->dbg, LEVEL_2, "Calculating ASAP of node %+F ... ", irn));
502
503         block    = get_nodes_block(irn);
504         node     = get_ilpsched_irn(env, irn);
505         na       = get_ilpsched_node_attr(node);
506         na->asap = 1;
507
508         for (i = get_irn_ins_or_deps(irn) - 1; i >= 0; --i) {
509                 ir_node *pred = skip_normal_Proj(env->arch_env->isa, get_irn_in_or_dep(irn, i));
510
511                 /* check for greatest distance to top */
512                 if (! is_Phi(pred) && ! is_NoMem(pred) && get_nodes_block(pred) == block) {
513                         be_ilpsched_irn_t    *pred_node = get_ilpsched_irn(env, pred);
514                         ilpsched_node_attr_t *pna       = get_ilpsched_node_attr(pred_node);
515                         unsigned             lat;
516
517                         lat         = fixed_latency(env->sel, pred, env->block_env);
518                         na->latency = lat;
519                         na->asap    = MAX(na->asap, pna->asap + lat);
520                 }
521         }
522
523         /* add node to ILP node list and update max_steps */
524         block_node = get_ilpsched_irn(env, block);
525         ba         = get_ilpsched_block_attr(block_node);
526
527         set_irn_link(irn, ba->head_ilp_nodes);
528         ba->head_ilp_nodes = irn;
529         ba->max_steps     += fixed_latency(env->sel, irn, env->block_env);
530
531         DB((env->dbg, LEVEL_2, "%u\n", na->asap));
532 }
533
534 /**
535  * Calculate the ALAP scheduling step of all irns in current block.
536  * Depends on max_steps being calculated.
537  */
538 static void calculate_block_alap(ir_node *block, void *walk_env) {
539         be_ilpsched_env_t     *env        = walk_env;
540         be_ilpsched_irn_t     *block_node = get_ilpsched_irn(env, block);
541         ilpsched_block_attr_t *ba         = get_ilpsched_block_attr(block_node);
542         waitq                 *cur_queue  = new_waitq();
543         plist_element_t       *el;
544
545         assert(is_Block(block));
546
547         DBG((env->dbg, LEVEL_2, "Calculating ALAP for nodes in %+F (%u nodes, %u max steps)\n",
548                 block, ba->n_interesting_nodes, ba->max_steps));
549
550         /* TODO: Might be faster to use out edges and call phase_reinit_single_irn_data */
551         //phase_reinit_block_irn_data(&env->ph, block);
552
553         /* init start queue */
554         foreach_plist(ba->root_nodes, el) {
555                 waitq_put(cur_queue, plist_element_get_value(el));
556         }
557
558         /* repeat until all nodes are processed */
559         while (! waitq_empty(cur_queue)) {
560                 waitq *next_queue = new_waitq();
561
562                 /* process all nodes in current step */
563                 while (! waitq_empty(cur_queue)) {
564                         ir_node              *cur_irn = waitq_get(cur_queue);
565                         be_ilpsched_irn_t    *node    = get_ilpsched_irn(env, cur_irn);
566                         ilpsched_node_attr_t *na      = get_ilpsched_node_attr(node);
567                         int                  i;
568
569                         /* cur_node has no alap set -> it's a root node, set to max alap */
570                         if (na->alap == 0) {
571                                 na->alap = ba->max_steps;
572                                 DBG((env->dbg, LEVEL_2, "setting ALAP of node %+F to %u, handling preds:\n",
573                                         cur_irn, na->alap));
574                         }
575                         else {
576                                 DBG((env->dbg, LEVEL_2, "ALAP of node %+F is %u, handling preds:\n",
577                                         cur_irn, na->alap));
578                         }
579
580                         /* set the alap's of all predecessors */
581                         for (i = get_irn_ins_or_deps(cur_irn) - 1; i >= 0; --i) {
582                                 ir_node *pred = skip_normal_Proj(env->arch_env->isa, get_irn_in_or_dep(cur_irn, i));
583
584                                 /* check for greatest distance to bottom */
585                                 if (! is_Phi(pred) && ! is_NoMem(pred) && get_nodes_block(pred) == block) {
586                                         be_ilpsched_irn_t    *pred_node = get_ilpsched_irn(env, pred);
587                                         ilpsched_node_attr_t *pna       = get_ilpsched_node_attr(pred_node);
588                                         unsigned             lat;
589
590                                         /* mark the predecessor as visited by current irn */
591                                         if (pna->visit_idx == get_irn_idx(cur_irn) && ! na->alap_changed)
592                                                 continue;
593                                         pna->visit_idx = get_irn_idx(cur_irn);
594
595                                         lat = fixed_latency(env->sel, pred, env->block_env);
596
597                                         /* set ALAP of current pred */
598                                         if (pna->alap == 0) {
599                                                 /* current ALAP is 0: node has not yet been visited */
600                                                 pna->alap_changed = 1;
601                                                 pna->alap         = na->alap - lat;
602                                         }
603                                         else if (pna->alap > na->alap - lat) {
604                                                 /* we found a longer path to root node: change ALAP */
605                                                 pna->alap         = na->alap - lat;
606                                                 pna->alap_changed = 1;
607                                         }
608                                         else {
609                                                 /* current ALAP is best found so far: keep it */
610                                                 pna->alap_changed = 0;
611                                         }
612
613                                         DBG((env->dbg, LEVEL_2, "\tsetting ALAP of node %+F to %u\n", pred, pna->alap));
614
615                                         /* enqueue node for next iteration */
616                                         if (get_irn_ins_or_deps(pred) > 0)
617                                                 waitq_put(next_queue, pred);
618                                 }
619                         }
620                 }
621
622                 /* prepare for next iteration */
623                 del_waitq(cur_queue);
624                 cur_queue = next_queue;
625         }
626 }
627
628 /**
629  * Free list of root nodes and the set of live-in nodes.
630  */
631 static void clear_unwanted_data(ir_node *block, void *walk_env) {
632         be_ilpsched_env_t     *env        = walk_env;
633         be_ilpsched_irn_t     *block_node = get_ilpsched_irn(env, block);
634         ilpsched_block_attr_t *ba         = get_ilpsched_block_attr(block_node);
635
636         plist_free(ba->root_nodes);
637         ba->root_nodes = NULL;
638         del_pset(ba->livein_nodes);
639         ba->livein_nodes = NULL;
640 }
641
642 /**
643  * Refine the {ASAP(n), ALAP(n)} interval for the nodes.
644  * Set the ASAP/ALAP times of Projs and Keeps to their ancestor ones.
645  */
646 static void refine_asap_alap_times(ir_node *irn, void *walk_env) {
647         be_ilpsched_env_t    *env  = walk_env;
648         ir_node              *pred = irn;
649         be_ilpsched_irn_t    *node, *pred_node;
650         ilpsched_node_attr_t *na, *pna;
651
652         if (! consider_for_sched(env->arch_env->isa, irn))
653                 return;
654
655         if (! is_Proj(irn) && ! be_is_Keep(irn))
656                 return;
657
658         /* go to the ancestor */
659         if (be_is_Keep(irn))
660                 pred = get_irn_n(irn, 0);
661         pred = skip_Proj(pred);
662
663         node      = get_ilpsched_irn(env, irn);
664         pred_node = get_ilpsched_irn(env, pred);
665         na        = get_ilpsched_node_attr(node);
666         pna       = get_ilpsched_node_attr(pred_node);
667
668         na->asap = pna->asap;
669         na->alap = pna->alap;
670
671         /* record all Projs and Keeps for this node */
672         if (! pna->projkeeps)
673                 pna->projkeeps = new_waitq();
674         waitq_put(pna->projkeeps, irn);
675
676         DBG((env->dbg, LEVEL_2, "fixing ASAP/ALAP of %+F to %u/%u\n", irn, na->asap, na->alap));
677 }
678
679 /*******************************************
680  *           _              _       _
681  *          | |            | |     | |
682  *  ___  ___| |__   ___  __| |_   _| | ___
683  * / __|/ __| '_ \ / _ \/ _` | | | | |/ _ \
684  * \__ \ (__| | | |  __/ (_| | |_| | |  __/
685  * |___/\___|_| |_|\___|\__,_|\__,_|_|\___|
686  *
687  *******************************************/
688
689 static INLINE void check_for_keeps(waitq *keeps, ir_node *block, ir_node *irn) {
690         const ir_edge_t *edge;
691         (void) block;
692
693         foreach_out_edge(irn, edge) {
694                 ir_node *user = get_edge_src_irn(edge);
695
696                 if (be_is_Keep(user)) {
697                         assert(get_nodes_block(user) == block && "Keep must not be in different block.");
698                         waitq_put(keeps, user);
699                 }
700         }
701 }
702
703 /**
704  * Inserts @p irn before @p before into schedule and notifies backend.
705  */
706 static INLINE void notified_sched_add_before(be_ilpsched_env_t *env,
707         ir_node *before, ir_node *irn, unsigned cycle)
708 {
709         be_ilp_sched_node_scheduled(env->sel, irn, cycle, env->block_env);
710         sched_add_before(before, irn);
711 }
712
713 /**
714  * Adds a node, it's Projs (in case of mode_T nodes) and
715  * it's Keeps to schedule.
716  */
717 static void add_to_sched(be_ilpsched_env_t *env, ir_node *block, ir_node *irn, unsigned cycle) {
718         const ir_edge_t *edge;
719         waitq           *keeps = new_waitq();
720
721         /* mode_M nodes are not scheduled */
722         if (get_irn_mode(irn) == mode_M)
723                 return;
724
725         if (! sched_is_scheduled(irn))
726                 notified_sched_add_before(env, block, irn, cycle);
727
728         /* add Projs */
729         if (get_irn_mode(irn) == mode_T) {
730                 foreach_out_edge(irn, edge) {
731                         ir_node *user = get_edge_src_irn(edge);
732
733                         if ((to_appear_in_schedule(user) || get_irn_mode(user) == mode_b) &&
734                                 get_irn_n_edges(user) > 0)
735                         {
736                                 notified_sched_add_before(env, block, user, cycle);
737                         }
738
739                         check_for_keeps(keeps, block, user);
740                 }
741         }
742         else {
743                 check_for_keeps(keeps, block, irn);
744         }
745
746         /* add Keeps */
747         while (! waitq_empty(keeps)) {
748                 ir_node *keep = waitq_get(keeps);
749                 if (! sched_is_scheduled(keep))
750                         notified_sched_add_before(env, block, keep, cycle);
751         }
752
753         del_waitq(keeps);
754 }
755
756 /**
757  * Schedule all nodes in the given block, according to the ILP solution.
758  */
759 static void apply_solution(be_ilpsched_env_t *env, lpp_t *lpp, ir_node *block) {
760         be_ilpsched_irn_t     *block_node = get_ilpsched_irn(env, block);
761         ilpsched_block_attr_t *ba         = get_ilpsched_block_attr(block_node);
762         sched_info_t          *info       = get_irn_sched_info(block);
763         be_ilpsched_irn_t     **sched_nodes;
764         unsigned              i, l;
765         ir_node               *cfop, *irn;
766         const ir_edge_t       *edge;
767
768         /* init block schedule list */
769         INIT_LIST_HEAD(&info->list);
770         info->scheduled = 1;
771
772         /* collect nodes and their scheduling time step */
773         sched_nodes = NEW_ARR_F(be_ilpsched_irn_t *, 0);
774         if (ba->n_interesting_nodes == 0) {
775                 /* ignore */
776         }
777         else if (ba->n_interesting_nodes == 1) {
778                 be_ilpsched_irn_t *node = get_ilpsched_irn(env, ba->head_ilp_nodes);
779
780                 /* add the single node */
781                 ARR_APP1(be_ilpsched_irn_t *, sched_nodes, node);
782         }
783         else {
784                 /* check all nodes for their positive solution */
785                 foreach_linked_irns(ba->head_ilp_nodes, irn) {
786                         be_ilpsched_irn_t    *node;
787                         ilpsched_node_attr_t *na;
788                         int                  tp_idx, found;
789                         unsigned             cur_var, t;
790
791                         node    = get_ilpsched_irn(env, irn);
792                         na      = get_ilpsched_node_attr(node);
793                         cur_var = 0;
794                         found   = 0;
795
796                         if (! na->is_dummy_node) {
797                                 for (tp_idx = na->n_unit_types - 1; ! found && tp_idx >= 0; --tp_idx) {
798                                         for (t = na->asap - 1; ! found && t <= na->alap - 1; ++t) {
799                                                 double cost = lpp_get_var_sol(lpp, na->ilp_vars.y[cur_var]);
800
801                                                 if (! LPP_VALUE_IS_0(cost)) {
802                                                         DBG((env->dbg, LEVEL_3, "%+F has additional regpressure costs of %f\n", irn, cost));
803                                                         found = 1;
804                                                 }
805                                         }
806                                 }
807                         }
808
809                         found = 0;
810                         /* go over all variables of a node until the non-zero one is found */
811                         for (tp_idx = na->n_unit_types - 1; ! found && tp_idx >= 0; --tp_idx) {
812                                 for (t = na->asap - 1; ! found && t <= na->alap - 1; ++t) {
813                                         double val = lpp_get_var_sol(lpp, na->ilp_vars.x[cur_var++]);
814
815                                         /* check, if variable is set to one (it's not zero then :) */
816                                         if (! LPP_VALUE_IS_0(val)) {
817                                                 na->sched_point = t;
818                                                 ARR_APP1(be_ilpsched_irn_t *, sched_nodes, node);
819                                                 DBG((env->dbg, LEVEL_2, "Schedpoint of %+F is %u at unit type %s\n",
820                                                         irn, t, na->type_info[tp_idx].tp->name));
821                                                 found = 1;
822                                         }
823                                 }
824                         }
825                 }
826
827                 glob_heights = env->height;
828                 /* sort nodes ascending by scheduling time step */
829                 qsort(sched_nodes, ARR_LEN(sched_nodes), sizeof(sched_nodes[0]), cmp_ilpsched_irn);
830         }
831
832         /* make all Phis ready and remember the single cf op */
833         cfop = NULL;
834         foreach_out_edge(block, edge) {
835                 irn = get_edge_src_irn(edge);
836
837                 switch (get_irn_opcode(irn)) {
838                         case iro_Phi:
839                                 add_to_sched(env, block, irn, 0);
840                                 break;
841                         case iro_Start:
842                         case iro_End:
843                         case iro_Proj:
844                         case iro_Bad:
845                         case iro_Unknown:
846                                 break;
847                         default:
848                                 if (is_cfop(irn)) {
849                                         assert(cfop == NULL && "Highlander - there can be only one");
850                                         cfop = irn;
851                                 }
852                         break;
853                 }
854         }
855
856         /* add all nodes from list */
857         for (i = 0, l = ARR_LEN(sched_nodes); i < l; ++i) {
858                 ilpsched_node_attr_t *na = get_ilpsched_node_attr(sched_nodes[i]);
859                 if (sched_nodes[i]->irn != cfop)
860                         add_to_sched(env, block, sched_nodes[i]->irn, na->sched_point);
861         }
862
863         /* schedule control flow node if not already done */
864         if (cfop && ! sched_is_scheduled(cfop))
865                 add_to_sched(env, block, cfop, 0);
866
867         DEL_ARR_F(sched_nodes);
868 }
869
870 /***************************************************************
871  *   _____ _      _____     _____           _   _
872  *  |_   _| |    |  __ \   / ____|         | | (_)
873  *    | | | |    | |__) | | (___   ___  ___| |_ _  ___  _ __
874  *    | | | |    |  ___/   \___ \ / _ \/ __| __| |/ _ \| '_ \
875  *   _| |_| |____| |       ____) |  __/ (__| |_| | (_) | | | |
876  *  |_____|______|_|      |_____/ \___|\___|\__|_|\___/|_| |_|
877  *
878  ***************************************************************/
879
880 /**
881  * Check if node can be executed on given unit type.
882  */
883 static INLINE int is_valid_unit_type_for_node(const be_execution_unit_type_t *tp, be_ilpsched_irn_t *node) {
884         int                  i;
885         ilpsched_node_attr_t *na = get_ilpsched_node_attr(node);
886
887         for (i = na->n_unit_types - 1; i >= 0; --i) {
888                 if (na->type_info[i].tp == tp)
889                         return i;
890         }
891
892         return -1;
893 }
894
895 /************************************************
896  *                   _       _     _
897  *                  (_)     | |   | |
898  *  __   ____ _ _ __ _  __ _| |__ | | ___  ___
899  *  \ \ / / _` | '__| |/ _` | '_ \| |/ _ \/ __|
900  *   \ V / (_| | |  | | (_| | |_) | |  __/\__ \
901  *    \_/ \__,_|_|  |_|\__,_|_.__/|_|\___||___/
902  *
903  ************************************************/
904
905 static int be_ilpsched_set_type_info(be_ilpsched_env_t *env, ir_node *irn, struct obstack *obst) {
906         const be_execution_unit_t ***execunits = arch_isa_get_allowed_execution_units(env->arch_env->isa, irn);
907         unsigned                  n_unit_types = 0;
908         be_ilpsched_irn_t         *node;
909         ilpsched_node_attr_t      *na;
910         unsigned                  unit_idx, tp_idx;
911
912         /* count number of available unit types for this node */
913         for (n_unit_types = 0; execunits[n_unit_types]; ++n_unit_types)
914                 /* just count */ ;
915
916         node = get_ilpsched_irn(env, irn);
917         na   = get_ilpsched_node_attr(node);
918
919         if (! na->type_info) {
920                 na->n_unit_types = n_unit_types;
921                 na->type_info    = NEW_ARR_D(unit_type_info_t, obst, n_unit_types);
922
923                 /* fill the type info array */
924                 for (tp_idx = 0; tp_idx < n_unit_types; ++tp_idx) {
925                         for (unit_idx = 0; execunits[tp_idx][unit_idx]; ++unit_idx) {
926                                 /* beware: we also count number of available units here */
927                                 if (be_machine_is_dummy_unit(execunits[tp_idx][unit_idx]))
928                                         na->is_dummy_node = 1;
929                         }
930
931                         na->type_info[tp_idx].tp      = execunits[tp_idx][0]->tp;
932                         na->type_info[tp_idx].n_units = unit_idx;
933                 }
934         }
935
936         return n_unit_types;
937 }
938
939 /**
940  * Returns the largest alap time of a user of @p irn.
941  * The user must be in block @p block.
942  */
943 static unsigned be_ilpsched_get_max_alap_user(be_ilpsched_env_t *env, ir_node *irn, ir_node *block) {
944         const ir_edge_t *edge;
945         unsigned        max_alap = 0;
946
947         foreach_out_edge(irn, edge) {
948                 ir_node *user = get_edge_src_irn(edge);
949
950                 if (get_nodes_block(user) == block) {
951                         be_ilpsched_irn_t    *node = get_ilpsched_irn(env, user);
952                         ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
953
954                         max_alap = MAX(max_alap, na->alap);
955                 }
956         }
957
958         assert(max_alap > 0);
959         return max_alap;
960 }
961
962 /**
963  * Create the following variables:
964  * - x_{nt}^k    binary     weigthed with: t
965  *      node n is scheduled at time step t to unit type k
966  * ==>> These variables represent the schedule
967  *
968  * - a_{nt}^k    binary     weighted with num_nodes
969  *      node n is alive at time step t on unit type k
970  *
971  * - y_{nt}^k    continuous  weighted with: num_nodes^2
972  *      register pressure over limit for unit type k
973  * ==>> These variables represent the register pressure
974  *
975  */
976 static void create_variables(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node, struct obstack *var_obst) {
977         char                  buf[1024];
978         ir_node               *irn;
979         unsigned              num_block_var, num_nodes;
980         ilp_livein_node_t     *livein;
981         ilpsched_block_attr_t *ba      = get_ilpsched_block_attr(block_node);
982         unsigned              weigth_y = ba->n_interesting_nodes * ba->n_interesting_nodes;
983         lc_timer_t            *t_var   = lc_timer_register("beilpsched_var", "create ilp variables");
984
985         ilp_timer_push(t_var);
986         num_block_var = num_nodes = 0;
987         foreach_linked_irns(ba->head_ilp_nodes, irn) {
988                 be_ilpsched_irn_t    *node;
989                 ilpsched_node_attr_t *na;
990                 unsigned             n_unit_types, tp_idx, n_var, cur_unit;
991                 unsigned             cur_var_ad, cur_var_x, cur_var_y, num_ad;
992                 int                  i;
993
994                 node         = get_ilpsched_irn(env, irn);
995                 na           = get_ilpsched_node_attr(node);
996                 n_unit_types = be_ilpsched_set_type_info(env, irn, var_obst);
997
998                 /* allocate space for ilp variables */
999                 na->ilp_vars.x = NEW_ARR_D(int, var_obst, n_unit_types * VALID_SCHED_INTERVAL(na));
1000                 memset(na->ilp_vars.x, -1, ARR_LEN(na->ilp_vars.x) * sizeof(na->ilp_vars.x[0]));
1001
1002                 /* we need these variables only for "real" nodes */
1003                 if (! na->is_dummy_node) {
1004                         na->ilp_vars.y = NEW_ARR_D(int, var_obst, n_unit_types * VALID_SCHED_INTERVAL(na));
1005                         memset(na->ilp_vars.y, -1, ARR_LEN(na->ilp_vars.y) * sizeof(na->ilp_vars.y[0]));
1006
1007                         num_ad         = ba->max_steps - na->asap + 1;
1008                         na->ilp_vars.a = NEW_ARR_D(int, var_obst, n_unit_types * num_ad);
1009                         memset(na->ilp_vars.a, -1, ARR_LEN(na->ilp_vars.a) * sizeof(na->ilp_vars.a[0]));
1010                 }
1011
1012                 DBG((env->dbg, LEVEL_3, "\thandling %+F (asap %u, alap %u, unit types %u):\n",
1013                         irn, na->asap, na->alap, na->n_unit_types));
1014
1015                 cur_var_x = cur_var_ad = cur_var_y = cur_unit = n_var = 0;
1016                 /* create variables */
1017                 for (tp_idx = 0; tp_idx < n_unit_types; ++tp_idx) {
1018                         unsigned t;
1019
1020                         for (t = na->asap - 1; t <= na->alap - 1; ++t) {
1021                                 /* x_{nt}^k variables */
1022                                 snprintf(buf, sizeof(buf), "x_n%u_%s_%u",
1023                                         get_irn_idx(irn), na->type_info[tp_idx].tp->name, t);
1024                                 na->ilp_vars.x[cur_var_x++] = lpp_add_var(lpp, buf, lpp_binary, (double)(t + 1));
1025                                 DBG((env->dbg, LEVEL_4, "\t\tcreated ILP variable %s\n", buf));
1026                                 /* variable counter */
1027                                 n_var++;
1028                                 num_block_var++;
1029
1030                                 if (! na->is_dummy_node) {
1031                                         /* y_{nt}^k variables */
1032                                         snprintf(buf, sizeof(buf), "y_n%u_%s_%u",
1033                                                 get_irn_idx(irn), na->type_info[tp_idx].tp->name, t);
1034                                         na->ilp_vars.y[cur_var_y++] = lpp_add_var(lpp, buf, lpp_continous, (double)(weigth_y));
1035                                         DBG((env->dbg, LEVEL_4, "\t\tcreated ILP variable %s\n", buf));
1036
1037                                         /* variable counter */
1038                                         n_var++;
1039                                         num_block_var++;
1040                                 }
1041                         }
1042
1043                         /* a node can die at any step t: asap(n) <= t <= U */
1044                         if (! na->is_dummy_node) {
1045                                 for (t = na->asap - 1; t <= ba->max_steps; ++t) {
1046
1047                                         /* a_{nt}^k variables */
1048                                         snprintf(buf, sizeof(buf), "a_n%u_%s_%u",
1049                                                 get_irn_idx(irn), na->type_info[tp_idx].tp->name, t);
1050                                         na->ilp_vars.a[cur_var_ad++] = lpp_add_var(lpp, buf, lpp_binary, (double)(ba->n_interesting_nodes));
1051                                         DBG((env->dbg, LEVEL_4, "\t\tcreated ILP variable %s\n", buf));
1052
1053                                         /* variable counter */
1054                                         n_var++;
1055                                         num_block_var++;
1056                                 }
1057                         }
1058
1059                         /* collect live-in nodes */
1060                         for (i = get_irn_arity(irn) - 1; i >= 0; --i) {
1061                                 ir_node *pred = get_irn_n(irn, i);
1062
1063                                 if (get_nodes_block(pred) != block_node->irn && consider_for_sched(env->arch_env->isa, pred)) {
1064                                         be_ilpsched_set_type_info(env, pred, var_obst);
1065                                         if (! na->is_dummy_node) {
1066                                                 ilp_livein_node_t *entry = obstack_alloc(var_obst, sizeof(*entry));
1067                                                 entry->irn = pred;
1068                                                 entry->a   = NULL;
1069                                                 pset_insert(ba->livein_nodes, entry, (unsigned)get_irn_idx(pred));
1070                                         }
1071                                 }
1072                         }
1073                 }
1074
1075                 DB((env->dbg, LEVEL_3, "%u variables created\n", n_var));
1076                 num_nodes++;
1077         }
1078
1079         /* create alive variables a_{nt}^k for live-ins */
1080         foreach_pset(ba->livein_nodes, livein) {
1081                 be_ilpsched_irn_t    *node;
1082                 ilpsched_node_attr_t *na;
1083                 unsigned             tp_idx, var_idx;
1084                 ir_node              *irn;
1085
1086                 irn  = livein->irn;
1087                 node = get_ilpsched_irn(env, irn);
1088                 na   = get_ilpsched_node_attr(node);
1089
1090                 livein->max_alive_steps = be_ilpsched_get_max_alap_user(env, irn, block_node->irn);
1091
1092                 livein->a = NEW_ARR_D(int, var_obst, na->n_unit_types * livein->max_alive_steps);
1093                 var_idx   = 0;
1094
1095                 /* create variables */
1096                 for (tp_idx = 0; tp_idx < na->n_unit_types; ++tp_idx) {
1097                         unsigned t;
1098
1099                         for (t = 0; t < livein->max_alive_steps; ++t) {
1100                                 /* a_{nt}^k variables */
1101                                 snprintf(buf, sizeof(buf), "al_n%u_%s_%u",
1102                                         get_irn_idx(irn), na->type_info[tp_idx].tp->name, t);
1103                                 livein->a[var_idx++] = lpp_add_var(lpp, buf, lpp_binary, (double)(ba->n_interesting_nodes));
1104                                 DBG((env->dbg, LEVEL_4, "\t\tcreated ILP variable %s\n", buf));
1105                                 num_block_var++;
1106                         }
1107                 }
1108         }
1109
1110         ilp_timer_pop();
1111         DBG((env->dbg, LEVEL_1, "... %u variables for %u nodes created (%g sec)\n",
1112                 num_block_var, num_nodes, ilp_timer_elapsed_usec(t_var) / 1000000.0));
1113 }
1114
1115 /*******************************************************
1116  *                      _             _       _
1117  *                     | |           (_)     | |
1118  *   ___ ___  _ __  ___| |_ _ __ __ _ _ _ __ | |_ ___
1119  *  / __/ _ \| '_ \/ __| __| '__/ _` | | '_ \| __/ __|
1120  * | (_| (_) | | | \__ \ |_| | | (_| | | | | | |_\__ \
1121  *  \___\___/|_| |_|___/\__|_|  \__,_|_|_| |_|\__|___/
1122  *
1123  *******************************************************/
1124
1125 /**
1126  * Collect all operands and nodes @p irn depends on.
1127  * If there is a Proj within the dependencies, all other Projs of the parent node are added as well.
1128  */
1129 static void sta_collect_in_deps(ir_node *irn, ir_nodeset_t *deps) {
1130         int i;
1131
1132         for (i = get_irn_ins_or_deps(irn) - 1; i >= 0; --i) {
1133                 ir_node *p = get_irn_in_or_dep(irn, i);
1134
1135                 if (is_Proj(p)) {
1136                         const ir_edge_t *edge;
1137
1138                         p = get_Proj_pred(p);
1139                         foreach_out_edge(p, edge) {
1140                                 ir_node *src = get_edge_src_irn(edge);
1141                                 ir_nodeset_insert(deps, src);
1142                         }
1143                 }
1144                 else {
1145                         ir_nodeset_insert(deps, p);
1146                 }
1147         }
1148 }
1149
1150 /**
1151  * Create following ILP constraints:
1152  * - the assignment constraints:
1153  *     assure each node is executed once by exactly one (allowed) execution unit
1154  * - the dead node assignment constraints:
1155  *     assure a node can only die at most once
1156  * - the precedence constraints:
1157  *     assure that no data dependencies are violated
1158  */
1159 static void create_assignment_and_precedence_constraints(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1160         unsigned              num_cst_assign, num_cst_prec, num_cst_dead;
1161         char                  buf[1024];
1162         ir_node               *irn;
1163         ilpsched_block_attr_t *ba            = get_ilpsched_block_attr(block_node);
1164         bitset_t              *bs_block_irns = bitset_alloca(ba->block_last_idx);
1165         lc_timer_t            *t_cst_assign  = lc_timer_register("beilpsched_cst_assign", "create assignment constraints");
1166         lc_timer_t            *t_cst_prec    = lc_timer_register("beilpsched_cst_prec",   "create precedence constraints");
1167
1168         num_cst_assign = num_cst_prec = num_cst_dead = 0;
1169         foreach_linked_irns(ba->head_ilp_nodes, irn) {
1170                 int                  cst, tp_idx;
1171                 unsigned             cur_var;
1172                 be_ilpsched_irn_t    *node;
1173                 ilpsched_node_attr_t *na;
1174                 ir_node              *pred;
1175                 ir_nodeset_t          deps;
1176                 ir_nodeset_iterator_t iter;
1177
1178                 ir_nodeset_init(&deps);
1179
1180                 node    = get_ilpsched_irn(env, irn);
1181                 na      = get_ilpsched_node_attr(node);
1182                 cur_var = 0;
1183
1184                 /* the assignment constraint */
1185                 ilp_timer_push(t_cst_assign);
1186                 snprintf(buf, sizeof(buf), "assignment_cst_n%u", get_irn_idx(irn));
1187                 cst = lpp_add_cst_uniq(lpp, buf, lpp_equal, 1.0);
1188                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1189                 num_cst_assign++;
1190
1191                 lpp_set_factor_fast_bulk(lpp, cst, na->ilp_vars.x, ARR_LEN(na->ilp_vars.x), 1.0);
1192                 ilp_timer_pop();
1193
1194                 /* We have separate constraints for Projs and Keeps */
1195                 // ILP becomes infeasible ?!?
1196 //              if (is_Proj(irn) || be_is_Keep(irn))
1197 //                      continue;
1198
1199                 /* the precedence constraints */
1200                 ilp_timer_push(t_cst_prec);
1201                 bs_block_irns = bitset_clear_all(bs_block_irns);
1202
1203                 sta_collect_in_deps(irn, &deps);
1204                 foreach_ir_nodeset(&deps, pred, iter) {
1205                         unsigned             t_low, t_high, t;
1206                         be_ilpsched_irn_t    *pred_node;
1207                         ilpsched_node_attr_t *pna;
1208                         unsigned             delay;
1209
1210                         pred = skip_normal_Proj(env->arch_env->isa, pred);
1211                         if (is_Phi(pred) || block_node->irn != get_nodes_block(pred) || is_NoMem(pred))
1212                                 continue;
1213
1214                         pred_node = get_ilpsched_irn(env, pred);
1215                         pna       = get_ilpsched_node_attr(pred_node);
1216
1217                         assert(pna->asap > 0 && pna->alap >= pna->asap && "Invalid scheduling interval.");
1218
1219                         if (! bitset_is_set(bs_block_irns, pna->block_idx))
1220                                 bitset_set(bs_block_irns, pna->block_idx);
1221                         else
1222                                 continue;
1223
1224                         /* irn = n, pred = m */
1225                         delay  = fixed_latency(env->sel, pred, env->block_env);
1226                         t_low  = MAX(na->asap, pna->asap + delay - 1);
1227                         t_high = MIN(na->alap, pna->alap + delay - 1);
1228                         for (t = t_low - 1; t <= t_high - 1; ++t) {
1229                                 unsigned tn, tm;
1230                                 int      *tmp_var_idx = NEW_ARR_F(int, 0);
1231
1232                                 snprintf(buf, sizeof(buf), "precedence_n%u_n%u_%u", get_irn_idx(pred), get_irn_idx(irn), t);
1233                                 cst = lpp_add_cst_uniq(lpp, buf, lpp_less, 1.0);
1234                                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1235                                 num_cst_prec++;
1236
1237                                 /* lpp_set_factor_fast_bulk needs variables sorted ascending by index */
1238                                 if (na->ilp_vars.x[0] < pna->ilp_vars.x[0]) {
1239                                         /* node variables have smaller index than pred variables */
1240                                         for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1241                                                 for (tn = na->asap - 1; tn <= t; ++tn) {
1242                                                         unsigned idx = ILPVAR_IDX(na, tp_idx, tn);
1243                                                         ARR_APP1(int, tmp_var_idx, na->ilp_vars.x[idx]);
1244                                                 }
1245                                         }
1246
1247                                         for (tp_idx = pna->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1248                                                 for (tm = t - delay + 1; tm < pna->alap; ++tm) {
1249                                                         unsigned idx = ILPVAR_IDX(pna, tp_idx, tm);
1250                                                         ARR_APP1(int, tmp_var_idx, pna->ilp_vars.x[idx]);
1251                                                 }
1252                                         }
1253                                 }
1254                                 else {
1255                                         /* pred variables have smaller index than node variables */
1256                                         for (tp_idx = pna->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1257                                                 for (tm = t - delay + 1; tm < pna->alap; ++tm) {
1258                                                         unsigned idx = ILPVAR_IDX(pna, tp_idx, tm);
1259                                                         ARR_APP1(int, tmp_var_idx, pna->ilp_vars.x[idx]);
1260                                                 }
1261                                         }
1262
1263                                         for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1264                                                 for (tn = na->asap - 1; tn <= t; ++tn) {
1265                                                         unsigned idx = ILPVAR_IDX(na, tp_idx, tn);
1266                                                         ARR_APP1(int, tmp_var_idx, na->ilp_vars.x[idx]);
1267                                                 }
1268                                         }
1269                                 }
1270
1271                                 if (ARR_LEN(tmp_var_idx) > 0)
1272                                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx, ARR_LEN(tmp_var_idx), 1.0);
1273
1274                                 DEL_ARR_F(tmp_var_idx);
1275                         }
1276                 }
1277                 ir_nodeset_destroy(&deps);
1278                 ilp_timer_pop();
1279         }
1280         DBG((env->dbg, LEVEL_1, "\t%u assignement constraints (%g sec)\n",
1281                 num_cst_assign, ilp_timer_elapsed_usec(t_cst_assign) / 1000000.0));
1282         DBG((env->dbg, LEVEL_1, "\t%u precedence constraints (%g sec)\n",
1283                 num_cst_prec, ilp_timer_elapsed_usec(t_cst_prec) / 1000000.0));
1284 }
1285
1286 /**
1287  * Create ILP resource constraints:
1288  * - assure that for each time step not more instructions are scheduled
1289  *   to the same unit types as units of this type are available
1290  */
1291 static void create_ressource_constraints(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1292         int                   glob_type_idx;
1293         char                  buf[1024];
1294         unsigned              num_cst_resrc = 0;
1295         ilpsched_block_attr_t *ba           = get_ilpsched_block_attr(block_node);
1296         lc_timer_t            *t_cst_rsrc   = lc_timer_register("beilpsched_cst_rsrc",   "create resource constraints");
1297
1298         ilp_timer_push(t_cst_rsrc);
1299         for (glob_type_idx = env->cpu->n_unit_types - 1; glob_type_idx >= 0; --glob_type_idx) {
1300                 unsigned                 t;
1301                 be_execution_unit_type_t *cur_tp = &env->cpu->unit_types[glob_type_idx];
1302
1303                 /* BEWARE: the DUMMY unit type is not in CPU, so it's skipped automatically */
1304
1305                 /* check each time step */
1306                 for (t = 0; t < ba->max_steps; ++t) {
1307                         ir_node *irn;
1308                         int     cst;
1309                         int     *tmp_var_idx = NEW_ARR_F(int, 0);
1310
1311                         snprintf(buf, sizeof(buf), "resource_cst_%s_%u", cur_tp->name, t);
1312                         cst = lpp_add_cst_uniq(lpp, buf, lpp_less, (double)cur_tp->n_units);
1313                         DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1314                         num_cst_resrc++;
1315
1316                         foreach_linked_irns(ba->head_ilp_nodes, irn) {
1317                                 be_ilpsched_irn_t    *node = get_ilpsched_irn(env, irn);
1318                                 ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
1319                                 int                  tp_idx;
1320
1321                                 tp_idx = is_valid_unit_type_for_node(cur_tp, node);
1322
1323                                 if (tp_idx >= 0 && t >= na->asap - 1 && t <= na->alap - 1) {
1324                                         int cur_var = ILPVAR_IDX(na, tp_idx, t);
1325                                         ARR_APP1(int, tmp_var_idx, na->ilp_vars.x[cur_var]);
1326                                 }
1327                         }
1328
1329                         /* set constraints if we have some */
1330                         if (ARR_LEN(tmp_var_idx) > 0)
1331                                 lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx, ARR_LEN(tmp_var_idx), 1.0);
1332
1333                         DEL_ARR_F(tmp_var_idx);
1334                 }
1335         }
1336         ilp_timer_pop();
1337         DBG((env->dbg, LEVEL_1, "\t%u resource constraints (%g sec)\n",
1338                 num_cst_resrc, ilp_timer_elapsed_usec(t_cst_rsrc) / 1000000.0));
1339 }
1340
1341 /**
1342  * Create ILP bundle constraints:
1343  * - assure, at most bundle_size * bundles_per_cycle instructions
1344  *   can be started at a certain point.
1345  */
1346 static void create_bundle_constraints(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1347         char                  buf[1024];
1348         unsigned              t;
1349         unsigned              num_cst_bundle = 0;
1350         unsigned              n_instr_max    = env->cpu->bundle_size * env->cpu->bundels_per_cycle;
1351         ilpsched_block_attr_t *ba            = get_ilpsched_block_attr(block_node);
1352         lc_timer_t            *t_cst_bundle  = lc_timer_register("beilpsched_cst_bundle", "create bundle constraints");
1353
1354         ilp_timer_push(t_cst_bundle);
1355         for (t = 0; t < ba->max_steps; ++t) {
1356                 ir_node *irn;
1357                 int     cst;
1358                 int     *tmp_var_idx = NEW_ARR_F(int, 0);
1359
1360                 snprintf(buf, sizeof(buf), "bundle_cst_%u", t);
1361                 cst = lpp_add_cst_uniq(lpp, buf, lpp_less, (double)n_instr_max);
1362                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1363                 num_cst_bundle++;
1364
1365                 foreach_linked_irns(ba->head_ilp_nodes, irn) {
1366                         be_ilpsched_irn_t    *node;
1367                         ilpsched_node_attr_t *na;
1368                         int                  tp_idx;
1369
1370                         /* Projs and Keeps do not contribute to bundle size */
1371                         if (is_Proj(irn) || be_is_Keep(irn))
1372                                 continue;
1373
1374                         node = get_ilpsched_irn(env, irn);
1375                         na   = get_ilpsched_node_attr(node);
1376
1377                         /* nodes assigned to DUMMY unit do not contribute to bundle size */
1378                         if (na->is_dummy_node)
1379                                 continue;
1380
1381                         if (t >= na->asap - 1 && t <= na->alap - 1) {
1382                                 for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1383                                         int idx = ILPVAR_IDX(na, tp_idx, t);
1384                                         ARR_APP1(int, tmp_var_idx, na->ilp_vars.x[idx]);
1385                                 }
1386                         }
1387                 }
1388
1389                 if (ARR_LEN(tmp_var_idx) > 0)
1390                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx, ARR_LEN(tmp_var_idx), 1.0);
1391
1392                 DEL_ARR_F(tmp_var_idx);
1393         }
1394         ilp_timer_pop();
1395         DBG((env->dbg, LEVEL_1, "\t%u bundle constraints (%g sec)\n",
1396                 num_cst_bundle, ilp_timer_elapsed_usec(t_cst_bundle) / 1000000.0));
1397 }
1398
1399 /**
1400  * Create ILP alive nodes constraints:
1401  * - set variable a_{nt}^k to 1 if nodes n is alive at step t on unit k
1402  */
1403 static void create_alive_nodes_constraint(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1404         char                  buf[1024];
1405         ir_node               *irn;
1406         unsigned              num_cst = 0;
1407         ilpsched_block_attr_t *ba     = get_ilpsched_block_attr(block_node);
1408         lc_timer_t            *t_cst  = lc_timer_register("beilpsched_cst_alive_nodes", "create alive nodes constraints");
1409
1410         ilp_timer_push(t_cst);
1411         /* for each node */
1412         foreach_linked_irns(ba->head_ilp_nodes, irn) {
1413                 be_ilpsched_irn_t    *node = get_ilpsched_irn(env, irn);
1414                 ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
1415                 unsigned             t;
1416
1417                 /* we ignore nodes assigned to dummy unit here */
1418                 if (na->is_dummy_node)
1419                         continue;
1420
1421                 /* check check all time steps: asap(n) <= t <= U */
1422                 for (t = na->asap - 1; t < ba->max_steps; ++t) {
1423                         int node_tp_idx;
1424
1425                         /* for all unit types available for this node */
1426                         for (node_tp_idx = na->n_unit_types - 1; node_tp_idx >= 0; --node_tp_idx) {
1427                                 unsigned tn, tn_max, idx;
1428                                 int      cst, i;
1429                                 int      *tmp_var_idx_n = NEW_ARR_F(int, 0);
1430                                 int      *tmp_var_idx_m = NEW_ARR_F(int, 0);
1431
1432                                 snprintf(buf, sizeof(buf), "alive_node_cst_%u_n%u_%s",
1433                                         t, get_irn_idx(irn), na->type_info[node_tp_idx].tp->name);
1434                                 cst = lpp_add_cst_uniq(lpp, buf, lpp_less, 0.0);
1435                                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1436                                 num_cst++;
1437
1438                                 tn_max = MIN(na->alap - 1, t);
1439                                 /* check if the node has been scheduled so far */
1440                                 for (tn = na->asap - 1; tn <= tn_max; ++tn) {
1441                                         int idx = ILPVAR_IDX(na, node_tp_idx, tn);
1442                                         ARR_APP1(int, tmp_var_idx_n, na->ilp_vars.x[idx]);
1443                                 }
1444
1445                                 if (ARR_LEN(tmp_var_idx_n) > 0)
1446                                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx_n, ARR_LEN(tmp_var_idx_n), (double)(na->n_consumer));
1447                                 DEL_ARR_F(tmp_var_idx_n);
1448
1449                                 /* subtract the number of consumer scheduled so far */
1450                                 for (i = ARR_LEN(na->block_consumer) - 1; i >= 0; --i) {
1451                                         be_ilpsched_irn_t    *cons = get_ilpsched_irn(env, na->block_consumer[i]);
1452                                         ilpsched_node_attr_t *ca   = get_ilpsched_node_attr(cons);
1453                                         int                  tp_idx;
1454                                         unsigned             tm, tm_max;
1455
1456                                         tm_max = MIN(ca->alap - 1, t);
1457                                         for (tp_idx = ca->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1458                                                 for (tm = ca->asap - 1; tm <= tm_max; ++tm) {
1459                                                         int idx = ILPVAR_IDX(ca, tp_idx, tm);
1460                                                         ARR_APP1(int, tmp_var_idx_m, ca->ilp_vars.x[idx]);
1461                                                 }
1462                                         }
1463                                 }
1464
1465                                 if (ARR_LEN(tmp_var_idx_m) > 0)
1466                                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx_m, ARR_LEN(tmp_var_idx_m), -1.0);
1467                                 DEL_ARR_F(tmp_var_idx_m);
1468
1469                                 /* -c * a_{nt}^k */
1470                                 idx = ILPVAR_IDX_DEAD(ba, na, node_tp_idx, t);
1471                                 lpp_set_factor_fast(lpp, cst, na->ilp_vars.a[idx], 0.0 - (double)(na->n_consumer));
1472
1473                         }
1474                 }
1475         }
1476         ilp_timer_pop();
1477         DBG((env->dbg, LEVEL_1, "\t%u alive nodes constraints (%g sec)\n",
1478                 num_cst, ilp_timer_elapsed_usec(t_cst) / 1000000.0));
1479 }
1480
1481 /**
1482  * Create ILP alive nodes constraints for live-in nodes:
1483  * - set variable a_{nt}^k to 1 if nodes n is alive at step t on unit k
1484  */
1485 static void create_alive_livein_nodes_constraint(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1486         char                  buf[1024];
1487         ilp_livein_node_t     *livein;
1488         unsigned              num_cst = 0;
1489         ilpsched_block_attr_t *ba     = get_ilpsched_block_attr(block_node);
1490         lc_timer_t            *t_cst  = lc_timer_register("beilpsched_cst_alive_livein_nodes", "create alive livein nodes constraints");
1491
1492         ilp_timer_push(t_cst);
1493         /* for each node */
1494         foreach_pset(ba->livein_nodes, livein) {
1495                 ir_node              *irn  = livein->irn;
1496                 be_ilpsched_irn_t    *node = get_ilpsched_irn(env, irn);
1497                 ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
1498                 unsigned             t;
1499
1500                 /* check check all time steps: 0 <= t < max_alive_steps */
1501                 for (t = 0; t < livein->max_alive_steps; ++t) {
1502                         int node_tp_idx;
1503
1504                         /* for all unit types available for this node */
1505                         for (node_tp_idx = na->n_unit_types - 1; node_tp_idx >= 0; --node_tp_idx) {
1506                                 const ir_edge_t *edge;
1507                                 unsigned idx;
1508                                 int      cst, num_block_user;
1509                                 int      *tmp_var_idx_m = NEW_ARR_F(int, 0);
1510
1511                                 /* check the number of consumer scheduled so far */
1512                                 num_block_user = 0;
1513                                 foreach_out_edge(irn, edge) {
1514                                         ir_node              *user = get_edge_src_irn(edge);
1515                                         be_ilpsched_irn_t    *cons;
1516                                         ilpsched_node_attr_t *ca;
1517                                         int                  tp_idx;
1518                                         unsigned             tm, tm_max;
1519
1520                                         /* check only users within current block */
1521                                         if (get_nodes_block(user) != block_node->irn)
1522                                                 continue;
1523
1524                                         num_block_user++;
1525                                         cons = get_ilpsched_irn(env, user);
1526                                         ca   = get_ilpsched_node_attr(cons);
1527
1528                                         tm_max = MIN(ca->alap - 1, t);
1529                                         for (tp_idx = ca->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1530                                                 for (tm = ca->asap - 1; tm <= tm_max; ++tm) {
1531                                                         int idx = ILPVAR_IDX(ca, tp_idx, tm);
1532                                                         ARR_APP1(int, tmp_var_idx_m, ca->ilp_vars.x[idx]);
1533                                                 }
1534                                         }
1535                                 }
1536
1537                                 snprintf(buf, sizeof(buf), "alive_livein_node_cst_%u_n%u_%s",
1538                                         t, get_irn_idx(irn), na->type_info[node_tp_idx].tp->name);
1539                                 cst = lpp_add_cst_uniq(lpp, buf, lpp_greater, (double)num_block_user);
1540                                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1541                                 num_cst++;
1542
1543                                 /* sum(scheduled users) */
1544                                 if (ARR_LEN(tmp_var_idx_m) > 0)
1545                                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx_m, ARR_LEN(tmp_var_idx_m), 1.0);
1546                                 DEL_ARR_F(tmp_var_idx_m);
1547
1548                                 /* + c * a_{nt}^k */
1549                                 idx = node_tp_idx * livein->max_alive_steps + t;
1550                                 lpp_set_factor_fast(lpp, cst, livein->a[idx], (double)(num_block_user));
1551                         }
1552                 }
1553         }
1554         ilp_timer_pop();
1555         DBG((env->dbg, LEVEL_1, "\t%u alive livein nodes constraints (%g sec)\n",
1556                 num_cst, ilp_timer_elapsed_usec(t_cst) / 1000000.0));
1557 }
1558
1559 /**
1560  * Create ILP pressure constraints, based on alive nodes:
1561  * - add additional costs to objective function if a node is scheduled
1562  *   on a unit although all units of this type are currently occupied
1563  */
1564 static void create_pressure_alive_constraint(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1565         char                  buf[1024];
1566         ir_node               *cur_irn;
1567         unsigned              num_cst = 0;
1568         ilpsched_block_attr_t *ba     = get_ilpsched_block_attr(block_node);
1569         lc_timer_t            *t_cst  = lc_timer_register("beilpsched_cst_pressure", "create pressure constraints");
1570
1571         ilp_timer_push(t_cst);
1572         /* y_{nt}^k is set for each node and timestep and unit type */
1573         foreach_linked_irns(ba->head_ilp_nodes, cur_irn) {
1574                 unsigned             cur_idx   = get_irn_idx(cur_irn);
1575                 be_ilpsched_irn_t    *cur_node = get_ilpsched_irn(env, cur_irn);
1576                 ilpsched_node_attr_t *cur_na   = get_ilpsched_node_attr(cur_node);
1577                 int                  glob_type_idx;
1578
1579                 /* we ignore nodes assigned to DUMMY unit here */
1580                 if (cur_na->is_dummy_node)
1581                         continue;
1582
1583                 /* for all types */
1584                 for (glob_type_idx = env->cpu->n_unit_types - 1; glob_type_idx >= 0; --glob_type_idx) {
1585                         be_execution_unit_type_t *cur_tp   = &env->cpu->unit_types[glob_type_idx];
1586                         int                      cur_tp_idx;
1587                         unsigned                 t;
1588
1589                         /* BEWARE: the DUMMY unit types is not in CPU, so it's skipped automatically */
1590
1591                         /* check if node can be executed on this unit type */
1592                         cur_tp_idx = is_valid_unit_type_for_node(cur_tp, cur_node);
1593                         if (cur_tp_idx < 0)
1594                                 continue;
1595
1596                         /* check all time_steps at which the current node can be scheduled */
1597                         for (t = cur_na->asap - 1; t <= cur_na->alap - 1; ++t) {
1598                                 int     cst, y_idx;
1599                                 ir_node *irn;
1600                                 int     *tmp_var_idx = NEW_ARR_F(int, 0);
1601                                 ilp_livein_node_t *livein;
1602
1603                                 snprintf(buf, sizeof(buf), "pressure_cst_n%u_%u_%s", cur_idx, t, cur_tp->name);
1604                                 cst = lpp_add_cst_uniq(lpp, buf, lpp_less, (double)(cur_tp->n_units - 1));
1605                                 DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1606                                 num_cst++;
1607
1608                                 /* - accumulate all nodes alive at point t on unit type k */
1609                                 foreach_linked_irns(ba->head_ilp_nodes, irn) {
1610                                         be_ilpsched_irn_t    *node = get_ilpsched_irn(env, irn);
1611                                         ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
1612                                         int                  a_idx, tp_idx;
1613
1614                                         /* check if node can be alive here */
1615                                         if (t < na->asap - 1)
1616                                                 continue;
1617
1618                                         tp_idx = is_valid_unit_type_for_node(cur_tp, node);
1619
1620                                         /* current type is not suitable */
1621                                         if (tp_idx < 0)
1622                                                 continue;
1623
1624                                         a_idx = ILPVAR_IDX_DEAD(ba, na, tp_idx, t);
1625                                         ARR_APP1(int, tmp_var_idx, na->ilp_vars.a[a_idx]);
1626                                 }
1627                                 /* do the same for livein nodes */
1628                                 foreach_pset(ba->livein_nodes, livein) {
1629                                         ir_node              *irn  = livein->irn;
1630                                         be_ilpsched_irn_t    *node = get_ilpsched_irn(env, irn);
1631                                         int                  a_idx, tp_idx;
1632
1633                                         /* check if node can be alive here */
1634                                         if (t >= livein->max_alive_steps)
1635                                                 continue;
1636
1637                                         tp_idx = is_valid_unit_type_for_node(cur_tp, node);
1638
1639                                         /* current type is not suitable */
1640                                         if (tp_idx < 0)
1641                                                 continue;
1642
1643                                         a_idx = tp_idx * livein->max_alive_steps + t;
1644                                         ARR_APP1(int, tmp_var_idx, livein->a[a_idx]);
1645                                 }
1646
1647                                 if (ARR_LEN(tmp_var_idx) > 0)
1648                                         lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx, ARR_LEN(tmp_var_idx), 1.0);
1649                                 DEL_ARR_F(tmp_var_idx);
1650
1651                                 /* - num_nodes * y_{nt}^k */
1652                                 y_idx = ILPVAR_IDX(cur_na, cur_tp_idx, t);
1653                                 lpp_set_factor_fast(lpp, cst, cur_na->ilp_vars.y[y_idx], -1.0);
1654                         }
1655                 }
1656         }
1657         ilp_timer_pop();
1658         DBG((env->dbg, LEVEL_1, "\t%u pressure constraints (%g sec)\n",
1659                 num_cst, ilp_timer_elapsed_usec(t_cst) / 1000000.0));
1660 }
1661
1662 /**
1663  * Create ILP branch constraints:
1664  * Assure, alle nodes are scheduled prior to cfg op.
1665  */
1666 static void create_branch_constraint(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1667         char                  buf[1024];
1668         ir_node               *cur_irn, *cfop;
1669         unsigned              num_cst          = 0;
1670         unsigned              num_non_branches = 0;
1671         ilpsched_block_attr_t *ba     = get_ilpsched_block_attr(block_node);
1672         lc_timer_t            *t_cst  = lc_timer_register("beilpsched_cst_branch", "create branch constraints");
1673
1674         ilp_timer_push(t_cst);
1675         cfop = NULL;
1676         /* determine number of non-branch nodes and the one and only branch node */
1677         foreach_linked_irns(ba->head_ilp_nodes, cur_irn) {
1678                 switch (get_irn_opcode(cur_irn)) {
1679                         case iro_Phi:
1680                         case iro_Start:
1681                         case iro_End:
1682                         case iro_Proj:
1683                         case iro_Bad:
1684                         case iro_Unknown:
1685                                 num_non_branches++;
1686                                 break;
1687                         default:
1688                                 if (is_cfop(cur_irn)) {
1689                                         assert(cfop == NULL && "Highlander - there can be only one to be constrained");
1690                                         cfop = cur_irn;
1691                                 }
1692                                 else {
1693                                         num_non_branches++;
1694                                 }
1695                                 break;
1696                 }
1697         }
1698
1699         if (cfop) {
1700                 be_ilpsched_irn_t    *cf_node = get_ilpsched_irn(env, cfop);
1701                 ilpsched_node_attr_t *cf_na   = get_ilpsched_node_attr(cf_node);
1702                 unsigned t;
1703
1704                 /* for each time step */
1705                 for (t = cf_na->asap - 1; t <= cf_na->alap - 1; ++t) {
1706                         int *non_branch_vars, *branch_vars;
1707                         int cst;
1708
1709                         snprintf(buf, sizeof(buf), "branch_cst_%u_n%u", t, get_irn_idx(cfop));
1710                         cst = lpp_add_cst_uniq(lpp, buf, lpp_greater, 0.0);
1711                         DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1712                         num_cst++;
1713
1714                         /* sum(overall non branches: n)x_{nt}^k - sum(overall branches: b)(num_non_branches * x_{bt}^k >= 0) */
1715                         non_branch_vars = NEW_ARR_F(int, 0);
1716                         branch_vars     = NEW_ARR_F(int, 0);
1717                         foreach_linked_irns(ba->head_ilp_nodes, cur_irn) {
1718                                 be_ilpsched_irn_t    *node = get_ilpsched_irn(env, cur_irn);
1719                                 ilpsched_node_attr_t *na   = get_ilpsched_node_attr(node);
1720                                 int                  tp_idx;
1721
1722                                 if (cur_irn == cfop) {
1723                                         /* for all unit types available for this node */
1724                                         for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1725                                                 unsigned idx = ILPVAR_IDX(na, tp_idx, t);
1726                                                 ARR_APP1(int, branch_vars, na->ilp_vars.x[idx]);
1727                                         }
1728                                 }
1729                                 else {
1730                                         /* sum up all possible schedule points for this node upto current timestep */
1731                                         for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1732                                                 unsigned tn;
1733                                                 unsigned tmax = MIN(t, na->alap - 1);
1734
1735                                                 for (tn = na->asap - 1; tn <= tmax; ++tn) {
1736                                                         unsigned idx = ILPVAR_IDX(na, tp_idx, tn);
1737                                                         ARR_APP1(int, non_branch_vars, na->ilp_vars.x[idx]);
1738                                                 }
1739                                         }
1740                                 }
1741
1742                         }
1743
1744                         if (ARR_LEN(non_branch_vars) > 0)
1745                                 lpp_set_factor_fast_bulk(lpp, cst, non_branch_vars, ARR_LEN(non_branch_vars), 1.0);
1746                         if (ARR_LEN(branch_vars) > 0)
1747                                 lpp_set_factor_fast_bulk(lpp, cst, branch_vars, ARR_LEN(branch_vars), 0.0 - (double)num_non_branches);
1748
1749                         DEL_ARR_F(branch_vars);
1750                         DEL_ARR_F(non_branch_vars);
1751                 }
1752         }
1753         ilp_timer_pop();
1754         DBG((env->dbg, LEVEL_1, "\t%u branch constraints (%g sec)\n",
1755                 num_cst, ilp_timer_elapsed_usec(t_cst) / 1000000.0));
1756 }
1757
1758 #if 0
1759 static void create_proj_keep_constraints(be_ilpsched_env_t *env, lpp_t *lpp, be_ilpsched_irn_t *block_node) {
1760         char                  buf[1024];
1761         ir_node               *irn;
1762         unsigned              num_cst = 0;
1763         ilpsched_block_attr_t *ba     = get_ilpsched_block_attr(block_node);
1764         lc_timer_t            *t_cst  = lc_timer_register("beilpsched_cst_projkeep", "create proj and keep constraints");
1765
1766         ilp_timer_push(t_cst);
1767         /* check all nodes */
1768         foreach_linked_irns(ba->head_ilp_nodes, irn) {
1769                 be_ilpsched_irn_t    *node;
1770                 ilpsched_node_attr_t *na;
1771                 unsigned             t;
1772                 ir_node              **pk;
1773
1774                 /* only mode_T nodes can have Projs and Keeps assigned */
1775                 if (get_irn_mode(irn) != mode_T)
1776                         continue;
1777
1778                 node = get_ilpsched_irn(env, irn);
1779                 na   = get_ilpsched_node_attr(node);
1780
1781                 /* check if has some Projs and Keeps assigned */
1782                 if (! na->projkeeps)
1783                         continue;
1784
1785                 /* we can run only once over the queue, so preserve the nodes */
1786                 pk = NEW_ARR_F(ir_node *, 0);
1787                 while (! waitq_empty(na->projkeeps))
1788                         ARR_APP1(ir_node *, pk, waitq_get(na->projkeeps));
1789                 del_waitq(na->projkeeps);
1790                 na->projkeeps = NULL;
1791
1792                 /* for all time steps at which this node can be scheduled */
1793                 for (t = na->asap - 1; t <= na->alap - 1; ++t) {
1794                         int cst, tp_idx, i;
1795                         int *tmp_var_idx_n = NEW_ARR_F(int, 0);
1796
1797                         /* add the constraint, assure, that a node is always scheduled along with it's Projs and Keeps */
1798                         snprintf(buf, sizeof(buf), "projkeep_cst_n%u_%u", get_irn_idx(irn), t);
1799                         cst = lpp_add_cst_uniq(lpp, buf, lpp_equal, 0.0);
1800                         DBG((env->dbg, LEVEL_2, "added constraint %s\n", buf));
1801                         num_cst++;
1802
1803                         /* sum up scheduling variables for this time step */
1804                         for (tp_idx = na->n_unit_types - 1; tp_idx >= 0; --tp_idx) {
1805                                 int idx = ILPVAR_IDX(na, tp_idx, t);
1806                                 ARR_APP1(int, tmp_var_idx_n, na->ilp_vars.x[idx]);
1807                         }
1808
1809                         if (ARR_LEN(tmp_var_idx_n) > 0)
1810                                 lpp_set_factor_fast_bulk(lpp, cst, tmp_var_idx_n, ARR_LEN(tmp_var_idx_n), (double)(ARR_LEN(pk)));
1811                         DEL_ARR_F(tmp_var_idx_n);
1812
1813                         /* subtract all Proj and Keep variables for this step */
1814                         for (i = ARR_LEN(pk) - 1; i >= 0; --i) {
1815                                 be_ilpsched_irn_t    *pk_node = get_ilpsched_irn(env, pk[i]);
1816                                 ilpsched_node_attr_t *pk_na   = get_ilpsched_node_attr(pk_node);
1817                                 int                  pk_tp_idx;
1818
1819                                 for (pk_tp_idx = pk_na->n_unit_types - 1; pk_tp_idx >= 0; --pk_tp_idx) {
1820                                         int idx = ILPVAR_IDX(pk_na, pk_tp_idx, t);
1821                                         lpp_set_factor_fast(lpp, cst, pk_na->ilp_vars.x[idx], -1.0);
1822                                 }
1823                         }
1824                 }
1825         }
1826         ilp_timer_pop();
1827         DBG((env->dbg, LEVEL_1, "\t%u Proj and Keep constraints (%g sec)\n",
1828                 num_cst, ilp_timer_elapsed_usec(t_cst) / 1000000.0));
1829 }
1830 #endif /* if 0 */
1831
1832 /***************************************************
1833  *  _____ _      _____                    _
1834  * |_   _| |    |  __ \                  (_)
1835  *   | | | |    | |__) |  _ __ ___   __ _ _ _ __
1836  *   | | | |    |  ___/  | '_ ` _ \ / _` | | '_ \
1837  *  _| |_| |____| |      | | | | | | (_| | | | | |
1838  * |_____|______|_|      |_| |_| |_|\__,_|_|_| |_|
1839  *
1840  ***************************************************/
1841
1842 /**
1843  * Create the ilp (add variables, build constraints, solve, build schedule from solution).
1844  */
1845 static void create_ilp(ir_node *block, void *walk_env) {
1846         be_ilpsched_env_t     *env           = walk_env;
1847         be_ilpsched_irn_t     *block_node    = get_ilpsched_irn(env, block);
1848         ilpsched_block_attr_t *ba            = get_ilpsched_block_attr(block_node);
1849         FILE                  *logfile       = NULL;
1850         lpp_t                 *lpp           = NULL;
1851         int                   need_heur      = 0;
1852         struct obstack        var_obst;
1853         char                  name[1024];
1854
1855         DBG((env->dbg, 255, "\n\n\n=========================================\n"));
1856         DBG((env->dbg, 255, "  ILP Scheduling for %+F\n", block));
1857         DBG((env->dbg, 255, "=========================================\n\n"));
1858
1859         DBG((env->dbg, LEVEL_1, "Creating ILP Variables for nodes in %+F (%u interesting nodes, %u max steps)\n",
1860                 block, ba->n_interesting_nodes, ba->max_steps));
1861
1862         /* notify backend and get block environment */
1863         env->block_env = be_ilp_sched_init_block_ilp_schedule(env->sel, block);
1864
1865         /* if we have less than two interesting nodes, there is no need to create the ILP */
1866         if (ba->n_interesting_nodes > 1) {
1867                 double fact_var        = ba->n_interesting_nodes > 25 ? 2.3 : 3;
1868                 double fact_cst        = ba->n_interesting_nodes > 25 ? 3   : 4.5;
1869                 int    base_num        = ba->n_interesting_nodes * ba->n_interesting_nodes;
1870                 int    estimated_n_var = (int)((double)base_num * fact_var);
1871                 int    estimated_n_cst = (int)((double)base_num * fact_cst);
1872
1873                 DBG((env->dbg, LEVEL_1, "Creating LPP with estimated numbers: %d vars, %d cst\n",
1874                         estimated_n_var, estimated_n_cst));
1875                 (void) estimated_n_var;
1876
1877                 /* set up the LPP object */
1878                 snprintf(name, sizeof(name), "ilp scheduling IRG %s", get_entity_ld_name(get_irg_entity(env->irg)));
1879
1880                 lpp = new_lpp_userdef(
1881                         (const char *)name,
1882                         lpp_minimize,
1883                         estimated_n_cst,     /* num vars */
1884                         estimated_n_cst + 1, /* num cst */
1885                         1.3);                /* grow factor */
1886                 obstack_init(&var_obst);
1887
1888                 /* create ILP variables */
1889                 create_variables(env, lpp, block_node, &var_obst);
1890
1891                 /* create ILP constraints */
1892                 DBG((env->dbg, LEVEL_1, "Creating constraints for nodes in %+F:\n", block));
1893                 create_assignment_and_precedence_constraints(env, lpp, block_node);
1894                 create_ressource_constraints(env, lpp, block_node);
1895                 create_bundle_constraints(env, lpp, block_node);
1896                 create_branch_constraint(env, lpp, block_node);
1897                 //create_proj_keep_constraints(env, lpp, block_node);
1898
1899                 if (env->opts->regpress) {
1900                         create_alive_nodes_constraint(env, lpp, block_node);
1901                         create_alive_livein_nodes_constraint(env, lpp, block_node);
1902                         create_pressure_alive_constraint(env, lpp, block_node);
1903                 }
1904
1905                 DBG((env->dbg, LEVEL_1, "ILP to solve: %u variables, %u constraints\n", lpp->var_next, lpp->cst_next));
1906
1907                 /* debug stuff, dump lpp when debugging is on  */
1908                 DEBUG_ONLY({
1909                         if (firm_dbg_get_mask(env->dbg) > 1) {
1910                                 char buf[1024];
1911                                 FILE *f;
1912
1913                                 snprintf(buf, sizeof(buf), "lpp_block_%lu.txt", get_irn_node_nr(block));
1914                                 f = fopen(buf, "w");
1915                                 lpp_dump_plain(lpp, f);
1916                                 fclose(f);
1917                                 snprintf(buf, sizeof(buf), "lpp_block_%lu.mps", get_irn_node_nr(block));
1918                                 lpp_dump(lpp, buf);
1919                         }
1920                 })
1921
1922                 /* set solve time limit */
1923                 lpp_set_time_limit(lpp, env->opts->time_limit);
1924
1925                 /* set logfile if requested */
1926                 if (strlen(env->opts->log_file) > 0) {
1927                         if (strcasecmp(env->opts->log_file, "stdout") == 0)
1928                                 lpp_set_log(lpp, stdout);
1929                         else if (strcasecmp(env->opts->log_file, "stderr") == 0)
1930                                 lpp_set_log(lpp, stderr);
1931                         else {
1932                                 logfile = fopen(env->opts->log_file, "w");
1933                                 if (! logfile)
1934                                         fprintf(stderr, "Could not open logfile '%s'! Logging disabled.\n", env->opts->log_file);
1935                                 else
1936                                         lpp_set_log(lpp, logfile);
1937                         }
1938                 }
1939
1940                 /* solve the ILP */
1941                 lpp_solve_net(lpp, env->main_env->options->ilp_server, env->main_env->options->ilp_solver);
1942
1943                 if (logfile)
1944                         fclose(logfile);
1945
1946                 /* check for valid solution */
1947                 if (! lpp_is_sol_valid(lpp)) {
1948                         DEBUG_ONLY({
1949                             char buf[1024];
1950                             FILE *f;
1951
1952                             if (firm_dbg_get_mask(env->dbg) >= 2) {
1953                               snprintf(buf, sizeof(buf), "lpp_block_%lu.infeasible.txt", get_irn_node_nr(block));
1954                               f = fopen(buf, "w");
1955                               lpp_dump_plain(lpp, f);
1956                               fclose(f);
1957                               snprintf(buf, sizeof(buf), "lpp_block_%lu.infeasible.mps", get_irn_node_nr(block));
1958                               lpp_dump(lpp, buf);
1959                               dump_ir_block_graph(env->irg, "-infeasible");
1960                            }
1961                         })
1962
1963                         ir_fprintf(stderr, "ILP found no solution within time (%+F, %+F), falling back to heuristics.\n", block, env->irg);
1964                         need_heur = 1;
1965                 }
1966
1967                 DBG((env->dbg, LEVEL_1, "\nSolution:\n"));
1968                 DBG((env->dbg, LEVEL_1, "\tsend time: %g sec\n", lpp->send_time / 1000000.0));
1969                 DBG((env->dbg, LEVEL_1, "\treceive time: %g sec\n", lpp->recv_time / 1000000.0));
1970                 DBG((env->dbg, LEVEL_1, "\tmatrix: %u elements, density %.2f%%, size %.2fMB\n", lpp->n_elems, lpp->density, (double)lpp->matrix_mem / 1024.0 / 1024.0));
1971                 DBG((env->dbg, LEVEL_1, "\titerations: %d\n", lpp->iterations));
1972                 DBG((env->dbg, LEVEL_1, "\tsolution time: %g\n", lpp->sol_time));
1973                 DBG((env->dbg, LEVEL_1, "\tobjective function: %g\n", LPP_VALUE_IS_0(lpp->objval) ? 0.0 : lpp->objval));
1974                 DBG((env->dbg, LEVEL_1, "\tbest bound: %g\n", LPP_VALUE_IS_0(lpp->best_bound) ? 0.0 : lpp->best_bound));
1975
1976                 DBG((env->dbg, LEVEL_1, "variables used %u bytes\n", obstack_memory_used(&var_obst)));
1977         }
1978
1979         /* apply solution */
1980         be_stat_ev("nodes", ba->block_last_idx);
1981         be_stat_ev("vars", lpp ? lpp->var_next : 0);
1982         be_stat_ev("csts", lpp ? lpp->cst_next : 0);
1983         if (need_heur) {
1984                         be_stat_ev("time", -1);
1985                         be_stat_ev_dbl("opt", 0.0);
1986                 list_sched_single_block(env->birg, block, env->be_opts);
1987         }
1988         else {
1989                 if (lpp) {
1990                         double opt = lpp->sol_state == lpp_optimal ? 100.0 : 100.0 * lpp->best_bound / lpp->objval;
1991                         be_stat_ev_dbl("time", lpp->sol_time);
1992                         be_stat_ev_dbl("opt", opt);
1993                 }
1994                 else {
1995                         be_stat_ev_dbl("time", 0.0);
1996                         be_stat_ev_dbl("opt", 100.0);
1997                 }
1998                 apply_solution(env, lpp, block);
1999         }
2000
2001         if (lpp)
2002                 free_lpp(lpp);
2003
2004         /* notify backend */
2005         be_ilp_sched_finish_block_ilp_schedule(env->sel, block, env->block_env);
2006 }
2007
2008 /**
2009  * Perform ILP scheduling on the given irg.
2010  */
2011 void be_ilp_sched(const be_irg_t *birg, be_options_t *be_opts) {
2012         be_ilpsched_env_t          env;
2013         const char                 *name     = "be ilp scheduling";
2014         ir_graph                   *irg      = be_get_birg_irg(birg);
2015         const arch_env_t           *arch_env = be_get_birg_arch_env(birg);
2016         const arch_isa_t           *isa      = arch_env->isa;
2017         const ilp_sched_selector_t *sel      = isa->impl->get_ilp_sched_selector(isa);
2018
2019         FIRM_DBG_REGISTER(env.dbg, "firm.be.sched.ilp");
2020
2021         stat_ev_ctx_push("ilpsched");
2022
2023 //      firm_dbg_set_mask(env.dbg, 1);
2024
2025         env.irg_env    = be_ilp_sched_init_irg_ilp_schedule(sel, irg);
2026         env.sel        = sel;
2027         env.irg        = irg;
2028         env.height     = heights_new(irg);
2029         env.main_env   = birg->main_env;
2030         env.arch_env   = arch_env;
2031         env.cpu        = arch_isa_get_machine(arch_env->isa);
2032         env.opts       = &ilp_opts;
2033         env.birg       = birg;
2034         env.be_opts    = be_opts;
2035         phase_init(&env.ph, name, env.irg, PHASE_DEFAULT_GROWTH, init_ilpsched_irn, NULL);
2036
2037         /* assign a unique per block number to all interesting nodes */
2038         irg_walk_in_or_dep_graph(env.irg, NULL, build_block_idx, &env);
2039
2040         /*
2041                 The block indices are completely build after the walk,
2042                 now we can allocate the bitsets (size depends on block indices)
2043                 for all nodes.
2044         */
2045         phase_reinit_irn_data(&env.ph);
2046
2047         /* Collect all root nodes (having no user in their block) and calculate ASAP. */
2048         irg_walk_in_or_dep_blkwise_graph(env.irg, collect_alap_root_nodes, calculate_irn_asap, &env);
2049
2050         /* Calculate ALAP of all irns */
2051         irg_block_walk_graph(env.irg, NULL, calculate_block_alap, &env);
2052
2053         /* We refine the {ASAP(n), ALAP(n)} interval and fix the time steps for Projs and Keeps */
2054         irg_walk_in_or_dep_blkwise_graph(env.irg, NULL, refine_asap_alap_times, &env);
2055
2056         /* perform ILP scheduling */
2057         irg_block_walk_graph(env.irg, NULL, create_ilp, &env);
2058
2059         DEBUG_ONLY(
2060                 if (firm_dbg_get_mask(env.dbg)) {
2061                         phase_stat_t stat;
2062                         phase_stat_t *stat_ptr = phase_stat(&env.ph, &stat);
2063
2064                         fprintf(stderr, "Phase used: %u bytes\n", stat_ptr->overall_bytes);
2065                 }
2066         );
2067
2068         /* free data allocated dynamically */
2069         irg_block_walk_graph(env.irg, NULL, clear_unwanted_data, &env);
2070
2071         /* free all allocated object */
2072         phase_free(&env.ph);
2073         heights_free(env.height);
2074
2075         /* notify backend */
2076         be_ilp_sched_finish_irg_ilp_schedule(sel, birg->irg, env.irg_env);
2077
2078         stat_ev_ctx_pop("ilpsched");
2079 }
2080
2081 /**
2082  * Register ILP scheduler options.
2083  */
2084 void be_init_ilpsched(void)
2085 {
2086         lc_opt_entry_t *be_grp = lc_opt_get_grp(firm_opt_get_root(), "be");
2087         lc_opt_entry_t *sched_grp = lc_opt_get_grp(be_grp, "ilpsched");
2088
2089         lc_opt_add_table(sched_grp, ilpsched_option_table);
2090 }
2091
2092 #else /* WITH_ILP */
2093
2094 static INLINE void some_picky_compiler_do_not_allow_empty_files(void)
2095 {}
2096
2097 #endif /* WITH_ILP */