pzhdivpressurebound.h

Go to the documentation of this file.

//
//  pzhdivpressurebound.h
//  PZ
//
//  Created by Agnaldo Farias on 25/06/13.
//
//

#ifndef __PZ__pzhdivpressurebound__
#define __PZ__pzhdivpressurebound__

#include <iostream>

#include "pzelchdivbound2.h"

template<class TSHAPE>
class TPZCompElHDivPressureBound : public TPZCompElHDivBound2<TSHAPE> {

                int fPressureOrder;
    
public:
    
    TPZCompElHDivPressureBound(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
    
    
                TPZCompElHDivPressureBound();
    
                TPZCompElHDivPressureBound(TPZCompMesh &mesh, const TPZCompElHDivPressureBound<TSHAPE> &copy, std::map<int64_t,int64_t> & gl2lcConMap, std::map<int64_t,int64_t> & gl2lcElMap);
    
    
    TPZCompElHDivPressureBound(TPZCompMesh &mesh, const TPZCompElHDivPressureBound<TSHAPE> &copy);
   
    
                virtual ~TPZCompElHDivPressureBound();
                
                virtual TPZCompEl *Clone(TPZCompMesh &mesh) const  override {
                                return new TPZCompElHDivPressureBound<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 TPZCompElHDivPressureBound<TSHAPE> (mesh, *this, gl2lcConMap, gl2lcElMap);
                }

    
    
    virtual int NConnects() const override;
    
                virtual void SetPressureOrder(int ord);
    
    virtual void SetConnectIndex(int i, int64_t connectindex) override;
    
    virtual int NConnectShapeF(int connect, int order) const override;
    
    virtual int SideConnectLocId(int node, int side) const override;
    
                virtual void SetCreateFunctions(TPZCompMesh *mesh) override;
                
    //          virtual MElementType Type();
                
                
                virtual int Dimension() const  override {
                                return TSHAPE::Dimension;
                }
                
                virtual int NCornerConnects() const  override {
                                return 0;
                }
                
//              virtual int NSideConnects(int side) const;
                
                virtual void GetInterpolationOrder(TPZVec<int> &ord) override;
                
    //          virtual int PreferredSideOrder(int iside);
                
                virtual void SetSideOrder(int side, int order) override;
                
                virtual int EffectiveSideOrder(int side) const override;
                
                virtual int ConnectOrder(int connect) const override;
    
                virtual void InitMaterialData(TPZMaterialData &data) override;
                
                void ComputeShapeIndex(TPZVec<int> &sides, TPZVec<int64_t> &shapeindex);
                
                void FirstShapeIndex(TPZVec<int64_t> &Index);
                
                virtual void SideShapeFunction(int side,TPZVec<REAL> &point,TPZFMatrix<REAL> &phi,TPZFMatrix<REAL> &dphi) override;
                
                void Shape(TPZVec<REAL> &pt, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi) override;
                
                void IndexShapeToVec(TPZVec<int> &fVectorSide,TPZVec<std::pair<int,int64_t> > & IndexVecShape);
    
                public:
int ClassId() const override;

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

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

#endif /* defined(__PZ__pzhdivpressurebound__) */