pzl2projection.h

Go to the documentation of this file.

#ifndef PZL2PROJECTION_H
#define PZL2PROJECTION_H

#include "TPZMaterial.h"
#include "pzdiscgal.h"

class TPZL2Projection : public TPZDiscontinuousGalerkin{
                
protected:
                
                int fDim;
                
                int fNStateVars;
                
                TPZVec<STATE> fSol;
                
                bool fIsReferred;
                
                int fIntegrationOrder;
    
    REAL fScale;
                
public:
                
                TPZL2Projection(int id, int dim, int nstate, TPZVec<STATE> &sol,
                                                                                int IntegrationOrder = -1);
                
                ~TPZL2Projection();
                
                TPZL2Projection(const TPZL2Projection &cp);
                
                enum ESolutionVars { ENone = 0, ESolution = 1 , EDerivative = 2};
                
    virtual int IntegrationRuleOrder(int elPMaxOrder) const override ;
                
    void SetScaleFactor(REAL scale)
    {
        fScale = scale;
    }
    
                virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override ;
    virtual void ContributeVecShape(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef);
                
                virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ef) override
                {
                                TPZDiscontinuousGalerkin::Contribute(data,weight,ef);
                }
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override {
                                //NOTHING TO BE DONE HERE
                }
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {
                                //NOTHING TO BE DONE HERE
                }
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ef) override {
                                //NOTHING TO BE DONE HERE
                }
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, REAL weight, TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {
                                //NOTHING TO BE DONE HERE
                }
                
                virtual int Dimension() const  override { return this->fDim; }
    
    virtual void SetDimension(int dim) { this->fDim = dim; }
                
                virtual int NStateVariables() const override { return this->fNStateVars; }
    
                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) ;
                }
                
                void SetIsReferred(bool val);
                
                virtual TPZMaterial * NewMaterial() override ;
                
                virtual int VariableIndex(const std::string &name) override ;
                
                virtual int NSolutionVariables(int var) override ;
                
protected:
                virtual void Solution(TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol,
                                                                                                  TPZFMatrix<REAL> &axes, int var, TPZVec<STATE> &Solout) override ;
public:
                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override 
                {
                                TPZMaterial::Solution(data,var,Solout);
                }
                
                virtual void SolutionDisc(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, int var, TPZVec<STATE> &Solout)
                {
                                TPZDiscontinuousGalerkin::SolutionDisc(data,dataleft,dataright,var,Solout);
                }
    
    virtual 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 ;
    public:
virtual int ClassId() const override ;

};

#endif