find_neighbour_elements_of_conditions_process.hpp

Go to the documentation of this file.

// KRATOS___
//     //   ) )
//    //         ___      ___
//   //  ____  //___) ) //   ) )
//  //    / / //       //   / /
// ((____/ / ((____   ((___/ /  MECHANICS
//
//  License:         geo_mechanics_application/license.txt
//
//  Main authors:    Vahid Galavi
//
 
 
#pragma once
 
// Project includes
#include "includes/model_part.h"
#include "processes/process.h"
 
namespace Kratos
{
    using hashmap = std::unordered_multimap<DenseVector<int>, std::vector<Condition::Pointer>, KeyHasherRange<DenseVector<int>>, KeyComparorRange<DenseVector<int>>>;
 
 
class KRATOS_API(GEO_MECHANICS_APPLICATION) FindNeighbourElementsOfConditionsProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(FindNeighbourElementsOfConditionsProcess);
 
    using IndexType = std::size_t;
 
    using NodeType = Node;
 
    // Definition of the geometry
    using GeometryType = Geometry<NodeType>;
 
 
 
 
    explicit FindNeighbourElementsOfConditionsProcess( ModelPart& rModelPart ): Process(),
            mrModelPart(rModelPart)
    {
    }
 
    FindNeighbourElementsOfConditionsProcess& operator=(const FindNeighbourElementsOfConditionsProcess&) = delete;
    FindNeighbourElementsOfConditionsProcess(const FindNeighbourElementsOfConditionsProcess&) = delete;
    ~FindNeighbourElementsOfConditionsProcess() override = default;
 
 
    void operator()()
    {
        Execute();
    }
 
 
 
    void Execute() override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "FindNeighbourElementsOfConditionsProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        this->PrintInfo(rOStream);
    }
 
 
private:
    hashmap::iterator FindFaceReorderingTetrahedra3D10(DenseVector<int> FaceIds, hashmap &FacesMap) const;
    hashmap::iterator FindFaceReorderingTetrahedra3D4(DenseVector<int> FaceIds, hashmap &FacesMap) const;
    hashmap::iterator FindFaceReorderingHexahedra3D8(DenseVector<int> FaceIds, hashmap &FacesMap) const;
    hashmap::iterator FindFaceReorderingHexahedra3D20(DenseVector<int> FaceIds, hashmap &FacesMap) const;
 
 
 
    ModelPart& mrModelPart;
 
 
    bool CheckIfAllConditionsAreVisited() const;
 
    void CheckIf1DElementIsNeighbour(hashmap& rFacesMap);
 
    static void CheckForMultipleConditionsOnElement(hashmap& rFacesMap, hashmap::iterator& rItFace,
        PointerVector<Element>::iterator pItElem);
 
}; // Class Process
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  FindNeighbourElementsOfConditionsProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const FindNeighbourElementsOfConditionsProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
}