summaryrefslogtreecommitdiffstats
path: root/bench/libdimension/prtree.c
blob: 1f1c5ddfa5d244707d0926347e1e4f90d8e9c1ce (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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
/*************************************************************************
 * Copyright (C) 2010 Tavian Barnes <tavianator@gmail.com>               *
 *                                                                       *
 * This file is part of The Dimension Test Suite.                        *
 *                                                                       *
 * The Dimension Test Suite is free software; you can redistribute it    *
 * and/or modify it under the terms of the GNU General Public License as *
 * published by the Free Software Foundation; either version 3 of the    *
 * License, or (at your option) any later version.                       *
 *                                                                       *
 * The Dimension Test Suite is distributed in the hope that it will be   *
 * useful, but WITHOUT ANY WARRANTY; without even the implied warranty   *
 * of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU  *
 * General Public License for more details.                              *
 *                                                                       *
 * You should have received a copy of the GNU General Public License     *
 * along with this program.  If not, see <http://www.gnu.org/licenses/>. *
 *************************************************************************/

#include "dimension_impl.h"
#include <sandglass.h>
#include <stdlib.h>

bool
dmnsn_fake_intersection_fn(const dmnsn_object *object, dmnsn_line line,
                           dmnsn_intersection *intersection)
{
  intersection->t = (object->bounding_box.min.z - line.x0.z)/line.n.z;
  return true;
}

void
dmnsn_randomize_bounding_box(dmnsn_object *object)
{
  dmnsn_vector a, b;

  a.x = 2.0*((double)rand())/RAND_MAX - 1.0;
  a.y = 2.0*((double)rand())/RAND_MAX - 1.0;
  a.z = 2.0*((double)rand())/RAND_MAX - 1.0;

  b.x = 2.0*((double)rand())/RAND_MAX - 1.0;
  b.y = 2.0*((double)rand())/RAND_MAX - 1.0;
  b.z = 2.0*((double)rand())/RAND_MAX - 1.0;

  object->bounding_box.min = dmnsn_vector_min(a, b);
  object->bounding_box.max = dmnsn_vector_max(a, b);
}

int
main()
{
  dmnsn_prtree *tree;
  dmnsn_intersection intersection;
  dmnsn_line ray;
  const size_t nobjects = 128;
  dmnsn_array *objects;

  sandglass_t sandglass;
  if (sandglass_init_monotonic(&sandglass, SANDGLASS_CPUTIME) != 0) {
    perror("sandglass_create()");
    return EXIT_FAILURE;
  }

  objects = dmnsn_new_array(sizeof(dmnsn_object *));
  for (size_t i = 0; i < nobjects; ++i) {
    dmnsn_object *object = dmnsn_new_object();
    /* Generate a bounding box in  (-1, -1, -1), (1, 1, 1) */
    dmnsn_randomize_bounding_box(object);
    object->intersection_fn = &dmnsn_fake_intersection_fn;
    dmnsn_array_push(objects, &object);
  }

  sandglass_bench_noprecache(&sandglass, {
    tree = dmnsn_new_prtree(objects);
  });
  printf("dmnsn_new_prtree(): %ld\n", sandglass.grains);

  /* dmnsn_prtree_intersection() */
  ray.x0 = dmnsn_new_vector(0.0, 0.0, -2.0);
  ray.n  = dmnsn_new_vector(0.0, 0.0, 1.0);

  sandglass_bench_fine(&sandglass, {
    dmnsn_prtree_intersection(tree, ray, &intersection);
  });
  printf("dmnsn_prtree_intersection(): %ld\n", sandglass.grains);

  /* Cleanup */
  dmnsn_delete_prtree(tree);
  for (size_t i = 0; i < nobjects; ++i) {
    dmnsn_object *object;
    dmnsn_array_get(objects, i, &object);
    dmnsn_delete_object(object);
  }

  dmnsn_delete_array(objects);
  return EXIT_SUCCESS;
}