0
|
1 /* Graph representation and manipulation functions.
|
|
2 Copyright (C) 2007
|
|
3 Free Software Foundation, Inc.
|
|
4
|
|
5 This file is part of GCC.
|
|
6
|
|
7 GCC is free software; you can redistribute it and/or modify it under
|
|
8 the terms of the GNU General Public License as published by the Free
|
|
9 Software Foundation; either version 3, or (at your option) any later
|
|
10 version.
|
|
11
|
|
12 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
|
|
13 WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
15 for more details.
|
|
16
|
|
17 You should have received a copy of the GNU General Public License
|
|
18 along with GCC; see the file COPYING3. If not see
|
|
19 <http://www.gnu.org/licenses/>. */
|
|
20
|
|
21 #include "config.h"
|
|
22 #include "system.h"
|
|
23 #include "coretypes.h"
|
|
24 #include "obstack.h"
|
|
25 #include "bitmap.h"
|
|
26 #include "vec.h"
|
|
27 #include "vecprim.h"
|
|
28 #include "graphds.h"
|
|
29
|
|
30 /* Dumps graph G into F. */
|
|
31
|
|
32 void
|
|
33 dump_graph (FILE *f, struct graph *g)
|
|
34 {
|
|
35 int i;
|
|
36 struct graph_edge *e;
|
|
37
|
|
38 for (i = 0; i < g->n_vertices; i++)
|
|
39 {
|
|
40 if (!g->vertices[i].pred
|
|
41 && !g->vertices[i].succ)
|
|
42 continue;
|
|
43
|
|
44 fprintf (f, "%d (%d)\t<-", i, g->vertices[i].component);
|
|
45 for (e = g->vertices[i].pred; e; e = e->pred_next)
|
|
46 fprintf (f, " %d", e->src);
|
|
47 fprintf (f, "\n");
|
|
48
|
|
49 fprintf (f, "\t->");
|
|
50 for (e = g->vertices[i].succ; e; e = e->succ_next)
|
|
51 fprintf (f, " %d", e->dest);
|
|
52 fprintf (f, "\n");
|
|
53 }
|
|
54 }
|
|
55
|
|
56 /* Creates a new graph with N_VERTICES vertices. */
|
|
57
|
|
58 struct graph *
|
|
59 new_graph (int n_vertices)
|
|
60 {
|
|
61 struct graph *g = XNEW (struct graph);
|
|
62
|
|
63 g->n_vertices = n_vertices;
|
|
64 g->vertices = XCNEWVEC (struct vertex, n_vertices);
|
|
65
|
|
66 return g;
|
|
67 }
|
|
68
|
|
69 /* Adds an edge from F to T to graph G. The new edge is returned. */
|
|
70
|
|
71 struct graph_edge *
|
|
72 add_edge (struct graph *g, int f, int t)
|
|
73 {
|
|
74 struct graph_edge *e = XNEW (struct graph_edge);
|
|
75 struct vertex *vf = &g->vertices[f], *vt = &g->vertices[t];
|
|
76
|
|
77
|
|
78 e->src = f;
|
|
79 e->dest = t;
|
|
80
|
|
81 e->pred_next = vt->pred;
|
|
82 vt->pred = e;
|
|
83
|
|
84 e->succ_next = vf->succ;
|
|
85 vf->succ = e;
|
|
86
|
|
87 return e;
|
|
88 }
|
|
89
|
|
90 /* Moves all the edges incident with U to V. */
|
|
91
|
|
92 void
|
|
93 identify_vertices (struct graph *g, int v, int u)
|
|
94 {
|
|
95 struct vertex *vv = &g->vertices[v];
|
|
96 struct vertex *uu = &g->vertices[u];
|
|
97 struct graph_edge *e, *next;
|
|
98
|
|
99 for (e = uu->succ; e; e = next)
|
|
100 {
|
|
101 next = e->succ_next;
|
|
102
|
|
103 e->src = v;
|
|
104 e->succ_next = vv->succ;
|
|
105 vv->succ = e;
|
|
106 }
|
|
107 uu->succ = NULL;
|
|
108
|
|
109 for (e = uu->pred; e; e = next)
|
|
110 {
|
|
111 next = e->pred_next;
|
|
112
|
|
113 e->dest = v;
|
|
114 e->pred_next = vv->pred;
|
|
115 vv->pred = e;
|
|
116 }
|
|
117 uu->pred = NULL;
|
|
118 }
|
|
119
|
|
120 /* Helper function for graphds_dfs. Returns the source vertex of E, in the
|
|
121 direction given by FORWARD. */
|
|
122
|
|
123 static inline int
|
|
124 dfs_edge_src (struct graph_edge *e, bool forward)
|
|
125 {
|
|
126 return forward ? e->src : e->dest;
|
|
127 }
|
|
128
|
|
129 /* Helper function for graphds_dfs. Returns the destination vertex of E, in
|
|
130 the direction given by FORWARD. */
|
|
131
|
|
132 static inline int
|
|
133 dfs_edge_dest (struct graph_edge *e, bool forward)
|
|
134 {
|
|
135 return forward ? e->dest : e->src;
|
|
136 }
|
|
137
|
|
138 /* Helper function for graphds_dfs. Returns the first edge after E (including
|
|
139 E), in the graph direction given by FORWARD, that belongs to SUBGRAPH. */
|
|
140
|
|
141 static inline struct graph_edge *
|
|
142 foll_in_subgraph (struct graph_edge *e, bool forward, bitmap subgraph)
|
|
143 {
|
|
144 int d;
|
|
145
|
|
146 if (!subgraph)
|
|
147 return e;
|
|
148
|
|
149 while (e)
|
|
150 {
|
|
151 d = dfs_edge_dest (e, forward);
|
|
152 if (bitmap_bit_p (subgraph, d))
|
|
153 return e;
|
|
154
|
|
155 e = forward ? e->succ_next : e->pred_next;
|
|
156 }
|
|
157
|
|
158 return e;
|
|
159 }
|
|
160
|
|
161 /* Helper function for graphds_dfs. Select the first edge from V in G, in the
|
|
162 direction given by FORWARD, that belongs to SUBGRAPH. */
|
|
163
|
|
164 static inline struct graph_edge *
|
|
165 dfs_fst_edge (struct graph *g, int v, bool forward, bitmap subgraph)
|
|
166 {
|
|
167 struct graph_edge *e;
|
|
168
|
|
169 e = (forward ? g->vertices[v].succ : g->vertices[v].pred);
|
|
170 return foll_in_subgraph (e, forward, subgraph);
|
|
171 }
|
|
172
|
|
173 /* Helper function for graphds_dfs. Returns the next edge after E, in the
|
|
174 graph direction given by FORWARD, that belongs to SUBGRAPH. */
|
|
175
|
|
176 static inline struct graph_edge *
|
|
177 dfs_next_edge (struct graph_edge *e, bool forward, bitmap subgraph)
|
|
178 {
|
|
179 return foll_in_subgraph (forward ? e->succ_next : e->pred_next,
|
|
180 forward, subgraph);
|
|
181 }
|
|
182
|
|
183 /* Runs dfs search over vertices of G, from NQ vertices in queue QS.
|
|
184 The vertices in postorder are stored into QT. If FORWARD is false,
|
|
185 backward dfs is run. If SUBGRAPH is not NULL, it specifies the
|
|
186 subgraph of G to run DFS on. Returns the number of the components
|
|
187 of the graph (number of the restarts of DFS). */
|
|
188
|
|
189 int
|
|
190 graphds_dfs (struct graph *g, int *qs, int nq, VEC (int, heap) **qt,
|
|
191 bool forward, bitmap subgraph)
|
|
192 {
|
|
193 int i, tick = 0, v, comp = 0, top;
|
|
194 struct graph_edge *e;
|
|
195 struct graph_edge **stack = XNEWVEC (struct graph_edge *, g->n_vertices);
|
|
196 bitmap_iterator bi;
|
|
197 unsigned av;
|
|
198
|
|
199 if (subgraph)
|
|
200 {
|
|
201 EXECUTE_IF_SET_IN_BITMAP (subgraph, 0, av, bi)
|
|
202 {
|
|
203 g->vertices[av].component = -1;
|
|
204 g->vertices[av].post = -1;
|
|
205 }
|
|
206 }
|
|
207 else
|
|
208 {
|
|
209 for (i = 0; i < g->n_vertices; i++)
|
|
210 {
|
|
211 g->vertices[i].component = -1;
|
|
212 g->vertices[i].post = -1;
|
|
213 }
|
|
214 }
|
|
215
|
|
216 for (i = 0; i < nq; i++)
|
|
217 {
|
|
218 v = qs[i];
|
|
219 if (g->vertices[v].post != -1)
|
|
220 continue;
|
|
221
|
|
222 g->vertices[v].component = comp++;
|
|
223 e = dfs_fst_edge (g, v, forward, subgraph);
|
|
224 top = 0;
|
|
225
|
|
226 while (1)
|
|
227 {
|
|
228 while (e)
|
|
229 {
|
|
230 if (g->vertices[dfs_edge_dest (e, forward)].component
|
|
231 == -1)
|
|
232 break;
|
|
233 e = dfs_next_edge (e, forward, subgraph);
|
|
234 }
|
|
235
|
|
236 if (!e)
|
|
237 {
|
|
238 if (qt)
|
|
239 VEC_safe_push (int, heap, *qt, v);
|
|
240 g->vertices[v].post = tick++;
|
|
241
|
|
242 if (!top)
|
|
243 break;
|
|
244
|
|
245 e = stack[--top];
|
|
246 v = dfs_edge_src (e, forward);
|
|
247 e = dfs_next_edge (e, forward, subgraph);
|
|
248 continue;
|
|
249 }
|
|
250
|
|
251 stack[top++] = e;
|
|
252 v = dfs_edge_dest (e, forward);
|
|
253 e = dfs_fst_edge (g, v, forward, subgraph);
|
|
254 g->vertices[v].component = comp - 1;
|
|
255 }
|
|
256 }
|
|
257
|
|
258 free (stack);
|
|
259
|
|
260 return comp;
|
|
261 }
|
|
262
|
|
263 /* Determines the strongly connected components of G, using the algorithm of
|
|
264 Tarjan -- first determine the postorder dfs numbering in reversed graph,
|
|
265 then run the dfs on the original graph in the order given by decreasing
|
|
266 numbers assigned by the previous pass. If SUBGRAPH is not NULL, it
|
|
267 specifies the subgraph of G whose strongly connected components we want
|
|
268 to determine.
|
|
269
|
|
270 After running this function, v->component is the number of the strongly
|
|
271 connected component for each vertex of G. Returns the number of the
|
|
272 sccs of G. */
|
|
273
|
|
274 int
|
|
275 graphds_scc (struct graph *g, bitmap subgraph)
|
|
276 {
|
|
277 int *queue = XNEWVEC (int, g->n_vertices);
|
|
278 VEC (int, heap) *postorder = NULL;
|
|
279 int nq, i, comp;
|
|
280 unsigned v;
|
|
281 bitmap_iterator bi;
|
|
282
|
|
283 if (subgraph)
|
|
284 {
|
|
285 nq = 0;
|
|
286 EXECUTE_IF_SET_IN_BITMAP (subgraph, 0, v, bi)
|
|
287 {
|
|
288 queue[nq++] = v;
|
|
289 }
|
|
290 }
|
|
291 else
|
|
292 {
|
|
293 for (i = 0; i < g->n_vertices; i++)
|
|
294 queue[i] = i;
|
|
295 nq = g->n_vertices;
|
|
296 }
|
|
297
|
|
298 graphds_dfs (g, queue, nq, &postorder, false, subgraph);
|
|
299 gcc_assert (VEC_length (int, postorder) == (unsigned) nq);
|
|
300
|
|
301 for (i = 0; i < nq; i++)
|
|
302 queue[i] = VEC_index (int, postorder, nq - i - 1);
|
|
303 comp = graphds_dfs (g, queue, nq, NULL, true, subgraph);
|
|
304
|
|
305 free (queue);
|
|
306 VEC_free (int, heap, postorder);
|
|
307
|
|
308 return comp;
|
|
309 }
|
|
310
|
|
311 /* Runs CALLBACK for all edges in G. */
|
|
312
|
|
313 void
|
|
314 for_each_edge (struct graph *g, graphds_edge_callback callback)
|
|
315 {
|
|
316 struct graph_edge *e;
|
|
317 int i;
|
|
318
|
|
319 for (i = 0; i < g->n_vertices; i++)
|
|
320 for (e = g->vertices[i].succ; e; e = e->succ_next)
|
|
321 callback (g, e);
|
|
322 }
|
|
323
|
|
324 /* Releases the memory occupied by G. */
|
|
325
|
|
326 void
|
|
327 free_graph (struct graph *g)
|
|
328 {
|
|
329 struct graph_edge *e, *n;
|
|
330 struct vertex *v;
|
|
331 int i;
|
|
332
|
|
333 for (i = 0; i < g->n_vertices; i++)
|
|
334 {
|
|
335 v = &g->vertices[i];
|
|
336 for (e = v->succ; e; e = n)
|
|
337 {
|
|
338 n = e->succ_next;
|
|
339 free (e);
|
|
340 }
|
|
341 }
|
|
342 free (g->vertices);
|
|
343 free (g);
|
|
344 }
|
|
345
|
|
346 /* Returns the nearest common ancestor of X and Y in tree whose parent
|
|
347 links are given by PARENT. MARKS is the array used to mark the
|
|
348 vertices of the tree, and MARK is the number currently used as a mark. */
|
|
349
|
|
350 static int
|
|
351 tree_nca (int x, int y, int *parent, int *marks, int mark)
|
|
352 {
|
|
353 if (x == -1 || x == y)
|
|
354 return y;
|
|
355
|
|
356 /* We climb with X and Y up the tree, marking the visited nodes. When
|
|
357 we first arrive to a marked node, it is the common ancestor. */
|
|
358 marks[x] = mark;
|
|
359 marks[y] = mark;
|
|
360
|
|
361 while (1)
|
|
362 {
|
|
363 x = parent[x];
|
|
364 if (x == -1)
|
|
365 break;
|
|
366 if (marks[x] == mark)
|
|
367 return x;
|
|
368 marks[x] = mark;
|
|
369
|
|
370 y = parent[y];
|
|
371 if (y == -1)
|
|
372 break;
|
|
373 if (marks[y] == mark)
|
|
374 return y;
|
|
375 marks[y] = mark;
|
|
376 }
|
|
377
|
|
378 /* If we reached the root with one of the vertices, continue
|
|
379 with the other one till we reach the marked part of the
|
|
380 tree. */
|
|
381 if (x == -1)
|
|
382 {
|
|
383 for (y = parent[y]; marks[y] != mark; y = parent[y])
|
|
384 continue;
|
|
385
|
|
386 return y;
|
|
387 }
|
|
388 else
|
|
389 {
|
|
390 for (x = parent[x]; marks[x] != mark; x = parent[x])
|
|
391 continue;
|
|
392
|
|
393 return x;
|
|
394 }
|
|
395 }
|
|
396
|
|
397 /* Determines the dominance tree of G (stored in the PARENT, SON and BROTHER
|
|
398 arrays), where the entry node is ENTRY. */
|
|
399
|
|
400 void
|
|
401 graphds_domtree (struct graph *g, int entry,
|
|
402 int *parent, int *son, int *brother)
|
|
403 {
|
|
404 VEC (int, heap) *postorder = NULL;
|
|
405 int *marks = XCNEWVEC (int, g->n_vertices);
|
|
406 int mark = 1, i, v, idom;
|
|
407 bool changed = true;
|
|
408 struct graph_edge *e;
|
|
409
|
|
410 /* We use a slight modification of the standard iterative algorithm, as
|
|
411 described in
|
|
412
|
|
413 K. D. Cooper, T. J. Harvey and K. Kennedy: A Simple, Fast Dominance
|
|
414 Algorithm
|
|
415
|
|
416 sort vertices in reverse postorder
|
|
417 foreach v
|
|
418 dom(v) = everything
|
|
419 dom(entry) = entry;
|
|
420
|
|
421 while (anything changes)
|
|
422 foreach v
|
|
423 dom(v) = {v} union (intersection of dom(p) over all predecessors of v)
|
|
424
|
|
425 The sets dom(v) are represented by the parent links in the current version
|
|
426 of the dominance tree. */
|
|
427
|
|
428 for (i = 0; i < g->n_vertices; i++)
|
|
429 {
|
|
430 parent[i] = -1;
|
|
431 son[i] = -1;
|
|
432 brother[i] = -1;
|
|
433 }
|
|
434 graphds_dfs (g, &entry, 1, &postorder, true, NULL);
|
|
435 gcc_assert (VEC_length (int, postorder) == (unsigned) g->n_vertices);
|
|
436 gcc_assert (VEC_index (int, postorder, g->n_vertices - 1) == entry);
|
|
437
|
|
438 while (changed)
|
|
439 {
|
|
440 changed = false;
|
|
441
|
|
442 for (i = g->n_vertices - 2; i >= 0; i--)
|
|
443 {
|
|
444 v = VEC_index (int, postorder, i);
|
|
445 idom = -1;
|
|
446 for (e = g->vertices[v].pred; e; e = e->pred_next)
|
|
447 {
|
|
448 if (e->src != entry
|
|
449 && parent[e->src] == -1)
|
|
450 continue;
|
|
451
|
|
452 idom = tree_nca (idom, e->src, parent, marks, mark++);
|
|
453 }
|
|
454
|
|
455 if (idom != parent[v])
|
|
456 {
|
|
457 parent[v] = idom;
|
|
458 changed = true;
|
|
459 }
|
|
460 }
|
|
461 }
|
|
462
|
|
463 free (marks);
|
|
464 VEC_free (int, heap, postorder);
|
|
465
|
|
466 for (i = 0; i < g->n_vertices; i++)
|
|
467 if (parent[i] != -1)
|
|
468 {
|
|
469 brother[i] = son[parent[i]];
|
|
470 son[parent[i]] = i;
|
|
471 }
|
|
472 }
|