reorderer.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics 
//
//  License:         BSD License 
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Pooyan Dadvand
//                    
//
 
#if !defined(KRATOS_REORDERER_H_INCLUDED )
#define  KRATOS_REORDERER_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
#include <vector>
 
// External includes
 
// Project includes
#include "includes/define.h"
 
namespace Kratos
{
 
 
 
 
 
 
template<class TSparseSpaceType, class TDenseSpaceType>
class Reorderer
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(Reorderer);
 
    typedef typename TSparseSpaceType::MatrixType SparseMatrixType;
 
    typedef typename TSparseSpaceType::VectorType VectorType;
 
    typedef typename TDenseSpaceType::MatrixType DenseMatrixType;
 
    typedef unsigned int IndexType;
 
    typedef unsigned int SizeType;
 
    typedef std::vector<IndexType> IndexVectorType;
 
 
    Reorderer() {}
 
    Reorderer(Reorderer const& rOther) {}
 
    virtual ~Reorderer() {}
 
 
 
    Reorderer& operator=(Reorderer const& rOther)
    {
        return *this;
    }
 
 
 
    virtual void Initialize(SparseMatrixType& rA, VectorType& rX, VectorType& rB)
    {
        CalculateIndexPermutation(rA);
    }
 
    virtual void Reorder(SparseMatrixType& rA, VectorType& rX, VectorType& rB)
    {
    }
 
    virtual void InverseReorder(SparseMatrixType& rA, VectorType& rX, VectorType& rB)
    {
    }
 
    virtual IndexVectorType& CalculateIndexPermutation(SparseMatrixType& rA, IndexType InitialIndex = IndexType())
    {
        SizeType size = TSparseSpaceType::Size1(rA);
 
        if(mIndexPermutation.size() != size)
            mIndexPermutation.resize(size);
 
        for(SizeType i = 0 ; i < size ; i++)
            mIndexPermutation[i] = i;
 
        return mIndexPermutation;
    }
 
 
    virtual IndexVectorType& GetIndexPermutation()
    {
        return mIndexPermutation;
    }
 
 
 
 
    virtual std::string Info() const
    {
        return "Reorderer";
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const
    {
        rOStream << "Reorderer";
    }
 
    virtual void PrintData(std::ostream& rOStream) const
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
    IndexVectorType mIndexPermutation;
 
 
 
 
 
 
 
 
 
 
 
 
 
}; // Class Reorderer
 
 
 
 
 
 
template<class TSparseSpaceType, class TDenseSpaceType>
inline std::istream& operator >> (std::istream& rIStream,
                                  Reorderer<TSparseSpaceType, TDenseSpaceType>& rThis);
 
template<class TSparseSpaceType, class TDenseSpaceType>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const Reorderer<TSparseSpaceType, TDenseSpaceType>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_REORDERER_H_INCLUDED  defined