IvfMat4d.h

//
// Copyright 1999-2006 by Structural Mechanics, Lund University.
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Library 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
// Library General Public License for more details.
//
// You should have received a copy of the GNU Library 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.
//
// Please report all bugs and problems to "ivf@byggmek.lth.se".
//

#ifndef _CIvfMat4d_h_
#define _CIvfMat4d_h_

#include <ivf/ivfconfig.h>
#include <ivfmath/IvfMathBase.h>
#include <ivfmath/IvfVec4d.h>

class IVFMATH_API CIvfMat4d : public CIvfMathBase {
private:
        double m_matrix[4][4];
        int LU_Decompose  (double matrix[4][4], int* index);
        void LU_Backsub  (double matrix[4][4], int* index, double* B);
public:
        CIvfMat4d();

        CIvfMat4d(CIvfMat4d& M);

        CIvfMat4d(double scalar);

        virtual ~CIvfMat4d();

        void one();

        void identity();

        void zero();

        void translate(double tx, double ty, double tz);

        void rotateX(double angle);

        void rotateY(double angle);

        void rotateZ(double angle);

        void rotate(double vx, double vy, double vz, double theta);

        void scaling(double sx, double sy, double sz);

#ifdef IVF_HAVE_IOSTREAM

        void print(std::ostream &out);
#else

        void print(ostream &out);
#endif

        void setComponent(int row, int col, double value);

        double getComponent(int row, int col);

        void setColumn(int col, double v1, double v2, double v3, double v4);

        void setRow(int row, double v1, double v2, double v3, double v4);

        CIvfMat4d& t();

        CIvfMat4d& inv();

        CIvfMat4d& operator=(CIvfMat4d& A);

        CIvfMat4d& operator=(double a);
        
        CIvfMat4d& operator+=(CIvfMat4d& A);
        
        CIvfMat4d& operator-=(CIvfMat4d& A);
        
        CIvfMat4d& operator*(CIvfMat4d& A);
        
        CIvfVec4d& operator*(CIvfVec4d& v);
        
        //CIvfVec4d operator*(CIvfVec4d& v);

        CIvfMat4d& operator*(double scalar);
        
        CIvfMat4d& operator/(double scalar);
        
        CIvfMat4d& operator+(CIvfMat4d& A);
};

IVFMATH_API CIvfMat4d& ivfGetTempMatrix4d();

#endif 

---