TAPsSparseSymmetricMatrix_Matrix3x3.hpp

Go to the documentation of this file.

/******************************************************************************
TAPsSparseSymmetricMatrix_Matrix3x3.hpp
******************************************************************************/
/******************************************************************************
SUKITTI PUNAK   (12/30/2009)
UPDATE          (04/10/2010)
******************************************************************************/
#ifndef TAPs_SPARSE_SYMMETRIC_MATRIX_MATRIX3X3_HPP
#define TAPs_SPARSE_SYMMETRIC_MATRIX_MATRIX3X3_HPP

#include "TAPsSparseVector_Matrix3x3.hpp"

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

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

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

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

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

    inline void ChangeNumOfRows ( unsigned int numOfRows );

    inline void DiagonalizeRowCol ( unsigned int ii, T diagonalValue = 1 );

    inline void DiagonalizeRowsCols ( 
        std::vector< bool > const & ii, 
        T diagonalValue = 1             
    );

    inline void DiagonalizeRowsCols ( 
        std::vector< unsigned int > const & ii, 
        T diagonalValue = 1                     
    );

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

    void SetDirectAccessors ();

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

    static void MPCGSolver ( 
        SparseSymmetricMatrix_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-10                
    );

    static void MPCGSolverSpecificForFEM (
        std::vector< bool > const &fixedNodes,  
        SparseSymmetricMatrix_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-10                
    );
    //-------------------------------------------------------------------------

    // 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 SparseSymmetricMatrix_Matrix3x3
//=============================================================================
END_NAMESPACE_TAPs
//-----------------------------------------------------------------------------
#include "TAPsSparseSymmetricMatrix_Matrix3x3.cpp"

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