cvmath 0.0.8

Computer Graphics Vector Math Library
Documentation 
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

use super::*;

fn leaf(index: u32) -> BvhLeaf {
	let mut leaf = BvhLeaf::EMPTY;
	leaf.children[0] = index;
	leaf
}

#[test]
fn trace_returns_far_hit_when_near_child_misses() {

	let bvh = Bvh3 {
		nodes: vec![
			BvhNode::Branch(BvhBranch {
				left_bounds: Bounds3(Point3(-5.0, -1.0, -1.0), Point3(-1.0, 1.0, 1.0)),
				right_bounds: Bounds3(Point3(-4.0, -1.0, -1.0), Point3(2.0, 1.0, 1.0)),
				left_child: 1,
				right_child: 2,
			}),
			BvhNode::Leaf(leaf(0)),
			BvhNode::Leaf(leaf(1)),
		],
	};

	let ray = Ray3(Point3(-10.0, 0.0, 0.0), Vec3(1.0, 0.0, 0.0), Interval(0.0, 20.0));
	let hit = bvh.trace(&ray, |index, _| match index {
		0 => None,
		1 => Some(Hit3 {
			point: Point3(-2.0, 0.0, 0.0),
			distance: 8.0,
			normal: Vec3(-1.0, 0.0, 0.0),
			index,
			side: HitSide::Entry,
		}),
		_ => None,
	}).unwrap();

	assert_eq!(hit.index, 1);
	assert_eq!(hit.distance, 8.0);
}

#[test]
fn trace_prefers_closest_hit_when_ray_starts_inside_overlapping_bounds() {
	let bvh = Bvh3 {
		nodes: vec![
			BvhNode::Branch(BvhBranch {
				left_bounds: Bounds3(Point3(-10.0, -10.0, -10.0), Point3(10.0, 10.0, 10.0)),
				right_bounds: Bounds3(Point3(-5.0, -5.0, -5.0), Point3(20.0, 20.0, 20.0)),
				left_child: 1,
				right_child: 2,
			}),
			BvhNode::Leaf(leaf(0)),
			BvhNode::Leaf(leaf(1)),
		],
	};

	let ray = Ray3(Point3(0.0, 0.0, 0.0), Vec3(1.0, 0.0, 0.0), Interval(0.0, 100.0));
	let hit = bvh.trace(&ray, |index, _| match index {
		0 => Some(Hit3 {
			point: Point3(5.0, 0.0, 0.0),
			distance: 5.0,
			normal: Vec3(-1.0, 0.0, 0.0),
			index,
			side: HitSide::Entry,
		}),
		1 => Some(Hit3 {
			point: Point3(1.0, 0.0, 0.0),
			distance: 1.0,
			normal: Vec3(-1.0, 0.0, 0.0),
			index,
			side: HitSide::Entry,
		}),
		_ => None,
	}).unwrap();

	assert_eq!(hit.index, 1);
	assert_eq!(hit.distance, 1.0);
}

#[test]
fn inspect_visits_bounds_depth_first() {
	let bvh = Bvh3 {
		nodes: vec![
			BvhNode::Branch(BvhBranch {
				left_bounds: Bounds3(Point3(-4.0, -4.0, -4.0), Point3(0.0, 4.0, 4.0)),
				right_bounds: Bounds3(Point3(0.0, -2.0, -2.0), Point3(4.0, 2.0, 2.0)),
				left_child: 1,
				right_child: 2,
			}),
			BvhNode::Branch(BvhBranch {
				left_bounds: Bounds3(Point3(-4.0, 0.0, 0.0), Point3(-2.0, 4.0, 4.0)),
				right_bounds: Bounds3(Point3(-2.0, -4.0, -4.0), Point3(0.0, 0.0, 0.0)),
				left_child: 3,
				right_child: 4,
			}),
			BvhNode::Leaf(leaf(0)),
			BvhNode::Leaf(leaf(1)),
			BvhNode::Leaf(leaf(2)),
		],
	};

	let mut visited = Vec::new();
	bvh.visit(|depth, bounds| {
		visited.push((depth, *bounds));
		true
	});

	assert_eq!(visited, vec![
		(1, Bounds3(Point3(-4.0, -4.0, -4.0), Point3(0.0, 4.0, 4.0))),
		(2, Bounds3(Point3(-4.0, 0.0, 0.0), Point3(-2.0, 4.0, 4.0))),
		(2, Bounds3(Point3(-2.0, -4.0, -4.0), Point3(0.0, 0.0, 0.0))),
		(1, Bounds3(Point3(0.0, -2.0, -2.0), Point3(4.0, 2.0, 2.0))),
	]);
}