Matrix3D.h

Go to the documentation of this file.

#ifndef CLASSIC_Matrix3D_HH
#define CLASSIC_Matrix3D_HH

// ------------------------------------------------------------------------
// $RCSfile: Matrix3D.h,v $
// ------------------------------------------------------------------------
// $Revision: 1.1.1.1 $
// ------------------------------------------------------------------------
// Copyright: see Copyright.readme
// ------------------------------------------------------------------------
//
// Class: Matrix3D
//
// ------------------------------------------------------------------------
// Class category: BeamlineGeometry
// ------------------------------------------------------------------------
//
// $Date: 2000/03/27 09:32:34 $
// $Author: fci $
//
// ------------------------------------------------------------------------

class Vector3D;


// Class Matrix3D
// ------------------------------------------------------------------------
//  The copy constructor, destructor, and assignment operator generated
//  by the compiler perform the correct operation.  For speed reasons they
//  are not implemented.

class Matrix3D {

public:

    //  Constructs identity.
    Matrix3D();

    //  Use the three vectors (a,b,c) as column vectors.
    Matrix3D(const Vector3D &a, const Vector3D &b, const Vector3D &c);

    //  Use the elements as matrix elements by rows.
    Matrix3D(double x11, double x12, double x13,
             double x21, double x22, double x23,
             double x31, double x32, double x33);

    bool operator==(const Matrix3D &) const;
    bool operator!=(const Matrix3D &) const;

    //  Return reference to matrix element (i,k).
    double &operator()(int i, int k)      ;

    //  Return value of matrix element (i,k).
    double operator()(int i, int k) const;

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

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

    Matrix3D &operator*=(const Matrix3D &rhs);

    Matrix3D &operator*=(double factor)      ;

    static Matrix3D Identity();

    Matrix3D inverse() const;

    //  Return true, if [b]this[/b] is an identity matrix.
    bool isIdentity() const;

    Matrix3D transpose() const;

private:

    // matrix elements
    double m[3][3];
};


// Global operators.
// ------------------------------------------------------------------------

Matrix3D operator+(const Matrix3D &lhs, const Matrix3D &rhs);

Matrix3D operator-(const Matrix3D &lhs, const Matrix3D &rhs);

Matrix3D operator*(const Matrix3D &lhs, const Matrix3D &rhs);

Vector3D operator*(const Matrix3D &lhs, const Vector3D &rhs);


// Inline functions.
// ------------------------------------------------------------------------

inline Matrix3D::Matrix3D() {
    m[0][0] = 1.0;
    m[0][1] = 0.0;
    m[0][2] = 0.0;
    m[1][0] = 0.0;
    m[1][1] = 1.0;
    m[1][2] = 0.0;
    m[2][0] = 0.0;
    m[2][1] = 0.0;
    m[2][2] = 1.0;
}


inline double &Matrix3D::operator()(int i, int k)
{ return m[i][k]; }


inline double Matrix3D::operator()(int i, int k) const
{ return m[i][k]; }

#endif // __Matrix3D_HH