TPZMatDualHybridPoisson.h

Go to the documentation of this file.

//---------------------------------------------------------------------------

#ifndef TPZMatDualHybridPoissonH
#define TPZMatDualHybridPoissonH


#include <iostream>
#include "pzdiscgal.h"
#include "pzfmatrix.h"

class TPZMatDualHybridPoisson : public TPZDiscontinuousGalerkin {
    
    protected :
    
    REAL fXf;
    
    REAL fBetaZero;

    int mydim;
    
public:
    

    
    TPZMatDualHybridPoisson(int nummat, REAL f, REAL betaZero);
    TPZMatDualHybridPoisson(int matid);
    
    TPZMatDualHybridPoisson();
    
    TPZMatDualHybridPoisson(const TPZMatDualHybridPoisson &copy);
    
    void SetDimension(int dim){
        mydim = dim;
    }
    
    virtual ~TPZMatDualHybridPoisson() override;
    
    REAL Beta(int p, REAL size) const{
        return p*p*this->fBetaZero/size;
    }
    
    virtual TPZMaterial * NewMaterial() override {
        return new TPZMatDualHybridPoisson(*this);
    }
    
    virtual int Dimension() const override { return mydim;}
    
    virtual int NStateVariables() const override{
        return 1;
    }
    
    virtual void Print(std::ostream & out) override;
    
    virtual std::string Name()  override { return "TPZMatDualHybridPoisson"; }
    
    virtual void Contribute(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) override;
    
    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 int VariableIndex(const std::string &name) override;
    
    virtual int NSolutionVariables(int var) override;
    
    virtual int NFluxes() override { return 2;}
    
public:
    
    virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override;
    
    virtual void Flux(TPZVec<REAL> &x, TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol, TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux) override
    {
        DebugStop();
    }
    
    void Errors(TPZVec<REAL> &x,TPZVec<STATE> &u,
                TPZFMatrix<STATE> &dudx, TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux,
                TPZVec<STATE> &u_exact,TPZFMatrix<STATE> &du_exact,TPZVec<REAL> &values) override;
    
    virtual int NEvalErrors()  override {return 3;}
    
    public:
virtual int ClassId() const override;
    
    virtual void Write(TPZStream &buf, int withclassid) const override {
        DebugStop();
    }
    
    virtual void Read(TPZStream &buf, void *context) override {
        DebugStop();
    }
    
};

#endif