Mesh_decl.hpp

#ifndef Mesh_decl_hpp
#define Mesh_decl_hpp
 
#define FULL_Mesh_TIMER
 
#include "feddlib/core/FEDDCore.hpp"
#include "feddlib/core/General/DefaultTypeDefs.hpp"
#include "feddlib/core/LinearAlgebra/MultiVector.hpp"
#include "feddlib/core/FE/Elements.hpp"
#include "feddlib/core/FE/FiniteElement.hpp"
#include "feddlib/core/FE/Helper.hpp"
#include "feddlib/core/Mesh/AABBTree.hpp"

 
namespace FEDD {
template <class SC = default_sc, class LO = default_lo, class GO = default_go, class NO = default_no>
class Mesh {
    
public:
    typedef Elements Elements_Type;
    typedef FiniteElement FiniteElement_Type;
    typedef Teuchos::RCP<FiniteElement_Type> FiniteElementPtr_Type;
    typedef Teuchos::RCP<Elements_Type> ElementsPtr_Type;
    
    typedef Teuchos::RCP<Mesh> Mesh_ptr_Type;
    
    typedef Teuchos::RCP<Teuchos::Comm<int> > CommPtr_Type;
    typedef Teuchos::RCP<const Teuchos::Comm<int> > CommConstPtr_Type;
    typedef const CommConstPtr_Type CommConstPtrConst_Type;
    
    typedef Map<LO,GO,NO> Map_Type;
    typedef Teuchos::RCP<Map_Type> MapPtr_Type;
    typedef Teuchos::RCP<const Map_Type> MapConstPtr_Type;
    typedef Teuchos::RCP<const Map_Type> MapConstPtrConst_Type;
 
    typedef MultiVector<SC,LO,GO,NO> MultiVector_Type;
    typedef Teuchos::RCP<MultiVector_Type> MultiVectorPtr_Type;
 
    typedef AABBTree<SC,LO,GO,NO> AABBTree_Type;
    typedef Teuchos::RCP<AABBTree_Type > AABBTreePtr_Type;
    
    /* ###################################################################### */
    //
    Mesh();
 
    Mesh(CommConstPtrConst_Type &comm);
 
    virtual ~Mesh() = default;

    void deleteData();
    
    void setParameterList( ParameterListPtr_Type& pL );

    ParameterListConstPtr_Type getParameterList( ) const;
    
    vec_int_ptr_Type getElementsFlag() const;
    
    MapConstPtr_Type getMapUnique() const;
    
    MapConstPtr_Type getMapRepeated() const;

    MapConstPtr_Type getElementMap() const;
    
    MapConstPtr_Type getEdgeMap(); // Edge Map
    
    vec2D_dbl_ptr_Type getPointsRepeated() const;

    vec2D_dbl_ptr_Type getPointsUnique() const;
    
    vec_int_ptr_Type getBCFlagRepeated() const;
    
    vec_int_ptr_Type getBCFlagUnique() const;

    ElementsPtr_Type getElementsC();

    ElementsPtr_Type getSurfaceElements();
    
    int getDimension();
    
    GO getNumElementsGlobal();
    
    LO getNumElements();
    
    LO getNumPoints(std::string type="Unique");
    
    int getOrderElement();

    CommConstPtrConst_Type getComm(){return comm_;};
    
    int setStructuredMeshFlags(int flags){return 0;};
    
    void setElementFlags(std::string type="");
    
    void setReferenceConfiguration();
    
    void moveMesh( MultiVectorPtr_Type displacementUnique, MultiVectorPtr_Type displacementRepeated );
    
    int getSurfaceElementOrder() { return surfaceElementOrder_; };

    int getEdgeElementOrder() { return edgesElementOrder_; };
 
 
    // Creates an AABBTree from own vertice- and elementlist.
    void create_AABBTree();
    
    vec_int_ptr_Type  findElemsForPoints(vec2D_dbl_ptr_Type query_points);
    
    vec_dbl_Type getBaryCoords(vec_dbl_Type point, int element);
    
    bool isPointInElem(vec_dbl_Type point, int element);

    tuple_intint_Type getRankRange() const {return rankRange_;};
    
    void deleteSurfaceElements(){ surfaceElements_.reset(); };

    void correctNormalDirections();

    void correctElementOrientation();    

    vec2D_int_ptr_Type getElements();
    
    /* ###################################################################### */
    
    int                     dim_;
    long long               numElementsGlob_;
 
    std::string             FEType_;
    MapPtr_Type             mapUnique_;
    MapPtr_Type             mapRepeated_;
    vec2D_dbl_ptr_Type      pointsRep_;
    vec2D_dbl_ptr_Type      pointsUni_;
    vec_int_ptr_Type        bcFlagRep_;
    vec_int_ptr_Type        bcFlagUni_;
 
    ElementsPtr_Type        surfaceElements_;
 
    ElementsPtr_Type        elementsC_;
    MapPtr_Type             elementMap_;
    MapPtr_Type             edgeMap_;
 
 
    CommConstPtrConst_Type  comm_;
    
    vec2D_int_ptr_Type  elementsVec_;
    
    vec2D_dbl_ptr_Type      pointsRepRef_; // Repeated Referenzkonfiguration
    vec2D_dbl_ptr_Type      pointsUniRef_; // Unique Referenzkonfiguration
    
    //vec_int_ptr_Type      elementFlag_;
 
    int elementOrder_;
    int surfaceElementOrder_;
    int edgesElementOrder_;
 
    AABBTreePtr_Type AABBTree_;
    
    tuple_intint_Type rankRange_;
    
    /* ###################################################################### */
private:
 
    void flipSurface(ElementsPtr_Type subEl, int surfaceNumber);
    void flipElement(ElementsPtr_Type elements, int elementNumber);
};
}
 
#endif