GRASS GIS 7 Programmer's Manual  7.9.dev(2021)-e5379bbd7
avl.c
Go to the documentation of this file.
1 /* Produced by texiweb from libavl.w. */
2 
3 /* libavl - library for manipulation of binary trees.
4  Copyright (C) 1998, 1999, 2000, 2001, 2002, 2004 Free Software
5  Foundation, Inc.
6 
7  This library is free software; you can redistribute it and/or
8  modify it under the terms of the GNU Lesser General Public
9  License as published by the Free Software Foundation; either
10  version 3 of the License, or (at your option) any later version.
11 
12  This library is distributed in the hope that it will be useful,
13  but WITHOUT ANY WARRANTY; without even the implied warranty of
14  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15  Lesser General Public License for more details.
16 
17  You should have received a copy of the GNU Lesser General Public
18  License along with this library; if not, write to the Free Software
19  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
20  02110-1301 USA.
21  */
22 
23 /* Nov 2016, Markus Metz
24  * from libavl-2.0.3
25  * added safety checks and speed optimizations
26  */
27 
28 #include <assert.h>
29 #include <stdio.h>
30 #include <stdlib.h>
31 #include <string.h>
32 #include "avl.h"
33 
34 /* Creates and returns a new table
35  with comparison function |compare| using parameter |param|
36  and memory allocator |allocator|.
37  Returns |NULL| if memory allocation failed. */
39  struct libavl_allocator *allocator)
40 {
41  struct avl_table *tree;
42 
43  assert(compare != NULL);
44 
45  if (allocator == NULL)
46  allocator = &avl_allocator_default;
47 
48  tree = allocator->libavl_malloc(allocator, sizeof *tree);
49  if (tree == NULL)
50  return NULL;
51 
52  tree->avl_root = NULL;
53  tree->avl_compare = compare;
54  tree->avl_param = param;
55  tree->avl_alloc = allocator;
56  tree->avl_count = 0;
57  tree->avl_generation = 0;
58 
59  return tree;
60 }
61 
62 /* Search |tree| for an item matching |item|, and return it if found.
63  Otherwise return |NULL|. */
64 void *avl_find(const struct avl_table *tree, const void *item)
65 {
66  const struct avl_node *p;
67 
68  assert(tree != NULL && item != NULL);
69 
70  p = tree->avl_root;
71  while (p != NULL) {
72  int cmp = tree->avl_compare(item, p->avl_data, tree->avl_param);
73 
74  if (cmp == 0)
75  return p->avl_data;
76 
77  p = p->avl_link[cmp > 0];
78  }
79 
80  return NULL;
81 }
82 
83 /* Inserts |item| into |tree| and returns a pointer to |item|'s address.
84  If a duplicate item is found in the tree,
85  returns a pointer to the duplicate without inserting |item|.
86  Returns |NULL| in case of memory allocation failure. */
87 void **avl_probe(struct avl_table *tree, void *item)
88 {
89  struct avl_node *y, *z; /* Top node to update balance factor, and parent. */
90  struct avl_node *p, *q; /* Iterator, and parent. */
91  struct avl_node *n; /* Newly inserted node. */
92  struct avl_node *w; /* New root of rebalanced subtree. */
93  int dir; /* Direction to descend. */
94 
95  unsigned char da[AVL_MAX_HEIGHT]; /* Cached comparison results. */
96  int k = 0; /* Number of cached results. */
97 
98  assert(tree != NULL && item != NULL);
99 
100  z = (struct avl_node *)&tree->avl_root;
101  y = tree->avl_root;
102  dir = 0;
103  q = z, p = y;
104  while (p != NULL) {
105  int cmp = tree->avl_compare(item, p->avl_data, tree->avl_param);
106 
107  if (cmp == 0)
108  return &p->avl_data;
109 
110  if (p->avl_balance != 0)
111  z = q, y = p, k = 0;
112  da[k++] = dir = cmp > 0;
113  q = p, p = p->avl_link[dir];
114  }
115 
116  n = tree->avl_alloc->libavl_malloc(tree->avl_alloc, sizeof *n);
117  if (n == NULL)
118  return NULL;
119 
120  tree->avl_count++;
121  tree->avl_generation++;
122  n->avl_link[0] = n->avl_link[1] = NULL;
123  n->avl_data = item;
124  n->avl_balance = 0;
125  if (y == NULL) {
126  tree->avl_root = n;
127 
128  return &n->avl_data;
129  }
130  q->avl_link[dir] = n;
131 
132  p = y, k = 0;
133  while (p != n) {
134  if (da[k] == 0)
135  p->avl_balance--;
136  else
137  p->avl_balance++;
138  p = p->avl_link[da[k]], k++;
139  }
140 
141  if (y->avl_balance == -2) {
142  struct avl_node *x = y->avl_link[0];
143 
144  if (x->avl_balance == -1) {
145  w = x;
146  y->avl_link[0] = x->avl_link[1];
147  x->avl_link[1] = y;
148  x->avl_balance = y->avl_balance = 0;
149  }
150  else {
151  assert(x->avl_balance == +1);
152  w = x->avl_link[1];
153  x->avl_link[1] = w->avl_link[0];
154  w->avl_link[0] = x;
155  y->avl_link[0] = w->avl_link[1];
156  w->avl_link[1] = y;
157  if (w->avl_balance == -1)
158  x->avl_balance = 0, y->avl_balance = +1;
159  else if (w->avl_balance == 0)
160  x->avl_balance = y->avl_balance = 0;
161  else /* |w->avl_balance == +1| */
162  x->avl_balance = -1, y->avl_balance = 0;
163  w->avl_balance = 0;
164  }
165  }
166  else if (y->avl_balance == +2) {
167  struct avl_node *x = y->avl_link[1];
168 
169  if (x->avl_balance == +1) {
170  w = x;
171  y->avl_link[1] = x->avl_link[0];
172  x->avl_link[0] = y;
173  x->avl_balance = y->avl_balance = 0;
174  }
175  else {
176  assert(x->avl_balance == -1);
177  w = x->avl_link[0];
178  x->avl_link[0] = w->avl_link[1];
179  w->avl_link[1] = x;
180  y->avl_link[1] = w->avl_link[0];
181  w->avl_link[0] = y;
182  if (w->avl_balance == +1)
183  x->avl_balance = 0, y->avl_balance = -1;
184  else if (w->avl_balance == 0)
185  x->avl_balance = y->avl_balance = 0;
186  else /* |w->avl_balance == -1| */
187  x->avl_balance = +1, y->avl_balance = 0;
188  w->avl_balance = 0;
189  }
190  }
191  else
192  return &n->avl_data;
193  z->avl_link[y != z->avl_link[0]] = w;
194 
195  return &n->avl_data;
196 }
197 
198 /* Inserts |item| into |table|.
199  Returns |NULL| if |item| was successfully inserted
200  or if a memory allocation error occurred.
201  Otherwise, returns the duplicate item. */
202 void *avl_insert(struct avl_table *table, void *item)
203 {
204  void **p = avl_probe(table, item);
205 
206  return p == NULL || *p == item ? NULL : *p;
207 }
208 
209 /* Inserts |item| into |table|, replacing any duplicate item.
210  Returns |NULL| if |item| was inserted without replacing a duplicate,
211  or if a memory allocation error occurred.
212  Otherwise, returns the item that was replaced. */
213 void *avl_replace(struct avl_table *table, void *item)
214 {
215  void **p = avl_probe(table, item);
216 
217  if (p == NULL || *p == item)
218  return NULL;
219  else {
220  void *r = *p;
221 
222  *p = item;
223 
224  return r;
225  }
226 }
227 
228 /* Deletes from |tree| and returns an item matching |item|.
229  Returns a null pointer if no matching item found. */
230 void *avl_delete(struct avl_table *tree, const void *item)
231 {
232  /* Stack of nodes. */
233  struct avl_node *pa[AVL_MAX_HEIGHT]; /* Nodes. */
234  unsigned char da[AVL_MAX_HEIGHT]; /* |avl_link[]| indexes. */
235  int k; /* Stack pointer. */
236 
237  struct avl_node *p; /* Traverses tree to find node to delete. */
238  int dir; /* Side of |pa[k]| on which |p| is linked. */
239  int cmp; /* Result of comparison between |item| and |p|. */
240 
241  assert(tree != NULL && item != NULL);
242 
243  pa[0] = (struct avl_node *)&tree->avl_root;
244  da[0] = 0;
245  k = 1;
246  p = tree->avl_root;
247  while (p != NULL) {
248  cmp = tree->avl_compare(item, p->avl_data, tree->avl_param);
249 
250  if (cmp == 0)
251  break;
252 
253  dir = cmp > 0;
254 
255  pa[k] = p;
256  da[k++] = dir;
257 
258  p = p->avl_link[dir];
259  }
260  if (p == NULL)
261  return NULL;
262 
263  item = p->avl_data;
264 
265  if (p->avl_link[1] == NULL)
266  pa[k - 1]->avl_link[da[k - 1]] = p->avl_link[0];
267  else {
268  struct avl_node *r = p->avl_link[1];
269 
270  if (r->avl_link[0] == NULL) {
271  r->avl_link[0] = p->avl_link[0];
272  r->avl_balance = p->avl_balance;
273  pa[k - 1]->avl_link[da[k - 1]] = r;
274  da[k] = 1;
275  pa[k++] = r;
276  }
277  else {
278  struct avl_node *s;
279  int j = k++;
280 
281  while (r != NULL) {
282  da[k] = 0;
283  pa[k++] = r;
284  s = r->avl_link[0];
285  if (s->avl_link[0] == NULL)
286  break;
287 
288  r = s;
289  }
290 
291  s->avl_link[0] = p->avl_link[0];
292  r->avl_link[0] = s->avl_link[1];
293  s->avl_link[1] = p->avl_link[1];
294  s->avl_balance = p->avl_balance;
295 
296  pa[j - 1]->avl_link[da[j - 1]] = s;
297  da[j] = 1;
298  pa[j] = s;
299  }
300  }
301 
302  tree->avl_alloc->libavl_free(tree->avl_alloc, p);
303 
304  assert(k > 0);
305  while (--k > 0) {
306  struct avl_node *y = pa[k];
307 
308  if (da[k] == 0) {
309  y->avl_balance++;
310  if (y->avl_balance == +1)
311  break;
312  else if (y->avl_balance == +2) {
313  struct avl_node *x = y->avl_link[1];
314 
315  if (x->avl_balance == -1) {
316  struct avl_node *w;
317 
318  assert(x->avl_balance == -1);
319  w = x->avl_link[0];
320  x->avl_link[0] = w->avl_link[1];
321  w->avl_link[1] = x;
322  y->avl_link[1] = w->avl_link[0];
323  w->avl_link[0] = y;
324  if (w->avl_balance == +1)
325  x->avl_balance = 0, y->avl_balance = -1;
326  else if (w->avl_balance == 0)
327  x->avl_balance = y->avl_balance = 0;
328  else /* |w->avl_balance == -1| */
329  x->avl_balance = +1, y->avl_balance = 0;
330  w->avl_balance = 0;
331  pa[k - 1]->avl_link[da[k - 1]] = w;
332  }
333  else {
334  y->avl_link[1] = x->avl_link[0];
335  x->avl_link[0] = y;
336  pa[k - 1]->avl_link[da[k - 1]] = x;
337  if (x->avl_balance == 0) {
338  x->avl_balance = -1;
339  y->avl_balance = +1;
340  break;
341  }
342  else
343  x->avl_balance = y->avl_balance = 0;
344  }
345  }
346  }
347  else {
348  y->avl_balance--;
349  if (y->avl_balance == -1)
350  break;
351  else if (y->avl_balance == -2) {
352  struct avl_node *x = y->avl_link[0];
353 
354  if (x->avl_balance == +1) {
355  struct avl_node *w;
356 
357  assert(x->avl_balance == +1);
358  w = x->avl_link[1];
359  x->avl_link[1] = w->avl_link[0];
360  w->avl_link[0] = x;
361  y->avl_link[0] = w->avl_link[1];
362  w->avl_link[1] = y;
363  if (w->avl_balance == -1)
364  x->avl_balance = 0, y->avl_balance = +1;
365  else if (w->avl_balance == 0)
366  x->avl_balance = y->avl_balance = 0;
367  else /* |w->avl_balance == +1| */
368  x->avl_balance = -1, y->avl_balance = 0;
369  w->avl_balance = 0;
370  pa[k - 1]->avl_link[da[k - 1]] = w;
371  }
372  else {
373  y->avl_link[0] = x->avl_link[1];
374  x->avl_link[1] = y;
375  pa[k - 1]->avl_link[da[k - 1]] = x;
376  if (x->avl_balance == 0) {
377  x->avl_balance = +1;
378  y->avl_balance = -1;
379  break;
380  }
381  else
382  x->avl_balance = y->avl_balance = 0;
383  }
384  }
385  }
386  }
387 
388  tree->avl_count--;
389  tree->avl_generation++;
390 
391  return (void *)item;
392 }
393 
394 /* Refreshes the stack of parent pointers in |trav|
395  and updates its generation number. */
396 static void trav_refresh(struct avl_traverser *trav)
397 {
398  assert(trav != NULL);
399 
400  trav->avl_generation = trav->avl_table->avl_generation;
401 
402  if (trav->avl_node != NULL) {
404  void *param = trav->avl_table->avl_param;
405  struct avl_node *node = trav->avl_node;
406  struct avl_node *i;
407 
408  trav->avl_height = 0;
409  for (i = trav->avl_table->avl_root; i != node;) {
411  assert(i != NULL);
412 
413  trav->avl_stack[trav->avl_height++] = i;
414  i = i->avl_link[cmp(node->avl_data, i->avl_data, param) > 0];
415  }
416  }
417 }
418 
419 /* Initializes |trav| for use with |tree|
420  and selects the null node. */
421 void avl_t_init(struct avl_traverser *trav, struct avl_table *tree)
422 {
423  trav->avl_table = tree;
424  trav->avl_node = NULL;
425  trav->avl_height = 0;
426  trav->avl_generation = tree->avl_generation;
427 }
428 
429 /* Initializes |trav| for |tree|
430  and selects and returns a pointer to its least-valued item.
431  Returns |NULL| if |tree| contains no nodes. */
432 void *avl_t_first(struct avl_traverser *trav, struct avl_table *tree)
433 {
434  struct avl_node *x;
435 
436  assert(tree != NULL && trav != NULL);
437 
438  trav->avl_table = tree;
439  trav->avl_height = 0;
440  trav->avl_generation = tree->avl_generation;
441 
442  x = tree->avl_root;
443  if (x != NULL)
444  while (x->avl_link[0] != NULL) {
446  trav->avl_stack[trav->avl_height++] = x;
447  x = x->avl_link[0];
448  }
449  trav->avl_node = x;
450 
451  return x != NULL ? x->avl_data : NULL;
452 }
453 
454 /* Initializes |trav| for |tree|
455  and selects and returns a pointer to its greatest-valued item.
456  Returns |NULL| if |tree| contains no nodes. */
457 void *avl_t_last(struct avl_traverser *trav, struct avl_table *tree)
458 {
459  struct avl_node *x;
460 
461  assert(tree != NULL && trav != NULL);
462 
463  trav->avl_table = tree;
464  trav->avl_height = 0;
465  trav->avl_generation = tree->avl_generation;
466 
467  x = tree->avl_root;
468  if (x != NULL)
469  while (x->avl_link[1] != NULL) {
471  trav->avl_stack[trav->avl_height++] = x;
472  x = x->avl_link[1];
473  }
474  trav->avl_node = x;
475 
476  return x != NULL ? x->avl_data : NULL;
477 }
478 
479 /* Searches for |item| in |tree|.
480  If found, initializes |trav| to the item found and returns the item
481  as well.
482  If there is no matching item, initializes |trav| to the null item
483  and returns |NULL|. */
484 void *avl_t_find(struct avl_traverser *trav, struct avl_table *tree,
485  void *item)
486 {
487  struct avl_node *p;
488 
489  assert(trav != NULL && tree != NULL && item != NULL);
490  trav->avl_table = tree;
491  trav->avl_height = 0;
492  trav->avl_generation = tree->avl_generation;
493  p = tree->avl_root;
494  while (p != NULL) {
495  int cmp = tree->avl_compare(item, p->avl_data, tree->avl_param);
496 
497  if (cmp == 0) {
498  trav->avl_node = p;
499 
500  return p->avl_data;
501  }
502 
504  trav->avl_stack[trav->avl_height++] = p;
505  p = p->avl_link[cmp > 0];
506  }
507 
508  trav->avl_height = 0;
509  trav->avl_node = NULL;
510 
511  return NULL;
512 }
513 
514 /* Attempts to insert |item| into |tree|.
515  If |item| is inserted successfully, it is returned and |trav| is
516  initialized to its location.
517  If a duplicate is found, it is returned and |trav| is initialized to
518  its location. No replacement of the item occurs.
519  If a memory allocation failure occurs, |NULL| is returned and |trav|
520  is initialized to the null item. */
521 void *avl_t_insert(struct avl_traverser *trav, struct avl_table *tree,
522  void *item)
523 {
524  void **p;
525 
526  assert(trav != NULL && tree != NULL && item != NULL);
527 
528  p = avl_probe(tree, item);
529  if (p != NULL) {
530  trav->avl_table = tree;
531  trav->avl_node = ((struct avl_node *)
532  ((char *)p - offsetof(struct avl_node, avl_data)));
533 
534  trav->avl_generation = tree->avl_generation - 1;
535 
536  return *p;
537  }
538  else {
539  avl_t_init(trav, tree);
540 
541  return NULL;
542  }
543 }
544 
545 /* Initializes |trav| to have the same current node as |src|. */
546 void *avl_t_copy(struct avl_traverser *trav, const struct avl_traverser *src)
547 {
548  assert(trav != NULL && src != NULL);
549 
550  if (trav != src) {
551  trav->avl_table = src->avl_table;
552  trav->avl_node = src->avl_node;
553  trav->avl_generation = src->avl_generation;
554  if (trav->avl_generation == trav->avl_table->avl_generation) {
555  trav->avl_height = src->avl_height;
556  memcpy(trav->avl_stack, (const void *)src->avl_stack,
557  sizeof *trav->avl_stack * trav->avl_height);
558  }
559  }
560 
561  return trav->avl_node != NULL ? trav->avl_node->avl_data : NULL;
562 }
563 
564 /* Returns the next data item in inorder
565  within the tree being traversed with |trav|,
566  or if there are no more data items returns |NULL|. */
567 void *avl_t_next(struct avl_traverser *trav)
568 {
569  struct avl_node *x;
570 
571  assert(trav != NULL);
572 
573  if (trav->avl_generation != trav->avl_table->avl_generation)
574  trav_refresh(trav);
575 
576  x = trav->avl_node;
577  if (x == NULL) {
578  return avl_t_first(trav, trav->avl_table);
579  }
580  else if (x->avl_link[1] != NULL) {
582  trav->avl_stack[trav->avl_height++] = x;
583  x = x->avl_link[1];
584 
585  while (x->avl_link[0] != NULL) {
587  trav->avl_stack[trav->avl_height++] = x;
588  x = x->avl_link[0];
589  }
590  }
591  else {
592  struct avl_node *y;
593 
594  do {
595  if (trav->avl_height == 0) {
596  trav->avl_node = NULL;
597  return NULL;
598  }
599 
600  y = x;
601  x = trav->avl_stack[--trav->avl_height];
602  }
603  while (y == x->avl_link[1]);
604  }
605  trav->avl_node = x;
606 
607  return x->avl_data;
608 }
609 
610 /* Returns the previous data item in inorder
611  within the tree being traversed with |trav|,
612  or if there are no more data items returns |NULL|. */
613 void *avl_t_prev(struct avl_traverser *trav)
614 {
615  struct avl_node *x;
616 
617  assert(trav != NULL);
618 
619  if (trav->avl_generation != trav->avl_table->avl_generation)
620  trav_refresh(trav);
621 
622  x = trav->avl_node;
623  if (x == NULL) {
624  return avl_t_last(trav, trav->avl_table);
625  }
626  else if (x->avl_link[0] != NULL) {
628  trav->avl_stack[trav->avl_height++] = x;
629  x = x->avl_link[0];
630 
631  while (x->avl_link[1] != NULL) {
633  trav->avl_stack[trav->avl_height++] = x;
634  x = x->avl_link[1];
635  }
636  }
637  else {
638  struct avl_node *y;
639 
640  do {
641  if (trav->avl_height == 0) {
642  trav->avl_node = NULL;
643  return NULL;
644  }
645 
646  y = x;
647  x = trav->avl_stack[--trav->avl_height];
648  }
649  while (y == x->avl_link[0]);
650  }
651  trav->avl_node = x;
652 
653  return x->avl_data;
654 }
655 
656 /* Returns |trav|'s current item. */
657 void *avl_t_cur(struct avl_traverser *trav)
658 {
659  assert(trav != NULL);
660 
661  return trav->avl_node != NULL ? trav->avl_node->avl_data : NULL;
662 }
663 
664 /* Replaces the current item in |trav| by |new| and returns the item replaced.
665  |trav| must not have the null item selected.
666  The new item must not upset the ordering of the tree. */
667 void *avl_t_replace(struct avl_traverser *trav, void *new)
668 {
669  void *old;
670 
671  assert(trav != NULL && trav->avl_node != NULL && new != NULL);
672  old = trav->avl_node->avl_data;
673  trav->avl_node->avl_data = new;
674 
675  return old;
676 }
677 
678 /* Destroys |new| with |avl_destroy (new, destroy)|,
679  first setting right links of nodes in |stack| within |new|
680  to null pointers to avoid touching uninitialized data. */
681 static void
682 copy_error_recovery(struct avl_node **stack, int height,
683  struct avl_table *new, avl_item_func * destroy)
684 {
685  assert(stack != NULL && height >= 0 && new != NULL);
686 
687  for (; height > 2; height -= 2)
688  stack[height - 1]->avl_link[1] = NULL;
689  avl_destroy(new, destroy);
690 }
691 
692 /* Copies |org| to a newly created tree, which is returned.
693  If |copy != NULL|, each data item in |org| is first passed to |copy|,
694  and the return values are inserted into the tree,
695  with |NULL| return values taken as indications of failure.
696  On failure, destroys the partially created new tree,
697  applying |destroy|, if non-null, to each item in the new tree so far,
698  and returns |NULL|.
699  If |allocator != NULL|, it is used for allocation in the new tree.
700  Otherwise, the same allocator used for |org| is used. */
701 struct avl_table *avl_copy(const struct avl_table *org, avl_copy_func * copy,
702  avl_item_func * destroy,
703  struct libavl_allocator *allocator)
704 {
705  struct avl_node *stack[2 * (AVL_MAX_HEIGHT + 1)];
706  int height = 0;
707 
708  struct avl_table *new;
709  const struct avl_node *x;
710  struct avl_node *y;
711 
712  assert(org != NULL);
713  new = avl_create(org->avl_compare, org->avl_param,
714  allocator != NULL ? allocator : org->avl_alloc);
715  if (new == NULL)
716  return NULL;
717 
718  new->avl_count = org->avl_count;
719  if (new->avl_count == 0)
720  return new;
721 
722  x = (const struct avl_node *)&org->avl_root;
723  y = (struct avl_node *)&new->avl_root;
724  while (x != NULL) {
725  while (x->avl_link[0] != NULL) {
726  assert(height < 2 * (AVL_MAX_HEIGHT + 1));
727 
728  y->avl_link[0] =
729  new->avl_alloc->libavl_malloc(new->avl_alloc,
730  sizeof *y->avl_link[0]);
731  if (y->avl_link[0] == NULL) {
732  if (y != (struct avl_node *)&new->avl_root) {
733  y->avl_data = NULL;
734  y->avl_link[1] = NULL;
735  }
736 
737  copy_error_recovery(stack, height, new, destroy);
738  return NULL;
739  }
740 
741  stack[height++] = (struct avl_node *)x;
742  stack[height++] = y;
743  x = x->avl_link[0];
744  y = y->avl_link[0];
745  }
746  y->avl_link[0] = NULL;
747 
748  for (;;) {
749  y->avl_balance = x->avl_balance;
750  if (copy == NULL)
751  y->avl_data = x->avl_data;
752  else {
753  y->avl_data = copy(x->avl_data, org->avl_param);
754  if (y->avl_data == NULL) {
755  y->avl_link[1] = NULL;
756  copy_error_recovery(stack, height, new, destroy);
757  return NULL;
758  }
759  }
760 
761  if (x->avl_link[1] != NULL) {
762  y->avl_link[1] =
763  new->avl_alloc->libavl_malloc(new->avl_alloc,
764  sizeof *y->avl_link[1]);
765  if (y->avl_link[1] == NULL) {
766  copy_error_recovery(stack, height, new, destroy);
767  return NULL;
768  }
769 
770  x = x->avl_link[1];
771  y = y->avl_link[1];
772  break;
773  }
774  else
775  y->avl_link[1] = NULL;
776 
777  if (height <= 2)
778  return new;
779 
780  y = stack[--height];
781  x = stack[--height];
782  }
783  }
784 
785  return new;
786 }
787 
788 /* Frees storage allocated for |tree|.
789  If |destroy != NULL|, applies it to each data item in inorder. */
790 void avl_destroy(struct avl_table *tree, avl_item_func * destroy)
791 {
792  struct avl_node *p, *q;
793 
794  assert(tree != NULL);
795 
796  p = tree->avl_root;
797  while (p != NULL) {
798  if (p->avl_link[0] == NULL) {
799  q = p->avl_link[1];
800  if (destroy != NULL && p->avl_data != NULL)
801  destroy(p->avl_data, tree->avl_param);
802  tree->avl_alloc->libavl_free(tree->avl_alloc, p);
803  }
804  else {
805  q = p->avl_link[0];
806  p->avl_link[0] = q->avl_link[1];
807  q->avl_link[1] = p;
808  }
809  p = q;
810  }
811 
812  tree->avl_alloc->libavl_free(tree->avl_alloc, tree);
813 }
814 
815 /* Allocates |size| bytes of space using |malloc()|.
816  Returns a null pointer if allocation fails. */
817 void *avl_malloc(struct libavl_allocator *allocator, size_t size)
818 {
819  assert(allocator != NULL && size > 0);
820 
821  return malloc(size);
822 }
823 
824 /* Frees |block|. */
825 void avl_free(struct libavl_allocator *allocator, void *block)
826 {
827  assert(allocator != NULL && block != NULL);
828  free(block);
829 }
830 
831 /* Default memory allocator that uses |malloc()| and |free()|. */
833  avl_malloc,
834  avl_free
835 };
836 
837 #undef NDEBUG
838 #include <assert.h>
839 
840 /* Asserts that |avl_insert()| succeeds at inserting |item| into |table|. */
841 void (avl_assert_insert) (struct avl_table * table, void *item)
842 {
843  void **p = avl_probe(table, item);
844 
845  assert(p != NULL && *p == item);
846 }
847 
848 /* Asserts that |avl_delete()| really removes |item| from |table|,
849  and returns the removed item. */
850 void *(avl_assert_delete) (struct avl_table * table, void *item)
851 {
852  void *p = avl_delete(table, item);
853 
854  assert(p != NULL);
855 
856  return p;
857 }
struct avl_table * avl_copy(const struct avl_table *org, avl_copy_func *copy, avl_item_func *destroy, struct libavl_allocator *allocator)
Definition: avl.c:701
void avl_destroy(struct avl_table *tree, avl_item_func *destroy)
Definition: avl.c:790
void avl_item_func(void *avl_item, void *avl_param)
Definition: avl.h:31
void * avl_delete(struct avl_table *tree, const void *item)
Definition: avl.c:230
void *(* libavl_malloc)(struct libavl_allocator *, size_t libavl_size)
Definition: avl.h:39
void * avl_t_next(struct avl_traverser *trav)
Definition: avl.c:567
struct avl_node * avl_node
Definition: avl.h:77
void *() avl_assert_delete(struct avl_table *table, void *item)
Definition: avl.c:850
size_t avl_count
Definition: avl.h:61
unsigned long avl_generation
Definition: avl.h:62
struct avl_table * avl_create(avl_comparison_func *compare, void *param, struct libavl_allocator *allocator)
Definition: avl.c:38
void * avl_t_insert(struct avl_traverser *trav, struct avl_table *tree, void *item)
Definition: avl.c:521
struct avl_node * avl_link[2]
Definition: avl.h:68
void() avl_assert_insert(struct avl_table *table, void *item)
Definition: avl.c:841
Definition: avl.h:66
unsigned long avl_generation
Definition: avl.h:81
void(* libavl_free)(struct libavl_allocator *, void *libavl_block)
Definition: avl.h:40
void avl_free(struct libavl_allocator *allocator, void *block)
Definition: avl.c:825
void free(void *)
void * avl_copy_func(void *avl_item, void *avl_param)
Definition: avl.h:32
#define NULL
Definition: ccmath.h:32
void * avl_insert(struct avl_table *table, void *item)
Definition: avl.c:202
#define x
void * avl_t_find(struct avl_traverser *trav, struct avl_table *tree, void *item)
Definition: avl.c:484
void * avl_t_prev(struct avl_traverser *trav)
Definition: avl.c:613
void * avl_find(const struct avl_table *tree, const void *item)
Definition: avl.c:64
void * malloc(YYSIZE_T)
void * avl_t_cur(struct avl_traverser *trav)
Definition: avl.c:657
struct avl_node * avl_stack[AVL_MAX_HEIGHT]
Definition: avl.h:78
#define assert(condition)
Definition: lz4.c:324
void * avl_t_replace(struct avl_traverser *trav, void *new)
Definition: avl.c:667
int avl_comparison_func(const void *avl_a, const void *avl_b, void *avl_param)
Definition: avl.h:29
void ** avl_probe(struct avl_table *tree, void *item)
Definition: avl.c:87
struct avl_node * avl_root
Definition: avl.h:57
struct libavl_allocator avl_allocator_default
Definition: avl.c:832
#define AVL_MAX_HEIGHT
Definition: avl.h:51
Definition: avl.h:55
void * avl_t_copy(struct avl_traverser *trav, const struct avl_traverser *src)
Definition: avl.c:546
signed char avl_balance
Definition: avl.h:70
void * avl_replace(struct avl_table *table, void *item)
Definition: avl.c:213
void * avl_data
Definition: avl.h:69
void * avl_t_first(struct avl_traverser *trav, struct avl_table *tree)
Definition: avl.c:432
int compare(const void *a, const void *b)
Definition: dgraph.c:177
avl_comparison_func * avl_compare
Definition: avl.h:58
void * avl_param
Definition: avl.h:59
size_t avl_height
Definition: avl.h:80
struct libavl_allocator * avl_alloc
Definition: avl.h:60
void * avl_malloc(struct libavl_allocator *allocator, size_t size)
Definition: avl.c:817
void avl_t_init(struct avl_traverser *trav, struct avl_table *tree)
Definition: avl.c:421
struct avl_table * avl_table
Definition: avl.h:76
double r
Definition: r_raster.c:39
void * avl_t_last(struct avl_traverser *trav, struct avl_table *tree)
Definition: avl.c:457