pznonlinbiharmonic.h

Go to the documentation of this file.

#ifndef TPZNONLINBIHARMONICHPP
#define TPZNONLINBIHARMONICHPP

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

class TPZNonLinBiharmonic : public TPZDiscontinuousGalerkin {
                
private:
                STATE  fXf;
                
                public :
                
                static STATE gLambda1, gLambda2, gSigmaA,gSigmaB, gL_alpha, gM_alpha, gL_betta,
                gM_betta, g_teta, Re;
                static int NorP;
                
                TPZNonLinBiharmonic(int nummat, STATE f);
                
                virtual ~TPZNonLinBiharmonic();
                
                virtual int NEvalErrors()  override {return 8;}
                
                void SetMaterial(STATE &xfin){
                                fXf = xfin;
                }
                
                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 "TPZBiharmonic"; }

                virtual void Contribute(TPZMaterialData &data,
                            REAL weight,
                            TPZFMatrix<STATE> &ek,
                            TPZFMatrix<STATE> &ef) override;
                virtual void Contribute(TPZMaterialData &data,
                            REAL weight,
                                                                                                                TPZFMatrix<STATE> &ef) override
                {
                                TPZDiscontinuousGalerkin::Contribute(data,weight,ef);
                }
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override;
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override
                {
                                TPZDiscontinuousGalerkin::ContributeBC(data,weight,ef,bc);
                }

                virtual int VariableIndex(const std::string &name) override;
                
                virtual int NSolutionVariables(int var) override;
                
                virtual int NFluxes() override { return 0;}
                
protected:
                virtual void Solution(TPZVec<STATE> &Sol,TPZFMatrix<STATE> &DSol,TPZFMatrix<REAL> &axes,int var,TPZVec<STATE> &Solout) override;
public:
                
                virtual void SolutionDisc(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, int var, TPZVec<STATE> &Solout)
                {
                                TPZDiscontinuousGalerkin::SolutionDisc(data,dataleft,dataright,var,Solout);
                }
                
                virtual void Flux(TPZVec<REAL> &x, TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol, TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux) 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 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 ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                                                                                                                                 REAL weight,
                                                                                                                                                 TPZFMatrix<STATE> &ef) override
                {
                                TPZDiscontinuousGalerkin::ContributeInterface(data,dataleft,dataright,weight,ef);
                }
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                                                                                                                                   REAL weight,
                                                                                                                                                   TPZFMatrix<STATE> &ef,
                                                                                                                                                   TPZBndCond &bc) override
                {
                                TPZDiscontinuousGalerkin::ContributeBCInterface(data,dataleft,weight,ef,bc);
                }
    public:
virtual int ClassId() const override;
 
};

#endif