tetrahedral_mesh_orientation_check.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:           BSD License
//                          Kratos default license: kratos/license.txt
//
//  Main authors:    Pooyan Dadvand
//                   Riccardo Rossi
//
//
 
#ifndef KRATOS_TETRAHEDRAL_MESH_ORIENTATION_CHECK_H
#define KRATOS_TETRAHEDRAL_MESH_ORIENTATION_CHECK_H
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "processes/process.h"
#include "geometries/geometry.h"
 
namespace Kratos
{
 
typedef Node NodeType;
 
typedef Geometry<NodeType> GeometryType;
 
 
class KRATOS_API(KRATOS_CORE) TetrahedralMeshOrientationCheck
    : public Process
{
public:
 
    //DEFINITION OF FLAGS TO CONTROL THE BEHAVIOUR
    KRATOS_DEFINE_LOCAL_FLAG(ASSIGN_NEIGHBOUR_ELEMENTS_TO_CONDITIONS);
    KRATOS_DEFINE_LOCAL_FLAG(COMPUTE_NODAL_NORMALS);
    KRATOS_DEFINE_LOCAL_FLAG(COMPUTE_CONDITION_NORMALS);
    KRATOS_DEFINE_LOCAL_FLAG(MAKE_VOLUMES_POSITIVE);
    KRATOS_DEFINE_LOCAL_FLAG(ALLOW_REPEATED_CONDITIONS);
 
    KRATOS_CLASS_POINTER_DEFINITION(TetrahedralMeshOrientationCheck);
 
    typedef std::size_t IndexType;
 
    typedef std::size_t SizeType;
 
    typedef Node NodeType;
 
    // Definition of the geometry
    typedef Geometry<NodeType> GeometryType;
 
 
 
    TetrahedralMeshOrientationCheck(
        ModelPart& rModelPart,
        bool ThrowErrors,
        const Flags options = COMPUTE_NODAL_NORMALS.AsFalse() | COMPUTE_CONDITION_NORMALS.AsFalse() | ASSIGN_NEIGHBOUR_ELEMENTS_TO_CONDITIONS.AsFalse() | ALLOW_REPEATED_CONDITIONS.AsFalse()
        ):  Process(),
            mrModelPart(rModelPart),
            mThrowErrors(ThrowErrors), //to be changed to a flag
            mrOptions(options)
 
    {
    }
 
    TetrahedralMeshOrientationCheck(
        ModelPart& rModelPart,
        const Flags options = COMPUTE_NODAL_NORMALS.AsFalse() | COMPUTE_CONDITION_NORMALS.AsFalse() | ASSIGN_NEIGHBOUR_ELEMENTS_TO_CONDITIONS.AsFalse() | ALLOW_REPEATED_CONDITIONS.AsFalse()
        ):  Process(),
            mrModelPart(rModelPart),
            mThrowErrors(false),
            mrOptions(options)
    {
    }
 
    ~TetrahedralMeshOrientationCheck() override {}
 
 
    void operator()()
    {
        Execute();
    }
 
 
 
    void Execute() override;
 
    void SwapAll();
 
    void SwapNegativeElements();
 
 
 
 
 
 
    std::string Info() const override
    {
        return "TetrahedralMeshOrientationCheck";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "TetrahedralMeshOrientationCheck";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        this->PrintInfo(rOStream);
    }
 
 
 
 
 
 
private:
 
 
 
    ModelPart& mrModelPart;
    const bool mThrowErrors;
    Flags mrOptions;
 
 
    bool Orient(GeometryType& rGeom);
 
 
    TetrahedralMeshOrientationCheck& operator=(TetrahedralMeshOrientationCheck const& rOther);
 
    TetrahedralMeshOrientationCheck(TetrahedralMeshOrientationCheck const& rOther);
 
 
}; // Class Process
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  TetrahedralMeshOrientationCheck& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const TetrahedralMeshOrientationCheck& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
 
} // namespace Kratos
 
 
#endif // KRATOS_TETRAHEDRAL_MESH_ORIENTATION_PROCESS_H