EMM_CanonicalMassMatrix.h

Go to the documentation of this file.

#ifndef EMM_CanonicalMassMatrix_H
#define EMM_CanonicalMassMatrix_H

//interface class
#include "EMM_MassMatrix.h"

//enable the mask vector product
#include "MATH_MaskMatrix.h"
#include "MATH_ArrayVector.h"


DEFINE_SPTR(EMM_CanonicalMassMatrix);
DEFINE_SVPTR(EMM_CanonicalMassMatrix);
class EMM_CanonicalMassMatrix : public  EMM_MassMatrix, public  virtual MATH_MaskMatrix {
  
    SP_OBJECT(EMM_CanonicalMassMatrix);
    // ATTRIBUTES
public:
    
  
    

private:

    //for solving A.Y=B with work array of W
    SP::MATH_MaskVector mY;
    SP::MATH_MaskVector mB;
    SP::MATH_MaskVector mW;

  
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_CanonicalMassMatrix(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_CanonicalMassMatrix(void);

public:
   
    virtual SP::MATH_Matrix NewInstance() const {
        return New();
    }
    
    inline static SP::EMM_CanonicalMassMatrix New() {
        SP::EMM_CanonicalMassMatrix p(new EMM_CanonicalMassMatrix(),
                                      EMM_CanonicalMassMatrix::Delete());
        p->setThis(p);
        return p;
    };
    
    // SET methods 
 
public:
    
    // GET methods
    virtual tUIndex getRowsNumber() const {
        return EMM_MassMatrix::getRowsNumber();
    }
    virtual tUIndex getColumnsNumber() const {
        return EMM_MassMatrix::getColumnsNumber();
    }
public:
    virtual tULLInt getMemorySize() const {
        return MATH_MaskMatrix::getMemorySize()+mW->getMemorySize()+mB->getMemorySize()+mY->getMemorySize();
    }
 
    virtual void setSize(const tUIndex& nRows,const tUIndex& nCols) {
        EMM_MassMatrix::setSize(nRows,nCols);
    }
    
public:
    // OTHERS methods
    virtual tBoolean isSymmetric() const {
        SP::MATH_ArrayVector W=MATH_ArrayVector::New();
        return MATH_Matrix::isSymmetric(1.e-12,*W.get());
    }
    
protected:
    virtual void product(const tUIndex& n,
                         const tReal& alpha,
                         const tUSInt& incX,const tReal* x,
                         const tReal& beta,
                         const tUSInt& incY,tReal* y) const;
    

protected:
    virtual tReal symmetricDot(const tUSInt& incX,const tReal* x,
                               const tUSInt& incY,const tReal* y) const;
    virtual tReal symmetricDot(const tReal *D,
                               const tUSInt& incX,const tReal* X,
                               const tUSInt& incY,const tReal* Y) const;
    

public:

    virtual tBoolean solve(MATH_Solver& solver,EMM_RealField& B);
    
    virtual void product(const tReal& alpha,const EMM_RealField& X,
                         const tReal& beta,EMM_RealField& Y) const;
  
    
    virtual tReal symmetricDot(const EMM_RealField& X) const; 
    virtual tReal symmetricDot(const EMM_RealArray& W,const EMM_RealField& X) const;

    virtual tString toString() const;
    
    
    
  
};

#endif