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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
/*************************************************************************
* Copyright (C) 2009-2014 Tavian Barnes <tavianator@tavianator.com> *
* *
* This file is part of The Dimension Benchmark Suite. *
* *
* The Dimension Benchmark 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 Benchmark 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 "../platform/platform.c"
#include "../concurrency/threads.c"
#include "../concurrency/future.c"
#include "../bvh/bvh.c"
#include "../bvh/prtree.c"
#include <sandglass.h>
#include <stdlib.h>
static bool
dmnsn_fake_intersection_fn(const dmnsn_object *object, dmnsn_ray ray,
dmnsn_intersection *intersection)
{
intersection->t = (object->aabb.min.Z - ray.x0.Z)/ray.n.Z;
intersection->normal = dmnsn_x;
return true;
}
static bool
dmnsn_fake_inside_fn(const dmnsn_object *object, dmnsn_vector point)
{
return true;
}
static dmnsn_aabb
dmnsn_fake_bounding_fn(const dmnsn_object *object, dmnsn_matrix trans)
{
dmnsn_vector a, b;
for (unsigned int i = 0; i < 3; ++i) {
a.n[i] = 2.0*((double)rand())/RAND_MAX - 1.0;
b.n[i] = 2.0*((double)rand())/RAND_MAX - 1.0;
}
return dmnsn_new_aabb(dmnsn_vector_min(a, b), dmnsn_vector_max(a, b));
}
static dmnsn_object_vtable dmnsn_fake_vtable = {
.intersection_fn = dmnsn_fake_intersection_fn,
.inside_fn = dmnsn_fake_inside_fn,
.bounding_fn = dmnsn_fake_bounding_fn,
};
static dmnsn_object *
dmnsn_new_fake_object(dmnsn_pool *pool)
{
dmnsn_object *object = dmnsn_new_object(pool);
object->vtable = &dmnsn_fake_vtable;
return object;
}
int
main(void)
{
const size_t nobjects = 10000;
sandglass_t sandglass;
if (sandglass_init_monotonic(&sandglass, SANDGLASS_CPUTIME) != 0) {
perror("sandglass_create()");
return EXIT_FAILURE;
}
dmnsn_pool *pool = dmnsn_new_pool();
dmnsn_array *objects = DMNSN_PALLOC_ARRAY(pool, dmnsn_object *);
dmnsn_texture *texture = dmnsn_new_texture(pool);
texture->pigment = dmnsn_new_pigment(pool);
for (size_t i = 0; i < nobjects; ++i) {
dmnsn_object *object = dmnsn_new_fake_object(pool);
object->texture = texture;
dmnsn_object_precompute(object);
dmnsn_array_push(objects, &object);
}
dmnsn_bvh *bvh;
sandglass_bench_noprecache(&sandglass, {
bvh = dmnsn_new_bvh(objects, DMNSN_BVH_PRTREE);
});
printf("dmnsn_new_bvh(DMNSN_BVH_PRTREE): %ld\n", sandglass.grains);
// dmnsn_bvh_intersection()
dmnsn_ray ray = dmnsn_new_ray(
dmnsn_new_vector( 1.0, 1.0, -2.0),
dmnsn_new_vector(-0.5, -0.5, 1.0)
);
dmnsn_intersection intersection;
sandglass_bench_fine(&sandglass, {
dmnsn_bvh_intersection(bvh, ray, &intersection, true);
});
printf("dmnsn_bvh_intersection(): %ld\n", sandglass.grains);
sandglass_bench_fine(&sandglass, {
dmnsn_bvh_intersection(bvh, ray, &intersection, false);
});
printf("dmnsn_bvh_intersection(nocache): %ld\n", sandglass.grains);
// dmnsn_bvh_inside()
sandglass_bench_fine(&sandglass, {
dmnsn_bvh_inside(bvh, dmnsn_zero);
});
printf("dmnsn_bvh_inside(): %ld\n", sandglass.grains);
// Cleanup
dmnsn_delete_bvh(bvh);
dmnsn_delete_pool(pool);
return EXIT_SUCCESS;
}
|
