Matrix3D.h

Go to the documentation of this file.

#ifndef OPALX_Matrix3D_HH
#define OPALX_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