3ce0f1e3387238c91996ffdc5ccc8393c4b7fe50
[libfirm] / ir / adt / pqueue.c
1 /*
2  * Copyright (C) 1995-2008 University of Karlsruhe.  All right reserved.
3  *
4  * This file is part of libFirm.
5  *
6  * This file may be distributed and/or modified under the terms of the
7  * GNU General Public License version 2 as published by the Free Software
8  * Foundation and appearing in the file LICENSE.GPL included in the
9  * packaging of this file.
10  *
11  * Licensees holding valid libFirm Professional Edition licenses may use
12  * this file in accordance with the libFirm Commercial License.
13  * Agreement provided with the Software.
14  *
15  * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
16  * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR
17  * PURPOSE.
18  */
19
20 /**
21  * @file
22  * @author  Christian Wuerdig, Matthias Braun
23  * @brief   Priority Queue implementation based on the heap datastructure
24  * @version $Id$
25  */
26 #include "array.h"
27 #include "pqueue.h"
28
29 /*
30  * Implements a heap.
31  *
32  * Implementation note: It might seem strange that we start indexing at 0
33  * but use 2*i and 2*i+1 to find the left and right sucessor of an index.
34  * The trick is that for index 0 the left successor is 0 again, and the
35  * right successor is 1 in this scheme. For the right successor 1 everything
36  * works like usual. We simply took care in the algorithms that they still
37  * work with the left child of 0 being 0 again. This was possible without
38  * any extra ifs or arithmetic.
39  * Thus we can save the wastage of 1 array position you can see in other
40  * implementations or the ugly (i+1)*2 - 1 and (i+1)*2 for calculating the
41  * left and right child. (At the expense that stuff easily breaks when you make
42  * changes and don't think that the left child of 0 is 0 :-/)
43  *
44  */
45
46 typedef struct _pqueue_el_t {
47         void *data;
48         int  key;
49 } pqueue_el_t;
50
51 struct _pqueue_t {
52         pqueue_el_t *elems;
53 };
54
55 /**
56  * Enforces the heap characteristics if the queue
57  * starting from element at position @p pos.
58  */
59 static void pqueue_heapify(pqueue_t *q, unsigned pos) {
60         unsigned len = ARR_LEN(q->elems);
61
62         while (pos * 2 < len) {
63                 pqueue_el_t tmp;
64                 unsigned    exchange = pos;
65
66                 if (q->elems[exchange].key < q->elems[pos * 2].key) {
67                         exchange = pos * 2;
68                 }
69
70                 if ((pos * 2 + 1) < len && q->elems[exchange].key < q->elems[pos * 2 + 1].key) {
71                         exchange = pos * 2 + 1;
72                 }
73
74                 if (exchange == pos)
75                         break;
76
77                 tmp                = q->elems[pos];
78                 q->elems[pos]      = q->elems[exchange];
79                 q->elems[exchange] = tmp;
80
81                 pos = exchange;
82         }
83 }
84
85 /**
86  * Sifts up a newly inserted element at position @p pos.
87  */
88 static void pqueue_sift_up(pqueue_t *q, unsigned pos) {
89         while(q->elems[pos].key > q->elems[pos / 2].key) {
90                 pqueue_el_t tmp;
91
92                 tmp               = q->elems[pos];
93                 q->elems[pos]     = q->elems[pos / 2];
94                 q->elems[pos / 2] = tmp;
95
96                 pos /= 2;
97         }
98 }
99
100 pqueue_t *new_pqueue(void) {
101         pqueue_t *res = xmalloc(sizeof(*res));
102         res->elems = NEW_ARR_F(pqueue_el_t, 0);
103         return res;
104 }
105
106 void del_pqueue(pqueue_t *q) {
107         DEL_ARR_F(q->elems);
108         free(q);
109 }
110
111 void pqueue_put(pqueue_t *q, void *data, int key) {
112         pqueue_el_t el;
113
114         el.data = data;
115         el.key  = key;
116
117         ARR_APP1(pqueue_el_t, q->elems, el);
118
119         pqueue_sift_up(q, ARR_LEN(q->elems) - 1);
120 }
121
122 void *pqueue_pop_front(pqueue_t *q) {
123         switch(ARR_LEN(q->elems)) {
124                 case 0:
125                         assert(0 && "Attempt to retrieve element from empty priority queue.");
126                         return NULL;
127                         break;
128                 case 1:
129                         ARR_SHRINKLEN(q->elems, 0);
130                         return q->elems[0].data;
131                         break;
132                 default: {
133                         void *data = q->elems[0].data;
134                         int  len   = ARR_LEN(q->elems) - 1;
135
136                         q->elems[0] = q->elems[len];
137                         ARR_SHRINKLEN(q->elems, len);
138                         pqueue_heapify(q, 0);
139
140                         return data;
141                 }
142         }
143 }
144
145 int pqueue_length(const pqueue_t *q) {
146         return ARR_LEN(q->elems);
147 }
148
149 int pqueue_empty(const pqueue_t *q) {
150         return ARR_LEN(q->elems) == 0;
151 }