pzspblockdiagpivot.h

Go to the documentation of this file.

#ifndef PZSPBLOCKDIAGPIVOT_H
#define PZSPBLOCKDIAGPIVOT_H

#include "tpzsparseblockdiagonal.h"

template<class TVar>
class TPZSpBlockDiagPivot : public TPZSparseBlockDiagonal<TVar>
{
public:
    TPZSpBlockDiagPivot();
                
    ~TPZSpBlockDiagPivot();
    
    virtual int Decompose_LU();
    
                virtual int Decompose_LU(std::list<int64_t> &singular)
                {
        return TPZBlockDiagonal<TVar>::Decompose_LU(singular);
    }
                
                virtual int Substitution( TPZFMatrix<TVar> * B ) const;
                
private:
                TPZVec<int> fPivotIndices;
                
                int Substitution2( TPZFMatrix<TVar> * B ) const;
                
};

#endif