TPZHybridizeHDiv.h

Go to the documentation of this file.

//
//  TPZHybridizeHDiv.hpp
//  ErrorEstimate
//
//  Created by Philippe Devloo on 16/05/18.
//

#ifndef TPZHybridizeHDiv_hpp
#define TPZHybridizeHDiv_hpp

#include <stdio.h>
#include <map>

template<class T>
class TPZVec;

class TPZCompMesh;
class TPZMaterial;
class TPZCompElSide;
class TPZInterpolatedElement;
class TPZMultiphysicsCompMesh;

struct TPZHybridizeHDiv {
    
    // boundary condition for flux with zero contribution
    int fHDivWrapMatid = -10;
    // material id of the lagrange interface (either pressure or displacement)
    int fLagrangeInterface = -9;
    // material id of the interface elements
    int fInterfaceMatid = -8;
    // number of state variables
    int fNState = 1;
    
    TPZHybridizeHDiv() = default;
    
    TPZHybridizeHDiv(TPZVec<TPZCompMesh *> &meshvec_Hybrid);
    
    void ComputeNState(TPZVec<TPZCompMesh*>& meshvec_Hybrid);

    void ComputePeriferalMaterialIds(TPZVec<TPZCompMesh *> &meshvec_Hybrid);
    
    void SetPeriferalMaterialIds(int HDivWrapMatid, int LagrangeInterface, int InterfaceMatid)
    {
        fHDivWrapMatid = HDivWrapMatid;
        fLagrangeInterface = LagrangeInterface;
        fInterfaceMatid = InterfaceMatid;
        
    }
    
    bool IsPeriferalMaterialId(int matid)
    {
        return matid == fHDivWrapMatid || matid == fLagrangeInterface || matid == fInterfaceMatid;
    }
    void HybridizeInternalSides(TPZVec<TPZCompMesh *> &meshvec_Hybrid);

    void CreateInterfaceElements(TPZCompMesh *cmesh_Hybrid, TPZVec<TPZCompMesh *> &meshvec_Hybrid);

    void CreateInterfaceElements(TPZMultiphysicsCompMesh *cmesh_Hybrid);
    
    TPZCompMesh * CreateMultiphysicsMesh(TPZCompMesh *cmesh_HDiv, TPZVec<TPZCompMesh *> &meshvec_Hybrid, double Lagrange_term_multiplier = 1.);
    
    TPZCompMesh * CreateMultiphysicsMesh(TPZMultiphysicsCompMesh *cmesh_HDiv, double Lagrange_term_multiplier = 1.);

    void ReCreateMultiphysicsMesh(TPZMultiphysicsCompMesh *cmesh_HDiv, double Lagrange_term_multiplier = 1.);
    
    // elementgroup[el] = index of the element with which the element should be grouped
    // this method only gives effective result for hybridized hdiv meshes
    static void AssociateElements(TPZCompMesh *cmesh, TPZVec<int64_t> &elementgroup);
    
    static void GroupandCondenseElements(TPZCompMesh *cmesh_Hybrid);
    
    void InsertPeriferalMaterialObjects(TPZVec<TPZCompMesh *> &meshvec_Hybrid);
    
    void InsertPeriferalMaterialObjects(TPZCompMesh *cmesh_Hybrid, double Lagrange_term_multiplier = 1.);
    
    std::tuple<TPZCompMesh *, TPZVec<TPZCompMesh *> > Hybridize(TPZCompMesh *cmesh_Multiphysics, TPZVec<TPZCompMesh *> &meshvec_HDiv, bool group_elements=true, double Lagrange_term_multiplier = 1.);
    
    TPZMultiphysicsCompMesh *Hybridize(TPZMultiphysicsCompMesh *multiphysics, bool group_elements=true, double Lagrange_term_multiplier = 1.);
    
    void HybridizeGivenMesh(TPZMultiphysicsCompMesh &multiphysics, bool group_elements=true, double Lagrange_term_multiplier = 1.);
    
    static void VerifySolutionConsistency(TPZCompMesh *fluxmesh, std::ostream &out);
private:
    
    std::tuple<int64_t,int> SplitConnects(const TPZCompElSide &left, const TPZCompElSide &right, TPZVec<TPZCompMesh *> &meshvec_Hybrid);

public:
    
    static TPZCompElSide RightElement(TPZInterpolatedElement *intel, int side);

};

#endif /* TPZHybridizeHDiv_hpp */