pzmaterialdata.h

Go to the documentation of this file.

#ifndef PZMATERIALDATA_H
#define PZMATERIALDATA_H

#include "pzmanvector.h"
#include "pzfmatrix.h"

#ifdef _AUTODIFF
#include "fad.h"
#endif


typedef TPZManVector<STATE, 10> TPZFemSol;
typedef TPZFNMatrix<15, STATE> TPZFemGradSol;
typedef TPZManVector<TPZFemSol,20> TPZSolVec;
typedef TPZManVector<TPZFemGradSol,20> TPZGradSolVec;


class TPZMaterialData : public TPZSavable {
    
public:
    
    
    enum MShapeFunctionType {EEmpty, EScalarShape, EVecandShape, EVecShape};
    // EScalarShape : regular shape functions: one shape function used for all state variables (default)
    // EVecandShape : HDiv type shape function and a scalar function
    // EVecShape : a vector valued shape function
    
    MShapeFunctionType fShapeType;
    bool fNeedsSol = false, fNeedsNeighborSol = false, fNeedsHSize = false, fNeedsNeighborCenter = false, fNeedsNormalVecFad = false;
    bool fNeedsNormal = false;
    bool fActiveApproxSpace = true;
    TPZFNMatrix<220, REAL> phi;
    TPZFNMatrix<660, REAL> dphi;
    TPZFNMatrix<660, REAL> dphix;
    TPZFNMatrix<220, REAL> divphi;
    TPZFNMatrix<9,REAL> axes;
    TPZFNMatrix<9,REAL> jacobian;
    TPZFNMatrix<9,REAL> jacinv;
    TPZManVector<REAL,3> normal;
    TPZManVector<REAL,3> x;
    TPZManVector<REAL,3> xParametric;
    int p;
    TPZSolVec sol;
    TPZGradSolVec dsol;
    TPZSolVec divsol;
    REAL HSize;
    REAL detjac;
    TPZManVector<REAL,3> XCenter;
    TPZFNMatrix<180> fDirectionsOnMaster;
    
    
    int numberdualfunctions;
    
    //Id of associated geo element
    int gelElId;
    
    TPZManVector<std::pair<int,int64_t> > fVecShapeIndex;
    TPZFNMatrix<180> fNormalVec;
#ifdef _AUTODIFF
    TPZFNMatrix<180,Fad<REAL>> fNormalVecFad;
#endif
    
    int intLocPtIndex;
    
    int intGlobPtIndex;
    
    int NintPts;
    
    TPZMaterialData();
    
    TPZMaterialData( const TPZMaterialData &cp );
    
    ~TPZMaterialData();
    
    void SetAllRequirements(bool set);
    
    //void InvertLeftRightData();
    
    TPZMaterialData &operator= (const TPZMaterialData &cp );
    
    void Print(std::ostream &out) const;
    void PrintMathematica(std::ostream &out) const;
    void Write(TPZStream &buf, int withclassid) const override;
    
    void Read(TPZStream &buf, void *context) override;
    
    virtual bool Compare(TPZSavable *copy, bool override = false) override;
    
    virtual bool Compare(TPZSavable *copy, bool override = false) const override;
    
    void ComputeFluxValues(TPZFMatrix<REAL> & fluxes);
    
    void ComputeFunctionDivergence();
    
public:
    int ClassId() const override;
    
};

#endif