/* dtkQuaternion.h ---
 *
 * Author: Thibaud Kloczko
 * Copyright (C) 2008 - Thibaud Kloczko, Inria.
 * Created: Tue Jun  8 15:00:12 2010 (+0200)
 * Version: $Id: d16d73941d5b84c0e9b9f72cb27f9cb850ad1391 $
 * Last-Updated: jeu. mars 26 11:22:43 2015 (+0100)
 *           By: Thibaud Kloczko
 *     Update #: 67
 */

/* Commentary:
 *
 */

/* Change log:
 *
 */

#ifndef DTKQUATERNION_H
#define DTKQUATERNION_H

#include <dtkMathSupportExport.h>
#include "dtkVector.h"
#include "dtkVector3D.h"

// /////////////////////////////////////////////////////////////////
// dtkQuaternion interface
// /////////////////////////////////////////////////////////////////

namespace dtkDeprecated {

    template <class T = double> class dtkQuaternion : public dtkVector<T>
    {
    public:
        dtkQuaternion(void) : dtkVector<T>(4U) {};
        dtkQuaternion(const dtkMatrix<T>& mat): dtkVector<T>(mat) {};
        dtkQuaternion(const dtkVector<T>& vec): dtkVector<T>(vec) {};
        dtkQuaternion(const dtkQuaternion& qtn): dtkVector<T>(qtn) {};
        dtkQuaternion(const dtkVector3D<T>&, const T& = dtkZero<T>());
        dtkQuaternion(const T&, const T&, const T&, const T&);
        dtkQuaternion(const T&);
        dtkQuaternion(const dtkMatrix<T>&, unsigned, unsigned);
        dtkQuaternion(const dtkVector<T>& vec, unsigned irowStart);
        dtkQuaternion(const T& r, const dtkVector3D<T>& vec, const T& angle);
        ~dtkQuaternion(void) {};

    public:
        QString identifier(void) const;

    public:
        void allocate(void) {
            dtkVector<T>::allocate(4);
        };

        void mapInto(const dtkMatrix<T>& mat,
                     unsigned irowStart,
                     unsigned icolStart);
        void mapInto(const dtkVector<T>& vec, unsigned irowStart);

        unsigned size(void) {
            return 4;
        };

    public:
        dtkQuaternion operator +(const dtkQuaternion&) const;
        dtkQuaternion operator -(const dtkQuaternion&) const;
        dtkQuaternion operator -(void) const;
        dtkQuaternion operator *(const T&) const;
        dtkQuaternion operator *(const dtkQuaternion&) const;
        dtkQuaternion operator /(const T& value) const {
            T tTmp = dtkUnity<T>();
            tTmp /= value;
            return (*this) * tTmp;
        }
        dtkQuaternion operator /(const dtkQuaternion&) const;

        dtkQuaternion& operator  =(const dtkQuaternion&);
        dtkQuaternion& operator  =(const dtkVector3D<T>&);
        dtkQuaternion& operator +=(const dtkQuaternion&);
        dtkQuaternion& operator -=(const dtkQuaternion&);
        dtkQuaternion& operator *=(const T&);
        dtkQuaternion& operator *=(const dtkQuaternion&);
        dtkQuaternion& operator /=(const T&);
        dtkQuaternion& operator /=(const dtkQuaternion&);

        void storeProduct(const dtkQuaternion&, const dtkQuaternion&);
        void storeConjugate(const dtkQuaternion&);

        void makeConjugate(void);

        dtkQuaternion unit(void) const {
            return (*this) / (this->norm());
        }
    };

// /////////////////////////////////////////////////////////////////
//
// /////////////////////////////////////////////////////////////////

    template <class T> class dtkUnity< dtkQuaternion<T> >: public dtkQuaternion<T>
    {
    public:
        dtkUnity(void) {
            (*this)[0] = (*this)[1] = (*this)[2] = dtkZero<T>();
            (*this)[3] = dtkUnity<T>();
        }
    };

// /////////////////////////////////////////////////////////////////
//
// /////////////////////////////////////////////////////////////////

    template <class T> class dtkZero< dtkQuaternion<T> >: public dtkQuaternion<T>
    {
    public:
        dtkZero(void) {
            fill(dtkZero<T>());
        }
    };

} // end of namespace

// /////////////////////////////////////////////////////////////////
// Implementation of the template class dtkQuaternion's methods
// /////////////////////////////////////////////////////////////////

#include "dtkQuaternion.tpp"

// /////////////////////////////////////////////////////////////////
//
// /////////////////////////////////////////////////////////////////

#include <QtCore>

typedef dtkDeprecated::dtkQuaternion<double> dtkQuaternionReal;

Q_DECLARE_METATYPE(dtkQuaternionReal);
Q_DECLARE_METATYPE(dtkQuaternionReal *);

#endif

// /////////////////////////////////////////////////////////////////
// Credits
// /////////////////////////////////////////////////////////////////

// $Id: d16d73941d5b84c0e9b9f72cb27f9cb850ad1391 $

// CwMtx matrix and vector math library
// Copyright (C) 1999-2001  Harry Kuiper
// Copyright (C) 2000  Will DeVore (template conversion)
// Copyright (C) 2000  Jiri Ecer (constructor from exponential form)

// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2 of the License, or (at your option) any later version.

// This library 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
// Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA