nuclear@4: #include <math.h>
nuclear@4: 
nuclear@4: namespace gph {
nuclear@4: 
nuclear@4: // ---- Vector3 ----
nuclear@4: 
nuclear@4: inline void Vector3::normalize()
nuclear@4: {
nuclear@4: 	float len = (float)sqrt(x * x + y * y + z * z);
nuclear@4: 	if(len != 0.0f) {
nuclear@4: 		x /= len;
nuclear@4: 		y /= len;
nuclear@4: 		z /= len;
nuclear@4: 	}
nuclear@4: }
nuclear@4: 
nuclear@4: inline float &Vector3::operator[] (int idx)
nuclear@4: {
nuclear@4: 	return idx == 0 ? x : (idx == 1 ? y : z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline const float &Vector3::operator[] (int idx) const
nuclear@4: {
nuclear@4: 	return idx == 0 ? x : (idx == 1 ? y : z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator +(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return Vector3(a.x + b.x, a.y + b.y, a.z + b.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator -(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return Vector3(a.x - b.x, a.y - b.y, a.z - b.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator *(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return Vector3(a.x * b.x, a.y * b.y, a.z * b.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator /(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return Vector3(a.x / b.x, a.y / b.y, a.z / b.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator *(const Vector3 &v, float s)
nuclear@4: {
nuclear@4: 	return Vector3(v.x * s, v.y * s, v.z * s);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator *(float s, const Vector3 &v)
nuclear@4: {
nuclear@4: 	return Vector3(s * v.x, s * v.y, s * v.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator /(const Vector3 &v, float s)
nuclear@4: {
nuclear@4: 	return Vector3(v.x / s, v.y / s, v.z / s);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 operator /(float s, const Vector3 &v)
nuclear@4: {
nuclear@4: 	return Vector3(s / v.x, s / v.y, s / v.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator +=(Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	a.x += b.x;
nuclear@4: 	a.y += b.y;
nuclear@4: 	a.z += b.z;
nuclear@4: 	return a;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator -=(Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	a.x -= b.x;
nuclear@4: 	a.y -= b.y;
nuclear@4: 	a.z -= b.z;
nuclear@4: 	return a;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator *=(Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	a.x *= b.x;
nuclear@4: 	a.y *= b.y;
nuclear@4: 	a.z *= b.z;
nuclear@4: 	return a;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator /=(Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	a.x /= b.x;
nuclear@4: 	a.y /= b.y;
nuclear@4: 	a.z /= b.z;
nuclear@4: 	return a;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator *=(Vector3 &v, float s)
nuclear@4: {
nuclear@4: 	v.x *= s;
nuclear@4: 	v.y *= s;
nuclear@4: 	v.z *= s;
nuclear@4: 	return v;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 &operator /=(Vector3 &v, float s)
nuclear@4: {
nuclear@4: 	v.x /= s;
nuclear@4: 	v.y /= s;
nuclear@4: 	v.z /= s;
nuclear@4: 	return v;
nuclear@4: }
nuclear@4: 
nuclear@4: inline bool operator ==(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return a.x == b.x && a.y == b.y && a.z == b.z;
nuclear@4: }
nuclear@4: 
nuclear@4: inline bool operator !=(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return !(a == b);
nuclear@4: }
nuclear@4: 
nuclear@4: inline float dot(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return a.x * b.x + a.y * b.y + a.z * b.z;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 cross(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return Vector3(a.y * b.z - a.z * b.y,
nuclear@4: 			a.z * b.x - a.x * b.z,
nuclear@4: 			a.x * b.y - a.y * b.x);
nuclear@4: }
nuclear@4: 
nuclear@4: inline float length(const Vector3 &v)
nuclear@4: {
nuclear@4: 	return (float)sqrt(v.x * v.x + v.y * v.y + v.z * v.z);
nuclear@4: }
nuclear@4: 
nuclear@4: inline float length_sq(const Vector3 &v)
nuclear@4: {
nuclear@4: 	return v.x * v.x + v.y * v.y + v.z * v.z;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 normalize(const Vector3 &v)
nuclear@4: {
nuclear@4: 	float len = length(v);
nuclear@4: 	if(len == 0.0f) {
nuclear@4: 		return v;
nuclear@4: 	}
nuclear@4: 
nuclear@4: 	return Vector3(v.x / len, v.y / len, v.z / len);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 reflect(const Vector3 &v, const Vector3 &n)
nuclear@4: {
nuclear@4: 	return v - n * dot(n, v) * 2.0;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 refract(const Vector3 &v, const Vector3 &n, float ior)
nuclear@4: {
nuclear@4: 	float ndotv = dot(n, v);
nuclear@4: 	float k = 1.0f - ior * ior * (1.0f - ndotv * ndotv);
nuclear@4: 	if(k < 0.0f) {
nuclear@4: 		return Vector3();
nuclear@4: 	}
nuclear@4: 	return ior * v - (ior * ndotv + sqrt(k)) * n;
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 refract(const Vector3 &v, const Vector3 &n, float from_ior, float to_ior)
nuclear@4: {
nuclear@4: 	if(to_ior == 0.0f) to_ior = 1.0f;
nuclear@4: 	return refract(v, n, from_ior / to_ior);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 distance(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return length(a - b);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 distance_sq(const Vector3 &a, const Vector3 &b)
nuclear@4: {
nuclear@4: 	return length_sq(a - b);
nuclear@4: }
nuclear@4: 
nuclear@4: inline Vector3 faceforward(const Vector3 &n, const Vector3 &vi, const Vector3 &ng)
nuclear@4: {
nuclear@4: 	return dot(ng, i) < 0.0f ? n : -n;
nuclear@4: }
nuclear@4: 
nuclear@4: // ---- Vector4 ----
nuclear@4: 
nuclear@4: 
nuclear@4: inline void Vector4::normalize()
nuclear@4: {
nuclear@4: 	float len = (float)sqrt(x * x + y * y + z * z + w * w);
nuclear@4: 	if(len != 0.0f) {
nuclear@4: 		x /= len;
nuclear@4: 		y /= len;
nuclear@4: 		z /= len;
nuclear@4: 		w /= len;
nuclear@4: 	}
nuclear@4: }
nuclear@4: 
nuclear@4: inline float &Vector4::operator[] (int idx)
nuclear@4: {
nuclear@4: 	return idx == 0 ? x : (idx == 1 ? y : (idx == 2 ? z : w));
nuclear@4: }
nuclear@4: 
nuclear@4: inline const float &Vector4::operator[] (int idx) const
nuclear@4: {
nuclear@4: 	return idx == 0 ? x : (idx == 1 ? y : (idx == 2 ? z : w));
nuclear@4: }
nuclear@4: 
nuclear@4: }	// namespace gph