TPZMHMixedMeshControl.h

Go to the documentation of this file.

//
//  TPZMHMixedMeshControl.hpp
//  PZ
//
//  Created by Philippe Devloo on 09/10/16.
//
//

#ifndef TPZMHMixedMeshControl_hpp
#define TPZMHMixedMeshControl_hpp

#include <stdio.h>

#include "TPZMHMeshControl.h"

class TPZMHMixedMeshControl : public TPZMHMeshControl
{
    
protected:
    // this mesh is the same as fCMesh if there are no lagrange multipliers assocated with the average pressure
    TPZAutoPointer<TPZCompMesh> fFluxMesh;
    
    TPZAutoPointer<TPZCompMesh> fRotationMesh;
    

public:
    
    TPZMHMixedMeshControl() : TPZMHMeshControl()
    {
        
    }
    
    TPZMHMixedMeshControl(int dimension);
    //    TPZMHMixedMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh, std::set<int64_t> &coarseindices);
    
    TPZMHMixedMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh, TPZVec<int64_t> &coarseindices);
    
    TPZMHMixedMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh);
    
    
    TPZMHMixedMeshControl(const TPZMHMixedMeshControl &copy) : TPZMHMeshControl(copy)
    {
        fFluxMesh = copy.fFluxMesh;
    }
    
    TPZMHMixedMeshControl &operator=(const TPZMHMixedMeshControl &cp)
    {
        fFluxMesh = cp.fFluxMesh;
        TPZMHMeshControl::operator=(cp);
        return *this;
    }
    
    virtual ~TPZMHMixedMeshControl();
    virtual void InsertPeriferalMaterialObjects();
    
    virtual void InsertPeriferalHdivMaterialObjects();
    
    virtual void InsertPeriferalPressureMaterialObjects();
    
    virtual void InsertPeriferalRotationMaterialObjects();
    
    virtual void BuildComputationalMesh(bool usersubstructure);
    
    TPZAutoPointer<TPZCompMesh> FluxMesh()
    {
        return fFluxMesh;
    }
    
    void SetLagrangeAveragePressure(bool flag)
    {
        if (flag == true) {
            DebugStop();
        }
    }
    
    virtual void Hybridize(bool flag)
    {
        DebugStop();
    }

    void GetMeshVec(TPZVec<TPZCompMesh *> &meshvec)
    {
        meshvec.Resize(2);
        meshvec[0] = fFluxMesh.operator->();
        meshvec[1] = fPressureFineMesh.operator->();
    }

    TPZVec<TPZAutoPointer<TPZCompMesh> > GetMeshes()
    {
        TPZManVector<TPZAutoPointer<TPZCompMesh>,3> result(2);
        result[0] = fFluxMesh;
        result[1] = fPressureFineMesh;
        return result;
    }
    

    void Print(std::ostream &out);

    virtual void PrintFriendly(std::ostream &out)
    {
        
    }
    
protected:
    
   

    void CreateHDivMHMMesh();
    
    virtual void CreatePressureMHMMesh();
    
    virtual void CreateRotationMesh();
    
    // create the elements domain per domain with approximation spaces disconnected from each other
    virtual void CreateInternalFluxElements();
    
    // create the approximation space associated with the skeleton and restrain the connects
    virtual void CreateSkeleton();
    

    void CreateHDivPressureMHMMesh();
    
    void BuildMultiPhysicsMesh();

    virtual void CreateMultiPhysicsInterfaceElements(int dim);
    
    virtual void CreateMultiPhysicsInterfaceElements(int dim, int pressmatid, std::pair<int,int> skelmatid);
    
    void HideTheElements();

    virtual void HybridizeSkeleton(int skeletonmatid, int pressurematid);
    
    void OptimizeInterfaceElements();
    
    virtual void GroupandCondenseElements();
    
    void DeletePressureElements();
};

#endif /* TPZMHMixedMeshControl_hpp */