nuclear@0: #ifndef GMATH_VEC_H_
nuclear@0: #define GMATH_VEC_H_
nuclear@0: 
nuclear@0: #include <math.h>
nuclear@0: 
nuclear@0: namespace gmath {
nuclear@0: 
nuclear@0: class Vec3 {
nuclear@0: public:
nuclear@0: 	float x, y, z;
nuclear@0: 
nuclear@0: 	Vec3() : x(0), y(0), z(0) {}
nuclear@0: 	Vec3(float x_, float y_, float z_) : x(x_), y(y_), z(z_) {}
nuclear@0: 
nuclear@0: 	void normalize()
nuclear@0: 	{
nuclear@0: 		float len = length(*this);
nuclear@0: 		if(len != 0.0f) {
nuclear@0: 			x /= len;
nuclear@0: 			y /= len;
nuclear@0: 			z /= len;
nuclear@0: 		}
nuclear@0: 	}
nuclear@0: };
nuclear@0: 
nuclear@0: inline Vec3 operator +(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return Vec3(a.x + b.x, a.y + b.y, a.z + b.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator -(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return Vec3(a.x - b.x, a.y - b.y, a.z - b.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator *(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return Vec3(a.x * b.x, a.y * b.y, a.z * b.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator /(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return Vec3(a.x / b.x, a.y / b.y, a.z / b.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator *(const Vec3 &v, float s)
nuclear@0: {
nuclear@0: 	return Vec3(v.x * s, v.y * s, v.z * s);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator *(float s, const Vec3 &v)
nuclear@0: {
nuclear@0: 	return Vec3(s * v.x, s * v.y, s * v.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator /(const Vec3 &v, float s)
nuclear@0: {
nuclear@0: 	return Vec3(v.x / s, v.y / s, v.z / s);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 operator /(float s, const Vec3 &v)
nuclear@0: {
nuclear@0: 	return Vec3(s / v.x, s / v.y, s / v.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator +=(Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	a.x += b.x;
nuclear@0: 	a.y += b.y;
nuclear@0: 	a.z += b.z;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator -=(Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	a.x -= b.x;
nuclear@0: 	a.y -= b.y;
nuclear@0: 	a.z -= b.z;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator *=(Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	a.x *= b.x;
nuclear@0: 	a.y *= b.y;
nuclear@0: 	a.z *= b.z;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator /=(Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	a.x /= b.x;
nuclear@0: 	a.y /= b.y;
nuclear@0: 	a.z /= b.z;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator *=(Vec3 &v, float s)
nuclear@0: {
nuclear@0: 	v.x *= s;
nuclear@0: 	v.y *= s;
nuclear@0: 	v.z *= s;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 &operator /=(Vec3 &v, float s)
nuclear@0: {
nuclear@0: 	v.x /= s;
nuclear@0: 	v.y /= s;
nuclear@0: 	v.z /= s;
nuclear@0: 	return *this;
nuclear@0: }
nuclear@0: 
nuclear@0: inline float dot(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return a.x * b.x + a.y * b.y + a.z * b.z;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 cross(const Vec3 &a, const Vec3 &b)
nuclear@0: {
nuclear@0: 	return Vec3(a.y * b.z - a.z * b.y, a.z * b.x - a.x * b.z, a.x * b.y - a.y * b.x);
nuclear@0: }
nuclear@0: 
nuclear@0: inline float length(const Vec3 &v)
nuclear@0: {
nuclear@0: 	return (float)sqrt(v.x * v.x + v.y * v.y + v.z * v.z);
nuclear@0: }
nuclear@0: 
nuclear@0: inline float length_sq(const Vec3 &v)
nuclear@0: {
nuclear@0: 	return v.x * v.x + v.y * v.y + v.z * v.z;
nuclear@0: }
nuclear@0: 
nuclear@0: inline Vec3 normalize(const Vec3 &v)
nuclear@0: {
nuclear@0: 	float len = length(v);
nuclear@0: 	if(len == 0.0f) {
nuclear@0: 		return v;
nuclear@0: 	}
nuclear@0: 
nuclear@0: 	return Vec3(v.x / len, v.y / len, v.z / len);
nuclear@0: }
nuclear@0: 
nuclear@0: }
nuclear@0: 
nuclear@0: #endif	/* GMATH_VEC_H_ */