summaryrefslogtreecommitdiffstats
path: root/kd-forest.h
blob: 4cf5750e673f24bb74ebbb16962556451b477f9d (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
/*********************************************************************
 * kd-forest                                                         *
 * Copyright (C) 2014 Tavian Barnes <tavianator@tavianator.com>      *
 *                                                                   *
 * This program is free software. It comes without any warranty, to  *
 * the extent permitted by applicable law. You can redistribute it   *
 * and/or modify it under the terms of the Do What The Fuck You Want *
 * To Public License, Version 2, as published by Sam Hocevar. See    *
 * the COPYING file or http://www.wtfpl.net/ for more details.       *
 *********************************************************************/

#ifndef KD_FOREST_H
#define KD_FOREST_H

#include <stdbool.h>
#include <stddef.h>
#include <stdint.h>

#define KD_DIMEN 3

// Single node in a k-d tree
typedef struct kd_node_t {
  // Node coordinates
  double coords[KD_DIMEN];
  // Sub-trees
  struct kd_node_t *left, *right;
  // Used to keep track of which sub-tree a node is in during construction
  bool is_left;
  // State flags
  bool added, removed;

  // Corresponding image position for this node
  unsigned int x, y;
} kd_node_t;

void kd_node_init(kd_node_t *node, unsigned int x, unsigned int y);

// A forest of k-d trees
typedef struct {
  // Array of k-d tree roots
  kd_node_t **roots;
  // Size and capacity of the roots array
  unsigned int roots_size, roots_capacity;
  // The actual size of this tree
  size_t size;
  // The size estimate for this tree, counting removed nodes
  size_t size_est;
} kd_forest_t;

void kdf_init(kd_forest_t *kdf);
void kdf_destroy(kd_forest_t *kdf);
void kdf_insert(kd_forest_t *kdf, kd_node_t *node);
void kdf_remove(kd_forest_t *kdf, kd_node_t *node);
kd_node_t *kdf_find_nearest(kd_forest_t *kdf, const kd_node_t *target);

#endif // KD_FOREST_H