src/tree.c

Wed, 20 Mar 2024 23:31:41 +0100

author
Mike Becker <universe@uap-core.de>
date
Wed, 20 Mar 2024 23:31:41 +0100
changeset 845
2615317311b7
parent 840
4f02995ce44e
child 848
6456036bbb37
permissions
-rw-r--r--

add cx_tree_visitor()

universe@816 1 /*
universe@816 2 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS HEADER.
universe@816 3 *
universe@816 4 * Copyright 2024 Mike Becker, Olaf Wintermann All rights reserved.
universe@816 5 *
universe@816 6 * Redistribution and use in source and binary forms, with or without
universe@816 7 * modification, are permitted provided that the following conditions are met:
universe@816 8 *
universe@816 9 * 1. Redistributions of source code must retain the above copyright
universe@816 10 * notice, this list of conditions and the following disclaimer.
universe@816 11 *
universe@816 12 * 2. Redistributions in binary form must reproduce the above copyright
universe@816 13 * notice, this list of conditions and the following disclaimer in the
universe@816 14 * documentation and/or other materials provided with the distribution.
universe@816 15 *
universe@816 16 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
universe@816 17 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
universe@816 18 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
universe@816 19 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
universe@816 20 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
universe@816 21 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
universe@816 22 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
universe@816 23 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
universe@816 24 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
universe@816 25 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
universe@816 26 * POSSIBILITY OF SUCH DAMAGE.
universe@816 27 */
universe@816 28
universe@816 29 #include "cx/tree.h"
universe@816 30
universe@826 31 #include "cx/array_list.h"
universe@826 32
universe@816 33 #include <assert.h>
universe@816 34
universe@816 35 #define CX_TREE_PTR(cur, off) (*(void**)(((char*)(cur))+(off)))
universe@816 36 #define CX_TREE_PTR(cur, off) (*(void**)(((char*)(cur))+(off)))
universe@816 37 #define tree_parent(node) CX_TREE_PTR(node, loc_parent)
universe@816 38 #define tree_children(node) CX_TREE_PTR(node, loc_children)
universe@816 39 #define tree_prev(node) CX_TREE_PTR(node, loc_prev)
universe@816 40 #define tree_next(node) CX_TREE_PTR(node, loc_next)
universe@816 41
universe@816 42 void cx_tree_link(
universe@816 43 void *restrict parent,
universe@816 44 void *restrict node,
universe@816 45 ptrdiff_t loc_parent,
universe@816 46 ptrdiff_t loc_children,
universe@816 47 ptrdiff_t loc_prev,
universe@816 48 ptrdiff_t loc_next
universe@816 49 ) {
universe@816 50 void *current_parent = tree_parent(node);
universe@816 51 if (current_parent == parent) return;
universe@816 52 if (current_parent != NULL) {
universe@816 53 cx_tree_unlink(node, loc_parent, loc_children,
universe@816 54 loc_prev, loc_next);
universe@816 55 }
universe@816 56
universe@816 57 if (tree_children(parent) == NULL) {
universe@816 58 tree_children(parent) = node;
universe@816 59 } else {
universe@816 60 void *children = tree_children(parent);
universe@816 61 tree_prev(children) = node;
universe@816 62 tree_next(node) = children;
universe@816 63 tree_children(parent) = node;
universe@816 64 }
universe@816 65 tree_parent(node) = parent;
universe@816 66 }
universe@816 67
universe@816 68 void cx_tree_unlink(
universe@816 69 void *node,
universe@816 70 ptrdiff_t loc_parent,
universe@816 71 ptrdiff_t loc_children,
universe@816 72 ptrdiff_t loc_prev,
universe@816 73 ptrdiff_t loc_next
universe@816 74 ) {
universe@816 75 if (tree_parent(node) == NULL) return;
universe@816 76
universe@816 77 void *left = tree_prev(node);
universe@816 78 void *right = tree_next(node);
universe@816 79 assert(left == NULL || tree_children(tree_parent(node)) != node);
universe@816 80 if (left == NULL) {
universe@816 81 tree_children(tree_parent(node)) = right;
universe@816 82 } else {
universe@816 83 tree_next(left) = right;
universe@816 84 }
universe@816 85 if (right != NULL) tree_prev(right) = left;
universe@816 86 tree_parent(node) = NULL;
universe@816 87 tree_prev(node) = NULL;
universe@816 88 tree_next(node) = NULL;
universe@816 89 }
universe@826 90
universe@826 91 int cx_tree_search(
universe@826 92 void const *root,
universe@826 93 void const *data,
universe@826 94 cx_tree_search_func sfunc,
universe@826 95 void **result,
universe@826 96 ptrdiff_t loc_children,
universe@826 97 ptrdiff_t loc_next
universe@826 98 ) {
universe@826 99 int ret;
universe@826 100 *result = NULL;
universe@826 101
universe@826 102 // shortcut: compare root before doing anything else
universe@826 103 ret = sfunc(root, data);
universe@826 104 if (ret < 0) {
universe@826 105 return ret;
universe@826 106 } else if (ret == 0 || tree_children(root) == NULL) {
universe@826 107 *result = (void*)root;
universe@826 108 return ret;
universe@826 109 }
universe@826 110
universe@826 111 // create a working stack
universe@834 112 CX_ARRAY_DECLARE(void const*, work);
universe@833 113 cx_array_initialize(work, 32);
universe@826 114
universe@826 115 // add the children of root to the working stack
universe@826 116 {
universe@826 117 void *c = tree_children(root);
universe@826 118 while (c != NULL) {
universe@833 119 cx_array_simple_add(work, c);
universe@826 120 c = tree_next(c);
universe@826 121 }
universe@826 122 }
universe@826 123
universe@826 124 // remember a candidate for adding the data
universe@826 125 // also remember the exact return code from sfunc
universe@826 126 void *candidate = NULL;
universe@826 127 int ret_candidate = -1;
universe@826 128
universe@826 129 // process the working stack
universe@826 130 while (work_size > 0) {
universe@826 131 // pop element
universe@826 132 void const *node = work[--work_size];
universe@826 133
universe@826 134 // apply the search function
universe@826 135 ret = sfunc(node, data);
universe@826 136
universe@826 137 if (ret == 0) {
universe@826 138 // if found, exit the search
universe@826 139 *result = (void*) node;
universe@826 140 work_size = 0;
universe@826 141 break;
universe@826 142 } else if (ret > 0) {
universe@826 143 // if children might contain the data, add them to the stack
universe@826 144 void *c = tree_children(node);
universe@826 145 while (c != NULL) {
universe@833 146 cx_array_simple_add(work, c);
universe@826 147 c = tree_next(c);
universe@826 148 }
universe@826 149
universe@826 150 // remember this node in case no child is suitable
universe@826 151 if (ret_candidate < 0 || ret < ret_candidate) {
universe@826 152 candidate = (void *) node;
universe@826 153 ret_candidate = ret;
universe@826 154 }
universe@826 155 }
universe@826 156 }
universe@826 157
universe@826 158 // not found, but was there a candidate?
universe@826 159 if (ret != 0 && candidate != NULL) {
universe@826 160 ret = ret_candidate;
universe@826 161 *result = candidate;
universe@826 162 }
universe@826 163
universe@826 164 // free the working queue and return
universe@826 165 free(work);
universe@826 166 return ret;
universe@826 167 }
universe@830 168
universe@830 169 static bool cx_tree_iter_valid(void const *it) {
universe@830 170 struct cx_tree_iterator_s const *iter = it;
universe@830 171 return iter->node != NULL;
universe@830 172 }
universe@830 173
universe@830 174 static void *cx_tree_iter_current(void const *it) {
universe@830 175 struct cx_tree_iterator_s const *iter = it;
universe@830 176 return iter->node;
universe@830 177 }
universe@830 178
universe@830 179 static void cx_tree_iter_next(void *it) {
universe@830 180 struct cx_tree_iterator_s *iter = it;
universe@845 181 ptrdiff_t const loc_next = iter->loc_next;
universe@845 182 ptrdiff_t const loc_children = iter->loc_children;
universe@830 183
universe@838 184 void *children;
universe@838 185
universe@838 186 // check if we are currently exiting or entering nodes
universe@838 187 if (iter->exiting) {
universe@838 188 children = NULL;
universe@838 189 } else {
universe@838 190 children = tree_children(iter->node);
universe@838 191 }
universe@838 192
universe@836 193 if (children == NULL) {
universe@836 194 // search for the next node
universe@838 195 void *next;
universe@838 196 cx_tree_iter_search_next:
universe@838 197 // check if there is a sibling
universe@840 198 if (iter->exiting) {
universe@840 199 next = iter->next;
universe@840 200 } else {
universe@840 201 next = tree_next(iter->node);
universe@840 202 iter->next = next;
universe@840 203 }
universe@838 204 if (next == NULL) {
universe@838 205 // no sibling, we are done with this node and exit
universe@838 206 if (iter->visit_on_exit && !iter->exiting) {
universe@838 207 // iter is supposed to visit the node again
universe@838 208 iter->exiting = true;
universe@838 209 } else {
universe@838 210 iter->exiting = false;
universe@838 211 if (iter->depth == 1) {
universe@836 212 // there is no parent - we have iterated the entire tree
universe@836 213 // invalidate the iterator and free the node stack
universe@840 214 iter->node = iter->next = NULL;
universe@838 215 iter->stack_capacity = iter->depth = 0;
universe@836 216 free(iter->stack);
universe@836 217 iter->stack = NULL;
universe@836 218 } else {
universe@836 219 // the parent node can be obtained from the top of stack
universe@836 220 // this way we can avoid the loc_parent in the iterator
universe@838 221 iter->depth--;
universe@838 222 iter->node = iter->stack[iter->depth - 1];
universe@838 223 // retry with the parent node to find a sibling
universe@838 224 goto cx_tree_iter_search_next;
universe@836 225 }
universe@838 226 }
universe@838 227 } else {
universe@838 228 if (iter->visit_on_exit && !iter->exiting) {
universe@838 229 // iter is supposed to visit the node again
universe@838 230 iter->exiting = true;
universe@836 231 } else {
universe@838 232 iter->exiting = false;
universe@838 233 // move to the sibling
universe@836 234 iter->counter++;
universe@836 235 iter->node = next;
universe@836 236 // new top of stack is the sibling
universe@836 237 iter->stack[iter->depth - 1] = next;
universe@836 238 }
universe@838 239 }
universe@836 240 } else {
universe@836 241 // node has children, push the first child onto the stack and enter it
universe@836 242 cx_array_simple_add(iter->stack, children);
universe@836 243 iter->node = children;
universe@836 244 iter->counter++;
universe@836 245 }
universe@830 246 }
universe@830 247
universe@830 248 CxTreeIterator cx_tree_iterator(
universe@830 249 void *root,
universe@833 250 bool visit_on_exit,
universe@830 251 ptrdiff_t loc_children,
universe@830 252 ptrdiff_t loc_next
universe@830 253 ) {
universe@830 254 CxTreeIterator iter;
universe@830 255 iter.loc_children = loc_children;
universe@830 256 iter.loc_next = loc_next;
universe@833 257 iter.visit_on_exit = visit_on_exit;
universe@830 258
universe@830 259 // allocate stack
universe@830 260 iter.stack_capacity = 16;
universe@830 261 iter.stack = malloc(sizeof(void *) * 16);
universe@830 262 iter.depth = 0;
universe@830 263
universe@833 264 // visit the root node
universe@833 265 iter.node = root;
universe@833 266 iter.counter = 1;
universe@833 267 iter.depth = 1;
universe@833 268 iter.stack[0] = root;
universe@833 269 iter.exiting = false;
universe@830 270
universe@830 271 // assign base iterator functions
universe@830 272 iter.base.mutating = false;
universe@830 273 iter.base.remove = false;
universe@830 274 iter.base.current_impl = NULL;
universe@830 275 iter.base.valid = cx_tree_iter_valid;
universe@830 276 iter.base.next = cx_tree_iter_next;
universe@830 277 iter.base.current = cx_tree_iter_current;
universe@830 278
universe@830 279 return iter;
universe@830 280 }
universe@845 281
universe@845 282 static bool cx_tree_visitor_valid(void const *it) {
universe@845 283 struct cx_tree_visitor_s const *iter = it;
universe@845 284 return iter->node != NULL;
universe@845 285 }
universe@845 286
universe@845 287 static void *cx_tree_visitor_current(void const *it) {
universe@845 288 struct cx_tree_visitor_s const *iter = it;
universe@845 289 return iter->node;
universe@845 290 }
universe@845 291
universe@845 292 __attribute__((__nonnull__))
universe@845 293 static void cx_tree_visitor_enqueue_siblings(
universe@845 294 struct cx_tree_visitor_s *iter, void *node, ptrdiff_t loc_next) {
universe@845 295 node = tree_next(node);
universe@845 296 while (node != NULL) {
universe@845 297 struct cx_tree_visitor_queue_s *q;
universe@845 298 q = malloc(sizeof(struct cx_tree_visitor_queue_s));
universe@845 299 q->depth = iter->queue_last->depth;
universe@845 300 q->node = node;
universe@845 301 iter->queue_last->next = q;
universe@845 302 iter->queue_last = q;
universe@845 303 node = tree_next(node);
universe@845 304 }
universe@845 305 iter->queue_last->next = NULL;
universe@845 306 }
universe@845 307
universe@845 308 static void cx_tree_visitor_next(void *it) {
universe@845 309 struct cx_tree_visitor_s *iter = it;
universe@845 310 ptrdiff_t const loc_next = iter->loc_next;
universe@845 311 ptrdiff_t const loc_children = iter->loc_children;
universe@845 312
universe@845 313 // check if there is a next node
universe@845 314 if (iter->queue_next == NULL) {
universe@845 315 iter->node = NULL;
universe@845 316 return;
universe@845 317 }
universe@845 318
universe@845 319 // dequeue the next node
universe@845 320 iter->node = iter->queue_next->node;
universe@845 321 iter->depth = iter->queue_next->depth;
universe@845 322 {
universe@845 323 struct cx_tree_visitor_queue_s *q = iter->queue_next;
universe@845 324 iter->queue_next = q->next;
universe@845 325 if (iter->queue_next == NULL) {
universe@845 326 assert(iter->queue_last == q);
universe@845 327 iter->queue_last = NULL;
universe@845 328 }
universe@845 329 free(q);
universe@845 330 }
universe@845 331
universe@845 332 // increment the node counter
universe@845 333 iter->counter++;
universe@845 334
universe@845 335 // add the children of the new node to the queue
universe@845 336 void *children = tree_children(iter->node);
universe@845 337 if (children != NULL) {
universe@845 338 struct cx_tree_visitor_queue_s *q;
universe@845 339 q = malloc(sizeof(struct cx_tree_visitor_queue_s));
universe@845 340 q->depth = iter->depth + 1;
universe@845 341 q->node = children;
universe@845 342 if (iter->queue_last == NULL) {
universe@845 343 assert(iter->queue_next == NULL);
universe@845 344 iter->queue_next = q;
universe@845 345 } else {
universe@845 346 iter->queue_last->next = q;
universe@845 347 }
universe@845 348 iter->queue_last = q;
universe@845 349 cx_tree_visitor_enqueue_siblings(iter, children, loc_next);
universe@845 350 }
universe@845 351 }
universe@845 352
universe@845 353 CxTreeVisitor cx_tree_visitor(
universe@845 354 void *root,
universe@845 355 ptrdiff_t loc_children,
universe@845 356 ptrdiff_t loc_next
universe@845 357 ) {
universe@845 358 CxTreeVisitor iter;
universe@845 359 iter.loc_children = loc_children;
universe@845 360 iter.loc_next = loc_next;
universe@845 361
universe@845 362 // allocate stack
universe@845 363 iter.depth = 0;
universe@845 364
universe@845 365 // visit the root node
universe@845 366 iter.node = root;
universe@845 367 iter.counter = 1;
universe@845 368 iter.depth = 1;
universe@845 369
universe@845 370 // put all children of root into the queue
universe@845 371 void *children = tree_children(root);
universe@845 372 if (children == NULL) {
universe@845 373 iter.queue_next = NULL;
universe@845 374 iter.queue_last = NULL;
universe@845 375 } else {
universe@845 376 iter.queue_next = malloc(sizeof(struct cx_tree_visitor_queue_s));
universe@845 377 iter.queue_next->depth = 2;
universe@845 378 iter.queue_next->node = children;
universe@845 379 iter.queue_last = iter.queue_next;
universe@845 380 cx_tree_visitor_enqueue_siblings(&iter, children, loc_next);
universe@845 381 }
universe@845 382
universe@845 383 // assign base iterator functions
universe@845 384 iter.base.mutating = false;
universe@845 385 iter.base.remove = false;
universe@845 386 iter.base.current_impl = NULL;
universe@845 387 iter.base.valid = cx_tree_visitor_valid;
universe@845 388 iter.base.next = cx_tree_visitor_next;
universe@845 389 iter.base.current = cx_tree_visitor_current;
universe@845 390
universe@845 391 return iter;
universe@845 392 }
universe@845 393

mercurial