TAPsSparseMatrix_Matrix3x3.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsSparseMatrix_Matrix3x3.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (12/22/2009)
UPDATE          (04/10/2010)
******************************************************************************/
#ifndef TAPs_SPARSE_MATRIX_MATRIX3X3_HPP
#define TAPs_SPARSE_MATRIX_MATRIX3X3_HPP

#include "TAPsSparseVector_Matrix3x3.hpp"
#include "TAPsSparseSymmetricMatrix_Matrix3x3.hpp"

BEGIN_NAMESPACE_TAPs
//=============================================================================
template <typename T>
class SparseMatrix_Matrix3x3 {
//=============================================================================
public:
    // Member Functions -------------------------------------------------------
    friend std::ostream & operator<< ( std::ostream &output, SparseMatrix_Matrix3x3<T> const &obj )
    {
        output << obj.StrInfo();
        return output;
    }
    //-------------------------------------------------------------------------
    SparseMatrix_Matrix3x3 ( unsigned int numOfRows = 0 );
    SparseMatrix_Matrix3x3 ( SparseMatrix_Matrix3x3<T> const &orig );
    virtual ~SparseMatrix_Matrix3x3 ();
    //-------------------------------------------------------------------------
    virtual std::string StrInfo () const;
    //-------------------------------------------------------------------------
    inline SparseMatrix_Matrix3x3<T> & operator= ( SparseMatrix_Matrix3x3<T> const &orig );

    //inline SparseMatrix_Matrix3x3<T> & operator= ( SparseSymmetricMatrix_Matrix3x3<T> &orig );

    inline void SetTo ( SparseMatrix_Matrix3x3<T> const &orig );

    //-------------------------------------------------------------------------
    
    //-------------------------------------------------------------------------
    // Operations

    inline void MulWith (
        std::vector< Vector3<T> > const &I, 
        std::vector< Vector3<T> > &O        
    ) const;

    inline void MulWith_wSkippedRows (
        std::vector< Vector3<T> > const &I,         
        std::vector< Vector3<T> > &O,               
        std::vector< bool > const &rowsFixStatus,   
        T diagval                                   
    ) const;

    inline void MulWith_wSkippedCols (
        std::vector< Vector3<T> > const &I,         
        std::vector< Vector3<T> > &O,               
        std::vector< bool > const &colsFixStatus,   
        T diagval                                   
    ) const;

    inline void MulWith_wSkippedRowsCols (
        std::vector< Vector3<T> > const &I,             
        std::vector< Vector3<T> > &O,                   
        std::vector< bool > const &rowscolsFixStatus,   
        T diagval                                       
    ) const;

    inline void ChangeNumOfRows ( unsigned int numOfRows );

    //===============================================================
    // Diagonalizes WITH A SPECIFIED VALUE
    //---------------------------------------------------------------
    inline void DiagonalizeRow ( unsigned int row, T diagonalValue );

    inline void DiagonalizeRows ( 
        std::vector< bool > const & rows,   
        T diagonalValue                     
    );

    inline void DiagonalizeRows ( 
        std::vector< unsigned int > const & rows,   
        T diagonalValue                             
    );

    inline void DiagonalizeCol ( unsigned int col, T diagonalValue );

    inline void DiagonalizeCols (
        std::vector< bool > const & cols,   
        T diagonalValue                     
    );

    inline void DiagonalizeRowCol ( unsigned int rowcol, T diagonalValue );

    inline void DiagonalizeRowsCols (
        std::vector< bool > const & rowscols,   
        T diagonalValue                         
    );
    //---------------------------------------------------------------

    //===============================================================
    // Diagonalizes
    //---------------------------------------------------------------
    inline void DiagonalizeRow ( unsigned int row );

    inline void DiagonalizeRows ( 
        std::vector< bool > const & rows    
    );

    inline void DiagonalizeRows ( 
        std::vector< unsigned int > const & rows    
    );

    inline void DiagonalizeCol ( unsigned int col );

    inline void DiagonalizeCols (
        std::vector< bool > const & cols    
    );

    inline void DiagonalizeRowCol ( unsigned int rowcol );

    inline void DiagonalizeRowsCols (
        std::vector< bool > const & rowscols    
    );
    //---------------------------------------------------------------

    inline unsigned int GetNumOfRows () const   { return m_Rows.size(); }

    void SetDirectAccessors ();

    inline SparseVector_Matrix3x3<T> * const DirectAccess ( unsigned int row );

    static void MPCGSolver ( 
        SparseMatrix_Matrix3x3<T> const &A,     
        std::vector< Vector3<T> > const &x,     
        std::vector< Vector3<T> > const &b,     
        std::vector< Vector3<T> > &x_out,       
        std::vector< Vector3<T> > const &ft,    
        std::vector< Vector3<T> > const &P,     
        int iterationLimit = 100,               
        T errorThreshold = 1E-3                 
    );

    //-------------------------------------------------------------------------

    // Data Members -----------------------------------------------------------
//=============================================================================
protected:
    // Member Functions -------------------------------------------------------

    // Data Members -----------------------------------------------------------
    std::list< SparseVector_Matrix3x3<T> >      m_Rows; 
    std::vector< SparseVector_Matrix3x3<T>* >   m_DirectAccessors;  
//=============================================================================
private:
    // Member Functions -------------------------------------------------------

    static inline std::string DEBUG_ToStr ( std::vector< Vector3<T> > const &V )
    {
        std::stringstream ss;
        ss << "Dense Vector of size "<<V.size()<<":\n";
        for ( unsigned int i = 0; i < V.size(); ++i ) {
            ss << V[i] << "\n";
        }
        return ss.str();
    }

    static inline std::string DEBUG_ToStr ( std::vector< Matrix3x3<T> > const &V )
    {
        std::stringstream ss;
        ss << "Dense Vector of size "<<V.size()<<":\n";
        for ( unsigned int i = 0; i < V.size(); ++i ) {
            ss << V[i] << "\n";
        }
        return ss.str();
    }
    // Data Members -----------------------------------------------------------
//=============================================================================

#if defined(__gl_h_) || defined(__GL_H__)
public:
    //virtual void Draw ();
#endif

//=============================================================================
}; // END CLASS SparseMatrix_Matrix3x3
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsSparseMatrix_Matrix3x3.cpp"

// Include definition if TAPs_USE_EXPORT is not defined
//#if !defined( TAPs_USE_EXPORT )
//  #include "TAPsSparseMatrix_Matrix3x3.cpp"
//#endif
//-----------------------------------------------------------------------------
#endif
//34567890123456789012345678901234567890123456789012345678901234567890123456789
//--+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----