nuclear@2: /*
nuclear@2: libvmath - a vector math library
nuclear@2: Copyright (C) 2004-2011 John Tsiombikas <nuclear@member.fsf.org>
nuclear@2: 
nuclear@2: This program is free software: you can redistribute it and/or modify
nuclear@2: it under the terms of the GNU Lesser General Public License as published
nuclear@2: by the Free Software Foundation, either version 3 of the License, or
nuclear@2: (at your option) any later version.
nuclear@2: 
nuclear@2: This program is distributed in the hope that it will be useful,
nuclear@2: but WITHOUT ANY WARRANTY; without even the implied warranty of
nuclear@2: MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
nuclear@2: GNU Lesser General Public License for more details.
nuclear@2: 
nuclear@2: You should have received a copy of the GNU Lesser General Public License
nuclear@2: along with this program.  If not, see <http://www.gnu.org/licenses/>.
nuclear@2: */
nuclear@2: 
nuclear@2: #ifndef VMATH_QUATERNION_H_
nuclear@2: #define VMATH_QUATERNION_H_
nuclear@2: 
nuclear@2: #include <stdio.h>
nuclear@2: #include "vmath_types.h"
nuclear@2: #include "vector.h"
nuclear@2: 
nuclear@2: #ifdef __cplusplus
nuclear@2: extern "C" {
nuclear@2: #endif	/* __cplusplus */
nuclear@2: 
nuclear@2: #define quat_cons(s, x, y, z)	v4_cons(x, y, z, s)
nuclear@2: #define quat_vec(q)				v3_cons((q).x, (q).y, (q).z)
nuclear@2: #define quat_s(q)				((q).w)
nuclear@2: #define quat_identity()			quat_cons(1.0, 0.0, 0.0, 0.0)
nuclear@2: void quat_print(FILE *fp, quat_t q);
nuclear@2: 
nuclear@2: #define quat_add		v4_add
nuclear@2: #define quat_sub		v4_sub
nuclear@2: #define quat_neg		v4_neg
nuclear@2: 
nuclear@2: static inline quat_t quat_mul(quat_t q1, quat_t q2);
nuclear@2: 
nuclear@2: static inline quat_t quat_conjugate(quat_t q);
nuclear@2: 
nuclear@2: #define quat_length		v4_length
nuclear@2: #define quat_length_sq	v4_length_sq
nuclear@2: 
nuclear@2: #define quat_normalize	v4_normalize
nuclear@2: static inline quat_t quat_inverse(quat_t q);
nuclear@2: 
nuclear@2: quat_t quat_rotate(quat_t q, scalar_t angle, scalar_t x, scalar_t y, scalar_t z);
nuclear@2: quat_t quat_rotate_quat(quat_t q, quat_t rotq);
nuclear@2: 
nuclear@2: static inline void quat_to_mat3(mat3_t res, quat_t q);
nuclear@2: static inline void quat_to_mat4(mat4_t res, quat_t q);
nuclear@2: 
nuclear@2: #define quat_lerp quat_slerp
nuclear@2: quat_t quat_slerp(quat_t q1, quat_t q2, scalar_t t);
nuclear@2: 
nuclear@2: 
nuclear@2: #ifdef __cplusplus
nuclear@2: }	/* extern "C" */
nuclear@2: 
nuclear@2: #include <iostream>
nuclear@2: 
nuclear@2: /* Quaternion */
nuclear@2: class Quaternion {
nuclear@2: public:
nuclear@2: 	scalar_t s;
nuclear@2: 	Vector3 v;
nuclear@2: 
nuclear@2: 	Quaternion();
nuclear@2: 	Quaternion(scalar_t s, const Vector3 &v);
nuclear@2: 	Quaternion(scalar_t s, scalar_t x, scalar_t y, scalar_t z);
nuclear@2: 	Quaternion(const Vector3 &axis, scalar_t angle);
nuclear@2: 	Quaternion(const quat_t &quat);
nuclear@2: 
nuclear@2: 	Quaternion operator +(const Quaternion &quat) const;
nuclear@2: 	Quaternion operator -(const Quaternion &quat) const;
nuclear@2: 	Quaternion operator -() const;
nuclear@2: 	Quaternion operator *(const Quaternion &quat) const;
nuclear@2: 
nuclear@2: 	void operator +=(const Quaternion &quat);
nuclear@2: 	void operator -=(const Quaternion &quat);
nuclear@2: 	void operator *=(const Quaternion &quat);
nuclear@2: 
nuclear@2: 	void reset_identity();
nuclear@2: 
nuclear@2: 	Quaternion conjugate() const;
nuclear@2: 
nuclear@2: 	scalar_t length() const;
nuclear@2: 	scalar_t length_sq() const;
nuclear@2: 
nuclear@2: 	void normalize();
nuclear@2: 	Quaternion normalized() const;
nuclear@2: 
nuclear@2: 	Quaternion inverse() const;
nuclear@2: 
nuclear@2: 	void set_rotation(const Vector3 &axis, scalar_t angle);
nuclear@2: 	void rotate(const Vector3 &axis, scalar_t angle);
nuclear@2: 	/* note: this is a totally different operation from the above
nuclear@2: 	 * this treats the quaternion as signifying direction and rotates
nuclear@2: 	 * it by a rotation quaternion by rot * q * rot'
nuclear@2: 	 */
nuclear@2: 	void rotate(const Quaternion &q);
nuclear@2: 
nuclear@2: 	Matrix3x3 get_rotation_matrix() const;
nuclear@2: 
nuclear@2: 	friend Quaternion slerp(const Quaternion &q1, const Quaternion &q2, scalar_t t);
nuclear@2: 
nuclear@2: 	friend std::ostream &operator <<(std::ostream &out, const Quaternion &q);
nuclear@2: };
nuclear@2: 
nuclear@2: Quaternion slerp(const Quaternion &q1, const Quaternion &q2, scalar_t t);
nuclear@2: inline Quaternion lerp(const Quaternion &q1, const Quaternion &q2, scalar_t t);
nuclear@2: #endif	/* __cplusplus */
nuclear@2: 
nuclear@2: #include "quat.inl"
nuclear@2: 
nuclear@2: #endif	/* VMATH_QUATERNION_H_ */