TPZLinearConvecDiff.h

Go to the documentation of this file.

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

#ifndef TPZLinearConvecDiffH
#define TPZLinearConvecDiffH
#include <iostream>
#include "pzfmatrix.h"
#include "TPZMaterial.h"


class TPZLinearConvecDiff : public TPZMaterial {

                protected :

                STATE fXf;

                STATE fK;

                REAL fConvDir[2];

                STATE fSD;

public:

                TPZLinearConvecDiff(int nummat, REAL k, const TPZVec<REAL> &conv, REAL f, REAL SD);

  TPZLinearConvecDiff(int matid);

                TPZLinearConvecDiff();
                
                TPZLinearConvecDiff(const TPZLinearConvecDiff &c);

                virtual ~TPZLinearConvecDiff();

                virtual TPZMaterial * NewMaterial() override {
                                return new TPZLinearConvecDiff(*this);
                }

                virtual int Dimension() const  override { return 2;}

                virtual int NStateVariables() const override { return 1; }

                virtual void Print(std::ostream & out) override;

                virtual std::string Name()  override { return "TPZLinearConvecDiff"; }

                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 int VariableIndex(const std::string &name) override;

                virtual int NSolutionVariables(int var) override;

                virtual int NFluxes() override { return 2;}

                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override;

                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