pzelchdiv.h

Go to the documentation of this file.

#ifndef PZELCHDIVHTT
#define PZELCHDIVHTT

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


template<class TSHAPE>
class TPZCompElHDiv : public TPZIntelGen<TSHAPE> {
                
    TPZManVector<int, TSHAPE::NFaces> fSideOrient;
    
    std::list<TPZOneShapeRestraint> fRestraints;

protected:
                void Append(TPZFMatrix<REAL> &u1, TPZFMatrix<REAL> &u2, TPZFMatrix<REAL> &u12);

public:
                    
                TPZCompElHDiv(TPZCompMesh &mesh, TPZGeoEl *gel, int64_t &index);
                
                TPZCompElHDiv(TPZCompMesh &mesh, const TPZCompElHDiv<TSHAPE> &copy);
                
                TPZCompElHDiv(TPZCompMesh &mesh,
                                                                  const TPZCompElHDiv<TSHAPE> &copy,
                                                                  std::map<int64_t,int64_t> & gl2lcConMap,
                                                                  std::map<int64_t,int64_t> & gl2lcElMap);
                
                TPZCompElHDiv();
                
                virtual ~TPZCompElHDiv();
                
                virtual TPZCompEl *Clone(TPZCompMesh &mesh) const  override {
                                return new TPZCompElHDiv<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 TPZCompElHDiv<TSHAPE> (mesh, *this, gl2lcConMap, gl2lcElMap);
                }
                
                virtual void SetCreateFunctions(TPZCompMesh *mesh) override;
                
                virtual void Print(std::ostream &out = std::cout) const override;
                

                
                virtual MElementType Type() override;
                
                virtual int NConnects() const override;
                
                virtual void SetConnectIndex(int i, int64_t connectindex) override;
                
    int RestrainedFace();
    
                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 NFluxShapeF() const;
                
                virtual int NSideConnects(int side) const override;
    
                virtual int SideConnectLocId(int node, int side) const override;
    
                virtual int ConnectSideLocId(int connect) const;
                
                virtual int64_t ConnectIndex(int con) const override;
    
    virtual void AddShapeRestraint(TPZOneShapeRestraint restraint) override
    {
        fRestraints.push_back(restraint);
    }
    
    virtual std::list<TPZOneShapeRestraint> GetShapeRestraints() const override
    {
        return fRestraints;
    }
    
    virtual void ResetShapeRestraints() override
    {
        fRestraints.clear();
    }

    virtual int GetSideOrient(int side) override;
    
    virtual void SetSideOrient(int side, int sideorient) override;
    
                
                virtual void SetIntegrationRule(int ord) override;
                
                virtual void GetInterpolationOrder(TPZVec<int> &ord) override;
                
                virtual int PreferredSideOrder(int iside) override;
                
                /*
     * @brief Sets the preferred interpolation order along a side \n
                 * This method only updates the datastructure of the element
                 * In order to change the interpolation order of an element, use the method PRefine
                 */
                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 override;
                int NShapeContinuous(TPZVec<int> &order);
    
    void FillOrder(TPZVec<int> &order) const ;
                
    virtual int MaxOrder() override;
    
    int SideOrient(int face)
    {
#ifdef PZDEBUG
        if (face < 0 || face >= TSHAPE::NFaces) {
            DebugStop();
        }
#endif
        return fSideOrient[face];
    }
                
                virtual void InitMaterialData(TPZMaterialData &data) override;
    
                virtual void ComputeRequiredData(TPZMaterialData &data,
                                                                                                                                                 TPZVec<REAL> &qsi) override;

                void ComputeShapeIndex(TPZVec<int> &sides, TPZVec<int64_t> &shapeindex);
                
                void FirstShapeIndex(TPZVec<int64_t> &Index) const;
    
                void IndexShapeToVec(TPZVec<int> &VectorSide,TPZVec<std::pair<int,int64_t> > & IndexVecShape, int pressureorder);
                
                void IndexShapeToVec(TPZVec<int> &VectorSide, TPZVec<int> &bilinear, TPZVec<int> &direction, TPZVec<std::pair<int,int64_t> > & IndexVecShape, int pressureorder);
    void IndexShapeToVec2(TPZVec<int> &VectorSide, TPZVec<int> &bilinear, TPZVec<int> &direction, TPZVec<std::pair<int,int64_t> > & IndexVecShape, int pressureorder);

                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;
    
                virtual void Solution( TPZVec<REAL> &qsi,int var,TPZVec<STATE> &sol) override;
                
public:
    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 ComputeSolutionHDiv(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;
    
    virtual void ComputeSolution(TPZVec<REAL> &qsi,
                                 TPZVec<REAL> &normal,
                                 TPZSolVec &leftsol, TPZGradSolVec &dleftsol,TPZFMatrix<REAL> &leftaxes,
                                 TPZSolVec &rightsol, TPZGradSolVec &drightsol,TPZFMatrix<REAL> &rightaxes) override
    {
        DebugStop();
    }

    
                void ComputeSolutionHDiv(TPZMaterialData &data);
                
                void CreateGraphicalElement(TPZGraphMesh &grafgrid, int dimension) override;
                
                
                TPZTransform<> TransformSideToElement(int side) override;
                
                public:
int ClassId() const override;

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

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

#include "pzcmesh.h"

template<class TSHAPE>
void TPZCompElHDiv<TSHAPE>::SetCreateFunctions(TPZCompMesh* mesh) {
    mesh->SetAllCreateFunctionsHDiv();
}


TPZCompEl *CreateHDivLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivCubeEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPrismEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivPyramEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivTetraEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivBoundPointEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivBoundLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivBoundQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl *CreateHDivBoundTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);

TPZCompEl * CreateRefHDivLinearEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivQuadEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivTriangleEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivCubeEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivPrismEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivPyramEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);
TPZCompEl * CreateRefHDivTetraEl(TPZGeoEl *gel,TPZCompMesh &mesh,int64_t &index);


#endif