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
126
127
128
129
130
131
132
|
//
// Created by Keuin on 2022/4/11.
//
#include <gtest/gtest.h>
#include "vec.h"
TEST(Vec, VecAdd) {
vec3i a{1, 1, 1}, b{2, 2, 2}, c{3, 3, 3};
ASSERT_EQ(a + b, c);
vec3d d{1.1, 2.2, 3.3}, e{4.4, 5.5, 6.6}, f{5.5, 7.7, 9.9};
ASSERT_EQ(d + e, f);
}
TEST(Vec, VecMin) {
vec3i a{1, 1, 1}, b{2, 2, 2}, c{-1, -1, -1};
ASSERT_EQ(a - b, c);
ASSERT_EQ(-a, c);
vec3d d{1.1, 2.2, 3.3}, e{-4.4, -5.5, -6.6}, f{5.5, 7.7, 9.9};
ASSERT_EQ(d - e, f);
ASSERT_EQ(d + (-e), f);
}
TEST(Vec, CrossProduct) {
vec3i a{1, 1, 1}, b{2, 2, 2}, c{3, 4, 5}, d{6, 7, 8}, e{-3, 6, -3};
ASSERT_EQ(cross(a, b), vec3i{});
ASSERT_EQ(a.cross(b), vec3i{});
ASSERT_EQ(cross(c, d), e);
ASSERT_EQ(c.cross(d), e);
vec3d f{3, 4, 5}, g{6, 7, 8}, h{-3, 6, -3};
ASSERT_EQ(cross(f, g), h);
ASSERT_EQ(f.cross(g), h);
vec3d u{2, -2, -2}, v{2.5, -5, -2.5}, w{-5, 0, -5};
ASSERT_EQ(cross(u, v), w);
vec3d x{1.6, 3.5, 11}, y{3, 7, 15}, z{-24.5, 9, 0.7};
ASSERT_EQ(cross(x, y), z);
}
TEST(Vec, DotProduct) {
vec3i c{3, 4, 5}, d{6, 7, 8};
int e = 18 + 28 + 40;
ASSERT_EQ(dot(c, d), e);
ASSERT_EQ(c.dot(d), e);
vec3d f{3, 4, 5}, g{6, 7, 8};
ASSERT_EQ(dot(f, g), e);
ASSERT_EQ(f.dot(g), e);
}
TEST(Vec, MulByScalar) {
vec3i a{1, 1, 1}, b{2, 2, 2};
vec3d c{7, 14, 21}, d{1, 2, 3};
ASSERT_EQ(a * 2, b);
ASSERT_EQ(a / 0.5, b);
ASSERT_EQ(2 * a, b);
ASSERT_EQ(c * (1.0 / 7), d);
ASSERT_EQ(c / 7, d);
ASSERT_EQ(c / 7.0, d);
}
TEST(Vec, Norm) {
vec3d a{1, 1, 1}, b{2, 2, 2};
ASSERT_EQ(a.norm(), sqrt(3));
ASSERT_EQ(b.norm(), sqrt(12));
}
TEST(Vec, UnitVec) {
vec3d a{1, 2, 2}, b{1.0 / 3, 2.0 / 3, 2.0 / 3};
ASSERT_EQ(a.unit_vec(), b);
}
TEST(Vec, Mod2) {
vec3i a{1, 2, 3};
vec3d b{2.5, 3, 1.2};
ASSERT_EQ(a.mod2(), 14);
ASSERT_LE(abs(b.mod2() - (2.5 * 2.5 + 3 * 3 + 1.2 * 1.2)), 1e-10);
}
TEST(Vec, IsZero) {
vec3i a{1, 0, 0}, b{0, -1, 0}, c{1, 2, 3};
ASSERT_FALSE(a.is_zero());
ASSERT_FALSE(b.is_zero());
ASSERT_FALSE(c.is_zero());
vec3d d{0.1, 0, 0}, e{0, -0.1, 0}, f{0.1, 0.1, 0.1};
ASSERT_FALSE(d.is_zero());
ASSERT_FALSE(e.is_zero());
ASSERT_FALSE(f.is_zero());
vec3i g{0, 0, 0};
vec3d h{0, 0, 0}, i{1e-10, 0, 0}, j{1e-10, 1e-10, 1e-10};
ASSERT_TRUE(g.is_zero());
ASSERT_TRUE(h.is_zero());
ASSERT_TRUE(i.is_zero());
ASSERT_TRUE(j.is_zero());
}
TEST(Vec, Reflect) {
vec3d n{1, 0, 0}, u{-1, 1.1, 0}, v{1, 1.1, 0};
ASSERT_EQ(v, n.reflect(u));
}
TEST(Vec, Refract) {
vec3d n{1, 0, 0}, u{-1, 0, -1}, v{-sqrt(14), 0, -sqrt(2)};
ASSERT_EQ(u.unit_vec(), n.refract<true>(u.unit_vec(), 1).unit_vec());
ASSERT_EQ(u.unit_vec(), n.refract<false>(u.unit_vec(), 1).unit_vec());
ASSERT_EQ(v.unit_vec(), n.refract<true>(u.unit_vec(), 0.5).unit_vec());
ASSERT_EQ(v.unit_vec(), n.refract<false>(u.unit_vec(), 0.5).unit_vec());
}
TEST(Vec, Refract_TIR) {
vec3d n{1, 0, 0}, u{-1, 0, -sqrt(3)}, v{1, 0, -sqrt(3)};
ASSERT_EQ(v.unit_vec(), n.refract<true>(u.unit_vec(), 2));
}
TEST(Vec, VecParallel) {
vec3d a{1, 1, 2}, b{1.1, 1.1, 2.2}, c{0,0,0}, d{1, 2, 1};
ASSERT_TRUE(a.parallel(b));
ASSERT_TRUE(b.parallel(a));
ASSERT_TRUE(a.parallel(c));
ASSERT_TRUE(c.parallel(a));
ASSERT_FALSE(a.parallel(d));
ASSERT_FALSE(d.parallel(a));
ASSERT_FALSE(b.parallel(d));
ASSERT_FALSE(d.parallel(b));
}
|
