matr3d.h

Go to the documentation of this file.

/*********************************************************************/ 
/*                                                                   */ 
/*             Config files library - (C) TGG 2019                   */ 
/*                                                                   */ 
/*********************************************************************/ 
/*                          Warszawa, 2019                           */ 
/*********************************************************************/ 
/*                                                                   */
/* File: matr3d.h                                                    */
/*                                                                   */
/* Author: F.A.Dul, modif. by T.Grabowski                            */
/*                                                                   */
/* Contents - matrix 3x3 class                                       */
/*                                                                   */ 
/*                                                                   */
/*********************************************************************/
/*                                                                   */
/*                                                                   */
 
#ifndef _MATR_3D_H_
#define _MATR_3D_H_
 
#include "vector3d.h"
 
 
class  MATRIX_3x3
{
    public:
//  row-column indices, e.g.,  xz -->  0,2
    double xx;  
    double xy;  
    double xz;  
    double yx;  
    double yy;  
    double yz;  
    double zx;  
    double zy;  
    double zz;  
 
    MATRIX_3x3( void ){ SetToZero(); };
    MATRIX_3x3( double xx_p , double xy_p , double xz_p ,
                double yx_p , double yy_p , double yz_p ,
                double zx_p , double zy_p , double zz_p ) :
                       xx(xx_p) , xy(xy_p) , xz(xz_p) ,
                       yx(yx_p) , yy(yy_p) , yz(yz_p) ,
                       zx(zx_p) , zy(zy_p) , zz(zz_p)
        {};
    MATRIX_3x3( double M[] ) :
        xx(M[0]) , xy(M[3]) , xz(M[6]) ,
        yx(M[1]) , yy(M[4]) , yz(M[7]) ,
        zx(M[2]) , zy(M[5]) , zz(M[8])
        {};
    MATRIX_3x3( double M[3][3] ) :
        xx(M[0][0]) , xy(M[0][1]) , xz(M[0][2]) ,
        yx(M[1][0]) , yy(M[1][1]) , yz(M[1][2]) ,
        zx(M[2][0]) , zy(M[2][1]) , zz(M[2][2])
        {};
    void GetFrom ( double xx_p , double xy_p , double xz_p ,
                  double yx_p , double yy_p , double yz_p ,
                  double zx_p , double zy_p , double zz_p )
        {
        xx = xx_p;   xy = xy_p;    xz = xz_p;
        yx = yx_p;   yy = yy_p;    yz = yz_p;
        zx = zx_p;   zy = zy_p;    zz = zz_p;
        };
    void  GetFrom( double A[] )
        {
        xx = A[0];   xy = A[1];    xz = A[2];
        yx = A[3];   yy = A[4];    yz = A[5];
        zx = A[6];   zy = A[7];    zz = A[8];
        };
    void  GetFrom( double A[3][3] )
        {
        xx = A[0][0];   xy = A[0][1];   xz = A[0][2];
        yx = A[1][0];   yy = A[1][1];   yz = A[1][2];
        zx = A[2][0];   zy = A[2][1];   zz = A[2][2];
        };
    void  SetToZero( void )
        {
        xx = 0.0;  xy = 0.0;  xz = 0.0;
        yx = 0.0;  yy = 0.0;  yz = 0.0;
        zx = 0.0;  zy = 0.0;  zz = 0.0;
        };
    void  Unit( void )
        {
        xx = 1.0;  xy = 0.0;  xz = 0.0;
        yx = 0.0;  yy = 1.0;  yz = 0.0;
        zx = 0.0;  zy = 0.0;  zz = 1.0;
        };
    double Det( void )
        {
        return xz*yy*zx + xx*yz*zy + xy*yx*zz;
        };
};
 
 
MATRIX_3x3  operator  *( const MATRIX_3x3 &A , const MATRIX_3x3 &B );
MATRIX_3x3  operator  &( const MATRIX_3x3 &A , const MATRIX_3x3 &B );
MATRIX_3x3  operator  !( const MATRIX_3x3 &A                       );
MATRIX_3x3  operator  ~( const MATRIX_3x3 &A                       );
MATRIX_3x3  operator  *( const MATRIX_3x3 &A , const double     &s );
MATRIX_3x3  operator  *( const double     &s , const MATRIX_3x3 &A );
MATRIX_3x3  operator  /( const MATRIX_3x3 &A , const double     &s );
MATRIX_3x3  operator  +( const MATRIX_3x3 &A , const MATRIX_3x3 &B );
MATRIX_3x3  operator  -( const MATRIX_3x3 &A , const MATRIX_3x3 &B );
MATRIX_3x3  operator  -( const MATRIX_3x3 &A                       );
void        operator /=(       MATRIX_3x3 &A , const double     &s );
 
VECTOR_3D   operator  *( const MATRIX_3x3 &A , const VECTOR_3D  &b );
VECTOR_3D   operator  *( const VECTOR_3D  &b , const MATRIX_3x3 &A );
VECTOR_3D   operator  /( const VECTOR_3D  &b ,       MATRIX_3x3 &A );
 
#endif /*_MATR_3D_H_*/