matrix3.h

Go to the documentation of this file.

#ifndef MatriX3_H
#define MatriX3_H
 
#ifdef _MCD_CHECK
#include "mcd.h"
#endif
 
#include "matrix.h"
#include "vector3.h"
 
namespace ModuleBase
{
 
class Matrix3
{
 
  public:
    double e11, e12, e13, e21, e22, e23, e31, e32, e33;
 
    Matrix3()
    {
        Identity();
    }
 
    Matrix3(const double &r11,
            const double &r12,
            const double &r13,
            const double &r21,
            const double &r22,
            const double &r23,
            const double &r31,
            const double &r32,
            const double &r33);
 
    void Identity(void);
 
    void Zero(void);
 
    double Det(void) const;
 
    Matrix3 Transpose(void) const;
 
    Matrix3 Inverse(void) const;
 
    Matrix3 &operator=(const Matrix3 &m);
 
    Matrix3 &operator+=(const Matrix3 &m);
 
    Matrix3 &operator-=(const Matrix3 &m);
 
    Matrix3 &operator*=(const double &s);
 
    Matrix3 &operator/=(const double &s);
 
    void print(void) const;
 
    ModuleBase::matrix to_matrix(void) const;
};
 
Matrix3 operator+(const Matrix3 &m1, const Matrix3 &m2);
 
Matrix3 operator-(const Matrix3 &m1, const Matrix3 &m2);
 
Matrix3 operator/(const Matrix3 &m, const double &s);
 
Matrix3 operator*(const Matrix3 &m1, const Matrix3 &m2);
 
Matrix3 operator*(const Matrix3 &m, const double &s);
 
Matrix3 operator*(const double &s, const Matrix3 &m);
 
template <typename T> ModuleBase::Vector3<double> operator*(const Matrix3 &m, const ModuleBase::Vector3<T> &u);
 
template <typename T> ModuleBase::Vector3<double> operator*(const ModuleBase::Vector3<T> &u, const Matrix3 &m);
 
bool operator==(const Matrix3 &m1, const Matrix3 &m2);
 
bool operator!=(const Matrix3 &m1, const Matrix3 &m2); // whethor m1 != m2
 
// m*u
template <typename T> ModuleBase::Vector3<double> operator*(const Matrix3 &m, const ModuleBase::Vector3<T> &u)
{
    return ModuleBase::Vector3<double>(m.e11 * u.x + m.e12 * u.y + m.e13 * u.z,
                                       m.e21 * u.x + m.e22 * u.y + m.e23 * u.z,
                                       m.e31 * u.x + m.e32 * u.y + m.e33 * u.z);
}
 
// u*m
template <typename T> ModuleBase::Vector3<double> operator*(const ModuleBase::Vector3<T> &u, const Matrix3 &m)
{
    return ModuleBase::Vector3<double>(u.x * m.e11 + u.y * m.e21 + u.z * m.e31,
                                       u.x * m.e12 + u.y * m.e22 + u.z * m.e32,
                                       u.x * m.e13 + u.y * m.e23 + u.z * m.e33);
}
 
} // namespace ModuleBase
 
#endif // MATRIX3_H