pzhdivpressure.h

Go to the documentation of this file.

#ifndef PZHDIVPRESSUREHTT
#define PZHDIVPRESSUREHTT

#include "pzelchdiv.h"


template<class TSHAPE>
class TPZCompElHDivPressure : public TPZCompElHDiv<TSHAPE> {
                
                int fPressureOrder;
                
                void Append(TPZFMatrix<REAL> &u1, TPZFMatrix<REAL> &u2, TPZFMatrix<REAL> &u12);
public:
                
                TPZCompElHDivPressure(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
                
                TPZCompElHDivPressure(TPZCompMesh &mesh, const TPZCompElHDivPressure<TSHAPE> &copy);
                
                TPZCompElHDivPressure(TPZCompMesh &mesh,
                                                                                                  const TPZCompElHDivPressure<TSHAPE> &copy,
                                                                                                  std::map<int64_t,int64_t> & gl2lcConMap,
                                                                                                  std::map<int64_t,int64_t> & gl2lcElMap);
                
                TPZCompElHDivPressure();
                
                virtual ~TPZCompElHDivPressure();
                
                virtual TPZCompEl *Clone(TPZCompMesh &mesh) const override {
                                return new TPZCompElHDivPressure<TSHAPE> (mesh, *this);
                }
                
                virtual TPZCompEl *ClonePatchEl(TPZCompMesh &mesh,std::map<int64_t,int64_t> & gl2lcConMap,std::map<int64_t,int64_t>&gl2lcElMap) const override
                {
                                return new TPZCompElHDivPressure<TSHAPE> (mesh, *this, gl2lcConMap, gl2lcElMap);
                }
                
#ifndef STATE_COMPLEX
  virtual void SetCreateFunctions(TPZCompMesh *mesh) override {
                                mesh->SetAllCreateFunctionsHDivPressure();
                }
#endif
                
                virtual MElementType Type() override;
                
                virtual int NConnects() const override;
                
                virtual void SetConnectIndex(int i, int64_t connectindex) override;
                
                virtual int NConnectShapeF(int connect, int order) const override;
                
                virtual int Dimension() const override {
                                return TSHAPE::Dimension;
                }
                
                virtual int NCornerConnects() const override {
                                return 0;
                }
                
                virtual int64_t ConnectIndex(int node) const override;
    
    virtual int PressureConnectIndex() const override
    {
        return NConnects()-1;
    }
                
                virtual void SetPressureOrder(int ord);
                int DualOrder();
                
                virtual void GetInterpolationOrder(TPZVec<int> &ord) override;
                
                virtual void SetPreferredOrder(int order) override;
                
                virtual int ConnectOrder(int connect) const override;
                
                virtual void InitMaterialData(TPZMaterialData &data) override;
                
                virtual void SideShapeFunction(int side,TPZVec<REAL> &point,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi) override;
                
                virtual void ShapeDual(TPZVec<REAL> &qsi, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi);
                
                void Shape(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi) override;
    
                virtual void Solution( TPZVec<REAL> &qsi,int var,TPZVec<STATE> &sol) override;
    
private:
                virtual         void ComputeSolution(TPZVec<REAL> &qsi, TPZSolVec &sol, TPZGradSolVec &dsol,TPZFMatrix<REAL> &axes) override;
    
public:
    
                virtual void ComputeSolution(TPZVec<REAL> &qsi, TPZMaterialData &data) override;
                
                void ComputeSolutionPressureHDiv(TPZVec<REAL> &qsi, TPZMaterialData &data);
                virtual void ComputeSolution(TPZVec<REAL> &qsi, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphix,
                                 const TPZFMatrix<REAL> &axes, TPZSolVec &sol, TPZGradSolVec &dsol) override;
                
                void ComputeSolutionPressureHDiv(TPZMaterialData &data);
                
                
                void CreateGraphicalElement(TPZGraphMesh &grafgrid, int dimension) override;
                
                TPZTransform<> TransformSideToElement(int side) override;
                
int ClassId() const override;

                void Write(TPZStream &buf, int withclassid) const override;
                
                void Read(TPZStream &buf, void *context) override;
                
};

template<class TSHAPE>
int TPZCompElHDivPressure<TSHAPE>::ClassId() const{
    return Hash("TPZCompElHDivPressure") ^ TPZCompElHDiv<TSHAPE>::ClassId() << 1;
}

TPZCompEl *CreateHDivPressurePointEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressureLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressureQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressureTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressureCubeEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressurePrismEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressurePyramEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPressureTetraEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);

#endif