1*e4b17023SJohn Marino /* ET-trees data structure implementation.
2*e4b17023SJohn Marino Contributed by Pavel Nejedly
3*e4b17023SJohn Marino Copyright (C) 2002, 2003, 2004, 2005, 2007, 2008, 2010 Free Software
4*e4b17023SJohn Marino Foundation, Inc.
5*e4b17023SJohn Marino
6*e4b17023SJohn Marino This file is part of the libiberty library.
7*e4b17023SJohn Marino Libiberty is free software; you can redistribute it and/or
8*e4b17023SJohn Marino modify it under the terms of the GNU Library General Public
9*e4b17023SJohn Marino License as published by the Free Software Foundation; either
10*e4b17023SJohn Marino version 3 of the License, or (at your option) any later version.
11*e4b17023SJohn Marino
12*e4b17023SJohn Marino Libiberty is distributed in the hope that it will be useful,
13*e4b17023SJohn Marino but WITHOUT ANY WARRANTY; without even the implied warranty of
14*e4b17023SJohn Marino MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15*e4b17023SJohn Marino Library General Public License for more details.
16*e4b17023SJohn Marino
17*e4b17023SJohn Marino You should have received a copy of the GNU Library General Public
18*e4b17023SJohn Marino License along with libiberty; see the file COPYING3. If not see
19*e4b17023SJohn Marino <http://www.gnu.org/licenses/>.
20*e4b17023SJohn Marino
21*e4b17023SJohn Marino The ET-forest structure is described in:
22*e4b17023SJohn Marino D. D. Sleator and R. E. Tarjan. A data structure for dynamic trees.
23*e4b17023SJohn Marino J. G'omput. System Sci., 26(3):362 381, 1983.
24*e4b17023SJohn Marino */
25*e4b17023SJohn Marino
26*e4b17023SJohn Marino #include "config.h"
27*e4b17023SJohn Marino #include "system.h"
28*e4b17023SJohn Marino #include "coretypes.h"
29*e4b17023SJohn Marino #include "et-forest.h"
30*e4b17023SJohn Marino #include "alloc-pool.h"
31*e4b17023SJohn Marino
32*e4b17023SJohn Marino /* We do not enable this with ENABLE_CHECKING, since it is awfully slow. */
33*e4b17023SJohn Marino #undef DEBUG_ET
34*e4b17023SJohn Marino
35*e4b17023SJohn Marino #ifdef DEBUG_ET
36*e4b17023SJohn Marino #include "basic-block.h"
37*e4b17023SJohn Marino #endif
38*e4b17023SJohn Marino
39*e4b17023SJohn Marino /* The occurrence of a node in the et tree. */
40*e4b17023SJohn Marino struct et_occ
41*e4b17023SJohn Marino {
42*e4b17023SJohn Marino struct et_node *of; /* The node. */
43*e4b17023SJohn Marino
44*e4b17023SJohn Marino struct et_occ *parent; /* Parent in the splay-tree. */
45*e4b17023SJohn Marino struct et_occ *prev; /* Left son in the splay-tree. */
46*e4b17023SJohn Marino struct et_occ *next; /* Right son in the splay-tree. */
47*e4b17023SJohn Marino
48*e4b17023SJohn Marino int depth; /* The depth of the node is the sum of depth
49*e4b17023SJohn Marino fields on the path to the root. */
50*e4b17023SJohn Marino int min; /* The minimum value of the depth in the subtree
51*e4b17023SJohn Marino is obtained by adding sum of depth fields
52*e4b17023SJohn Marino on the path to the root. */
53*e4b17023SJohn Marino struct et_occ *min_occ; /* The occurrence in the subtree with the minimal
54*e4b17023SJohn Marino depth. */
55*e4b17023SJohn Marino };
56*e4b17023SJohn Marino
57*e4b17023SJohn Marino static alloc_pool et_nodes;
58*e4b17023SJohn Marino static alloc_pool et_occurrences;
59*e4b17023SJohn Marino
60*e4b17023SJohn Marino /* Changes depth of OCC to D. */
61*e4b17023SJohn Marino
62*e4b17023SJohn Marino static inline void
set_depth(struct et_occ * occ,int d)63*e4b17023SJohn Marino set_depth (struct et_occ *occ, int d)
64*e4b17023SJohn Marino {
65*e4b17023SJohn Marino if (!occ)
66*e4b17023SJohn Marino return;
67*e4b17023SJohn Marino
68*e4b17023SJohn Marino occ->min += d - occ->depth;
69*e4b17023SJohn Marino occ->depth = d;
70*e4b17023SJohn Marino }
71*e4b17023SJohn Marino
72*e4b17023SJohn Marino /* Adds D to the depth of OCC. */
73*e4b17023SJohn Marino
74*e4b17023SJohn Marino static inline void
set_depth_add(struct et_occ * occ,int d)75*e4b17023SJohn Marino set_depth_add (struct et_occ *occ, int d)
76*e4b17023SJohn Marino {
77*e4b17023SJohn Marino if (!occ)
78*e4b17023SJohn Marino return;
79*e4b17023SJohn Marino
80*e4b17023SJohn Marino occ->min += d;
81*e4b17023SJohn Marino occ->depth += d;
82*e4b17023SJohn Marino }
83*e4b17023SJohn Marino
84*e4b17023SJohn Marino /* Sets prev field of OCC to P. */
85*e4b17023SJohn Marino
86*e4b17023SJohn Marino static inline void
set_prev(struct et_occ * occ,struct et_occ * t)87*e4b17023SJohn Marino set_prev (struct et_occ *occ, struct et_occ *t)
88*e4b17023SJohn Marino {
89*e4b17023SJohn Marino #ifdef DEBUG_ET
90*e4b17023SJohn Marino gcc_assert (occ != t);
91*e4b17023SJohn Marino #endif
92*e4b17023SJohn Marino
93*e4b17023SJohn Marino occ->prev = t;
94*e4b17023SJohn Marino if (t)
95*e4b17023SJohn Marino t->parent = occ;
96*e4b17023SJohn Marino }
97*e4b17023SJohn Marino
98*e4b17023SJohn Marino /* Sets next field of OCC to P. */
99*e4b17023SJohn Marino
100*e4b17023SJohn Marino static inline void
set_next(struct et_occ * occ,struct et_occ * t)101*e4b17023SJohn Marino set_next (struct et_occ *occ, struct et_occ *t)
102*e4b17023SJohn Marino {
103*e4b17023SJohn Marino #ifdef DEBUG_ET
104*e4b17023SJohn Marino gcc_assert (occ != t);
105*e4b17023SJohn Marino #endif
106*e4b17023SJohn Marino
107*e4b17023SJohn Marino occ->next = t;
108*e4b17023SJohn Marino if (t)
109*e4b17023SJohn Marino t->parent = occ;
110*e4b17023SJohn Marino }
111*e4b17023SJohn Marino
112*e4b17023SJohn Marino /* Recompute minimum for occurrence OCC. */
113*e4b17023SJohn Marino
114*e4b17023SJohn Marino static inline void
et_recomp_min(struct et_occ * occ)115*e4b17023SJohn Marino et_recomp_min (struct et_occ *occ)
116*e4b17023SJohn Marino {
117*e4b17023SJohn Marino struct et_occ *mson = occ->prev;
118*e4b17023SJohn Marino
119*e4b17023SJohn Marino if (!mson
120*e4b17023SJohn Marino || (occ->next
121*e4b17023SJohn Marino && mson->min > occ->next->min))
122*e4b17023SJohn Marino mson = occ->next;
123*e4b17023SJohn Marino
124*e4b17023SJohn Marino if (mson && mson->min < 0)
125*e4b17023SJohn Marino {
126*e4b17023SJohn Marino occ->min = mson->min + occ->depth;
127*e4b17023SJohn Marino occ->min_occ = mson->min_occ;
128*e4b17023SJohn Marino }
129*e4b17023SJohn Marino else
130*e4b17023SJohn Marino {
131*e4b17023SJohn Marino occ->min = occ->depth;
132*e4b17023SJohn Marino occ->min_occ = occ;
133*e4b17023SJohn Marino }
134*e4b17023SJohn Marino }
135*e4b17023SJohn Marino
136*e4b17023SJohn Marino #ifdef DEBUG_ET
137*e4b17023SJohn Marino /* Checks whether neighborhood of OCC seems sane. */
138*e4b17023SJohn Marino
139*e4b17023SJohn Marino static void
et_check_occ_sanity(struct et_occ * occ)140*e4b17023SJohn Marino et_check_occ_sanity (struct et_occ *occ)
141*e4b17023SJohn Marino {
142*e4b17023SJohn Marino if (!occ)
143*e4b17023SJohn Marino return;
144*e4b17023SJohn Marino
145*e4b17023SJohn Marino gcc_assert (occ->parent != occ);
146*e4b17023SJohn Marino gcc_assert (occ->prev != occ);
147*e4b17023SJohn Marino gcc_assert (occ->next != occ);
148*e4b17023SJohn Marino gcc_assert (!occ->next || occ->next != occ->prev);
149*e4b17023SJohn Marino
150*e4b17023SJohn Marino if (occ->next)
151*e4b17023SJohn Marino {
152*e4b17023SJohn Marino gcc_assert (occ->next != occ->parent);
153*e4b17023SJohn Marino gcc_assert (occ->next->parent == occ);
154*e4b17023SJohn Marino }
155*e4b17023SJohn Marino
156*e4b17023SJohn Marino if (occ->prev)
157*e4b17023SJohn Marino {
158*e4b17023SJohn Marino gcc_assert (occ->prev != occ->parent);
159*e4b17023SJohn Marino gcc_assert (occ->prev->parent == occ);
160*e4b17023SJohn Marino }
161*e4b17023SJohn Marino
162*e4b17023SJohn Marino gcc_assert (!occ->parent
163*e4b17023SJohn Marino || occ->parent->prev == occ
164*e4b17023SJohn Marino || occ->parent->next == occ);
165*e4b17023SJohn Marino }
166*e4b17023SJohn Marino
167*e4b17023SJohn Marino /* Checks whether tree rooted at OCC is sane. */
168*e4b17023SJohn Marino
169*e4b17023SJohn Marino static void
et_check_sanity(struct et_occ * occ)170*e4b17023SJohn Marino et_check_sanity (struct et_occ *occ)
171*e4b17023SJohn Marino {
172*e4b17023SJohn Marino et_check_occ_sanity (occ);
173*e4b17023SJohn Marino if (occ->prev)
174*e4b17023SJohn Marino et_check_sanity (occ->prev);
175*e4b17023SJohn Marino if (occ->next)
176*e4b17023SJohn Marino et_check_sanity (occ->next);
177*e4b17023SJohn Marino }
178*e4b17023SJohn Marino
179*e4b17023SJohn Marino /* Checks whether tree containing OCC is sane. */
180*e4b17023SJohn Marino
181*e4b17023SJohn Marino static void
et_check_tree_sanity(struct et_occ * occ)182*e4b17023SJohn Marino et_check_tree_sanity (struct et_occ *occ)
183*e4b17023SJohn Marino {
184*e4b17023SJohn Marino while (occ->parent)
185*e4b17023SJohn Marino occ = occ->parent;
186*e4b17023SJohn Marino
187*e4b17023SJohn Marino et_check_sanity (occ);
188*e4b17023SJohn Marino }
189*e4b17023SJohn Marino
190*e4b17023SJohn Marino /* For recording the paths. */
191*e4b17023SJohn Marino
192*e4b17023SJohn Marino /* An ad-hoc constant; if the function has more blocks, this won't work,
193*e4b17023SJohn Marino but since it is used for debugging only, it does not matter. */
194*e4b17023SJohn Marino #define MAX_NODES 100000
195*e4b17023SJohn Marino
196*e4b17023SJohn Marino static int len;
197*e4b17023SJohn Marino static void *datas[MAX_NODES];
198*e4b17023SJohn Marino static int depths[MAX_NODES];
199*e4b17023SJohn Marino
200*e4b17023SJohn Marino /* Records the path represented by OCC, with depth incremented by DEPTH. */
201*e4b17023SJohn Marino
202*e4b17023SJohn Marino static int
record_path_before_1(struct et_occ * occ,int depth)203*e4b17023SJohn Marino record_path_before_1 (struct et_occ *occ, int depth)
204*e4b17023SJohn Marino {
205*e4b17023SJohn Marino int mn, m;
206*e4b17023SJohn Marino
207*e4b17023SJohn Marino depth += occ->depth;
208*e4b17023SJohn Marino mn = depth;
209*e4b17023SJohn Marino
210*e4b17023SJohn Marino if (occ->prev)
211*e4b17023SJohn Marino {
212*e4b17023SJohn Marino m = record_path_before_1 (occ->prev, depth);
213*e4b17023SJohn Marino if (m < mn)
214*e4b17023SJohn Marino mn = m;
215*e4b17023SJohn Marino }
216*e4b17023SJohn Marino
217*e4b17023SJohn Marino fprintf (stderr, "%d (%d); ", ((basic_block) occ->of->data)->index, depth);
218*e4b17023SJohn Marino
219*e4b17023SJohn Marino gcc_assert (len < MAX_NODES);
220*e4b17023SJohn Marino
221*e4b17023SJohn Marino depths[len] = depth;
222*e4b17023SJohn Marino datas[len] = occ->of;
223*e4b17023SJohn Marino len++;
224*e4b17023SJohn Marino
225*e4b17023SJohn Marino if (occ->next)
226*e4b17023SJohn Marino {
227*e4b17023SJohn Marino m = record_path_before_1 (occ->next, depth);
228*e4b17023SJohn Marino if (m < mn)
229*e4b17023SJohn Marino mn = m;
230*e4b17023SJohn Marino }
231*e4b17023SJohn Marino
232*e4b17023SJohn Marino gcc_assert (mn == occ->min + depth - occ->depth);
233*e4b17023SJohn Marino
234*e4b17023SJohn Marino return mn;
235*e4b17023SJohn Marino }
236*e4b17023SJohn Marino
237*e4b17023SJohn Marino /* Records the path represented by a tree containing OCC. */
238*e4b17023SJohn Marino
239*e4b17023SJohn Marino static void
record_path_before(struct et_occ * occ)240*e4b17023SJohn Marino record_path_before (struct et_occ *occ)
241*e4b17023SJohn Marino {
242*e4b17023SJohn Marino while (occ->parent)
243*e4b17023SJohn Marino occ = occ->parent;
244*e4b17023SJohn Marino
245*e4b17023SJohn Marino len = 0;
246*e4b17023SJohn Marino record_path_before_1 (occ, 0);
247*e4b17023SJohn Marino fprintf (stderr, "\n");
248*e4b17023SJohn Marino }
249*e4b17023SJohn Marino
250*e4b17023SJohn Marino /* Checks whether the path represented by OCC, with depth incremented by DEPTH,
251*e4b17023SJohn Marino was not changed since the last recording. */
252*e4b17023SJohn Marino
253*e4b17023SJohn Marino static int
check_path_after_1(struct et_occ * occ,int depth)254*e4b17023SJohn Marino check_path_after_1 (struct et_occ *occ, int depth)
255*e4b17023SJohn Marino {
256*e4b17023SJohn Marino int mn, m;
257*e4b17023SJohn Marino
258*e4b17023SJohn Marino depth += occ->depth;
259*e4b17023SJohn Marino mn = depth;
260*e4b17023SJohn Marino
261*e4b17023SJohn Marino if (occ->next)
262*e4b17023SJohn Marino {
263*e4b17023SJohn Marino m = check_path_after_1 (occ->next, depth);
264*e4b17023SJohn Marino if (m < mn)
265*e4b17023SJohn Marino mn = m;
266*e4b17023SJohn Marino }
267*e4b17023SJohn Marino
268*e4b17023SJohn Marino len--;
269*e4b17023SJohn Marino gcc_assert (depths[len] == depth && datas[len] == occ->of);
270*e4b17023SJohn Marino
271*e4b17023SJohn Marino if (occ->prev)
272*e4b17023SJohn Marino {
273*e4b17023SJohn Marino m = check_path_after_1 (occ->prev, depth);
274*e4b17023SJohn Marino if (m < mn)
275*e4b17023SJohn Marino mn = m;
276*e4b17023SJohn Marino }
277*e4b17023SJohn Marino
278*e4b17023SJohn Marino gcc_assert (mn == occ->min + depth - occ->depth);
279*e4b17023SJohn Marino
280*e4b17023SJohn Marino return mn;
281*e4b17023SJohn Marino }
282*e4b17023SJohn Marino
283*e4b17023SJohn Marino /* Checks whether the path represented by a tree containing OCC was
284*e4b17023SJohn Marino not changed since the last recording. */
285*e4b17023SJohn Marino
286*e4b17023SJohn Marino static void
check_path_after(struct et_occ * occ)287*e4b17023SJohn Marino check_path_after (struct et_occ *occ)
288*e4b17023SJohn Marino {
289*e4b17023SJohn Marino while (occ->parent)
290*e4b17023SJohn Marino occ = occ->parent;
291*e4b17023SJohn Marino
292*e4b17023SJohn Marino check_path_after_1 (occ, 0);
293*e4b17023SJohn Marino gcc_assert (!len);
294*e4b17023SJohn Marino }
295*e4b17023SJohn Marino
296*e4b17023SJohn Marino #endif
297*e4b17023SJohn Marino
298*e4b17023SJohn Marino /* Splay the occurrence OCC to the root of the tree. */
299*e4b17023SJohn Marino
300*e4b17023SJohn Marino static void
et_splay(struct et_occ * occ)301*e4b17023SJohn Marino et_splay (struct et_occ *occ)
302*e4b17023SJohn Marino {
303*e4b17023SJohn Marino struct et_occ *f, *gf, *ggf;
304*e4b17023SJohn Marino int occ_depth, f_depth, gf_depth;
305*e4b17023SJohn Marino
306*e4b17023SJohn Marino #ifdef DEBUG_ET
307*e4b17023SJohn Marino record_path_before (occ);
308*e4b17023SJohn Marino et_check_tree_sanity (occ);
309*e4b17023SJohn Marino #endif
310*e4b17023SJohn Marino
311*e4b17023SJohn Marino while (occ->parent)
312*e4b17023SJohn Marino {
313*e4b17023SJohn Marino occ_depth = occ->depth;
314*e4b17023SJohn Marino
315*e4b17023SJohn Marino f = occ->parent;
316*e4b17023SJohn Marino f_depth = f->depth;
317*e4b17023SJohn Marino
318*e4b17023SJohn Marino gf = f->parent;
319*e4b17023SJohn Marino
320*e4b17023SJohn Marino if (!gf)
321*e4b17023SJohn Marino {
322*e4b17023SJohn Marino set_depth_add (occ, f_depth);
323*e4b17023SJohn Marino occ->min_occ = f->min_occ;
324*e4b17023SJohn Marino occ->min = f->min;
325*e4b17023SJohn Marino
326*e4b17023SJohn Marino if (f->prev == occ)
327*e4b17023SJohn Marino {
328*e4b17023SJohn Marino /* zig */
329*e4b17023SJohn Marino set_prev (f, occ->next);
330*e4b17023SJohn Marino set_next (occ, f);
331*e4b17023SJohn Marino set_depth_add (f->prev, occ_depth);
332*e4b17023SJohn Marino }
333*e4b17023SJohn Marino else
334*e4b17023SJohn Marino {
335*e4b17023SJohn Marino /* zag */
336*e4b17023SJohn Marino set_next (f, occ->prev);
337*e4b17023SJohn Marino set_prev (occ, f);
338*e4b17023SJohn Marino set_depth_add (f->next, occ_depth);
339*e4b17023SJohn Marino }
340*e4b17023SJohn Marino set_depth (f, -occ_depth);
341*e4b17023SJohn Marino occ->parent = NULL;
342*e4b17023SJohn Marino
343*e4b17023SJohn Marino et_recomp_min (f);
344*e4b17023SJohn Marino #ifdef DEBUG_ET
345*e4b17023SJohn Marino et_check_tree_sanity (occ);
346*e4b17023SJohn Marino check_path_after (occ);
347*e4b17023SJohn Marino #endif
348*e4b17023SJohn Marino return;
349*e4b17023SJohn Marino }
350*e4b17023SJohn Marino
351*e4b17023SJohn Marino gf_depth = gf->depth;
352*e4b17023SJohn Marino
353*e4b17023SJohn Marino set_depth_add (occ, f_depth + gf_depth);
354*e4b17023SJohn Marino occ->min_occ = gf->min_occ;
355*e4b17023SJohn Marino occ->min = gf->min;
356*e4b17023SJohn Marino
357*e4b17023SJohn Marino ggf = gf->parent;
358*e4b17023SJohn Marino
359*e4b17023SJohn Marino if (gf->prev == f)
360*e4b17023SJohn Marino {
361*e4b17023SJohn Marino if (f->prev == occ)
362*e4b17023SJohn Marino {
363*e4b17023SJohn Marino /* zig zig */
364*e4b17023SJohn Marino set_prev (gf, f->next);
365*e4b17023SJohn Marino set_prev (f, occ->next);
366*e4b17023SJohn Marino set_next (occ, f);
367*e4b17023SJohn Marino set_next (f, gf);
368*e4b17023SJohn Marino
369*e4b17023SJohn Marino set_depth (f, -occ_depth);
370*e4b17023SJohn Marino set_depth_add (f->prev, occ_depth);
371*e4b17023SJohn Marino set_depth (gf, -f_depth);
372*e4b17023SJohn Marino set_depth_add (gf->prev, f_depth);
373*e4b17023SJohn Marino }
374*e4b17023SJohn Marino else
375*e4b17023SJohn Marino {
376*e4b17023SJohn Marino /* zag zig */
377*e4b17023SJohn Marino set_prev (gf, occ->next);
378*e4b17023SJohn Marino set_next (f, occ->prev);
379*e4b17023SJohn Marino set_prev (occ, f);
380*e4b17023SJohn Marino set_next (occ, gf);
381*e4b17023SJohn Marino
382*e4b17023SJohn Marino set_depth (f, -occ_depth);
383*e4b17023SJohn Marino set_depth_add (f->next, occ_depth);
384*e4b17023SJohn Marino set_depth (gf, -occ_depth - f_depth);
385*e4b17023SJohn Marino set_depth_add (gf->prev, occ_depth + f_depth);
386*e4b17023SJohn Marino }
387*e4b17023SJohn Marino }
388*e4b17023SJohn Marino else
389*e4b17023SJohn Marino {
390*e4b17023SJohn Marino if (f->prev == occ)
391*e4b17023SJohn Marino {
392*e4b17023SJohn Marino /* zig zag */
393*e4b17023SJohn Marino set_next (gf, occ->prev);
394*e4b17023SJohn Marino set_prev (f, occ->next);
395*e4b17023SJohn Marino set_prev (occ, gf);
396*e4b17023SJohn Marino set_next (occ, f);
397*e4b17023SJohn Marino
398*e4b17023SJohn Marino set_depth (f, -occ_depth);
399*e4b17023SJohn Marino set_depth_add (f->prev, occ_depth);
400*e4b17023SJohn Marino set_depth (gf, -occ_depth - f_depth);
401*e4b17023SJohn Marino set_depth_add (gf->next, occ_depth + f_depth);
402*e4b17023SJohn Marino }
403*e4b17023SJohn Marino else
404*e4b17023SJohn Marino {
405*e4b17023SJohn Marino /* zag zag */
406*e4b17023SJohn Marino set_next (gf, f->prev);
407*e4b17023SJohn Marino set_next (f, occ->prev);
408*e4b17023SJohn Marino set_prev (occ, f);
409*e4b17023SJohn Marino set_prev (f, gf);
410*e4b17023SJohn Marino
411*e4b17023SJohn Marino set_depth (f, -occ_depth);
412*e4b17023SJohn Marino set_depth_add (f->next, occ_depth);
413*e4b17023SJohn Marino set_depth (gf, -f_depth);
414*e4b17023SJohn Marino set_depth_add (gf->next, f_depth);
415*e4b17023SJohn Marino }
416*e4b17023SJohn Marino }
417*e4b17023SJohn Marino
418*e4b17023SJohn Marino occ->parent = ggf;
419*e4b17023SJohn Marino if (ggf)
420*e4b17023SJohn Marino {
421*e4b17023SJohn Marino if (ggf->prev == gf)
422*e4b17023SJohn Marino ggf->prev = occ;
423*e4b17023SJohn Marino else
424*e4b17023SJohn Marino ggf->next = occ;
425*e4b17023SJohn Marino }
426*e4b17023SJohn Marino
427*e4b17023SJohn Marino et_recomp_min (gf);
428*e4b17023SJohn Marino et_recomp_min (f);
429*e4b17023SJohn Marino #ifdef DEBUG_ET
430*e4b17023SJohn Marino et_check_tree_sanity (occ);
431*e4b17023SJohn Marino #endif
432*e4b17023SJohn Marino }
433*e4b17023SJohn Marino
434*e4b17023SJohn Marino #ifdef DEBUG_ET
435*e4b17023SJohn Marino et_check_sanity (occ);
436*e4b17023SJohn Marino check_path_after (occ);
437*e4b17023SJohn Marino #endif
438*e4b17023SJohn Marino }
439*e4b17023SJohn Marino
440*e4b17023SJohn Marino /* Create a new et tree occurrence of NODE. */
441*e4b17023SJohn Marino
442*e4b17023SJohn Marino static struct et_occ *
et_new_occ(struct et_node * node)443*e4b17023SJohn Marino et_new_occ (struct et_node *node)
444*e4b17023SJohn Marino {
445*e4b17023SJohn Marino struct et_occ *nw;
446*e4b17023SJohn Marino
447*e4b17023SJohn Marino if (!et_occurrences)
448*e4b17023SJohn Marino et_occurrences = create_alloc_pool ("et_occ pool", sizeof (struct et_occ), 300);
449*e4b17023SJohn Marino nw = (struct et_occ *) pool_alloc (et_occurrences);
450*e4b17023SJohn Marino
451*e4b17023SJohn Marino nw->of = node;
452*e4b17023SJohn Marino nw->parent = NULL;
453*e4b17023SJohn Marino nw->prev = NULL;
454*e4b17023SJohn Marino nw->next = NULL;
455*e4b17023SJohn Marino
456*e4b17023SJohn Marino nw->depth = 0;
457*e4b17023SJohn Marino nw->min_occ = nw;
458*e4b17023SJohn Marino nw->min = 0;
459*e4b17023SJohn Marino
460*e4b17023SJohn Marino return nw;
461*e4b17023SJohn Marino }
462*e4b17023SJohn Marino
463*e4b17023SJohn Marino /* Create a new et tree containing DATA. */
464*e4b17023SJohn Marino
465*e4b17023SJohn Marino struct et_node *
et_new_tree(void * data)466*e4b17023SJohn Marino et_new_tree (void *data)
467*e4b17023SJohn Marino {
468*e4b17023SJohn Marino struct et_node *nw;
469*e4b17023SJohn Marino
470*e4b17023SJohn Marino if (!et_nodes)
471*e4b17023SJohn Marino et_nodes = create_alloc_pool ("et_node pool", sizeof (struct et_node), 300);
472*e4b17023SJohn Marino nw = (struct et_node *) pool_alloc (et_nodes);
473*e4b17023SJohn Marino
474*e4b17023SJohn Marino nw->data = data;
475*e4b17023SJohn Marino nw->father = NULL;
476*e4b17023SJohn Marino nw->left = NULL;
477*e4b17023SJohn Marino nw->right = NULL;
478*e4b17023SJohn Marino nw->son = NULL;
479*e4b17023SJohn Marino
480*e4b17023SJohn Marino nw->rightmost_occ = et_new_occ (nw);
481*e4b17023SJohn Marino nw->parent_occ = NULL;
482*e4b17023SJohn Marino
483*e4b17023SJohn Marino return nw;
484*e4b17023SJohn Marino }
485*e4b17023SJohn Marino
486*e4b17023SJohn Marino /* Releases et tree T. */
487*e4b17023SJohn Marino
488*e4b17023SJohn Marino void
et_free_tree(struct et_node * t)489*e4b17023SJohn Marino et_free_tree (struct et_node *t)
490*e4b17023SJohn Marino {
491*e4b17023SJohn Marino while (t->son)
492*e4b17023SJohn Marino et_split (t->son);
493*e4b17023SJohn Marino
494*e4b17023SJohn Marino if (t->father)
495*e4b17023SJohn Marino et_split (t);
496*e4b17023SJohn Marino
497*e4b17023SJohn Marino pool_free (et_occurrences, t->rightmost_occ);
498*e4b17023SJohn Marino pool_free (et_nodes, t);
499*e4b17023SJohn Marino }
500*e4b17023SJohn Marino
501*e4b17023SJohn Marino /* Releases et tree T without maintaining other nodes. */
502*e4b17023SJohn Marino
503*e4b17023SJohn Marino void
et_free_tree_force(struct et_node * t)504*e4b17023SJohn Marino et_free_tree_force (struct et_node *t)
505*e4b17023SJohn Marino {
506*e4b17023SJohn Marino pool_free (et_occurrences, t->rightmost_occ);
507*e4b17023SJohn Marino if (t->parent_occ)
508*e4b17023SJohn Marino pool_free (et_occurrences, t->parent_occ);
509*e4b17023SJohn Marino pool_free (et_nodes, t);
510*e4b17023SJohn Marino }
511*e4b17023SJohn Marino
512*e4b17023SJohn Marino /* Release the alloc pools, if they are empty. */
513*e4b17023SJohn Marino
514*e4b17023SJohn Marino void
et_free_pools(void)515*e4b17023SJohn Marino et_free_pools (void)
516*e4b17023SJohn Marino {
517*e4b17023SJohn Marino free_alloc_pool_if_empty (&et_occurrences);
518*e4b17023SJohn Marino free_alloc_pool_if_empty (&et_nodes);
519*e4b17023SJohn Marino }
520*e4b17023SJohn Marino
521*e4b17023SJohn Marino /* Sets father of et tree T to FATHER. */
522*e4b17023SJohn Marino
523*e4b17023SJohn Marino void
et_set_father(struct et_node * t,struct et_node * father)524*e4b17023SJohn Marino et_set_father (struct et_node *t, struct et_node *father)
525*e4b17023SJohn Marino {
526*e4b17023SJohn Marino struct et_node *left, *right;
527*e4b17023SJohn Marino struct et_occ *rmost, *left_part, *new_f_occ, *p;
528*e4b17023SJohn Marino
529*e4b17023SJohn Marino /* Update the path represented in the splay tree. */
530*e4b17023SJohn Marino new_f_occ = et_new_occ (father);
531*e4b17023SJohn Marino
532*e4b17023SJohn Marino rmost = father->rightmost_occ;
533*e4b17023SJohn Marino et_splay (rmost);
534*e4b17023SJohn Marino
535*e4b17023SJohn Marino left_part = rmost->prev;
536*e4b17023SJohn Marino
537*e4b17023SJohn Marino p = t->rightmost_occ;
538*e4b17023SJohn Marino et_splay (p);
539*e4b17023SJohn Marino
540*e4b17023SJohn Marino set_prev (new_f_occ, left_part);
541*e4b17023SJohn Marino set_next (new_f_occ, p);
542*e4b17023SJohn Marino
543*e4b17023SJohn Marino p->depth++;
544*e4b17023SJohn Marino p->min++;
545*e4b17023SJohn Marino et_recomp_min (new_f_occ);
546*e4b17023SJohn Marino
547*e4b17023SJohn Marino set_prev (rmost, new_f_occ);
548*e4b17023SJohn Marino
549*e4b17023SJohn Marino if (new_f_occ->min + rmost->depth < rmost->min)
550*e4b17023SJohn Marino {
551*e4b17023SJohn Marino rmost->min = new_f_occ->min + rmost->depth;
552*e4b17023SJohn Marino rmost->min_occ = new_f_occ->min_occ;
553*e4b17023SJohn Marino }
554*e4b17023SJohn Marino
555*e4b17023SJohn Marino t->parent_occ = new_f_occ;
556*e4b17023SJohn Marino
557*e4b17023SJohn Marino /* Update the tree. */
558*e4b17023SJohn Marino t->father = father;
559*e4b17023SJohn Marino right = father->son;
560*e4b17023SJohn Marino if (right)
561*e4b17023SJohn Marino left = right->left;
562*e4b17023SJohn Marino else
563*e4b17023SJohn Marino left = right = t;
564*e4b17023SJohn Marino
565*e4b17023SJohn Marino left->right = t;
566*e4b17023SJohn Marino right->left = t;
567*e4b17023SJohn Marino t->left = left;
568*e4b17023SJohn Marino t->right = right;
569*e4b17023SJohn Marino
570*e4b17023SJohn Marino father->son = t;
571*e4b17023SJohn Marino
572*e4b17023SJohn Marino #ifdef DEBUG_ET
573*e4b17023SJohn Marino et_check_tree_sanity (rmost);
574*e4b17023SJohn Marino record_path_before (rmost);
575*e4b17023SJohn Marino #endif
576*e4b17023SJohn Marino }
577*e4b17023SJohn Marino
578*e4b17023SJohn Marino /* Splits the edge from T to its father. */
579*e4b17023SJohn Marino
580*e4b17023SJohn Marino void
et_split(struct et_node * t)581*e4b17023SJohn Marino et_split (struct et_node *t)
582*e4b17023SJohn Marino {
583*e4b17023SJohn Marino struct et_node *father = t->father;
584*e4b17023SJohn Marino struct et_occ *r, *l, *rmost, *p_occ;
585*e4b17023SJohn Marino
586*e4b17023SJohn Marino /* Update the path represented by the splay tree. */
587*e4b17023SJohn Marino rmost = t->rightmost_occ;
588*e4b17023SJohn Marino et_splay (rmost);
589*e4b17023SJohn Marino
590*e4b17023SJohn Marino for (r = rmost->next; r->prev; r = r->prev)
591*e4b17023SJohn Marino continue;
592*e4b17023SJohn Marino et_splay (r);
593*e4b17023SJohn Marino
594*e4b17023SJohn Marino r->prev->parent = NULL;
595*e4b17023SJohn Marino p_occ = t->parent_occ;
596*e4b17023SJohn Marino et_splay (p_occ);
597*e4b17023SJohn Marino t->parent_occ = NULL;
598*e4b17023SJohn Marino
599*e4b17023SJohn Marino l = p_occ->prev;
600*e4b17023SJohn Marino p_occ->next->parent = NULL;
601*e4b17023SJohn Marino
602*e4b17023SJohn Marino set_prev (r, l);
603*e4b17023SJohn Marino
604*e4b17023SJohn Marino et_recomp_min (r);
605*e4b17023SJohn Marino
606*e4b17023SJohn Marino et_splay (rmost);
607*e4b17023SJohn Marino rmost->depth = 0;
608*e4b17023SJohn Marino rmost->min = 0;
609*e4b17023SJohn Marino
610*e4b17023SJohn Marino pool_free (et_occurrences, p_occ);
611*e4b17023SJohn Marino
612*e4b17023SJohn Marino /* Update the tree. */
613*e4b17023SJohn Marino if (father->son == t)
614*e4b17023SJohn Marino father->son = t->right;
615*e4b17023SJohn Marino if (father->son == t)
616*e4b17023SJohn Marino father->son = NULL;
617*e4b17023SJohn Marino else
618*e4b17023SJohn Marino {
619*e4b17023SJohn Marino t->left->right = t->right;
620*e4b17023SJohn Marino t->right->left = t->left;
621*e4b17023SJohn Marino }
622*e4b17023SJohn Marino t->left = t->right = NULL;
623*e4b17023SJohn Marino t->father = NULL;
624*e4b17023SJohn Marino
625*e4b17023SJohn Marino #ifdef DEBUG_ET
626*e4b17023SJohn Marino et_check_tree_sanity (rmost);
627*e4b17023SJohn Marino record_path_before (rmost);
628*e4b17023SJohn Marino
629*e4b17023SJohn Marino et_check_tree_sanity (r);
630*e4b17023SJohn Marino record_path_before (r);
631*e4b17023SJohn Marino #endif
632*e4b17023SJohn Marino }
633*e4b17023SJohn Marino
634*e4b17023SJohn Marino /* Finds the nearest common ancestor of the nodes N1 and N2. */
635*e4b17023SJohn Marino
636*e4b17023SJohn Marino struct et_node *
et_nca(struct et_node * n1,struct et_node * n2)637*e4b17023SJohn Marino et_nca (struct et_node *n1, struct et_node *n2)
638*e4b17023SJohn Marino {
639*e4b17023SJohn Marino struct et_occ *o1 = n1->rightmost_occ, *o2 = n2->rightmost_occ, *om;
640*e4b17023SJohn Marino struct et_occ *l, *r, *ret;
641*e4b17023SJohn Marino int mn;
642*e4b17023SJohn Marino
643*e4b17023SJohn Marino if (n1 == n2)
644*e4b17023SJohn Marino return n1;
645*e4b17023SJohn Marino
646*e4b17023SJohn Marino et_splay (o1);
647*e4b17023SJohn Marino l = o1->prev;
648*e4b17023SJohn Marino r = o1->next;
649*e4b17023SJohn Marino if (l)
650*e4b17023SJohn Marino l->parent = NULL;
651*e4b17023SJohn Marino if (r)
652*e4b17023SJohn Marino r->parent = NULL;
653*e4b17023SJohn Marino et_splay (o2);
654*e4b17023SJohn Marino
655*e4b17023SJohn Marino if (l == o2 || (l && l->parent != NULL))
656*e4b17023SJohn Marino {
657*e4b17023SJohn Marino ret = o2->next;
658*e4b17023SJohn Marino
659*e4b17023SJohn Marino set_prev (o1, o2);
660*e4b17023SJohn Marino if (r)
661*e4b17023SJohn Marino r->parent = o1;
662*e4b17023SJohn Marino }
663*e4b17023SJohn Marino else if (r == o2 || (r && r->parent != NULL))
664*e4b17023SJohn Marino {
665*e4b17023SJohn Marino ret = o2->prev;
666*e4b17023SJohn Marino
667*e4b17023SJohn Marino set_next (o1, o2);
668*e4b17023SJohn Marino if (l)
669*e4b17023SJohn Marino l->parent = o1;
670*e4b17023SJohn Marino }
671*e4b17023SJohn Marino else
672*e4b17023SJohn Marino {
673*e4b17023SJohn Marino /* O1 and O2 are in different components of the forest. */
674*e4b17023SJohn Marino if (l)
675*e4b17023SJohn Marino l->parent = o1;
676*e4b17023SJohn Marino if (r)
677*e4b17023SJohn Marino r->parent = o1;
678*e4b17023SJohn Marino return NULL;
679*e4b17023SJohn Marino }
680*e4b17023SJohn Marino
681*e4b17023SJohn Marino if (0 < o2->depth)
682*e4b17023SJohn Marino {
683*e4b17023SJohn Marino om = o1;
684*e4b17023SJohn Marino mn = o1->depth;
685*e4b17023SJohn Marino }
686*e4b17023SJohn Marino else
687*e4b17023SJohn Marino {
688*e4b17023SJohn Marino om = o2;
689*e4b17023SJohn Marino mn = o2->depth + o1->depth;
690*e4b17023SJohn Marino }
691*e4b17023SJohn Marino
692*e4b17023SJohn Marino #ifdef DEBUG_ET
693*e4b17023SJohn Marino et_check_tree_sanity (o2);
694*e4b17023SJohn Marino #endif
695*e4b17023SJohn Marino
696*e4b17023SJohn Marino if (ret && ret->min + o1->depth + o2->depth < mn)
697*e4b17023SJohn Marino return ret->min_occ->of;
698*e4b17023SJohn Marino else
699*e4b17023SJohn Marino return om->of;
700*e4b17023SJohn Marino }
701*e4b17023SJohn Marino
702*e4b17023SJohn Marino /* Checks whether the node UP is an ancestor of the node DOWN. */
703*e4b17023SJohn Marino
704*e4b17023SJohn Marino bool
et_below(struct et_node * down,struct et_node * up)705*e4b17023SJohn Marino et_below (struct et_node *down, struct et_node *up)
706*e4b17023SJohn Marino {
707*e4b17023SJohn Marino struct et_occ *u = up->rightmost_occ, *d = down->rightmost_occ;
708*e4b17023SJohn Marino struct et_occ *l, *r;
709*e4b17023SJohn Marino
710*e4b17023SJohn Marino if (up == down)
711*e4b17023SJohn Marino return true;
712*e4b17023SJohn Marino
713*e4b17023SJohn Marino et_splay (u);
714*e4b17023SJohn Marino l = u->prev;
715*e4b17023SJohn Marino r = u->next;
716*e4b17023SJohn Marino
717*e4b17023SJohn Marino if (!l)
718*e4b17023SJohn Marino return false;
719*e4b17023SJohn Marino
720*e4b17023SJohn Marino l->parent = NULL;
721*e4b17023SJohn Marino
722*e4b17023SJohn Marino if (r)
723*e4b17023SJohn Marino r->parent = NULL;
724*e4b17023SJohn Marino
725*e4b17023SJohn Marino et_splay (d);
726*e4b17023SJohn Marino
727*e4b17023SJohn Marino if (l == d || l->parent != NULL)
728*e4b17023SJohn Marino {
729*e4b17023SJohn Marino if (r)
730*e4b17023SJohn Marino r->parent = u;
731*e4b17023SJohn Marino set_prev (u, d);
732*e4b17023SJohn Marino #ifdef DEBUG_ET
733*e4b17023SJohn Marino et_check_tree_sanity (u);
734*e4b17023SJohn Marino #endif
735*e4b17023SJohn Marino }
736*e4b17023SJohn Marino else
737*e4b17023SJohn Marino {
738*e4b17023SJohn Marino l->parent = u;
739*e4b17023SJohn Marino
740*e4b17023SJohn Marino /* In case O1 and O2 are in two different trees, we must just restore the
741*e4b17023SJohn Marino original state. */
742*e4b17023SJohn Marino if (r && r->parent != NULL)
743*e4b17023SJohn Marino set_next (u, d);
744*e4b17023SJohn Marino else
745*e4b17023SJohn Marino set_next (u, r);
746*e4b17023SJohn Marino
747*e4b17023SJohn Marino #ifdef DEBUG_ET
748*e4b17023SJohn Marino et_check_tree_sanity (u);
749*e4b17023SJohn Marino #endif
750*e4b17023SJohn Marino return false;
751*e4b17023SJohn Marino }
752*e4b17023SJohn Marino
753*e4b17023SJohn Marino if (0 >= d->depth)
754*e4b17023SJohn Marino return false;
755*e4b17023SJohn Marino
756*e4b17023SJohn Marino return !d->next || d->next->min + d->depth >= 0;
757*e4b17023SJohn Marino }
758*e4b17023SJohn Marino
759*e4b17023SJohn Marino /* Returns the root of the tree that contains NODE. */
760*e4b17023SJohn Marino
761*e4b17023SJohn Marino struct et_node *
et_root(struct et_node * node)762*e4b17023SJohn Marino et_root (struct et_node *node)
763*e4b17023SJohn Marino {
764*e4b17023SJohn Marino struct et_occ *occ = node->rightmost_occ, *r;
765*e4b17023SJohn Marino
766*e4b17023SJohn Marino /* The root of the tree corresponds to the rightmost occurrence in the
767*e4b17023SJohn Marino represented path. */
768*e4b17023SJohn Marino et_splay (occ);
769*e4b17023SJohn Marino for (r = occ; r->next; r = r->next)
770*e4b17023SJohn Marino continue;
771*e4b17023SJohn Marino et_splay (r);
772*e4b17023SJohn Marino
773*e4b17023SJohn Marino return r->of;
774*e4b17023SJohn Marino }
775