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
133
134
135
136
137
138
139
140
141
142
143
144
145
146
|
/*
* (C) 2003-2006 Gabest
* (C) 2006-2012 see Authors.txt
*
* This file is part of MPC-HC.
*
* MPC-HC 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.
*
* MPC-HC 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/>.
*
*/
#pragma once
#include <math.h>
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif
#define EPSILON (1e-7)
#define BIGNUMBER (1e+9)
#define IsZero(d) (fabs(d) < EPSILON)
#define Sgn(d) (IsZero(d) ? 0 : (d) > 0 ? 1 : -1)
//#define SgnPow(d, p) (IsZero(d) ? 0 : (pow(fabs(d), p) * Sgn(d)))
class Vector
{
public:
float x, y, z;
Vector() { x = y = z = 0; }
Vector(float x, float y, float z);
void Set(float x, float y, float z);
Vector Normal(Vector& a, Vector& b);
float Angle(Vector& a, Vector& b);
float Angle(Vector& a);
void Angle(float& u, float& v); // returns spherical coords in radian, -M_PI_2 <= u <= M_PI_2, -M_PI <= v <= M_PI
Vector Angle(); // does like prev., returns 'u' in 'ret.x', and 'v' in 'ret.y'
Vector Unit();
Vector& Unitalize();
float Length();
float Sum(); // x + y + z
float CrossSum(); // xy + xz + yz
Vector Cross(); // xy, xz, yz
Vector Pow(float exp);
Vector& Min(Vector& a);
Vector& Max(Vector& a);
Vector Abs();
Vector Reflect(Vector& n);
Vector Refract(Vector& n, float nFront, float nBack, float* nOut = NULL);
Vector Refract2(Vector& n, float nFrom, float nTo, float* nOut = NULL);
Vector operator - ();
float& operator [](size_t i);
float operator | (Vector& v); // dot
Vector operator % (Vector& v); // cross
bool operator == (const Vector& v) const;
bool operator != (const Vector& v) const;
Vector operator + (float d);
Vector operator + (Vector& v);
Vector operator - (float d);
Vector operator - (Vector& v);
Vector operator * (float d);
Vector operator * (Vector& v);
Vector operator / (float d);
Vector operator / (Vector& v);
Vector& operator += (float d);
Vector& operator += (Vector& v);
Vector& operator -= (float d);
Vector& operator -= (Vector& v);
Vector& operator *= (float d);
Vector& operator *= (Vector& v);
Vector& operator /= (float d);
Vector& operator /= (Vector& v);
template<typename T> static float DegToRad(T angle) { return (float)(angle * M_PI / 180); }
};
class Ray
{
public:
Vector p, d;
Ray() {}
Ray(Vector& p, Vector& d);
void Set(Vector& p, Vector& d);
float GetDistanceFrom(Ray& r); // r = plane
float GetDistanceFrom(Vector& v); // v = point
Vector operator [](float t);
};
class XForm
{
class Matrix
{
public:
float mat[4][4];
Matrix();
void Initalize();
Matrix operator * (Matrix& m);
Matrix& operator *= (Matrix& m);
} m;
bool m_isWorldToLocal;
public:
XForm() : m_isWorldToLocal(true) {}
XForm(Ray& r, Vector& s, bool isWorldToLocal = true);
void Initalize();
void Initalize(Ray& r, Vector& s, bool isWorldToLocal = true);
void operator *= (Vector& s); // scale
void operator += (Vector& t); // translate
void operator <<= (Vector& r); // rotate
void operator /= (Vector& s); // scale
void operator -= (Vector& t); // translate
void operator >>= (Vector& r); // rotate
// transformations
Vector operator < (Vector& n); // normal
Vector operator << (Vector& v); // vector
Ray operator << (Ray& r); // ray
};
|