TPZFrontMatrix.h

Go to the documentation of this file.

#ifndef TPZFRONTMATRIX_H
#define TPZFRONTMATRIX_H


#include "TPZFront.h"
#include "TPZFrontNonSym.h"
#include "TPZFrontSym.h"


#include "TPZStackEqnStorage.h"
#include "pzvec.h"
#include "TPZFileEqnStorage.h"
#include "pzmatrix.h"
#include "pzfmatrix.h"
#include "pzlog.h"

template<class TVar>
class TPZAbstractFrontMatrix : public TPZMatrix<TVar>
{
public:
                
                TPZAbstractFrontMatrix() : TPZMatrix<TVar>()
                {
                }
                
                TPZAbstractFrontMatrix(int64_t ieq, int64_t jeq) : TPZMatrix<TVar>(ieq,jeq)
                {
                }
                
                virtual TPZFront<TVar> & GetFront() = 0;
                
int ClassId() const override;
};

template<class TVar>
int TPZAbstractFrontMatrix<TVar>::ClassId() const{
    return Hash("TPZAbstractFrontMatrix") ^ TPZMatrix<TVar>::ClassId() << 1;
}

template <class TVar, class store, class front>
class TPZFrontMatrix : public TPZAbstractFrontMatrix<TVar> {
protected:
                
                store fStorage;
                
                front fFront;
public:
                int Work();
                void FinishWriting();
                void ReOpen();
                
                virtual int Zero() override;
                void AllocData();
                
                TPZFront<TVar> &GetFront() override { return fFront;}
                void CheckCompress();
                static void main();
                void Print(const char * name, std::ostream & out ,const MatrixOutputFormat form = EFormatted) const override;
    ~TPZFrontMatrix();
int ClassId() const override;
    TPZFrontMatrix();
    TPZFrontMatrix(int64_t globalsize);
                
                TPZFrontMatrix(const TPZFrontMatrix &cp) : TPZRegisterClassId(&TPZFrontMatrix::ClassId),TPZAbstractFrontMatrix<TVar>(cp), fStorage(cp.fStorage),
                fFront(cp.fFront),fNumEq(cp.fNumEq),fLastDecomposed(cp.fLastDecomposed), fNumElConnected(cp.fNumElConnected),fNumElConnectedBackup(cp.fNumElConnectedBackup)
    {
    }
    
  //  CLONEDEF(TPZFrontMatrix)
                virtual TPZMatrix<TVar>*Clone() const  override { return new TPZFrontMatrix(*this); }
    void EquationsToDecompose(TPZVec<int64_t> &destinationindex, int64_t &lower_eq, int64_t &upper_eq);
                
                
                //    void AssembleMatrix(TPZVec < int > & eqnumbers, TPZFMatrix<REAL> & ek, TPZFMatrix<REAL> & ef);
                
    void SetNumElConnected(TPZVec < int > &numelconnected);
                
                virtual void AddKel(TPZFMatrix<TVar> & elmat, TPZVec < int64_t > & destinationindex) override;
                
    void SymbolicAddKel(TPZVec < int64_t > & destinationindex);
                
    virtual void AddKel(TPZFMatrix<TVar> & elmat, TPZVec < int64_t > & sourceindex, TPZVec < int64_t > & destinationindex) override ;
                
    virtual int SolveDirect(TPZFMatrix<TVar> &B ,const DecomposeType dt, std::list<int64_t> &singular) override;
                int Subst_Forward(TPZFMatrix<TVar> *b) const override;
                
                int Subst_Backward(TPZFMatrix<TVar> *b) const override;
                int Substitution(TPZFMatrix<TVar> *) const override;
    /*
                 void SetFileName(
                 const char *name = SetTempFileName() //! Name of the file
                 );
                 */
                
    void SetFileName(char option, const char *name);
                
protected:
                int64_t fNumEq;
                int64_t fLastDecomposed;
                
                //    TPZFront fFront;
                
    //TPZStackEqnStorage fStorage;
                
    TPZVec <int> fNumElConnected;
                
    TPZVec <int> fNumElConnectedBackup;
                
    //TPZFileEqnStorage fFileEqnStorage;
};

#endif //TPZFRONTMATRIX_H