pzelctemp.h

Go to the documentation of this file.

#ifndef PZELCTEMPH
#define PZELCTEMPH

#include "pzintel.h"
#include "pzquad.h"

template<class TSHAPE>
class TPZIntelGen : public TPZInterpolatedElement {
                
protected:
                
    TPZManVector<int64_t,TSHAPE::NSides> fConnectIndexes;
                
                typename TSHAPE::IntruleType fIntRule;
                
public:
                
                TPZIntelGen(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
                
                TPZIntelGen(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index, int nocreate);
                
                TPZIntelGen(TPZCompMesh &mesh, const TPZIntelGen<TSHAPE> &copy);
                
                TPZIntelGen(TPZCompMesh &mesh,
                                                                const TPZIntelGen<TSHAPE> &copy,
                                                                std::map<int64_t,int64_t> & gl2lcConMap,
                                                                std::map<int64_t,int64_t> & gl2lcElMap);
                
                TPZIntelGen();
                
                virtual ~TPZIntelGen();
                
                virtual TPZCompEl *Clone(TPZCompMesh &mesh) const override {
                                return new TPZIntelGen<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 TPZIntelGen<TSHAPE> (mesh, *this, gl2lcConMap, gl2lcElMap);
                }
                
                
                virtual MElementType Type() override;
                
                virtual int NConnects() const override{
                                return fConnectIndexes.size();
                }
                
                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 TSHAPE::NCornerNodes;
                }
                
                virtual int NSideConnects(int side) const override;
                
                virtual int SideConnectLocId(int node, int side) const override;
                
                virtual int64_t ConnectIndex(int node) const  override;
    
                virtual void SetIntegrationRule(int ord) override;
                
                virtual void SetInterpolationOrder(int order);
                
                virtual void GetInterpolationOrder(TPZVec<int> &ord) override;
                
                virtual int PreferredSideOrder(int iside) override;
                
                virtual void SetPreferredOrder(int order) override;
                
                virtual void SetSideOrder(int side, int order) override;
                
    virtual int EffectiveSideOrder(int side) const override;
                virtual int ConnectOrder(int connect) const;

                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 CreateGraphicalElement(TPZGraphMesh &grafgrid, int dimension) override;
                
                TPZTransform<> TransformSideToElement(int side) override;
                
                virtual const TPZIntPoints &GetIntegrationRule() const  override {
    if (this->fIntegrationRule) {
      return *fIntegrationRule;
    }
    else
    {
      return fIntRule;
    }
                }
                
                virtual TPZIntPoints &GetIntegrationRule()  override {
    if (fIntegrationRule) {
      return *fIntegrationRule;
    }
    else
    {
      return fIntRule;
    }
                }
                
                public:
    virtual int ClassId() const override;

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

template<class TSHAPE>
int TPZIntelGen<TSHAPE>::ClassId() const{
    return Hash("TPZIntelGen") ^ TPZInterpolatedElement::ClassId() << 1 ^ TSHAPE().ClassId() << 2;
}

#endif