glc_matrix4x4.cpp

Go to the documentation of this file.

/****************************************************************************

 This file is part of the GLC-lib library.
 Copyright (C) 2005-2008 Laurent Ribon (laumaya@users.sourceforge.net)
 http://glc-lib.sourceforge.net

 GLC-lib 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 3 of the License, or
 (at your option) any later version.

 GLC-lib 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 GLC-lib; if not, write to the Free Software
 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA

*****************************************************************************/

#include "glc_matrix4x4.h"

#include <QtDebug>

// Set Functions

GLC_Matrix4x4& GLC_Matrix4x4::fromEuler(const double angle_x, const double angle_y, const double angle_z)
{
    const double A= cos(angle_x);
    const double B= sin(angle_x);
    const double C= cos(angle_y);
    const double D= sin(angle_y);
    const double E= cos(angle_z);
    const double F= sin(angle_z);

    const double AD= A * D;
    const double BD= B * D;

    m_Matrix[0]  = C * E;
    m_Matrix[4]  = -C * F;
    m_Matrix[8]  = -D;
    m_Matrix[1]  = -BD * E + A * F;
    m_Matrix[5]  = BD * F + A * E;
    m_Matrix[9]  = -B * C;
    m_Matrix[2]  = AD * E + B * F;
    m_Matrix[6]  = -AD * F + B * E;
    m_Matrix[10] = A * C;

    m_Matrix[12]=  0.0; m_Matrix[13]= 0.0; m_Matrix[14]= 0.0; m_Matrix[3]= 0.0; m_Matrix[7]= 0.0; m_Matrix[11] = 0.0;
    m_Matrix[15] =  1.0;

        return *this;
}

GLC_Matrix4x4& GLC_Matrix4x4::setColumn(int index, const GLC_Vector3d& vector)
{
        Q_ASSERT(index < 4);
        index= index * 4;
        m_Matrix[index]= vector.x();
        m_Matrix[index + 1]= vector.y();
        m_Matrix[index + 2]= vector.z();

        m_Type= General;

        return *this;
}
// Private services function

QVector<double> GLC_Matrix4x4::toEuler(void) const
{
        double angle_x;
        double angle_y;
        double angle_z;
        double tracex, tracey;
        angle_y= -asin(m_Matrix[8]);
        double C= cos(angle_y);

        if (!qFuzzyCompare(C, 0.0)) // Gimball lock?
        {
                tracex= m_Matrix[10] / C;
                tracey= - m_Matrix[9] / C;
                angle_x= atan2( tracey, tracex);

                tracex= m_Matrix[0] / C;
                tracey= - m_Matrix[4] / C;
                angle_z= atan2( tracey, tracex);
        }
        else // Gimball lock?
        {
                angle_x= 0.0;
                tracex= m_Matrix[5] / C;
                tracey= m_Matrix[1] / C;
                angle_z= atan2( tracey, tracex);
        }
        QVector<double> result;
        result.append(fmod(angle_x, 2.0 * glc::PI));
        result.append(fmod(angle_y, 2.0 * glc::PI));
        result.append(fmod(angle_z, 2.0 * glc::PI));

        return result;
}

QString GLC_Matrix4x4::toString() const
{
        QString result;
        for (int i= 0; i < DIMMAT4X4; ++i)
        {
                result+= (QString::number(m_Matrix[0 + i])) + QString(" ");
                result+= (QString::number(m_Matrix[4 + i])) + QString(" ");
                result+= (QString::number(m_Matrix[8 + i])) + QString(" ");
                result+= (QString::number(m_Matrix[12 + i])) + QString("\n");
        }
        result.remove(result.size() - 1, 1);
        return result;
}