pzburger.h

Go to the documentation of this file.

#ifndef BURGERH
#define BURGERH

#include "pzpoisson3d.h"
#include "pzpoisson3dreferred.h"

class TPZBurger : public TPZMatPoisson3dReferred {
                
public:
                enum EStabilizationScheme{ESUPG = 1, EGRADIENT = 2};
                
                static int gStabilizationScheme;
                TPZBurger(int nummat, int dim);
                TPZBurger(const TPZBurger &cp);
                virtual ~TPZBurger();
                
                bool IsReferred(){ return this->fIsReferred;}
                
                void SetReferred(bool Is){ this->fIsReferred = Is; }
                
                REAL fSolRef;
                
                virtual void Contribute(TPZMaterialData &data,
                            REAL weight,
                            TPZFMatrix<STATE> &ek,
                            TPZFMatrix<STATE> &ef) override {
        int numbersol = data.sol.size();
        if(numbersol != 1) 
        {
            DebugStop();
        }
                                if (TPZBurger::gStabilizationScheme == ESUPG){
                                                this->ContributeSUPG(data.x,data.jacinv,data.sol[0],data.dsol[0],weight,data.axes,data.phi,data.dphix,ek,ef);
                                }
                                if (TPZBurger::gStabilizationScheme == EGRADIENT){
                                                this->ContributeGradStab(data.x,data.jacinv,data.sol[0],data.dsol[0],weight,data.axes,data.phi,data.dphix,ek,ef);
                                }
                }//Contribute
                
                virtual void Contribute(TPZMaterialData &data,
                            REAL weight,
                                                                                                                TPZFMatrix<STATE> &ef) override
                {
                                TPZMatPoisson3dReferred::Contribute(data,weight,ef);
                }
                
                void ContributeGradStab(TPZVec<REAL> &x,TPZFMatrix<REAL> &jacinv,TPZVec<STATE> &sol,TPZFMatrix<STATE> &dsol,REAL weight,
                                                                                                                TPZFMatrix<REAL> &axes,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef);
                void ContributeSUPG(TPZVec<REAL> &x,TPZFMatrix<REAL> &jacinv,TPZVec<STATE> &sol,TPZFMatrix<STATE> &dsol,REAL weight,
                                                                                                TPZFMatrix<REAL> &axes,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef);
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override;
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                                                                                                                                 REAL weight,
                                                                                                                                                 TPZFMatrix<STATE> &ek,
                                                                                                                                                 TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                                                                                                                                   REAL weight,
                                                                                                                                                   TPZFMatrix<STATE> &ek,
                                                                                                                                                   TPZFMatrix<STATE> &ef,
                                                                                                                                                   TPZBndCond &bc) override;
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override
                {
                                TPZMatPoisson3dReferred::ContributeBC(data,weight,ef,bc);
                }
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                                                                                                                                 REAL weight,
                                                                                                                                                 TPZFMatrix<STATE> &ef) override
                {
                                TPZMatPoisson3dReferred::ContributeInterface(data,dataleft,dataright,weight,ef);
                }
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                                                                                                                                   REAL weight,
                                                                                                                                                   TPZFMatrix<STATE> &ef,
                                                                                                                                                   TPZBndCond &bc) override
                {
                TPZMatPoisson3dReferred::ContributeBCInterface(data,dataleft,weight,ef,bc);
                }
                
public:
int ClassId() const override;

protected:
    bool fIsReferred;

                
};

#endif