mpi_normal_calculation_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main author:     Jordi Cotela
//
 
#ifndef KRATOS_MPI_NORMAL_CALCULATION_UTILITIES_H
#define KRATOS_MPI_NORMAL_CALCULATION_UTILITIES_H
 
// System includes
 
// External includes
 
 
// Project includes
#include "includes/model_part.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
class KRATOS_API(KRATOS_MPI_CORE) MPINormalCalculationUtils
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(MPINormalCalculationUtils);
 
 
    MPINormalCalculationUtils();
 
    virtual ~MPINormalCalculationUtils();
 
 
 
 
 
    int Check(ModelPart& rModelPart);
 
    void OrientFaces(ModelPart& rModelPart,
                     bool OutwardsPositive);
 
    void CalculateOnSimplex(ModelPart& rModelPart,
                            int Dimension,
                            const Variable<double>& rVariable,
                            const double rAlpha);
 
 
 
 
 
 
    std::string Info() const;
 
    void PrintInfo(std::ostream& rOStream) const;
 
    void PrintData(std::ostream& rOStream) const;
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
 
 
 
 
 
    void IdentifyFaces(ModelPart& rModelPart,
                       const Variable<double>& rVariable,
                       int& MaxNeigh,
                       int& NumNodes);
 
 
    void InitializeNormalData(ModelPart& rModelPart,
                              const Variable<double>& rVariable,
                              int MaxNeigh,
                              std::vector<double>& rNormals,
                              std::vector<int>& pActiveNeigh,
                              const int NumNodes);
 
    void DetectEdges(ModelPart& rModelPart,
                     unsigned int Dimension,
                     double MaxAngle,
                     const std::vector<double>& rNormals,
                     const std::vector<int>& rActiveNeigh,
                     const int NumNodes,
                     const int MaxNeighs);
 
 
    bool OrientElement(Geometry< Node >& rGeom);
 
 
    void NormalContribution(Geometry< Node >&rGeom);
 
 
    void FaceNormal2D(array_1d<double,3>& An,
                      Geometry<Node >& rGeometry);
 
 
    void FaceNormal3D(array_1d<double,3>& An,
                      Geometry<Node >& rGeometry);
 
    void UpdateNodeNormals(ModelPart& rModelPart,
                           const unsigned int Dimension,
                           const Variable<double>& rVariable);
 
 
 
 
 
 
 
    MPINormalCalculationUtils& operator=(MPINormalCalculationUtils const& rOther);
 
    MPINormalCalculationUtils(MPINormalCalculationUtils const& rOther);
 
 
 
}; // Class MPINormalCalculationUtils
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  MPINormalCalculationUtils& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const MPINormalCalculationUtils& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_MPI_NORMAL_CALCULATION_UTILITIES_H