pzelchdivbound2.h

Go to the documentation of this file.

#ifndef PZELCHDIVBOUND_H_
#define PZELCHDIVBOUND_H_

#include "pzelctemp.h"
#include "TPZOneShapeRestraint.h"

template<class TSHAPE>
class TPZCompElHDivBound2 : public TPZIntelGen<TSHAPE> {
    
    int fSideOrient;
                
                void Append(TPZFMatrix<REAL> &u1, TPZFMatrix<REAL> &u2, TPZFMatrix<REAL> &u12);
    
    TPZCompElSide fneighbour;
    
    TPZOneShapeRestraint fRestraint;
public:
                
                TPZCompElHDivBound2(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
                
                TPZCompElHDivBound2(TPZCompMesh &mesh, const TPZCompElHDivBound2<TSHAPE> &copy);

                TPZCompElHDivBound2(TPZCompMesh &mesh,
                                                                                                const TPZCompElHDivBound2<TSHAPE> &copy,
                                                                                                std::map<int64_t,int64_t> & gl2lcConMap,
                                                                                                std::map<int64_t,int64_t> & gl2lcElMap);
                
                TPZCompElHDivBound2();
                virtual ~TPZCompElHDivBound2();
                
                virtual TPZCompEl *Clone(TPZCompMesh &mesh) const override {
                                return new TPZCompElHDivBound2<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 TPZCompElHDivBound2<TSHAPE> (mesh, *this, gl2lcConMap, gl2lcElMap);
                }
                
                virtual void SetCreateFunctions(TPZCompMesh *mesh) override;
                
                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 int NSideConnects(int side) const override;
                
                virtual int SideConnectLocId(int node, int side) const override;
    
    
    virtual void SetSideOrient(int side, int sideorient) override;
    virtual int GetSideOrient(int side) override;
                
                
                virtual void GetInterpolationOrder(TPZVec<int> &ord) override;
                
                virtual int PreferredSideOrder(int iside) override;
                
                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 ComputeShape(TPZVec<REAL> &intpoint, TPZMaterialData &data) override;
    
    void ComputeShape(TPZVec<REAL> &intpoint, TPZVec<REAL> &X, TPZFMatrix<REAL> &jacobian, TPZFMatrix<REAL> &axes,
                      REAL &detjac, TPZFMatrix<REAL> &jacinv, TPZFMatrix<REAL> &phi, TPZFMatrix<REAL> &dphi, TPZFMatrix<REAL> &dphidx) 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);
    
    virtual void AddShapeRestraint(TPZOneShapeRestraint restraint) override
    {
        if (fRestraint.IsInitialized()) {
            std::cout << "****************** Overwriting a constraint\n";
        }
        fRestraint = restraint;
    }
    

    virtual std::list<TPZOneShapeRestraint> GetShapeRestraints() const override
    {
        std::list<TPZOneShapeRestraint> loc;
        if (fRestraint.IsInitialized()) {
            loc.push_back(fRestraint);
        }
        return loc;
    }

    virtual void ResetShapeRestraints() override
    {
        fRestraint = TPZOneShapeRestraint();
    }

                public:
int ClassId() const override;

                void Write(TPZStream &buf, int withclassid) const override;
                
                void Read(TPZStream &buf, void *context) override;
    
    virtual void Print(std::ostream &out) const override;

};

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

TPZCompEl *CreateRefHDivBoundPointEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateRefHDivBoundLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateRefHDivBoundQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateRefHDivBoundTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);

#endif /* PZELCHDIVBOUND_H_ */