TPZVecL2.h

Go to the documentation of this file.

#ifndef PZVECL2HPP
#define PZVECL2HPP

#include "TPZMaterial.h"

class  TPZVecL2 : public TPZMaterial
{
    
    
public:
    
    TPZVecL2(int id);
    
    TPZVecL2();
    
    TPZVecL2(const TPZVecL2 &mat);
    
    TPZVecL2 &operator=(const TPZVecL2 &mat);
    
    virtual ~TPZVecL2();
    
    virtual std::string Name()  override { return "TPZVecL2"; }
    
    virtual int Dimension() const  override
    {
        return fDim;
    }
    
    void SetDimension(int dim){
        fDim = dim;
    }
    
    virtual int NStateVariables() const override
    {
        return fNState;
    }
    
    void SetNStateVariables(int nstate)
    {
        fNState = nstate;
    }
    
    virtual int NFluxes()  override {return 0;}
    
                
    virtual void Print(std::ostream &out = std::cout) override;
    
    virtual int VariableIndex(const std::string &name) override;
    
    virtual int NSolutionVariables(int var) override;
    
    virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override;
                
    virtual void Solution(TPZVec<TPZMaterialData> &datavec, int var, TPZVec<STATE> &Solout) override;
                
    virtual void Solution(TPZMaterialData &data, TPZVec<TPZMaterialData> &dataleftvec, TPZVec<TPZMaterialData> &datarightvec, int var, TPZVec<STATE> &Solout) override;
                
    virtual void Solution(TPZMaterialData &data, TPZVec<TPZMaterialData> &dataleftvec, TPZVec<TPZMaterialData> &datarightvec, int var, TPZVec<STATE> &Solout, TPZCompEl * left, TPZCompEl * ritgh) override;
    
protected:
    virtual void Solution(TPZVec<STATE> &Sol,TPZFMatrix<STATE> &DSol,TPZFMatrix<REAL> &axes,int var,TPZVec<STATE> &Solout) override;
    
    int fDim;
    
    int fNState = 1;
    
public:
    
    virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
    
                
    virtual void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
                
    virtual void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ef) override
    {
        TPZFMatrix<STATE> ek(ef.Rows(),ef.Rows(),0.);
        Contribute(datavec, weight, ek, ef);
    }
    
    virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override;
  
                
    virtual void ContributeBC(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override
    {
        
    }
    
    virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ef) override;
    
    virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override;
                
    virtual TPZMaterial * NewMaterial() override;
    
    public:
virtual int ClassId() const override;

    
    virtual void Write(TPZStream &buf, int withclassid) const override;
    
    virtual void Read(TPZStream &buf, void *context) override;
    
    void ErrorsHdiv(TPZMaterialData &data,TPZVec<STATE> &u_exact,TPZFMatrix<STATE> &du_exact,TPZVec<REAL> &values) override;
    
};

#endif