calculate_distance_to_boundary_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Miguel Maso Sotomayor
//
 
#ifndef KRATOS_CALCULATE_DISTANCE_TO_BOUNDARY_PROCESS_H_INCLUDED
#define KRATOS_CALCULATE_DISTANCE_TO_BOUNDARY_PROCESS_H_INCLUDED
 
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "containers/model.h"
#include "processes/process.h"
#include "geometries/geometry.h"
#include "includes/kratos_parameters.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
class ModelPart;
 
class KRATOS_API(SHALLOW_WATER_APPLICATION) CalculateDistanceToBoundaryProcess : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(CalculateDistanceToBoundaryProcess);
 
    typedef Geometry<Point> GeometryType;
 
    typedef Node NodeType;
 
 
    CalculateDistanceToBoundaryProcess(
        Model& rModel,
        Parameters ThisParameters) :
            Process(),
            mrModelPart(rModel.GetModelPart(ThisParameters["computing_model_part_name"].GetString())),
            mrBoundaryPart(rModel.GetModelPart(ThisParameters["absorbing_boundary_name"].GetString()))
    {
        msInstancesCount++;
        mInitializeDistance = (msInstancesCount == 1);
        ThisParameters.ValidateAndAssignDefaults(GetDefaultParameters());
        mRSquaredThreshold = ThisParameters["r_squared_threshold"].GetDouble();
        FindApproximatingGeometry(mpBoundary, mrBoundaryPart);
    }
 
    CalculateDistanceToBoundaryProcess(
        ModelPart& rComputingModelPart,
        ModelPart& rBoundaryModelPart,
        double RSquaredThreshold = 0.99) :
            Process(),
            mrModelPart(rComputingModelPart),
            mrBoundaryPart(rBoundaryModelPart)
    {
        msInstancesCount++;
        mInitializeDistance = (msInstancesCount == 1);
        mRSquaredThreshold = RSquaredThreshold;
        FindApproximatingGeometry(mpBoundary, mrBoundaryPart);
    }
 
    virtual ~CalculateDistanceToBoundaryProcess() {
        msInstancesCount--;
    }
 
 
 
 
    int Check() override;
 
    void ExecuteBeforeSolutionLoop() override;
 
    const Parameters GetDefaultParameters() const override;
 
 
 
 
 
 
    virtual std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "CalculateDistanceToBoundaryProcess";
        return buffer.str();
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info();
    }
 
    virtual void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
private:
 
    static std::size_t msInstancesCount;
 
 
    ModelPart& mrModelPart;
    ModelPart& mrBoundaryPart;
    GeometryType::Pointer mpBoundary;
    double mRSquaredThreshold;
    bool mBruteForceSearch;
    bool mInitializeDistance;
 
 
 
 
    void FindApproximatingGeometry(GeometryType::Pointer& pEntity, const ModelPart& rBoundaryPart);
 
    double RSquared(const GeometryType& rLine, const ModelPart& rBoundaryPart);
 
    double SquaredDistance(const Point& rPointA, const Point& rPointB);
 
    double Distance(const Point& rPointA, const Point& rPointB);
 
 
 
 
 
 
    CalculateDistanceToBoundaryProcess& operator=(CalculateDistanceToBoundaryProcess const& rOther) = delete;
 
    CalculateDistanceToBoundaryProcess(CalculateDistanceToBoundaryProcess const& rOther) = delete;
 
 
 
}; // Class CalculateDistanceToBoundaryProcess
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                CalculateDistanceToBoundaryProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                const CalculateDistanceToBoundaryProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_CALCULATE_DISTANCE_TO_BOUNDARY_PROCESS_H_INCLUDED  defined