TPZFrontNonSym.h

Go to the documentation of this file.

/* Generated by Together */
template<class TVar>
class TPZEqnArray;

#ifndef TPZFRONTNONSYM_H
#define TPZFRONTNONSYM_H


#include "pzmatrix.h"
#include "pzstack.h"
#include "pzvec.h"
#include "TPZFront.h"

#include <math.h>
#include <iostream>
#include <fstream>
#include "TPZStackEqnStorage.h"
#include "TPZFileEqnStorage.h"

#ifdef USING_BLAS
#ifdef USING_MKL
#include <mkl.h>
#elif MACOSX
#include <Accelerate/Accelerate.h>
#else
#ifdef MACOSX
#include <Accelerate/Accelerate.h>
#else
extern "C"{
#include "cblas.h"
                //#include "g2c.h"
                //#include "fblaswr.h"
};
#endif
#endif
#endif

#ifdef USING_ATLAS
extern "C"{
#include <cblas.h>
                //#include "g2c.h"
                //#include "fblaswr.h"
};
#endif


template<class TVar>
class TPZFrontNonSym : public TPZFront<TVar> {
public:
                std::string GetMatrixType();
                
                static void main();
    ~TPZFrontNonSym();
    TPZFrontNonSym();
                TPZFrontNonSym(int64_t GlobalSize);
                
                TPZFrontNonSym(const TPZFrontNonSym &cp);
    
    virtual void SetDecomposeType(DecomposeType dectype) override
    {
        if (dectype == ELU) {
            this->fDecomposeType = dectype;
        }
        else
        {
            DebugStop();
        }
    }
    

    void DecomposeEquations(int64_t mineq, int64_t maxeq, TPZEqnArray<TVar> & result);
                
    void SymbolicDecomposeEquations(int64_t mineq, int64_t maxeq);
                
                void SymbolicAddKel(TPZVec < int64_t > & destinationindex);
                
    void Compress();
                
                void Expand(int largefrontsize);
                
                TVar & Element(int64_t i, int64_t j) {
                                return this->fData[this->fMaxFront*j + i];
                }
    
    const TVar & Element(int64_t i, int64_t j) const{
        return this->fData[this->fMaxFront*j + i];
    }
    
    void AddKel(TPZFMatrix<TVar> &elmat, TPZVec<int64_t> &destinationindex);
                
    virtual void AddKel(TPZFMatrix<TVar> &elmat, TPZVec<int64_t> &sourceindex,  TPZVec<int64_t> &destinationindex);
                
                virtual void ExtractFrontMatrix(TPZFMatrix<TVar> &front) override;
        
        public:
int ClassId() const override;
                
private:    
                
    void DecomposeOneEquation(int64_t ieq, TPZEqnArray<TVar> &eqnarray);
                
    void FreeGlobal(int64_t global);
    int Local(int64_t global);
                
public:
                virtual int64_t NFree() override;
                void Reset(int64_t GlobalSize=0);
                void AllocData();
                
                void Print(const char *name, std::ostream& out) const;
                void PrintGlobal(const char *name, std::ostream& out);
                
                
private:
                
    /*#  TPZStackEqnStorage lnkTPZStackEqnStorage; */
                
    /*#  TPZFileEqnStorage lnkTPZFileEqnStorage; */
                
public:
                
                virtual void TensorProductIJ(int ithread, typename TPZFront<TVar>::STensorProductMTData *data) override;
                
};

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


#endif //TPZFRONTNONSYM_H