metrics_levelset_process.h

Go to the documentation of this file.

// KRATOS  __  __ _____ ____  _   _ ___ _   _  ____
//        |  \/  | ____/ ___|| | | |_ _| \ | |/ ___|
//        | |\/| |  _| \___ \| |_| || ||  \| | |  _
//        | |  | | |___ ___) |  _  || || |\  | |_| |
//        |_|  |_|_____|____/|_| |_|___|_| \_|\____| APPLICATION
//
//  License:         BSD License
//                       license: MeshingApplication/license.txt
//
//  Main authors:    Vicente Mataix Ferrandiz
//
 
#if !defined(KRATOS_LEVELSET_METRICS_PROCESS)
#define KRATOS_LEVELSET_METRICS_PROCESS
 
// Project includes
#include "includes/kratos_parameters.h"
#include "processes/process.h"
#include "includes/model_part.h"
#include "meshing_application_variables.h"
 
namespace Kratos
{
 
 
    typedef std::size_t SizeType;
 
 
 
 
template<SizeType TDim>
class KRATOS_API(MESHING_APPLICATION) ComputeLevelSetSolMetricProcess
    : public Process
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(ComputeLevelSetSolMetricProcess);
 
    typedef Node                                                   NodeType;
 
    typedef ModelPart::NodesContainerType                        NodesArrayType;
    typedef ModelPart::ElementsContainerType                  ElementsArrayType;
    typedef ModelPart::ConditionsContainerType              ConditionsArrayType;
 
    typedef std::size_t                                               IndexType;
 
    typedef typename std::conditional<TDim == 2, array_1d<double, 3>, array_1d<double, 6>>::type TensorArrayType;
 
 
    enum class Interpolation {CONSTANT = 0, LINEAR = 1, EXPONENTIAL = 2, PIECEWISE_LINEAR = 3};
 
 
    ComputeLevelSetSolMetricProcess(
        ModelPart& rThisModelPart,
        const Variable<array_1d<double,3>>& rVariableGradient = DISTANCE_GRADIENT,
        Parameters ThisParameters = Parameters(R"({})")
        );
 
    ~ComputeLevelSetSolMetricProcess() override = default;
 
 
    void operator()()
    {
        Execute();
    }
 
    const Parameters GetDefaultParameters() const override;
 
 
    void Execute() override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "ComputeLevelSetSolMetricProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "ComputeLevelSetSolMetricProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    ModelPart& mThisModelPart;                             
    const Variable<array_1d<double,3>>& mVariableGradient; 
    std::string mRatioReferenceVariable = "DISTANCE";      
    std::string mSizeReferenceVariable = "DISTANCE";       
    double mMinSize;                                       
    double mMaxSize;                                       
    bool mEnforceCurrent;                                  
    double mAnisotropicRatio;                              
    double mBoundLayer;                                    
    double mSizeBoundLayer;                                
    Table<double> mSizeDistributionTable;                  
    Interpolation mInterpolation;                          
    Interpolation mSizeInterpolation;                      
 
 
 
    TensorArrayType ComputeLevelSetMetricTensor(
        const array_1d<double, 3>& GradientValue,
        const double Ratio,
        const double ElementSize
        );
 
    Interpolation ConvertInter(const std::string& Str)
    {
        if(Str == "Constant" || Str == "CONSTANT" || Str == "constant")
            return Interpolation::CONSTANT;
        else if(Str == "Linear" || Str == "LINEAR"  || Str == "linear")
            return Interpolation::LINEAR;
        else if(Str == "Exponential" || Str == "EXPONENTIAL"  || Str == "exponential")
            return Interpolation::EXPONENTIAL;
        else if(Str == "piecewise_linear" || Str == "PIECEWISE_LINEAR")
            return Interpolation::PIECEWISE_LINEAR;
        else
            return Interpolation::LINEAR;
    }
 
    double CalculateAnisotropicRatio(
        const double Distance
        );
 
    double CalculateElementSize(
        const double Distance,
        const double NodalH
        );
 
 
 
 
 
    ComputeLevelSetSolMetricProcess& operator=(ComputeLevelSetSolMetricProcess const& rOther)
    {
        return *this;
    };
 
    //ComputeLevelSetSolMetricProcess(ComputeLevelSetSolMetricProcess const& rOther);
 
};// class ComputeLevelSetSolMetricProcess
 
 
 
template<unsigned int TDim>
inline std::istream& operator >> (std::istream& rIStream,
                                  ComputeLevelSetSolMetricProcess<TDim>& rThis);
 
template<unsigned int TDim>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const ComputeLevelSetSolMetricProcess<TDim>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
};// namespace Kratos.
#endif /* KRATOS_LEVELSET_METRICS_PROCESS defined */