TPZGeoElement.h

Go to the documentation of this file.

#ifndef TPZGEOELEMENTH
#define TPZGEOELEMENTH

#include "pzgeoelrefless.h"

class TPZGeoElSide;
class TPZCompMesh;
class TPZCompEl;
template<class T,int N>
class TPZStack;

template <class TGeo, class TRef>
class TPZGeoElement : public TPZGeoElRefLess<TGeo> {
                
                int64_t fSubEl[TRef::NSubEl];
public:
                typedef TGeo TGeoLoc;
                
public:
                TPZGeoElement();
                TPZGeoElement(int64_t id,TPZVec<int64_t> &nodeindexes,int matind,TPZGeoMesh &mesh);
                TPZGeoElement(TPZVec<int64_t> &nodeindices,int matind,TPZGeoMesh &mesh);
                TPZGeoElement(TGeo &geo, int matind, TPZGeoMesh &mesh);
                TPZGeoElement(TPZVec<int64_t> &nodeindices,int matind,TPZGeoMesh &mesh,int64_t &index);
                
                TPZGeoElement(TPZGeoMesh &DestMesh, const TPZGeoElement &cp);
                
                TPZGeoElement ( TPZGeoMesh &DestMesh,
                                                                   const TPZGeoElement &cp,
                                                                   std::map<int64_t,int64_t> &gl2lcNdIdx,
                                                                   std::map<int64_t,int64_t> &gl2lcElIdx );
                virtual ~TPZGeoElement(){};
                
                public:
int ClassId() const override;

                
                void Read(TPZStream &str, void *context) override;
                
                void Write(TPZStream &str, int withclassid) const override;
                
                virtual TPZGeoEl * Clone(TPZGeoMesh &DestMesh) const override;
                
                virtual TPZGeoEl * ClonePatchEl(TPZGeoMesh &DestMesh,
                                                                                                                                                std::map<int64_t,int64_t> &gl2lcNdIdx,
                                                                                                                                                std::map<int64_t,int64_t> &gl2lcElIdx) const override;
                
                int HasSubElement() const override {return fSubEl[0]!=-1;}
                
                void SetSubElement(int id, TPZGeoEl *el) override;
                
                
                REAL RefElVolume() override;
                
                void MidSideNodeIndex(int side,int64_t &index) const override;
                
                int NSubElements() const override;
                
                int NSideSubElements(int side) const override;
                
                TPZGeoEl *SubElement(int is) const override;
                
                TPZGeoElSide SideSubElement(int side,int position);
                
                TPZTransform<> GetTransform(int side,int son) override;
                
                virtual int FatherSide(int side, int son)  override {
                                return TRef::FatherSide(side,son);
                }
                
                virtual void Divide(TPZVec<TPZGeoEl *> &pv) override;
                
                virtual void GetSubElements2(int side, TPZStack<TPZGeoElSide> &subel) const override;
                
                virtual void ResetSubElements() override;
    
    virtual TPZGeoEl *CreateGeoElement(MElementType type,
                                       TPZVec<int64_t>& nodeindexes,
                                       int matid,
                                       int64_t& index) override
    {
        return this->Mesh()->CreateGeoElement(type,nodeindexes,matid,index,0);
    }
    

                
};

template <class TGeo, class TRef>
int TPZGeoElement<TGeo, TRef> ::ClassId() const{
    return Hash("TPZGeoElement") ^ TPZGeoElRefLess<TGeo>::ClassId() << 1 ^ TRef().ClassId() << 2;
}
#endif