[ViewVC] Diff of: OpenMD/trunk/src/math/Quaternion.hpp

Comparing trunk/src/math/Quaternion.hpp (file contents):
Revision 92 by tim, Sat Oct 16 01:31:28 2004 UTC vs.
Revision 110 by tim, Tue Oct 19 21:28:55 2004 UTC

+#ifndef MATH_QUATERNION_HPP
+#define MATH_QUATERNION_HPP
-+
+#include "math/Vector.hpp"
-+
+#include "math/SquareMatrix.hpp"
-+
+namespace oopse{
+    /**
+     * @class Quaternion Quaternion.hpp "math/Quaternion.hpp"
-<
+     * @brief
->
+     * Quaternion is a sort of a higher-level complex number.
->
+     * It is defined as Q = w + x*i + y*j + z*k,
->
+     * where w, x, y, and z are numbers of type T (e.g. double), and
->
+     * i*i = -1; j*j = -1; k*k = -1;
->
+     * i*j = k; j*k = i; k*i = j;
+     */
+    template<typename Real>
+    class Quaternion : public Vector<Real, 4> {
-+
+        public:
-+
+            Quaternion() : Vector<Real, 4>() {}
-<
+    };
->
+            /** Constructs and initializes a Quaternion from w, x, y, z values */
->
+            Quaternion(Real w, Real x, Real y, Real z) {
->
+                data_[0] = w;
->
+                data_[1] = x;
->
+                data_[2] = y;
->
+                data_[3] = z;
->
+            }
->
->
+            /** Constructs and initializes a Quaternion from a  Vector<Real,4> */
->
+            Quaternion(const Vector<Real,4>& v)
->
+                : Vector<Real, 4>(v){
->
+            }
-+
+            /** copy assignment */
-+
+            Quaternion& operator =(const Vector<Real, 4>& v){
-+
+                if (this == & v)
-+
+                    return *this;
-+
-+
+                Vector<Real, 4>::operator=(v);
-+
-+
+                return *this;
-+
+            }
-+
-+
+            /**
-+
+             * Returns the value of the first element of this quaternion.
-+
+             * @return the value of the first element of this quaternion
-+
+             */
-+
+            Real w() const {
-+
+                return data_[0];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the reference of the first element of this quaternion.
-+
+             * @return the reference of the first element of this quaternion
-+
+             */
-+
+            Real& w() {
-+
+                return data_[0];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the value of the first element of this quaternion.
-+
+             * @return the value of the first element of this quaternion
-+
+             */
-+
+            Real x() const {
-+
+                return data_[1];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the reference of the second element of this quaternion.
-+
+             * @return the reference of the second element of this quaternion
-+
+             */
-+
+            Real& x() {
-+
+                return data_[1];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the value of the thirf element of this quaternion.
-+
+             * @return the value of the third element of this quaternion
-+
+             */
-+
+            Real y() const {
-+
+                return data_[2];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the reference of the third element of this quaternion.
-+
+             * @return the reference of the third element of this quaternion
-+
+             */
-+
+            Real& y() {
-+
+                return data_[2];
-+
+            }
-+
-+
+            /**
-+
+             * Returns the value of the fourth element of this quaternion.
-+
+             * @return the value of the fourth element of this quaternion
-+
+             */
-+
+            Real z() const {
-+
+                return data_[3];
-+
+            }
-+
+            /**
-+
+             * Returns the reference of the fourth element of this quaternion.
-+
+             * @return the reference of the fourth element of this quaternion
-+
+             */
-+
+            Real& z() {
-+
+                return data_[3];
-+
+            }
-+
-+
+            /**
-+
+             * Tests if this quaternion is equal to other quaternion
-+
+             * @return true if equal, otherwise return false
-+
+             * @param q quaternion to be compared
-+
+             */
-+
+             inline bool operator ==(const Quaternion<Real>& q) {
-+
-+
+                for (unsigned int i = 0; i < 4; i ++) {
-+
+                    if (!equal(data_[i], q[i])) {
-+
+                        return false;
-+
+                    }
-+
+                }
-+
-+
+                return true;
-+
+            }
-+
-+
+            /**
-+
+             * Returns the inverse of this quaternion
-+
+             * @return inverse
-+
+             * @note since quaternion is a complex number, the inverse of quaternion
-+
+             * q = w + xi + yj+ zk is inv_q = (w -xi - yj - zk)/(|q|^2)
-+
+             */
-+
+            Quaternion<Real> inverse() {
-+
+                Quaternion<Real> q;
-+
+                Real d = this->lengthSquare();
-+
-+
+                q.w() = w() / d;
-+
+                q.x() = -x() / d;
-+
+                q.y() = -y() / d;
-+
+                q.z() = -z() / d;
-+
-+
+                return q;
-+
+            }
-+
-+
+            /**
-+
+             * Sets the value to the multiplication of itself and another quaternion
-+
+             * @param q the other quaternion
-+
+             */
-+
+            void mul(const Quaternion<Real>& q) {
-+
+                Quaternion<Real> tmp(*this);
-+
-+
+                data_[0] = (tmp[0]*q[0]) -(tmp[1]*q[1]) - (tmp[2]*q[2]) - (tmp[3]*q[3]);
-+
+                data_[1] = (tmp[0]*q[1]) + (tmp[1]*q[0]) + (tmp[2]*q[3]) - (tmp[3]*q[2]);
-+
+                data_[2] = (tmp[0]*q[2]) + (tmp[2]*q[0]) + (tmp[3]*q[1]) - (tmp[1]*q[3]);
-+
+                data_[3] = (tmp[0]*q[3]) + (tmp[3]*q[0]) + (tmp[1]*q[2]) - (tmp[2]*q[1]);
-+
+            }
-+
-+
+            void mul(const Real& s) {
-+
+                data_[0] *= s;
-+
+                data_[1] *= s;
-+
+                data_[2] *= s;
-+
+                data_[3] *= s;
-+
+            }
-+
-+
+            /** Set the value of this quaternion to the division of itself by another quaternion */
-+
+            void div(Quaternion<Real>& q) {
-+
+                mul(q.inverse());
-+
+            }
-+
-+
+            void div(const Real& s) {
-+
+                data_[0] /= s;
-+
+                data_[1] /= s;
-+
+                data_[2] /= s;
-+
+                data_[3] /= s;
-+
+            }
-+
-+
+            Quaternion<Real>& operator *=(const Quaternion<Real>& q) {
-+
+                mul(q);
-+
+                return *this;
-+
+            }
-+
-+
+            Quaternion<Real>& operator *=(const Real& s) {
-+
+                mul(s);
-+
+                return *this;
-+
+            }
-+
-+
+            Quaternion<Real>& operator /=(Quaternion<Real>& q) {
-+
+                *this *= q.inverse();
-+
+                return *this;
-+
+            }
-+
-+
+            Quaternion<Real>& operator /=(const Real& s) {
-+
+                div(s);
-+
+                return *this;
-+
+            }
-+
+            /**
-+
+             * Returns the conjugate quaternion of this quaternion
-+
+             * @return the conjugate quaternion of this quaternion
-+
+             */
-+
+            Quaternion<Real> conjugate() {
-+
+                return Quaternion<Real>(w(), -x(), -y(), -z());
-+
+            }
-+
-+
+            /**
-+
+             * Returns the corresponding rotation matrix (3x3)
-+
+             * @return a 3x3 rotation matrix
-+
+             */
-+
+            SquareMatrix<Real, 3> toRotationMatrix3() {
-+
+                SquareMatrix<Real, 3> rotMat3;
-+
-+
+                Real w2;
-+
+                Real x2;
-+
+                Real y2;
-+
+                Real z2;
-+
-+
+                if (!isNormalized())
-+
+                    normalize();
-+
-+
+                w2 = w() * w();
-+
+                x2 = x() * x();
-+
+                y2 = y() * y();
-+
+                z2 = z() * z();
-+
-+
+                rotMat3(0, 0) = w2 + x2 - y2 - z2;
-+
+                rotMat3(0, 1) = 2.0 * ( x() * y() + w() * z() );
-+
+                rotMat3(0, 2) = 2.0 * ( x() * z() - w() * y() );
-+
-+
+                rotMat3(1, 0) = 2.0 * ( x() * y() - w() * z() );
-+
+                rotMat3(1, 1) = w2 - x2 + y2 - z2;
-+
+                rotMat3(1, 2) = 2.0 * ( y() * z() + w() * x() );
-+
-+
+                rotMat3(2, 0) = 2.0 * ( x() * z() + w() * y() );
-+
+                rotMat3(2, 1) = 2.0 * ( y() * z() - w() * x() );
-+
+                rotMat3(2, 2) = w2 - x2 -y2 +z2;
-+
-+
+                return rotMat3;
-+
+            }
-+
-+
+    };//end Quaternion
-+
-+
-+
+    /**
-+
+     * Returns the vaule of scalar multiplication of this quaterion q (q * s).
-+
+     * @return  the vaule of scalar multiplication of this vector
-+
+     * @param q the source quaternion
-+
+     * @param s the scalar value
-+
+     */
-+
+    template<typename Real, unsigned int Dim>
-+
+    Quaternion<Real> operator * ( const Quaternion<Real>& q, Real s) {
-+
+        Quaternion<Real> result(q);
-+
+        result.mul(s);
-+
+        return result;
-+
+    }
-+
-+
+    /**
-+
+     * Returns the vaule of scalar multiplication of this quaterion q (q * s).
-+
+     * @return  the vaule of scalar multiplication of this vector
-+
+     * @param s the scalar value
-+
+     * @param q the source quaternion
-+
+     */
-+
+    template<typename Real, unsigned int Dim>
-+
+    Quaternion<Real> operator * ( const Real& s, const Quaternion<Real>& q ) {
-+
+        Quaternion<Real> result(q);
-+
+        result.mul(s);
-+
+        return result;
-+
+    }
-+
-+
+    /**
-+
+     * Returns the multiplication of two quaternion
-+
+     * @return the multiplication of two quaternion
-+
+     * @param q1 the first quaternion
-+
+     * @param q2 the second quaternion
-+
+     */
-+
+    template<typename Real>
-+
+    inline Quaternion<Real> operator *(const Quaternion<Real>& q1, const Quaternion<Real>& q2) {
-+
+        Quaternion<Real> result(q1);
-+
+        result *= q2;
-+
+        return result;
-+
+    }
-+
-+
+    /**
-+
+     * Returns the division of two quaternion
-+
+     * @param q1 divisor
-+
+     * @param q2 dividen
-+
+     */
-+
-+
+    template<typename Real>
-+
+    inline Quaternion<Real> operator /( Quaternion<Real>& q1,  Quaternion<Real>& q2) {
-+
+        return q1 * q2.inverse();
-+
+    }
-+
-+
+    /**
-+
+     * Returns the value of the division of a scalar by a quaternion
-+
+     * @return the value of the division of a scalar by a quaternion
-+
+     * @param s scalar
-+
+     * @param q quaternion
-+
+     * @note for a quaternion q, 1/q = q.inverse()
-+
+     */
-+
+    template<typename Real>
-+
+    Quaternion<Real> operator /(const Real& s,  Quaternion<Real>& q) {
-+
-+
+        Quaternion<Real> x;
-+
+        x = q.inverse();
-+
+        x *= s;
-+
+        return x;
-+
+    }
-+
-+
+    template <class T>
-+
+    inline bool operator==(const Quaternion<T>& lhs, const Quaternion<T>& rhs) {
-+
+        return equal(lhs[0] ,rhs[0]) && equal(lhs[1] , rhs[1]) && equal(lhs[2], rhs[2]) && equal(lhs[3], rhs[3]);
-+
+    }
-+
-+
+    typedef Quaternion<double> Quat4d;
+#endif //MATH_QUATERNION_HPP

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing trunk/src/math/Quaternion.hpp (file contents): Revision 92 by tim, Sat Oct 16 01:31:28 2004 UTC vs. Revision 110 by tim, Tue Oct 19 21:28:55 2004 UTC

Diff Legend

Comparing trunk/src/math/Quaternion.hpp (file contents):
Revision 92 by tim, Sat Oct 16 01:31:28 2004 UTC vs.
Revision 110 by tim, Tue Oct 19 21:28:55 2004 UTC