rb.c revision 1.3 1 1.3 matt /* $NetBSD: rb.c,v 1.3 2008/06/30 20:54:19 matt Exp $ */
2 1.1 matt
3 1.1 matt /*-
4 1.1 matt * Copyright (c) 2001 The NetBSD Foundation, Inc.
5 1.1 matt * All rights reserved.
6 1.1 matt *
7 1.1 matt * This code is derived from software contributed to The NetBSD Foundation
8 1.1 matt * by Matt Thomas <matt (at) 3am-software.com>.
9 1.1 matt *
10 1.1 matt * Redistribution and use in source and binary forms, with or without
11 1.1 matt * modification, are permitted provided that the following conditions
12 1.1 matt * are met:
13 1.1 matt * 1. Redistributions of source code must retain the above copyright
14 1.1 matt * notice, this list of conditions and the following disclaimer.
15 1.1 matt * 2. Redistributions in binary form must reproduce the above copyright
16 1.1 matt * notice, this list of conditions and the following disclaimer in the
17 1.1 matt * documentation and/or other materials provided with the distribution.
18 1.1 matt *
19 1.1 matt * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20 1.1 matt * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21 1.1 matt * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22 1.1 matt * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23 1.1 matt * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24 1.1 matt * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25 1.1 matt * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26 1.1 matt * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27 1.1 matt * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28 1.1 matt * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29 1.1 matt * POSSIBILITY OF SUCH DAMAGE.
30 1.1 matt */
31 1.1 matt
32 1.1 matt #if !defined(_KERNEL) && !defined(_STANDALONE)
33 1.1 matt #include <sys/types.h>
34 1.1 matt #include <sys/types.h>
35 1.1 matt #include <stddef.h>
36 1.1 matt #include <assert.h>
37 1.1 matt #include <stdbool.h>
38 1.1 matt #ifdef RBDEBUG
39 1.1 matt #define KASSERT(s) assert(s)
40 1.1 matt #else
41 1.3 matt #define KASSERT(s) do { } while (/*CONSTCOND*/ 0)
42 1.1 matt #endif
43 1.1 matt #else
44 1.1 matt #include <lib/libkern/libkern.h>
45 1.1 matt #endif
46 1.1 matt
47 1.1 matt #ifdef _LIBC
48 1.1 matt __weak_alias(rb_tree_init, _rb_tree_init)
49 1.1 matt __weak_alias(rb_tree_find_node, _rb_tree_find_node)
50 1.1 matt __weak_alias(rb_tree_find_node_geq, _rb_tree_find_node_geq)
51 1.1 matt __weak_alias(rb_tree_find_node_leq, _rb_tree_find_node_leq)
52 1.1 matt __weak_alias(rb_tree_insert_node, _rb_tree_insert_node)
53 1.1 matt __weak_alias(rb_tree_remove_node, _rb_tree_remove_node)
54 1.1 matt __weak_alias(rb_tree_iterate, _rb_tree_iterate)
55 1.1 matt #ifdef RBDEBUG
56 1.1 matt __weak_alias(rb_tree_check, _rb_tree_check)
57 1.1 matt __weak_alias(rb_tree_depths, _rb_tree_depths)
58 1.1 matt #endif
59 1.1 matt
60 1.1 matt #define rb_tree_init _rb_tree_init
61 1.1 matt #define rb_tree_find_node _rb_tree_find_node
62 1.1 matt #define rb_tree_find_node_geq _rb_tree_find_node_geq
63 1.1 matt #define rb_tree_find_node_leq _rb_tree_find_node_leq
64 1.1 matt #define rb_tree_insert_node _rb_tree_insert_node
65 1.1 matt #define rb_tree_remove_node _rb_tree_remove_node
66 1.1 matt #define rb_tree_iterate _rb_tree_iterate
67 1.1 matt #ifdef RBDEBUG
68 1.1 matt #define rb_tree_check _rb_tree_check
69 1.1 matt #define rb_tree_depths _rb_tree_depths
70 1.1 matt #endif
71 1.1 matt #endif
72 1.1 matt
73 1.1 matt #ifdef RBTEST
74 1.1 matt #include "rb.h"
75 1.1 matt #else
76 1.1 matt #include <sys/rb.h>
77 1.1 matt #endif
78 1.1 matt
79 1.1 matt static void rb_tree_insert_rebalance(struct rb_tree *, struct rb_node *);
80 1.1 matt static void rb_tree_removal_rebalance(struct rb_tree *, struct rb_node *,
81 1.1 matt unsigned int);
82 1.1 matt #ifdef RBDEBUG
83 1.1 matt static const struct rb_node *rb_tree_iterate_const(const struct rb_tree *,
84 1.1 matt const struct rb_node *, const unsigned int);
85 1.1 matt static bool rb_tree_check_node(const struct rb_tree *, const struct rb_node *,
86 1.1 matt const struct rb_node *, bool);
87 1.1 matt #else
88 1.1 matt #define rb_tree_check_node(a, b, c, d) true
89 1.1 matt #endif
90 1.1 matt
91 1.1 matt #define RB_SENTINEL_NODE NULL
92 1.1 matt
93 1.1 matt void
94 1.1 matt rb_tree_init(struct rb_tree *rbt, const struct rb_tree_ops *ops)
95 1.1 matt {
96 1.1 matt rbt->rbt_ops = ops;
97 1.1 matt *((const struct rb_node **)&rbt->rbt_root) = RB_SENTINEL_NODE;
98 1.1 matt RB_TAILQ_INIT(&rbt->rbt_nodes);
99 1.1 matt #ifndef RBSMALL
100 1.1 matt rbt->rbt_minmax[RB_DIR_LEFT] = rbt->rbt_root; /* minimum node */
101 1.1 matt rbt->rbt_minmax[RB_DIR_RIGHT] = rbt->rbt_root; /* maximum node */
102 1.1 matt #endif
103 1.1 matt #ifdef RBSTATS
104 1.1 matt rbt->rbt_count = 0;
105 1.1 matt rbt->rbt_insertions = 0;
106 1.1 matt rbt->rbt_removals = 0;
107 1.1 matt rbt->rbt_insertion_rebalance_calls = 0;
108 1.1 matt rbt->rbt_insertion_rebalance_passes = 0;
109 1.1 matt rbt->rbt_removal_rebalance_calls = 0;
110 1.1 matt rbt->rbt_removal_rebalance_passes = 0;
111 1.1 matt #endif
112 1.1 matt }
113 1.1 matt
114 1.1 matt struct rb_node *
115 1.1 matt rb_tree_find_node(struct rb_tree *rbt, const void *key)
116 1.1 matt {
117 1.2 matt rbto_compare_key_fn compare_key = rbt->rbt_ops->rbto_compare_key;
118 1.1 matt struct rb_node *parent = rbt->rbt_root;
119 1.1 matt
120 1.1 matt while (!RB_SENTINEL_P(parent)) {
121 1.1 matt const signed int diff = (*compare_key)(parent, key);
122 1.1 matt if (diff == 0)
123 1.1 matt return parent;
124 1.1 matt parent = parent->rb_nodes[diff > 0];
125 1.1 matt }
126 1.1 matt
127 1.1 matt return NULL;
128 1.1 matt }
129 1.1 matt
130 1.1 matt struct rb_node *
131 1.1 matt rb_tree_find_node_geq(struct rb_tree *rbt, const void *key)
132 1.1 matt {
133 1.2 matt rbto_compare_key_fn compare_key = rbt->rbt_ops->rbto_compare_key;
134 1.1 matt struct rb_node *parent = rbt->rbt_root;
135 1.1 matt struct rb_node *last = NULL;
136 1.1 matt
137 1.1 matt while (!RB_SENTINEL_P(parent)) {
138 1.1 matt const signed int diff = (*compare_key)(parent, key);
139 1.1 matt if (diff == 0)
140 1.1 matt return parent;
141 1.1 matt if (diff < 0)
142 1.1 matt last = parent;
143 1.1 matt parent = parent->rb_nodes[diff > 0];
144 1.1 matt }
145 1.1 matt
146 1.1 matt return last;
147 1.1 matt }
148 1.1 matt
149 1.1 matt struct rb_node *
150 1.1 matt rb_tree_find_node_leq(struct rb_tree *rbt, const void *key)
151 1.1 matt {
152 1.2 matt rbto_compare_key_fn compare_key = rbt->rbt_ops->rbto_compare_key;
153 1.1 matt struct rb_node *parent = rbt->rbt_root;
154 1.1 matt struct rb_node *last = NULL;
155 1.1 matt
156 1.1 matt while (!RB_SENTINEL_P(parent)) {
157 1.1 matt const signed int diff = (*compare_key)(parent, key);
158 1.1 matt if (diff == 0)
159 1.1 matt return parent;
160 1.1 matt if (diff > 0)
161 1.1 matt last = parent;
162 1.1 matt parent = parent->rb_nodes[diff > 0];
163 1.1 matt }
164 1.1 matt
165 1.1 matt return last;
166 1.1 matt }
167 1.1 matt
168 1.1 matt bool
170 1.1 matt rb_tree_insert_node(struct rb_tree *rbt, struct rb_node *self)
171 1.2 matt {
172 1.1 matt rbto_compare_nodes_fn compare_nodes = rbt->rbt_ops->rbto_compare_nodes;
173 1.1 matt struct rb_node *parent, *tmp;
174 1.1 matt unsigned int position;
175 1.1 matt bool rebalance;
176 1.1 matt
177 1.1 matt RBSTAT_INC(rbt->rbt_insertions);
178 1.1 matt
179 1.1 matt tmp = rbt->rbt_root;
180 1.1 matt /*
181 1.1 matt * This is a hack. Because rbt->rbt_root is just a struct rb_node *,
182 1.1 matt * just like rb_node->rb_nodes[RB_DIR_LEFT], we can use this fact to
183 1.1 matt * avoid a lot of tests for root and know that even at root,
184 1.1 matt * updating RB_FATHER(rb_node)->rb_nodes[RB_POSITION(rb_node)] will
185 1.1 matt * update rbt->rbt_root.
186 1.3 matt */
187 1.1 matt parent = (struct rb_node *)(void *)&rbt->rbt_root;
188 1.1 matt position = RB_DIR_LEFT;
189 1.1 matt
190 1.1 matt /*
191 1.1 matt * Find out where to place this new leaf.
192 1.1 matt */
193 1.1 matt while (!RB_SENTINEL_P(tmp)) {
194 1.1 matt const signed int diff = (*compare_nodes)(tmp, self);
195 1.1 matt if (__predict_false(diff == 0)) {
196 1.1 matt /*
197 1.1 matt * Node already exists; don't insert.
198 1.1 matt */
199 1.1 matt return false;
200 1.1 matt }
201 1.1 matt parent = tmp;
202 1.1 matt position = (diff > 0);
203 1.1 matt tmp = parent->rb_nodes[position];
204 1.1 matt }
205 1.1 matt
206 1.1 matt #ifdef RBDEBUG
207 1.1 matt {
208 1.1 matt struct rb_node *prev = NULL, *next = NULL;
209 1.1 matt
210 1.1 matt if (position == RB_DIR_RIGHT)
211 1.1 matt prev = parent;
212 1.1 matt else if (tmp != rbt->rbt_root)
213 1.1 matt next = parent;
214 1.1 matt
215 1.1 matt /*
216 1.1 matt * Verify our sequential position
217 1.1 matt */
218 1.1 matt KASSERT(prev == NULL || !RB_SENTINEL_P(prev));
219 1.1 matt KASSERT(next == NULL || !RB_SENTINEL_P(next));
220 1.1 matt if (prev != NULL && next == NULL)
221 1.1 matt next = TAILQ_NEXT(prev, rb_link);
222 1.1 matt if (prev == NULL && next != NULL)
223 1.1 matt prev = TAILQ_PREV(next, rb_node_qh, rb_link);
224 1.1 matt KASSERT(prev == NULL || !RB_SENTINEL_P(prev));
225 1.1 matt KASSERT(next == NULL || !RB_SENTINEL_P(next));
226 1.1 matt KASSERT(prev == NULL || (*compare_nodes)(prev, self) > 0);
227 1.1 matt KASSERT(next == NULL || (*compare_nodes)(self, next) > 0);
228 1.1 matt }
229 1.1 matt #endif
230 1.1 matt
231 1.1 matt /*
232 1.1 matt * Initialize the node and insert as a leaf into the tree.
233 1.1 matt */
234 1.1 matt RB_SET_FATHER(self, parent);
235 1.3 matt RB_SET_POSITION(self, position);
236 1.1 matt if (__predict_false(parent == (struct rb_node *)(void *)&rbt->rbt_root)) {
237 1.1 matt RB_MARK_BLACK(self); /* root is always black */
238 1.1 matt #ifndef RBSMALL
239 1.1 matt rbt->rbt_minmax[RB_DIR_LEFT] = self;
240 1.1 matt rbt->rbt_minmax[RB_DIR_RIGHT] = self;
241 1.1 matt #endif
242 1.1 matt rebalance = false;
243 1.1 matt } else {
244 1.1 matt KASSERT(position == RB_DIR_LEFT || position == RB_DIR_RIGHT);
245 1.1 matt #ifndef RBSMALL
246 1.1 matt /*
247 1.1 matt * Keep track of the minimum and maximum nodes. If our
248 1.1 matt * parent is a minmax node and we on their min/max side,
249 1.1 matt * we must be the new min/max node.
250 1.1 matt */
251 1.1 matt if (parent == rbt->rbt_minmax[position])
252 1.1 matt rbt->rbt_minmax[position] = self;
253 1.1 matt #endif /* !RBSMALL */
254 1.1 matt /*
255 1.1 matt * All new nodes are colored red. We only need to rebalance
256 1.1 matt * if our parent is also red.
257 1.1 matt */
258 1.1 matt RB_MARK_RED(self);
259 1.1 matt rebalance = RB_RED_P(parent);
260 1.1 matt }
261 1.1 matt KASSERT(RB_SENTINEL_P(parent->rb_nodes[position]));
262 1.1 matt self->rb_left = parent->rb_nodes[position];
263 1.1 matt self->rb_right = parent->rb_nodes[position];
264 1.1 matt parent->rb_nodes[position] = self;
265 1.1 matt KASSERT(RB_CHILDLESS_P(self));
266 1.1 matt
267 1.1 matt /*
268 1.1 matt * Insert the new node into a sorted list for easy sequential access
269 1.1 matt */
270 1.1 matt RBSTAT_INC(rbt->rbt_count);
271 1.1 matt #ifdef RBDEBUG
272 1.1 matt if (RB_ROOT_P(rbt, self)) {
273 1.1 matt RB_TAILQ_INSERT_HEAD(&rbt->rbt_nodes, self, rb_link);
274 1.1 matt } else if (position == RB_DIR_LEFT) {
275 1.1 matt KASSERT((*compare_nodes)(self, RB_FATHER(self)) > 0);
276 1.1 matt RB_TAILQ_INSERT_BEFORE(RB_FATHER(self), self, rb_link);
277 1.1 matt } else {
278 1.1 matt KASSERT((*compare_nodes)(RB_FATHER(self), self) > 0);
279 1.1 matt RB_TAILQ_INSERT_AFTER(&rbt->rbt_nodes, RB_FATHER(self),
280 1.1 matt self, rb_link);
281 1.1 matt }
282 1.1 matt #endif
283 1.1 matt KASSERT(rb_tree_check_node(rbt, self, NULL, !rebalance));
284 1.1 matt
285 1.1 matt /*
286 1.1 matt * Rebalance tree after insertion
287 1.1 matt */
288 1.1 matt if (rebalance) {
289 1.1 matt rb_tree_insert_rebalance(rbt, self);
290 1.1 matt KASSERT(rb_tree_check_node(rbt, self, NULL, true));
291 1.1 matt }
292 1.1 matt
293 1.1 matt return true;
294 1.1 matt }
295 1.1 matt
296 1.1 matt /*
298 1.1 matt * Swap the location and colors of 'self' and its child @ which. The child
299 1.1 matt * can not be a sentinel node. This is our rotation function. However,
300 1.1 matt * since it preserves coloring, it great simplifies both insertion and
301 1.1 matt * removal since rotation almost always involves the exchanging of colors
302 1.3 matt * as a separate step.
303 1.1 matt */
304 1.1 matt /*ARGSUSED*/
305 1.1 matt static void
306 1.1 matt rb_tree_reparent_nodes(struct rb_tree *rbt, struct rb_node *old_father,
307 1.1 matt const unsigned int which)
308 1.1 matt {
309 1.1 matt const unsigned int other = which ^ RB_DIR_OTHER;
310 1.1 matt struct rb_node * const grandpa = RB_FATHER(old_father);
311 1.1 matt struct rb_node * const old_child = old_father->rb_nodes[which];
312 1.1 matt struct rb_node * const new_father = old_child;
313 1.1 matt struct rb_node * const new_child = old_father;
314 1.1 matt
315 1.1 matt KASSERT(which == RB_DIR_LEFT || which == RB_DIR_RIGHT);
316 1.1 matt
317 1.1 matt KASSERT(!RB_SENTINEL_P(old_child));
318 1.1 matt KASSERT(RB_FATHER(old_child) == old_father);
319 1.1 matt
320 1.1 matt KASSERT(rb_tree_check_node(rbt, old_father, NULL, false));
321 1.1 matt KASSERT(rb_tree_check_node(rbt, old_child, NULL, false));
322 1.1 matt KASSERT(RB_ROOT_P(rbt, old_father) || rb_tree_check_node(rbt, grandpa, NULL, false));
323 1.1 matt
324 1.1 matt /*
325 1.1 matt * Exchange descendant linkages.
326 1.1 matt */
327 1.1 matt grandpa->rb_nodes[RB_POSITION(old_father)] = new_father;
328 1.1 matt new_child->rb_nodes[which] = old_child->rb_nodes[other];
329 1.1 matt new_father->rb_nodes[other] = new_child;
330 1.1 matt
331 1.1 matt /*
332 1.1 matt * Update ancestor linkages
333 1.1 matt */
334 1.1 matt RB_SET_FATHER(new_father, grandpa);
335 1.1 matt RB_SET_FATHER(new_child, new_father);
336 1.1 matt
337 1.1 matt /*
338 1.1 matt * Exchange properties between new_father and new_child. The only
339 1.1 matt * change is that new_child's position is now on the other side.
340 1.1 matt */
341 1.1 matt #if 0
342 1.1 matt {
343 1.1 matt struct rb_node tmp;
344 1.1 matt tmp.rb_info = 0;
345 1.1 matt RB_COPY_PROPERTIES(&tmp, old_child);
346 1.1 matt RB_COPY_PROPERTIES(new_father, old_father);
347 1.1 matt RB_COPY_PROPERTIES(new_child, &tmp);
348 1.1 matt }
349 1.1 matt #else
350 1.1 matt RB_SWAP_PROPERTIES(new_father, new_child);
351 1.1 matt #endif
352 1.1 matt RB_SET_POSITION(new_child, other);
353 1.1 matt
354 1.1 matt /*
355 1.1 matt * Make sure to reparent the new child to ourself.
356 1.1 matt */
357 1.1 matt if (!RB_SENTINEL_P(new_child->rb_nodes[which])) {
358 1.1 matt RB_SET_FATHER(new_child->rb_nodes[which], new_child);
359 1.1 matt RB_SET_POSITION(new_child->rb_nodes[which], which);
360 1.1 matt }
361 1.1 matt
362 1.1 matt KASSERT(rb_tree_check_node(rbt, new_father, NULL, false));
363 1.1 matt KASSERT(rb_tree_check_node(rbt, new_child, NULL, false));
364 1.1 matt KASSERT(RB_ROOT_P(rbt, new_father) || rb_tree_check_node(rbt, grandpa, NULL, false));
365 1.1 matt }
366 1.1 matt
367 1.1 matt static void
369 1.1 matt rb_tree_insert_rebalance(struct rb_tree *rbt, struct rb_node *self)
370 1.1 matt {
371 1.1 matt struct rb_node * father = RB_FATHER(self);
372 1.1 matt struct rb_node * grandpa = RB_FATHER(father);
373 1.1 matt struct rb_node * uncle;
374 1.1 matt unsigned int which;
375 1.1 matt unsigned int other;
376 1.1 matt
377 1.1 matt KASSERT(!RB_ROOT_P(rbt, self));
378 1.1 matt KASSERT(RB_RED_P(self));
379 1.1 matt KASSERT(RB_RED_P(father));
380 1.1 matt RBSTAT_INC(rbt->rbt_insertion_rebalance_calls);
381 1.1 matt
382 1.1 matt for (;;) {
383 1.1 matt KASSERT(!RB_SENTINEL_P(self));
384 1.1 matt
385 1.1 matt KASSERT(RB_RED_P(self));
386 1.1 matt KASSERT(RB_RED_P(father));
387 1.1 matt /*
388 1.1 matt * We are red and our parent is red, therefore we must have a
389 1.1 matt * grandfather and he must be black.
390 1.1 matt */
391 1.1 matt grandpa = RB_FATHER(father);
392 1.1 matt KASSERT(RB_BLACK_P(grandpa));
393 1.1 matt KASSERT(RB_DIR_RIGHT == 1 && RB_DIR_LEFT == 0);
394 1.1 matt which = (father == grandpa->rb_right);
395 1.1 matt other = which ^ RB_DIR_OTHER;
396 1.1 matt uncle = grandpa->rb_nodes[other];
397 1.1 matt
398 1.1 matt if (RB_BLACK_P(uncle))
399 1.1 matt break;
400 1.1 matt
401 1.1 matt RBSTAT_INC(rbt->rbt_insertion_rebalance_passes);
402 1.1 matt /*
403 1.1 matt * Case 1: our uncle is red
404 1.1 matt * Simply invert the colors of our parent and
405 1.1 matt * uncle and make our grandparent red. And
406 1.1 matt * then solve the problem up at his level.
407 1.1 matt */
408 1.1 matt RB_MARK_BLACK(uncle);
409 1.1 matt RB_MARK_BLACK(father);
410 1.1 matt if (__predict_false(RB_ROOT_P(rbt, grandpa))) {
411 1.1 matt /*
412 1.1 matt * If our grandpa is root, don't bother
413 1.1 matt * setting him to red, just return.
414 1.1 matt */
415 1.1 matt KASSERT(RB_BLACK_P(grandpa));
416 1.1 matt return;
417 1.1 matt }
418 1.1 matt RB_MARK_RED(grandpa);
419 1.1 matt self = grandpa;
420 1.1 matt father = RB_FATHER(self);
421 1.1 matt KASSERT(RB_RED_P(self));
422 1.1 matt if (RB_BLACK_P(father)) {
423 1.1 matt /*
424 1.1 matt * If our greatgrandpa is black, we're done.
425 1.1 matt */
426 1.1 matt KASSERT(RB_BLACK_P(rbt->rbt_root));
427 1.1 matt return;
428 1.1 matt }
429 1.1 matt }
430 1.1 matt
431 1.1 matt KASSERT(!RB_ROOT_P(rbt, self));
432 1.1 matt KASSERT(RB_RED_P(self));
433 1.1 matt KASSERT(RB_RED_P(father));
434 1.1 matt KASSERT(RB_BLACK_P(uncle));
435 1.1 matt KASSERT(RB_BLACK_P(grandpa));
436 1.1 matt /*
437 1.1 matt * Case 2&3: our uncle is black.
438 1.1 matt */
439 1.1 matt if (self == father->rb_nodes[other]) {
440 1.1 matt /*
441 1.1 matt * Case 2: we are on the same side as our uncle
442 1.1 matt * Swap ourselves with our parent so this case
443 1.1 matt * becomes case 3. Basically our parent becomes our
444 1.1 matt * child.
445 1.1 matt */
446 1.1 matt rb_tree_reparent_nodes(rbt, father, other);
447 1.1 matt KASSERT(RB_FATHER(father) == self);
448 1.1 matt KASSERT(self->rb_nodes[which] == father);
449 1.1 matt KASSERT(RB_FATHER(self) == grandpa);
450 1.1 matt self = father;
451 1.1 matt father = RB_FATHER(self);
452 1.1 matt }
453 1.1 matt KASSERT(RB_RED_P(self) && RB_RED_P(father));
454 1.1 matt KASSERT(grandpa->rb_nodes[which] == father);
455 1.1 matt /*
456 1.1 matt * Case 3: we are opposite a child of a black uncle.
457 1.1 matt * Swap our parent and grandparent. Since our grandfather
458 1.1 matt * is black, our father will become black and our new sibling
459 1.1 matt * (former grandparent) will become red.
460 1.1 matt */
461 1.1 matt rb_tree_reparent_nodes(rbt, grandpa, which);
462 1.1 matt KASSERT(RB_FATHER(self) == father);
463 1.1 matt KASSERT(RB_FATHER(self)->rb_nodes[RB_POSITION(self) ^ RB_DIR_OTHER] == grandpa);
464 1.1 matt KASSERT(RB_RED_P(self));
465 1.1 matt KASSERT(RB_BLACK_P(father));
466 1.1 matt KASSERT(RB_RED_P(grandpa));
467 1.1 matt
468 1.1 matt /*
469 1.1 matt * Final step: Set the root to black.
470 1.1 matt */
471 1.1 matt RB_MARK_BLACK(rbt->rbt_root);
472 1.1 matt }
473 1.1 matt
474 1.1 matt static void
476 1.1 matt rb_tree_prune_node(struct rb_tree *rbt, struct rb_node *self, bool rebalance)
477 1.1 matt {
478 1.1 matt const unsigned int which = RB_POSITION(self);
479 1.1 matt struct rb_node *father = RB_FATHER(self);
480 1.1 matt const bool was_root = RB_ROOT_P(rbt, self);
481 1.1 matt
482 1.1 matt KASSERT(rebalance || (RB_ROOT_P(rbt, self) || RB_RED_P(self)));
483 1.1 matt KASSERT(!rebalance || RB_BLACK_P(self));
484 1.1 matt KASSERT(RB_CHILDLESS_P(self));
485 1.1 matt KASSERT(rb_tree_check_node(rbt, self, NULL, false));
486 1.1 matt
487 1.1 matt /*
488 1.1 matt * Since we are childless, we know that self->rb_left is pointing
489 1.1 matt * to the sentinel node.
490 1.1 matt */
491 1.1 matt father->rb_nodes[which] = self->rb_left;
492 1.1 matt
493 1.1 matt /*
494 1.1 matt * Remove ourselves from the node list, decrement the count,
495 1.1 matt * and update min/max.
496 1.1 matt */
497 1.1 matt RB_TAILQ_REMOVE(&rbt->rbt_nodes, self, rb_link);
498 1.1 matt RBSTAT_DEC(rbt->rbt_count);
499 1.1 matt #ifndef RBSMALL
500 1.1 matt if (__predict_false(rbt->rbt_minmax[RB_POSITION(self)] == self)) {
501 1.1 matt rbt->rbt_minmax[RB_POSITION(self)] = father;
502 1.1 matt /*
503 1.1 matt * When removing the root, rbt->rbt_minmax[RB_DIR_LEFT] is
504 1.1 matt * updated automatically, but we also need to update
505 1.1 matt * rbt->rbt_minmax[RB_DIR_RIGHT];
506 1.1 matt */
507 1.1 matt if (__predict_false(was_root)) {
508 1.1 matt rbt->rbt_minmax[RB_DIR_RIGHT] = father;
509 1.1 matt }
510 1.1 matt }
511 1.1 matt RB_SET_FATHER(self, NULL);
512 1.1 matt #endif
513 1.1 matt
514 1.1 matt /*
515 1.1 matt * Rebalance if requested.
516 1.1 matt */
517 1.1 matt if (rebalance)
518 1.1 matt rb_tree_removal_rebalance(rbt, father, which);
519 1.1 matt KASSERT(was_root || rb_tree_check_node(rbt, father, NULL, true));
520 1.1 matt }
521 1.1 matt
522 1.1 matt /*
524 1.1 matt * When deleting an interior node
525 1.1 matt */
526 1.1 matt static void
527 1.1 matt rb_tree_swap_prune_and_rebalance(struct rb_tree *rbt, struct rb_node *self,
528 1.1 matt struct rb_node *standin)
529 1.1 matt {
530 1.1 matt const unsigned int standin_which = RB_POSITION(standin);
531 1.1 matt unsigned int standin_other = standin_which ^ RB_DIR_OTHER;
532 1.1 matt struct rb_node *standin_son;
533 1.1 matt struct rb_node *standin_father = RB_FATHER(standin);
534 1.1 matt bool rebalance = RB_BLACK_P(standin);
535 1.1 matt
536 1.1 matt if (standin_father == self) {
537 1.1 matt /*
538 1.1 matt * As a child of self, any childen would be opposite of
539 1.1 matt * our parent.
540 1.1 matt */
541 1.1 matt KASSERT(RB_SENTINEL_P(standin->rb_nodes[standin_other]));
542 1.1 matt standin_son = standin->rb_nodes[standin_which];
543 1.1 matt } else {
544 1.1 matt /*
545 1.1 matt * Since we aren't a child of self, any childen would be
546 1.1 matt * on the same side as our parent.
547 1.1 matt */
548 1.1 matt KASSERT(RB_SENTINEL_P(standin->rb_nodes[standin_which]));
549 1.1 matt standin_son = standin->rb_nodes[standin_other];
550 1.1 matt }
551 1.1 matt
552 1.1 matt /*
553 1.1 matt * the node we are removing must have two children.
554 1.1 matt */
555 1.1 matt KASSERT(RB_TWOCHILDREN_P(self));
556 1.1 matt /*
557 1.1 matt * If standin has a child, it must be red.
558 1.1 matt */
559 1.1 matt KASSERT(RB_SENTINEL_P(standin_son) || RB_RED_P(standin_son));
560 1.1 matt
561 1.1 matt /*
562 1.1 matt * Verify things are sane.
563 1.1 matt */
564 1.1 matt KASSERT(rb_tree_check_node(rbt, self, NULL, false));
565 1.1 matt KASSERT(rb_tree_check_node(rbt, standin, NULL, false));
566 1.1 matt
567 1.1 matt if (__predict_false(RB_RED_P(standin_son))) {
568 1.1 matt /*
569 1.1 matt * We know we have a red child so if we flip it to black
570 1.1 matt * we don't have to rebalance.
571 1.1 matt */
572 1.1 matt KASSERT(rb_tree_check_node(rbt, standin_son, NULL, true));
573 1.1 matt RB_MARK_BLACK(standin_son);
574 1.1 matt rebalance = false;
575 1.1 matt
576 1.1 matt if (standin_father == self) {
577 1.1 matt KASSERT(RB_POSITION(standin_son) == standin_which);
578 1.1 matt } else {
579 1.1 matt KASSERT(RB_POSITION(standin_son) == standin_other);
580 1.1 matt /*
581 1.1 matt * Change the son's parentage to point to his grandpa.
582 1.1 matt */
583 1.1 matt RB_SET_FATHER(standin_son, standin_father);
584 1.1 matt RB_SET_POSITION(standin_son, standin_which);
585 1.1 matt }
586 1.1 matt }
587 1.1 matt
588 1.1 matt if (standin_father == self) {
589 1.1 matt /*
590 1.1 matt * If we are about to delete the standin's father, then when
591 1.1 matt * we call rebalance, we need to use ourselves as our father.
592 1.1 matt * Otherwise remember our original father. Also, sincef we are
593 1.1 matt * our standin's father we only need to reparent the standin's
594 1.1 matt * brother.
595 1.1 matt *
596 1.1 matt * | R --> S |
597 1.1 matt * | Q S --> Q T |
598 1.1 matt * | t --> |
599 1.1 matt */
600 1.1 matt KASSERT(RB_SENTINEL_P(standin->rb_nodes[standin_other]));
601 1.1 matt KASSERT(!RB_SENTINEL_P(self->rb_nodes[standin_other]));
602 1.1 matt KASSERT(self->rb_nodes[standin_which] == standin);
603 1.1 matt /*
604 1.1 matt * Have our son/standin adopt his brother as his new son.
605 1.1 matt */
606 1.1 matt standin_father = standin;
607 1.1 matt } else {
608 1.1 matt /*
609 1.1 matt * | R --> S . |
610 1.1 matt * | / \ | T --> / \ | / |
611 1.1 matt * | ..... | S --> ..... | T |
612 1.1 matt *
613 1.1 matt * Sever standin's connection to his father.
614 1.1 matt */
615 1.1 matt standin_father->rb_nodes[standin_which] = standin_son;
616 1.1 matt /*
617 1.1 matt * Adopt the far son.
618 1.1 matt */
619 1.1 matt standin->rb_nodes[standin_other] = self->rb_nodes[standin_other];
620 1.1 matt RB_SET_FATHER(standin->rb_nodes[standin_other], standin);
621 1.1 matt KASSERT(RB_POSITION(self->rb_nodes[standin_other]) == standin_other);
622 1.1 matt /*
623 1.1 matt * Use standin_other because we need to preserve standin_which
624 1.1 matt * for the removal_rebalance.
625 1.1 matt */
626 1.1 matt standin_other = standin_which;
627 1.1 matt }
628 1.1 matt
629 1.1 matt /*
630 1.1 matt * Move the only remaining son to our standin. If our standin is our
631 1.1 matt * son, this will be the only son needed to be moved.
632 1.1 matt */
633 1.1 matt KASSERT(standin->rb_nodes[standin_other] != self->rb_nodes[standin_other]);
634 1.1 matt standin->rb_nodes[standin_other] = self->rb_nodes[standin_other];
635 1.1 matt RB_SET_FATHER(standin->rb_nodes[standin_other], standin);
636 1.1 matt
637 1.1 matt /*
638 1.1 matt * Now copy the result of self to standin and then replace
639 1.1 matt * self with standin in the tree.
640 1.1 matt */
641 1.1 matt RB_COPY_PROPERTIES(standin, self);
642 1.1 matt RB_SET_FATHER(standin, RB_FATHER(self));
643 1.1 matt RB_FATHER(standin)->rb_nodes[RB_POSITION(standin)] = standin;
644 1.1 matt
645 1.1 matt /*
646 1.1 matt * Remove ourselves from the node list, decrement the count,
647 1.1 matt * and update min/max.
648 1.1 matt */
649 1.1 matt RB_TAILQ_REMOVE(&rbt->rbt_nodes, self, rb_link);
650 1.1 matt RBSTAT_DEC(rbt->rbt_count);
651 1.1 matt #ifndef RBSMALL
652 1.1 matt if (__predict_false(rbt->rbt_minmax[RB_POSITION(self)] == self))
653 1.1 matt rbt->rbt_minmax[RB_POSITION(self)] = RB_FATHER(self);
654 1.1 matt RB_SET_FATHER(self, NULL);
655 1.1 matt #endif
656 1.1 matt
657 1.1 matt KASSERT(rb_tree_check_node(rbt, standin, NULL, false));
658 1.1 matt KASSERT(RB_FATHER_SENTINEL_P(standin)
659 1.1 matt || rb_tree_check_node(rbt, standin_father, NULL, false));
660 1.1 matt KASSERT(RB_LEFT_SENTINEL_P(standin)
661 1.1 matt || rb_tree_check_node(rbt, standin->rb_left, NULL, false));
662 1.1 matt KASSERT(RB_RIGHT_SENTINEL_P(standin)
663 1.1 matt || rb_tree_check_node(rbt, standin->rb_right, NULL, false));
664 1.1 matt
665 1.1 matt if (!rebalance)
666 1.1 matt return;
667 1.1 matt
668 1.1 matt rb_tree_removal_rebalance(rbt, standin_father, standin_which);
669 1.1 matt KASSERT(rb_tree_check_node(rbt, standin, NULL, true));
670 1.1 matt }
671 1.1 matt
672 1.1 matt /*
673 1.1 matt * We could do this by doing
674 1.1 matt * rb_tree_node_swap(rbt, self, which);
675 1.1 matt * rb_tree_prune_node(rbt, self, false);
676 1.1 matt *
677 1.1 matt * But it's more efficient to just evalate and recolor the child.
678 1.1 matt */
679 1.1 matt static void
680 1.1 matt rb_tree_prune_blackred_branch(struct rb_tree *rbt, struct rb_node *self,
681 1.1 matt unsigned int which)
682 1.1 matt {
683 1.1 matt struct rb_node *father = RB_FATHER(self);
684 1.1 matt struct rb_node *son = self->rb_nodes[which];
685 1.1 matt const bool was_root = RB_ROOT_P(rbt, self);
686 1.1 matt
687 1.1 matt KASSERT(which == RB_DIR_LEFT || which == RB_DIR_RIGHT);
688 1.1 matt KASSERT(RB_BLACK_P(self) && RB_RED_P(son));
689 1.1 matt KASSERT(!RB_TWOCHILDREN_P(son));
690 1.1 matt KASSERT(RB_CHILDLESS_P(son));
691 1.1 matt KASSERT(rb_tree_check_node(rbt, self, NULL, false));
692 1.1 matt KASSERT(rb_tree_check_node(rbt, son, NULL, false));
693 1.1 matt
694 1.1 matt /*
695 1.1 matt * Remove ourselves from the tree and give our former child our
696 1.1 matt * properties (position, color, root).
697 1.1 matt */
698 1.1 matt RB_COPY_PROPERTIES(son, self);
699 1.1 matt father->rb_nodes[RB_POSITION(son)] = son;
700 1.1 matt RB_SET_FATHER(son, father);
701 1.1 matt
702 1.1 matt /*
703 1.1 matt * Remove ourselves from the node list, decrement the count,
704 1.1 matt * and update minmax.
705 1.1 matt */
706 1.1 matt RB_TAILQ_REMOVE(&rbt->rbt_nodes, self, rb_link);
707 1.1 matt RBSTAT_DEC(rbt->rbt_count);
708 1.1 matt #ifndef RBSMALL
709 1.1 matt if (__predict_false(was_root)) {
710 1.1 matt KASSERT(rbt->rbt_minmax[which] == son);
711 1.1 matt rbt->rbt_minmax[which ^ RB_DIR_OTHER] = son;
712 1.1 matt } else if (rbt->rbt_minmax[RB_POSITION(self)] == self) {
713 1.1 matt rbt->rbt_minmax[RB_POSITION(self)] = son;
714 1.1 matt }
715 1.1 matt RB_SET_FATHER(self, NULL);
716 1.1 matt #endif
717 1.1 matt
718 1.1 matt KASSERT(was_root || rb_tree_check_node(rbt, father, NULL, true));
719 1.1 matt KASSERT(rb_tree_check_node(rbt, son, NULL, true));
720 1.1 matt }
721 1.1 matt /*
722 1.1 matt *
723 1.1 matt */
724 1.1 matt void
725 1.1 matt rb_tree_remove_node(struct rb_tree *rbt, struct rb_node *self)
726 1.1 matt {
727 1.1 matt struct rb_node *standin;
728 1.1 matt unsigned int which;
729 1.1 matt
730 1.1 matt KASSERT(!RB_SENTINEL_P(self));
731 1.1 matt RBSTAT_INC(rbt->rbt_removals);
732 1.1 matt
733 1.1 matt /*
734 1.1 matt * In the following diagrams, we (the node to be removed) are S. Red
735 1.1 matt * nodes are lowercase. T could be either red or black.
736 1.1 matt *
737 1.1 matt * Remember the major axiom of the red-black tree: the number of
738 1.1 matt * black nodes from the root to each leaf is constant across all
739 1.1 matt * leaves, only the number of red nodes varies.
740 1.1 matt *
741 1.1 matt * Thus removing a red leaf doesn't require any other changes to a
742 1.1 matt * red-black tree. So if we must remove a node, attempt to rearrange
743 1.1 matt * the tree so we can remove a red node.
744 1.1 matt *
745 1.1 matt * The simpliest case is a childless red node or a childless root node:
746 1.1 matt *
747 1.1 matt * | T --> T | or | R --> * |
748 1.1 matt * | s --> * |
749 1.1 matt */
750 1.1 matt if (RB_CHILDLESS_P(self)) {
751 1.1 matt const bool rebalance = RB_BLACK_P(self) && !RB_ROOT_P(rbt, self);
752 1.1 matt rb_tree_prune_node(rbt, self, rebalance);
753 1.1 matt return;
754 1.1 matt }
755 1.1 matt KASSERT(!RB_CHILDLESS_P(self));
756 1.1 matt if (!RB_TWOCHILDREN_P(self)) {
757 1.1 matt /*
758 1.1 matt * The next simpliest case is the node we are deleting is
759 1.1 matt * black and has one red child.
760 1.1 matt *
761 1.1 matt * | T --> T --> T |
762 1.1 matt * | S --> R --> R |
763 1.1 matt * | r --> s --> * |
764 1.1 matt */
765 1.1 matt which = RB_LEFT_SENTINEL_P(self) ? RB_DIR_RIGHT : RB_DIR_LEFT;
766 1.1 matt KASSERT(RB_BLACK_P(self));
767 1.1 matt KASSERT(RB_RED_P(self->rb_nodes[which]));
768 1.1 matt KASSERT(RB_CHILDLESS_P(self->rb_nodes[which]));
769 1.1 matt rb_tree_prune_blackred_branch(rbt, self, which);
770 1.1 matt return;
771 1.1 matt }
772 1.1 matt KASSERT(RB_TWOCHILDREN_P(self));
773 1.1 matt
774 1.1 matt /*
775 1.1 matt * We invert these because we prefer to remove from the inside of
776 1.1 matt * the tree.
777 1.1 matt */
778 1.1 matt which = RB_POSITION(self) ^ RB_DIR_OTHER;
779 1.1 matt
780 1.1 matt /*
781 1.1 matt * Let's find the node closes to us opposite of our parent
782 1.1 matt * Now swap it with ourself, "prune" it, and rebalance, if needed.
783 1.1 matt */
784 1.1 matt standin = rb_tree_iterate(rbt, self, which);
785 1.1 matt rb_tree_swap_prune_and_rebalance(rbt, self, standin);
786 1.1 matt }
787 1.1 matt
788 1.1 matt static void
789 1.1 matt rb_tree_removal_rebalance(struct rb_tree *rbt, struct rb_node *parent,
790 1.1 matt unsigned int which)
791 1.1 matt {
792 1.1 matt KASSERT(!RB_SENTINEL_P(parent));
793 1.1 matt KASSERT(RB_SENTINEL_P(parent->rb_nodes[which]));
794 1.1 matt KASSERT(which == RB_DIR_LEFT || which == RB_DIR_RIGHT);
795 1.1 matt RBSTAT_INC(rbt->rbt_removal_rebalance_calls);
796 1.1 matt
797 1.1 matt while (RB_BLACK_P(parent->rb_nodes[which])) {
798 1.1 matt unsigned int other = which ^ RB_DIR_OTHER;
799 1.1 matt struct rb_node *brother = parent->rb_nodes[other];
800 1.1 matt
801 1.1 matt RBSTAT_INC(rbt->rbt_removal_rebalance_passes);
802 1.1 matt
803 1.1 matt KASSERT(!RB_SENTINEL_P(brother));
804 1.1 matt /*
805 1.1 matt * For cases 1, 2a, and 2b, our brother's children must
806 1.1 matt * be black and our father must be black
807 1.1 matt */
808 1.1 matt if (RB_BLACK_P(parent)
809 1.1 matt && RB_BLACK_P(brother->rb_left)
810 1.1 matt && RB_BLACK_P(brother->rb_right)) {
811 1.1 matt if (RB_RED_P(brother)) {
812 1.1 matt /*
813 1.1 matt * Case 1: Our brother is red, swap its
814 1.1 matt * position (and colors) with our parent.
815 1.1 matt * This should now be case 2b (unless C or E
816 1.1 matt * has a red child which is case 3; thus no
817 1.1 matt * explicit branch to case 2b).
818 1.1 matt *
819 1.1 matt * B -> D
820 1.1 matt * A d -> b E
821 1.1 matt * C E -> A C
822 1.1 matt */
823 1.1 matt KASSERT(RB_BLACK_P(parent));
824 1.1 matt rb_tree_reparent_nodes(rbt, parent, other);
825 1.1 matt brother = parent->rb_nodes[other];
826 1.1 matt KASSERT(!RB_SENTINEL_P(brother));
827 1.1 matt KASSERT(RB_RED_P(parent));
828 1.1 matt KASSERT(RB_BLACK_P(brother));
829 1.1 matt KASSERT(rb_tree_check_node(rbt, brother, NULL, false));
830 1.1 matt KASSERT(rb_tree_check_node(rbt, parent, NULL, false));
831 1.1 matt } else {
832 1.1 matt /*
833 1.1 matt * Both our parent and brother are black.
834 1.1 matt * Change our brother to red, advance up rank
835 1.1 matt * and go through the loop again.
836 1.1 matt *
837 1.1 matt * B -> *B
838 1.1 matt * *A D -> A d
839 1.1 matt * C E -> C E
840 1.1 matt */
841 1.1 matt RB_MARK_RED(brother);
842 1.1 matt KASSERT(RB_BLACK_P(brother->rb_left));
843 1.1 matt KASSERT(RB_BLACK_P(brother->rb_right));
844 1.1 matt if (RB_ROOT_P(rbt, parent))
845 1.1 matt return; /* root == parent == black */
846 1.1 matt KASSERT(rb_tree_check_node(rbt, brother, NULL, false));
847 1.1 matt KASSERT(rb_tree_check_node(rbt, parent, NULL, false));
848 1.1 matt which = RB_POSITION(parent);
849 1.1 matt parent = RB_FATHER(parent);
850 1.1 matt continue;
851 1.1 matt }
852 1.1 matt }
853 1.1 matt /*
854 1.1 matt * Avoid an else here so that case 2a above can hit either
855 1.1 matt * case 2b, 3, or 4.
856 1.1 matt */
857 1.1 matt if (RB_RED_P(parent)
858 1.1 matt && RB_BLACK_P(brother)
859 1.1 matt && RB_BLACK_P(brother->rb_left)
860 1.1 matt && RB_BLACK_P(brother->rb_right)) {
861 1.1 matt KASSERT(RB_RED_P(parent));
862 1.1 matt KASSERT(RB_BLACK_P(brother));
863 1.1 matt KASSERT(RB_BLACK_P(brother->rb_left));
864 1.1 matt KASSERT(RB_BLACK_P(brother->rb_right));
865 1.1 matt /*
866 1.1 matt * We are black, our father is red, our brother and
867 1.1 matt * both nephews are black. Simply invert/exchange the
868 1.1 matt * colors of our father and brother (to black and red
869 1.1 matt * respectively).
870 1.1 matt *
871 1.1 matt * | f --> F |
872 1.1 matt * | * B --> * b |
873 1.1 matt * | N N --> N N |
874 1.1 matt */
875 1.1 matt RB_MARK_BLACK(parent);
876 1.1 matt RB_MARK_RED(brother);
877 1.1 matt KASSERT(rb_tree_check_node(rbt, brother, NULL, true));
878 1.1 matt break; /* We're done! */
879 1.1 matt } else {
880 1.1 matt /*
881 1.1 matt * Our brother must be black and have at least one
882 1.1 matt * red child (it may have two).
883 1.1 matt */
884 1.1 matt KASSERT(RB_BLACK_P(brother));
885 1.1 matt KASSERT(RB_RED_P(brother->rb_nodes[which]) ||
886 1.1 matt RB_RED_P(brother->rb_nodes[other]));
887 1.1 matt if (RB_BLACK_P(brother->rb_nodes[other])) {
888 1.1 matt /*
889 1.1 matt * Case 3: our brother is black, our near
890 1.1 matt * nephew is red, and our far nephew is black.
891 1.1 matt * Swap our brother with our near nephew.
892 1.1 matt * This result in a tree that matches case 4.
893 1.1 matt * (Our father could be red or black).
894 1.1 matt *
895 1.1 matt * | F --> F |
896 1.1 matt * | x B --> x B |
897 1.1 matt * | n --> n |
898 1.1 matt */
899 1.1 matt KASSERT(RB_RED_P(brother->rb_nodes[which]));
900 1.1 matt rb_tree_reparent_nodes(rbt, brother, which);
901 1.1 matt KASSERT(RB_FATHER(brother) == parent->rb_nodes[other]);
902 1.1 matt brother = parent->rb_nodes[other];
903 1.1 matt KASSERT(RB_RED_P(brother->rb_nodes[other]));
904 1.1 matt }
905 1.1 matt /*
906 1.1 matt * Case 4: our brother is black and our far nephew
907 1.1 matt * is red. Swap our father and brother locations and
908 1.1 matt * change our far nephew to black. (these can be
909 1.1 matt * done in either order so we change the color first).
910 1.1 matt * The result is a valid red-black tree and is a
911 1.1 matt * terminal case. (again we don't care about the
912 1.1 matt * father's color)
913 1.1 matt *
914 1.1 matt * If the father is red, we will get a red-black-black
915 1.1 matt * tree:
916 1.1 matt * | f -> f --> b |
917 1.1 matt * | B -> B --> F N |
918 1.1 matt * | n -> N --> |
919 1.1 matt *
920 1.1 matt * If the father is black, we will get an all black
921 1.1 matt * tree:
922 1.1 matt * | F -> F --> B |
923 1.1 matt * | B -> B --> F N |
924 1.1 matt * | n -> N --> |
925 1.1 matt *
926 1.1 matt * If we had two red nephews, then after the swap,
927 1.1 matt * our former father would have a red grandson.
928 1.1 matt */
929 1.1 matt KASSERT(RB_BLACK_P(brother));
930 1.1 matt KASSERT(RB_RED_P(brother->rb_nodes[other]));
931 1.1 matt RB_MARK_BLACK(brother->rb_nodes[other]);
932 1.1 matt rb_tree_reparent_nodes(rbt, parent, other);
933 1.1 matt break; /* We're done! */
934 1.1 matt }
935 1.1 matt }
936 1.1 matt KASSERT(rb_tree_check_node(rbt, parent, NULL, true));
937 1.1 matt }
938 1.1 matt
939 1.1 matt struct rb_node *
940 1.1 matt rb_tree_iterate(struct rb_tree *rbt, struct rb_node *self,
941 1.1 matt const unsigned int direction)
942 1.1 matt {
943 1.1 matt const unsigned int other = direction ^ RB_DIR_OTHER;
944 1.1 matt KASSERT(direction == RB_DIR_LEFT || direction == RB_DIR_RIGHT);
945 1.1 matt
946 1.1 matt if (self == NULL) {
947 1.1 matt #ifndef RBSMALL
948 1.1 matt if (RB_SENTINEL_P(rbt->rbt_root))
949 1.1 matt return NULL;
950 1.1 matt return rbt->rbt_minmax[direction];
951 1.1 matt #else
952 1.1 matt self = rbt->rbt_root;
953 1.1 matt if (RB_SENTINEL_P(self))
954 1.1 matt return NULL;
955 1.1 matt while (!RB_SENTINEL_P(self->rb_nodes[other]))
956 1.1 matt self = self->rb_nodes[other];
957 1.1 matt return self;
958 1.1 matt #endif /* !RBSMALL */
959 1.1 matt }
960 1.1 matt KASSERT(!RB_SENTINEL_P(self));
961 1.1 matt /*
962 1.1 matt * We can't go any further in this direction. We proceed up in the
963 1.1 matt * opposite direction until our parent is in direction we want to go.
964 1.1 matt */
965 1.1 matt if (RB_SENTINEL_P(self->rb_nodes[direction])) {
966 1.1 matt while (!RB_ROOT_P(rbt, self)) {
967 1.1 matt if (other == RB_POSITION(self))
968 1.1 matt return RB_FATHER(self);
969 1.1 matt self = RB_FATHER(self);
970 1.1 matt }
971 1.1 matt return NULL;
972 1.1 matt }
973 1.1 matt
974 1.1 matt /*
975 1.1 matt * Advance down one in current direction and go down as far as possible
976 1.1 matt * in the opposite direction.
977 1.1 matt */
978 1.1 matt self = self->rb_nodes[direction];
979 1.1 matt KASSERT(!RB_SENTINEL_P(self));
980 1.1 matt while (!RB_SENTINEL_P(self->rb_nodes[other]))
981 1.1 matt self = self->rb_nodes[other];
982 1.1 matt return self;
983 1.1 matt }
984 1.1 matt
985 1.1 matt #ifdef RBDEBUG
986 1.1 matt static const struct rb_node *
987 1.1 matt rb_tree_iterate_const(const struct rb_tree *rbt, const struct rb_node *self,
988 1.1 matt const unsigned int direction)
989 1.1 matt {
990 1.1 matt const unsigned int other = direction ^ RB_DIR_OTHER;
991 1.1 matt KASSERT(direction == RB_DIR_LEFT || direction == RB_DIR_RIGHT);
992 1.1 matt
993 1.1 matt if (self == NULL) {
994 1.1 matt #ifndef RBSMALL
995 1.1 matt if (RB_SENTINEL_P(rbt->rbt_root))
996 1.1 matt return NULL;
997 1.1 matt return rbt->rbt_minmax[direction];
998 1.1 matt #else
999 1.1 matt self = rbt->rbt_root;
1000 1.1 matt if (RB_SENTINEL_P(self))
1001 1.1 matt return NULL;
1002 1.1 matt while (!RB_SENTINEL_P(self->rb_nodes[other]))
1003 1.1 matt self = self->rb_nodes[other];
1004 1.1 matt return self;
1005 1.1 matt #endif /* !RBSMALL */
1006 1.1 matt }
1007 1.1 matt KASSERT(!RB_SENTINEL_P(self));
1008 1.1 matt /*
1009 1.1 matt * We can't go any further in this direction. We proceed up in the
1010 1.1 matt * opposite direction until our parent is in direction we want to go.
1011 1.1 matt */
1012 1.1 matt if (RB_SENTINEL_P(self->rb_nodes[direction])) {
1013 1.1 matt while (!RB_ROOT_P(rbt, self)) {
1014 1.1 matt if (other == RB_POSITION(self))
1015 1.1 matt return RB_FATHER(self);
1016 1.1 matt self = RB_FATHER(self);
1017 1.1 matt }
1018 1.1 matt return NULL;
1019 1.1 matt }
1020 1.1 matt
1021 1.1 matt /*
1022 1.1 matt * Advance down one in current direction and go down as far as possible
1023 1.1 matt * in the opposite direction.
1024 1.1 matt */
1025 1.1 matt self = self->rb_nodes[direction];
1026 1.1 matt KASSERT(!RB_SENTINEL_P(self));
1027 1.1 matt while (!RB_SENTINEL_P(self->rb_nodes[other]))
1028 1.1 matt self = self->rb_nodes[other];
1029 1.1 matt return self;
1030 1.1 matt }
1031 1.1 matt
1032 1.1 matt static unsigned int
1033 1.1 matt rb_tree_count_black(const struct rb_node *self)
1034 1.1 matt {
1035 1.1 matt unsigned int left, right;
1036 1.1 matt
1037 1.1 matt if (RB_SENTINEL_P(self))
1038 1.1 matt return 0;
1039 1.1 matt
1040 1.1 matt left = rb_tree_count_black(self->rb_left);
1041 1.1 matt right = rb_tree_count_black(self->rb_right);
1042 1.1 matt
1043 1.1 matt KASSERT(left == right);
1044 1.1 matt
1045 1.1 matt return left + RB_BLACK_P(self);
1046 1.1 matt }
1047 1.2 matt
1048 1.1 matt static bool
1049 1.1 matt rb_tree_check_node(const struct rb_tree *rbt, const struct rb_node *self,
1050 1.1 matt const struct rb_node *prev, bool red_check)
1051 1.1 matt {
1052 1.1 matt rbto_compare_nodes_fn compare_nodes = rbt->rbt_ops->rbto_compare_nodes;
1053 1.1 matt
1054 1.1 matt KASSERT(!RB_SENTINEL_P(self));
1055 1.1 matt KASSERT(prev == NULL || (*compare_nodes)(prev, self) > 0);
1056 1.1 matt
1057 1.1 matt /*
1058 1.1 matt * Verify our relationship to our parent.
1059 1.1 matt */
1060 1.1 matt if (RB_ROOT_P(rbt, self)) {
1061 1.1 matt KASSERT(self == rbt->rbt_root);
1062 1.1 matt KASSERT(RB_POSITION(self) == RB_DIR_LEFT);
1063 1.1 matt KASSERT(RB_FATHER(self)->rb_nodes[RB_DIR_LEFT] == self);
1064 1.1 matt KASSERT(RB_FATHER(self) == (const struct rb_node *) &rbt->rbt_root);
1065 1.1 matt } else {
1066 1.1 matt KASSERT(self != rbt->rbt_root);
1067 1.1 matt KASSERT(!RB_FATHER_SENTINEL_P(self));
1068 1.1 matt if (RB_POSITION(self) == RB_DIR_LEFT) {
1069 1.1 matt KASSERT((*compare_nodes)(self, RB_FATHER(self)) > 0);
1070 1.1 matt KASSERT(RB_FATHER(self)->rb_nodes[RB_DIR_LEFT] == self);
1071 1.1 matt } else {
1072 1.1 matt KASSERT((*compare_nodes)(self, RB_FATHER(self)) < 0);
1073 1.1 matt KASSERT(RB_FATHER(self)->rb_nodes[RB_DIR_RIGHT] == self);
1074 1.1 matt }
1075 1.1 matt }
1076 1.1 matt
1077 1.1 matt /*
1078 1.1 matt * Verify our position in the linked list against the tree itself.
1079 1.1 matt */
1080 1.1 matt {
1081 1.1 matt const struct rb_node *prev0 = rb_tree_iterate_const(rbt, self, RB_DIR_LEFT);
1082 1.1 matt const struct rb_node *next0 = rb_tree_iterate_const(rbt, self, RB_DIR_RIGHT);
1083 1.1 matt KASSERT(prev0 == TAILQ_PREV(self, rb_node_qh, rb_link));
1084 1.1 matt KASSERT(next0 == TAILQ_NEXT(self, rb_link));
1085 1.1 matt #ifndef RBSMALL
1086 1.1 matt KASSERT(prev0 != NULL || self == rbt->rbt_minmax[RB_DIR_LEFT]);
1087 1.1 matt KASSERT(next0 != NULL || self == rbt->rbt_minmax[RB_DIR_RIGHT]);
1088 1.1 matt #endif
1089 1.1 matt }
1090 1.1 matt
1091 1.1 matt /*
1092 1.1 matt * The root must be black.
1093 1.1 matt * There can never be two adjacent red nodes.
1094 1.1 matt */
1095 1.1 matt if (red_check) {
1096 1.1 matt KASSERT(!RB_ROOT_P(rbt, self) || RB_BLACK_P(self));
1097 1.1 matt (void) rb_tree_count_black(self);
1098 1.1 matt if (RB_RED_P(self)) {
1099 1.1 matt const struct rb_node *brother;
1100 1.1 matt KASSERT(!RB_ROOT_P(rbt, self));
1101 1.1 matt brother = RB_FATHER(self)->rb_nodes[RB_POSITION(self) ^ RB_DIR_OTHER];
1102 1.1 matt KASSERT(RB_BLACK_P(RB_FATHER(self)));
1103 1.1 matt /*
1104 1.1 matt * I'm red and have no children, then I must either
1105 1.1 matt * have no brother or my brother also be red and
1106 1.1 matt * also have no children. (black count == 0)
1107 1.1 matt */
1108 1.1 matt KASSERT(!RB_CHILDLESS_P(self)
1109 1.1 matt || RB_SENTINEL_P(brother)
1110 1.1 matt || RB_RED_P(brother)
1111 1.1 matt || RB_CHILDLESS_P(brother));
1112 1.1 matt /*
1113 1.1 matt * If I'm not childless, I must have two children
1114 1.1 matt * and they must be both be black.
1115 1.1 matt */
1116 1.1 matt KASSERT(RB_CHILDLESS_P(self)
1117 1.1 matt || (RB_TWOCHILDREN_P(self)
1118 1.1 matt && RB_BLACK_P(self->rb_left)
1119 1.1 matt && RB_BLACK_P(self->rb_right)));
1120 1.1 matt /*
1121 1.1 matt * If I'm not childless, thus I have black children,
1122 1.1 matt * then my brother must either be black or have two
1123 1.1 matt * black children.
1124 1.1 matt */
1125 1.1 matt KASSERT(RB_CHILDLESS_P(self)
1126 1.1 matt || RB_BLACK_P(brother)
1127 1.1 matt || (RB_TWOCHILDREN_P(brother)
1128 1.1 matt && RB_BLACK_P(brother->rb_left)
1129 1.1 matt && RB_BLACK_P(brother->rb_right)));
1130 1.1 matt } else {
1131 1.1 matt /*
1132 1.1 matt * If I'm black and have one child, that child must
1133 1.1 matt * be red and childless.
1134 1.1 matt */
1135 1.1 matt KASSERT(RB_CHILDLESS_P(self)
1136 1.1 matt || RB_TWOCHILDREN_P(self)
1137 1.1 matt || (!RB_LEFT_SENTINEL_P(self)
1138 1.1 matt && RB_RIGHT_SENTINEL_P(self)
1139 1.1 matt && RB_RED_P(self->rb_left)
1140 1.1 matt && RB_CHILDLESS_P(self->rb_left))
1141 1.1 matt || (!RB_RIGHT_SENTINEL_P(self)
1142 1.1 matt && RB_LEFT_SENTINEL_P(self)
1143 1.1 matt && RB_RED_P(self->rb_right)
1144 1.1 matt && RB_CHILDLESS_P(self->rb_right)));
1145 1.1 matt
1146 1.1 matt /*
1147 1.1 matt * If I'm a childless black node and my parent is
1148 1.1 matt * black, my 2nd closet relative away from my parent
1149 1.1 matt * is either red or has a red parent or red children.
1150 1.1 matt */
1151 1.1 matt if (!RB_ROOT_P(rbt, self)
1152 1.1 matt && RB_CHILDLESS_P(self)
1153 1.1 matt && RB_BLACK_P(RB_FATHER(self))) {
1154 1.1 matt const unsigned int which = RB_POSITION(self);
1155 1.1 matt const unsigned int other = which ^ RB_DIR_OTHER;
1156 1.1 matt const struct rb_node *relative0, *relative;
1157 1.1 matt
1158 1.1 matt relative0 = rb_tree_iterate_const(rbt,
1159 1.1 matt self, other);
1160 1.1 matt KASSERT(relative0 != NULL);
1161 1.1 matt relative = rb_tree_iterate_const(rbt,
1162 1.1 matt relative0, other);
1163 1.1 matt KASSERT(relative != NULL);
1164 1.1 matt KASSERT(RB_SENTINEL_P(relative->rb_nodes[which]));
1165 1.1 matt #if 0
1166 1.1 matt KASSERT(RB_RED_P(relative)
1167 1.1 matt || RB_RED_P(relative->rb_left)
1168 1.1 matt || RB_RED_P(relative->rb_right)
1169 1.1 matt || RB_RED_P(RB_FATHER(relative)));
1170 1.1 matt #endif
1171 1.1 matt }
1172 1.1 matt }
1173 1.1 matt /*
1174 1.1 matt * A grandparent's children must be real nodes and not
1175 1.1 matt * sentinels. First check out grandparent.
1176 1.1 matt */
1177 1.1 matt KASSERT(RB_ROOT_P(rbt, self)
1178 1.1 matt || RB_ROOT_P(rbt, RB_FATHER(self))
1179 1.1 matt || RB_TWOCHILDREN_P(RB_FATHER(RB_FATHER(self))));
1180 1.1 matt /*
1181 1.1 matt * If we are have grandchildren on our left, then
1182 1.1 matt * we must have a child on our right.
1183 1.1 matt */
1184 1.1 matt KASSERT(RB_LEFT_SENTINEL_P(self)
1185 1.1 matt || RB_CHILDLESS_P(self->rb_left)
1186 1.1 matt || !RB_RIGHT_SENTINEL_P(self));
1187 1.1 matt /*
1188 1.1 matt * If we are have grandchildren on our right, then
1189 1.1 matt * we must have a child on our left.
1190 1.1 matt */
1191 1.1 matt KASSERT(RB_RIGHT_SENTINEL_P(self)
1192 1.1 matt || RB_CHILDLESS_P(self->rb_right)
1193 1.1 matt || !RB_LEFT_SENTINEL_P(self));
1194 1.1 matt
1195 1.1 matt /*
1196 1.1 matt * If we have a child on the left and it doesn't have two
1197 1.1 matt * children make sure we don't have great-great-grandchildren on
1198 1.1 matt * the right.
1199 1.1 matt */
1200 1.1 matt KASSERT(RB_TWOCHILDREN_P(self->rb_left)
1201 1.1 matt || RB_CHILDLESS_P(self->rb_right)
1202 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_left)
1203 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_left->rb_left)
1204 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_left->rb_right)
1205 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_right)
1206 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_right->rb_left)
1207 1.1 matt || RB_CHILDLESS_P(self->rb_right->rb_right->rb_right));
1208 1.1 matt
1209 1.1 matt /*
1210 1.1 matt * If we have a child on the right and it doesn't have two
1211 1.1 matt * children make sure we don't have great-great-grandchildren on
1212 1.1 matt * the left.
1213 1.1 matt */
1214 1.1 matt KASSERT(RB_TWOCHILDREN_P(self->rb_right)
1215 1.1 matt || RB_CHILDLESS_P(self->rb_left)
1216 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_left)
1217 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_left->rb_left)
1218 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_left->rb_right)
1219 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_right)
1220 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_right->rb_left)
1221 1.1 matt || RB_CHILDLESS_P(self->rb_left->rb_right->rb_right));
1222 1.1 matt
1223 1.1 matt /*
1224 1.1 matt * If we are fully interior node, then our predecessors and
1225 1.1 matt * successors must have no children in our direction.
1226 1.1 matt */
1227 1.1 matt if (RB_TWOCHILDREN_P(self)) {
1228 1.1 matt const struct rb_node *prev0;
1229 1.1 matt const struct rb_node *next0;
1230 1.1 matt
1231 1.1 matt prev0 = rb_tree_iterate_const(rbt, self, RB_DIR_LEFT);
1232 1.1 matt KASSERT(prev0 != NULL);
1233 1.1 matt KASSERT(RB_RIGHT_SENTINEL_P(prev0));
1234 1.1 matt
1235 1.1 matt next0 = rb_tree_iterate_const(rbt, self, RB_DIR_RIGHT);
1236 1.1 matt KASSERT(next0 != NULL);
1237 1.1 matt KASSERT(RB_LEFT_SENTINEL_P(next0));
1238 1.1 matt }
1239 1.1 matt }
1240 1.1 matt
1241 1.1 matt return true;
1242 1.1 matt }
1243 1.1 matt
1244 1.1 matt void
1245 1.1 matt rb_tree_check(const struct rb_tree *rbt, bool red_check)
1246 1.1 matt {
1247 1.1 matt const struct rb_node *self;
1248 1.1 matt const struct rb_node *prev;
1249 1.1 matt #ifdef RBSTATS
1250 1.1 matt unsigned int count = 0;
1251 1.1 matt #endif
1252 1.1 matt
1253 1.1 matt KASSERT(rbt->rbt_root != NULL);
1254 1.1 matt KASSERT(RB_LEFT_P(rbt->rbt_root));
1255 1.1 matt
1256 1.1 matt #if defined(RBSTATS) && !defined(RBSMALL)
1257 1.1 matt KASSERT(rbt->rbt_count > 1
1258 1.1 matt || rbt->rbt_minmax[RB_DIR_LEFT] == rbt->rbt_minmax[RB_DIR_RIGHT]);
1259 1.1 matt #endif
1260 1.1 matt
1261 1.1 matt prev = NULL;
1262 1.1 matt TAILQ_FOREACH(self, &rbt->rbt_nodes, rb_link) {
1263 1.1 matt rb_tree_check_node(rbt, self, prev, false);
1264 1.1 matt #ifdef RBSTATS
1265 1.1 matt count++;
1266 1.1 matt #endif
1267 1.1 matt }
1268 1.1 matt #ifdef RBSTATS
1269 1.1 matt KASSERT(rbt->rbt_count == count);
1270 1.1 matt #endif
1271 1.1 matt if (red_check) {
1272 1.1 matt KASSERT(RB_BLACK_P(rbt->rbt_root));
1273 1.1 matt KASSERT(RB_SENTINEL_P(rbt->rbt_root)
1274 1.1 matt || rb_tree_count_black(rbt->rbt_root));
1275 1.1 matt
1276 1.1 matt /*
1277 1.1 matt * The root must be black.
1278 1.1 matt * There can never be two adjacent red nodes.
1279 1.1 matt */
1280 1.1 matt TAILQ_FOREACH(self, &rbt->rbt_nodes, rb_link) {
1281 1.1 matt rb_tree_check_node(rbt, self, NULL, true);
1282 1.1 matt }
1283 1.1 matt }
1284 1.1 matt }
1285 1.1 matt #endif /* RBDEBUG */
1286 1.1 matt
1287 1.1 matt #ifdef RBSTATS
1288 1.1 matt static void
1289 1.1 matt rb_tree_mark_depth(const struct rb_tree *rbt, const struct rb_node *self,
1290 1.1 matt size_t *depths, size_t depth)
1291 1.1 matt {
1292 1.1 matt if (RB_SENTINEL_P(self))
1293 1.1 matt return;
1294 1.1 matt
1295 1.1 matt if (RB_TWOCHILDREN_P(self)) {
1296 1.1 matt rb_tree_mark_depth(rbt, self->rb_left, depths, depth + 1);
1297 1.1 matt rb_tree_mark_depth(rbt, self->rb_right, depths, depth + 1);
1298 1.1 matt return;
1299 1.1 matt }
1300 1.1 matt depths[depth]++;
1301 1.1 matt if (!RB_LEFT_SENTINEL_P(self)) {
1302 1.1 matt rb_tree_mark_depth(rbt, self->rb_left, depths, depth + 1);
1303 1.1 matt }
1304 1.1 matt if (!RB_RIGHT_SENTINEL_P(self)) {
1305 1.1 matt rb_tree_mark_depth(rbt, self->rb_right, depths, depth + 1);
1306 1.1 matt }
1307 1.1 matt }
1308 1.1 matt
1309 1.1 matt void
1310 rb_tree_depths(const struct rb_tree *rbt, size_t *depths)
1311 {
1312 rb_tree_mark_depth(rbt, rbt->rbt_root, depths, 1);
1313 }
1314 #endif /* RBSTATS */
1315