Matrix6x6.h

Go to the documentation of this file.

#ifndef MATRIX6X6_H_
#define MATRIX6X6_H_

#include "I2D.h"
#include "Point6D.h"
#include "Matrix3x3.h"

class Matrix6x6: public I2D<double>
{
public:
    Matrix6x6();

    Matrix6x6( const double& init );

    Matrix6x6( Matrix3x3& m3x3i11, Matrix3x3& m3x3i12,
               Matrix3x3& m3x3i21, Matrix3x3& m3x3i22 );

    Matrix6x6( Point6D pos );

    // as a Str
    Str toString() const;

    bool determinantZero();

    void invert();

    void decomposeLU();

    // Solve for x in Ax = b, where this is the LU decompostion of
    // M and b is supplied here.
    void solveLU( Point6D& b, Point6D& x ) const;

    //const Point6D operator*( const Point6D& rhs ) const;

    bool operator==( const Matrix6x6& rhs );

    Matrix6x6& operator=( const Matrix6x6& rhs );

    Matrix6x6& operator=( const double rhs );

    void set( int index1, int index2, Matrix3x3& value );

    void set( int index1, int index2, const double value );

    double& operator()( int index1, int index2 );

    void set( const double& v00, const double& v01, const double& v02, const double& v03, const double& v04, const double& v05,
              const double& v10, const double& v11, const double& v12, const double& v13, const double& v14, const double& v15,
              const double& v20, const double& v21, const double& v22, const double& v23, const double& v24, const double& v25,
              const double& v30, const double& v31, const double& v32, const double& v33, const double& v34, const double& v35,
              const double& v40, const double& v41, const double& v42, const double& v43, const double& v44, const double& v45,
              const double& v50, const double& v51, const double& v52, const double& v53, const double& v54, const double& v55 );

    virtual double** getPtr2d()
    {
        return value2d_;
    }

    virtual double* getPtr1d()
    {
        return value_;
    }

    virtual int getM() const
    {
        return 6;
    }

    virtual int getN() const
    {
        return 6;
    }

protected:

    // The actual data
    double value_[ 36 ];

    // For I2D operations
    double* value2d_[ 6 ];

};

#endif /*MATRIX6X6_H_*/