xref: /dflybsd-src/contrib/gcc-4.7/gcc/et-forest.c (revision 04febcfb30580676d3e95f58a16c5137ee478b32)
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