becopyilp: fix size_reduction code
[libfirm] / ir / be / becopyilp.c
1 /*
2  * Copyright (C) 1995-2008 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       Common stuff used by all ILP formulations.
23  * @author      Daniel Grund
24  * @date        28.02.2006
25  */
26 #include "config.h"
27
28 #include <stdbool.h>
29
30 #include "irtools.h"
31 #include "irprintf.h"
32
33 #include "bestatevent.h"
34 #include "beirg.h"
35 #include "bemodule.h"
36 #include "error.h"
37
38 #include "lpp.h"
39
40 #include "lc_opts.h"
41 #include "lc_opts_enum.h"
42
43 #define DUMP_ILP 1
44 #define DUMP_SOL 2
45
46 static int time_limit = 60;
47 static int solve_log  = 0;
48 static unsigned dump_flags = 0;
49
50 static const lc_opt_enum_mask_items_t dump_items[] = {
51         { "ilp",   DUMP_ILP },
52         { "sol",   DUMP_SOL },
53         { NULL,    0 }
54 };
55
56 static lc_opt_enum_mask_var_t dump_var = {
57         &dump_flags, dump_items
58 };
59
60 static const lc_opt_table_entry_t options[] = {
61         LC_OPT_ENT_INT      ("limit", "time limit for solving in seconds (0 for unlimited)", &time_limit),
62         LC_OPT_ENT_BOOL     ("log",   "show ilp solving log",              &solve_log),
63         LC_OPT_ENT_ENUM_MASK("dump",  "dump flags",             &dump_var),
64         LC_OPT_LAST
65 };
66
67 BE_REGISTER_MODULE_CONSTRUCTOR(be_init_copyilp)
68 void be_init_copyilp(void)
69 {
70         lc_opt_entry_t *be_grp = lc_opt_get_grp(firm_opt_get_root(), "be");
71         lc_opt_entry_t *ra_grp = lc_opt_get_grp(be_grp, "ra");
72         lc_opt_entry_t *chordal_grp = lc_opt_get_grp(ra_grp, "chordal");
73         lc_opt_entry_t *co_grp = lc_opt_get_grp(chordal_grp, "co");
74         lc_opt_entry_t *ilp_grp = lc_opt_get_grp(co_grp, "ilp");
75
76         lc_opt_add_table(ilp_grp, options);
77 }
78
79 #include "becopyilp_t.h"
80 #include "beifg.h"
81
82 /******************************************************************************
83     _____ _                        _            _   _
84    / ____(_)                      | |          | | (_)
85   | (___  _ _______   _ __ ___  __| |_   _  ___| |_ _  ___  _ __
86    \___ \| |_  / _ \ | '__/ _ \/ _` | | | |/ __| __| |/ _ \| '_ \
87    ____) | |/ /  __/ | | |  __/ (_| | |_| | (__| |_| | (_) | | | |
88   |_____/|_/___\___| |_|  \___|\__,_|\__,_|\___|\__|_|\___/|_| |_|
89
90  *****************************************************************************/
91
92
93 size_red_t *new_size_red(copy_opt_t *co)
94 {
95         size_red_t *res = XMALLOC(size_red_t);
96
97         res->co = co;
98         res->all_removed = pset_new_ptr_default();
99         res->col_suff = NULL;
100         obstack_init(&res->ob);
101
102         return res;
103 }
104
105 /**
106  * Checks if a node is simplicial in the graph heeding the already removed nodes.
107  */
108 static inline bool sr_is_simplicial(size_red_t *sr, const ir_node *ifn)
109 {
110         be_ifg_t *ifg  = sr->co->cenv->ifg;
111         neighbours_iter_t iter;
112         ir_node **all  = ALLOCAN(ir_node*, be_ifg_degree(ifg, ifn));
113         ir_node  *curr;
114         int       size = 0;
115         int       i;
116         int       o;
117
118         /* get all non-removed neighbors */
119         be_ifg_foreach_neighbour(ifg, &iter, ifn, curr)
120                 if (!sr_is_removed(sr, curr))
121                         all[size++] = curr;
122
123         /* check if these form a clique */
124         for (i=0; i<size; ++i)
125                 for (o=i+1; o<size; ++o)
126                         if (!be_ifg_connected(ifg, all[i], all[o]))
127                                 return false;
128
129         /* all edges exist so this is a clique */
130         return true;
131 }
132
133 void sr_remove(size_red_t *sr)
134 {
135         ir_node *irn;
136         bool redo = true;
137         const be_ifg_t *ifg = sr->co->cenv->ifg;
138         nodes_iter_t iter;
139
140         while (redo) {
141                 redo = false;
142                 be_ifg_foreach_node(ifg, &iter, irn) {
143                         const arch_register_req_t *req = arch_get_irn_register_req(irn);
144                         coloring_suffix_t *cs;
145
146                         if (arch_register_req_is(req, limited) || sr_is_removed(sr, irn))
147                                 continue;
148                         if (co_gs_is_optimizable(sr->co, irn))
149                                 continue;
150                         if (!sr_is_simplicial(sr, irn))
151                                 continue;
152
153                         cs = OALLOC(&sr->ob, coloring_suffix_t);
154                         cs->irn = irn;
155                         cs->next = sr->col_suff;
156                         sr->col_suff = cs;
157
158                         pset_insert_ptr(sr->all_removed, irn);
159
160                         redo = true;
161                 }
162         }
163 }
164
165 void sr_reinsert(size_red_t *sr)
166 {
167         coloring_suffix_t *cs;
168         ir_graph *irg        = sr->co->irg;
169         be_ifg_t *ifg        = sr->co->cenv->ifg;
170         unsigned  n_regs     = arch_register_class_n_regs(sr->co->cls);
171         unsigned *possible_cols;
172         unsigned *allocatable_cols;
173
174         rbitset_alloca(allocatable_cols, n_regs);
175         be_set_allocatable_regs(irg, sr->co->cls, allocatable_cols);
176
177         rbitset_alloca(possible_cols, n_regs);
178         neighbours_iter_t iter;
179
180         /* color the removed nodes in right order */
181         for (cs = sr->col_suff; cs; cs = cs->next) {
182                 unsigned free_col;
183                 ir_node *other;
184                 ir_node *irn = cs->irn;
185
186                 rbitset_copy(possible_cols, allocatable_cols, n_regs);
187
188                 /* get free color by inspecting all neighbors */
189                 be_ifg_foreach_neighbour(ifg, &iter, irn, other) {
190                         /* only inspect nodes which are in graph right now */
191                         if (sr_is_removed(sr, other))
192                                 continue;
193                         rbitset_clear(possible_cols, get_irn_col(other));
194                 }
195
196                 /* now all bits not set are possible colors */
197                 free_col = (unsigned)rbitset_next(possible_cols, 0, true);
198                 assert(!rbitset_is_empty(possible_cols, n_regs) && "No free color found. This can not be.");
199                 set_irn_col(sr->co, irn, free_col);
200                 pset_remove_ptr(sr->all_removed, irn); /* irn is back in graph again */
201         }
202 }
203
204 void free_size_red(size_red_t *sr)
205 {
206         del_pset(sr->all_removed);
207         obstack_free(&sr->ob, NULL);
208         free(sr);
209 }
210
211 /******************************************************************************
212     _____                      _        _____ _      _____
213    / ____|                    (_)      |_   _| |    |  __ \
214   | |  __  ___ _ __   ___ _ __ _  ___    | | | |    | |__) |
215   | | |_ |/ _ \ '_ \ / _ \ '__| |/ __|   | | | |    |  ___/
216   | |__| |  __/ | | |  __/ |  | | (__   _| |_| |____| |
217    \_____|\___|_| |_|\___|_|  |_|\___| |_____|______|_|
218
219  *****************************************************************************/
220
221 #include <stdio.h>
222
223 ilp_env_t *new_ilp_env(copy_opt_t *co, ilp_callback build, ilp_callback apply, void *env)
224 {
225         ilp_env_t *res = XMALLOC(ilp_env_t);
226
227         res->co         = co;
228         res->build      = build;
229         res->apply      = apply;
230         res->env        = env;
231         res->sr         = new_size_red(co);
232
233         return res;
234 }
235
236 lpp_sol_state_t ilp_go(ilp_env_t *ienv)
237 {
238         ir_graph     *irg     = ienv->co->irg;
239         be_options_t *options = be_get_irg_options(irg);
240
241         sr_remove(ienv->sr);
242
243         ienv->build(ienv);
244
245         if (dump_flags & DUMP_ILP) {
246                 char buf[128];
247                 FILE *f;
248
249                 ir_snprintf(buf, sizeof(buf), "%F_%s-co.ilp", irg,
250                             ienv->co->cenv->cls->name);
251                 f = fopen(buf, "wt");
252                 if (f == NULL) {
253                         panic("Couldn't open '%s' for writing", buf);
254                 }
255                 lpp_dump_plain(ienv->lp, f);
256                 fclose(f);
257         }
258
259         lpp_set_time_limit(ienv->lp, time_limit);
260         if (solve_log)
261                 lpp_set_log(ienv->lp, stdout);
262
263         lpp_solve(ienv->lp, options->ilp_server, options->ilp_solver);
264
265         //be_stat_ev_dbl("co_ilp_objval",     ienv->lp->objval);
266         //be_stat_ev_dbl("co_ilp_best_bound", ienv->lp->best_bound);
267         be_stat_ev    ("co_ilp_iter",       lpp_get_iter_cnt(ienv->lp));
268         be_stat_ev_dbl("co_ilp_sol_time",   lpp_get_sol_time(ienv->lp));
269
270         ienv->apply(ienv);
271
272         sr_reinsert(ienv->sr);
273
274         return lpp_get_sol_state(ienv->lp);
275 }
276
277 void free_ilp_env(ilp_env_t *ienv)
278 {
279         free_size_red(ienv->sr);
280         lpp_free(ienv->lp);
281         free(ienv);
282 }