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