LAP_DoubleMatrixStorage.h

Go to the documentation of this file.

#ifndef LAP_MatrixStorage_H
#define LAP_MatrixStorage_H

#include "LAP_Matrix.h"
#include "LAP_DoubleVector.h"


DEFINE_SPTR(LAP_DoubleMatrixStorage);

class LAP_DoubleMatrixStorage: public LAP_Object{
    
    SP_OBJECT(LAP_DoubleMatrixStorage);
    
public:
    // ATTRIBUTES
    
    
   
    // to define zero value 
   

private:
    // ASSOCIATIONS

  

    SP::LAP_DoubleVector mValues;
    // METHODS
    
  
protected:
    // CONSTRUCTORS
    LAP_DoubleMatrixStorage() {
        mValues=LAP_DoubleVector::New();
    }
  
    
    
    
    
    // DESTRUCTORS
    virtual ~LAP_DoubleMatrixStorage() {
    }
    
    
    
public:
 
    //-----------
    // OPERATORS
    // ----------
public:
    
    inline const double& operator[](const tLVectorIndex& index) const {
        return (*mValues.get())[index];
    };
    inline double& operator[](const tLVectorIndex& index)  {
        return (*mValues.get())[index];
    };
    virtual double& operator()(const tLVectorIndex& i,const tLVectorIndex& j)=0;
    
    virtual const double& operator()(const tLVectorIndex& i,const tLVectorIndex& j) const=0;
    
    //copy method
    //===========
    
    virtual tBoolean copy(const tBoolean& isSymmetric,
                          const tBoolean& isUpper,
                          const tBoolean& isTransposed,
                          const LAP_DoubleMatrixStorage& x) {
        try {
            getValues().copy(x.getValues());
        } catch(std::exception e) {
            cout << e.what()<<"\n";
            return false;
        }
        return true;
    }
    
   
  
    // ----------------
    // MATRIX SETTINGS
    // ----------------

public:
  
    
public:

   
    


    //values settings
    //===============
    void init(const double& v) {
        LAP_DoubleVector* vs=mValues.get();
        if (vs!=null) vs->init(v);
    }
    
    virtual void setValuesPointer(SP::LAP_DoubleVector v) {
        if (v.get()==null)
            throw CORE_Exception("math/linalg/core",
                                 "LAP_DoubleMatrixStorage::setValuesPointer()",
                                 "can not set a null values of a matrix");
        
        mValues=v;
    };
    
    virtual void setValues(const tLVectorIndex& n,
                           const double* values) {
        getValues().setValues(n,values);
    }

    // dimension of values settings
    // ===========================
    virtual void setValuesNumber(const tLVectorIndex& n) {
        getValues().setSize(n);
    };

    
 
    // GET
    // ===
    virtual tLVectorIndex getNullValuesNumber() const {
        tLVectorIndex i,n=mValues->getSize();
        const double *vs=&(*mValues.get())[0];
        tLVectorIndex s=0;
        for (i=0;i<n;i++) {
            if (fabs(*vs)<ZERO_EPSILON) s++;
            vs++;
        }
        return s;
    }
    
    inline tLVectorIndex getValuesNumber() const {
        return mValues->getSize();
    }
    inline tLVectorIndex getCapacity() const {
        return mValues->getCapacity();
    }
    
    inline double getFillRate() const {
        return  ((int) ((1.-((double)getNullValuesNumber())/((double)getValuesNumber()))*100000))/1000.;
    }
   
  
   
    inline const LAP_DoubleVector& getValues() const {
        return *mValues.get();
    }
    inline  LAP_DoubleVector& getValues()  {
        return *mValues.get();
    };

    /* \brief get the vector values of the matrix as a shared pointer
     */
    virtual void getValuesPointer(SP::LAP_DoubleVector& v) {
        v=mValues;
    }
  
   
  
    
    //computation method
    virtual double trace() const=0;

    virtual void addDiagonal(const double& alpha)=0;

    virtual pair<tLVectorIndex,tLVectorIndex> computeBandsNumber() const=0;
  
   
   
};

#endif