metrics_divergencefree_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Riccardo Tosi
//                   Brendan Keith
//
 
#if !defined(KRATOS_DIVERGENCEFREE_METRICS_PROCESS)
#define KRATOS_DIVERGENCEFREE_METRICS_PROCESS
 
// Project includes
#include "includes/kratos_parameters.h"
#include "includes/model_part.h"
#include "processes/process.h"
 
namespace Kratos
{
 
 
    typedef std::size_t SizeType;
 
 
 
 
template<SizeType TDim>
class KRATOS_API(MESHING_APPLICATION) MetricDivergenceFreeProcess
    : public Process
{
public:
 
    typedef ModelPart::NodesContainerType                                     NodesArrayType;
    typedef ModelPart::ElementsContainerType                               ElementsArrayType;
    typedef ModelPart::ConditionsContainerType                           ConditionsArrayType;
 
    typedef Node                                                                NodeType;
 
    typedef WeakPointerVector< Element >::iterator                         WeakElementItType;
    typedef NodesArrayType::iterator                                              NodeItType;
    typedef ElementsArrayType::iterator                                        ElementItType;
 
    typedef std::size_t                                                            IndexType;
 
    typedef BoundedMatrix<double, TDim, TDim> MatrixType;
 
    typedef typename std::conditional<TDim == 2, array_1d<double, 3>, array_1d<double, 6>>::type TensorArrayType;
 
    KRATOS_CLASS_POINTER_DEFINITION(MetricDivergenceFreeProcess);
 
 
    // Constructor
 
    MetricDivergenceFreeProcess(
        ModelPart& rThisModelPart,
        Parameters ThisParameters = Parameters(R"({})")
        );
 
    ~MetricDivergenceFreeProcess() override = default;
 
 
    void operator()()
    {
        Execute();
    }
 
 
    void Execute() override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "MetricDivergenceFreeProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "MetricDivergenceFreeProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    ModelPart& mrThisModelPart; 
 
    double mMinSize;                 
    double mMaxSize;                 
    enum RefinementStrategies {MaximumStrategy, MeanDistributionStrategy, GlobalToleranceStrategy};
    RefinementStrategies mRefinementStrategy; 
    SizeType mEchoLevel;             
    std::string mReferenceVariable;
 
    // Mean dsitribution strategy
    std::string mMeanStrategyReferenceNorm;
    double mMeanStrategyTargetRefinementCoefficient;
    double mMeanStrategyRefinementBound;
    double mMeanStrategyDivergenceFreeOverAllDomain;
 
    // Maximum strategy
    double mMaxStrategyTargetRefinementCoefficient;
    double mMaxStrategyRefinementCoefficient;
    double mDivergenceFreeMaxValue;
 
    // Interpolation error strategy
    double mGlobalErrorStrategyGlobalTolerance;
    double mGlobalErrorStrategyMeshConstant;
 
 
 
    void InitializeRefinementStrategy();
 
    void CalculateMetric();
 
 
 
 
 
    MetricDivergenceFreeProcess& operator=(MetricDivergenceFreeProcess const& rOther);
 
    //MetricDivergenceFreeProcess(MetricDivergenceFreeProcess const& rOther);
 
};// class MetricDivergenceFreeProcess
 
 
 
template<unsigned int TDim, class TVarType>
inline std::istream& operator >> (std::istream& rIStream,
                                  MetricDivergenceFreeProcess<TDim>& rThis);
 
template<unsigned int TDim, class TVarType>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const MetricDivergenceFreeProcess<TDim>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
};// namespace Kratos.
#endif /* KRATOS_DIVERGENCEFREE_METRICS_PROCESS defined */