Matrix3x3.h

Go to the documentation of this file.

#ifndef MATRIX3X3_H_
#define MATRIX3X3_H_

#include "I2D.h"
#include "Point3D.h"
#include "Point6D.h"

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

    Matrix3x3( const double& init );

    Matrix3x3( const double value11, const double value12, const double value13,
               const double value21, const double value22, const double value23,
               const double value31, const double value32, const double value33 );

    Matrix3x3( Point6D pos );

    Str toString() const;

    void transpose();

    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( Point3D& b,  Point3D& x ) const;

    const void mult( Point3D& lhs, Point3D& rhs ) const;

    bool operator==( const Matrix3x3& rhs );

    Matrix3x3& operator=( const Matrix3x3& rhs );

    Matrix3x3& operator=( const double rhs );

    Matrix3x3& operator+=( const Matrix3x3& rhs );

    Matrix3x3& operator-=( const Matrix3x3& rhs );

    Matrix3x3& operator*=( const double rhs );

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

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

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

    virtual double* getPtr1d()
    {
        return value_;
    }

    virtual int getM() const
    {
        return 3;
    }

    virtual int getN() const
    {
        return 3;
    }

protected:

    // The actual data
    double value_[9];
    double* value2d_[3];
};

#endif /*MATRIX3X3_H_*/