metrics_error_process.h

Go to the documentation of this file.

// KRATOS  __  __ _____ ____  _   _ ___ _   _  ____
//        |  \/  | ____/ ___|| | | |_ _| \ | |/ ___|
//        | |\/| |  _| \___ \| |_| || ||  \| | |  _
//        | |  | | |___ ___) |  _  || || |\  | |_| |
//        |_|  |_|_____|____/|_| |_|___|_| \_|\____| APPLICATION
//
//  License:         BSD License
//                       license: MeshingApplication/license.txt
//
//  Main authors:    Vicente Mataix Ferrandiz 
//                   Anna Rehr
//
 
#if !defined(KRATOS_ERROR_METRICS_PROCESS)
#define KRATOS_ERROR_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) MetricErrorProcess
    : public Process
{
public:
 
    typedef ModelPart::NodesContainerType                                     NodesArrayType;
    typedef ModelPart::ElementsContainerType                               ElementsArrayType;
    typedef ModelPart::ConditionsContainerType                           ConditionsArrayType;
 
    typedef Node                                                                NodeType;
 
    typedef GlobalPointersVector< 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(MetricErrorProcess);
 
 
    // Constructor
 
    MetricErrorProcess(
        ModelPart& rThisModelPart,
        Parameters ThisParameters = Parameters(R"({})")
        );
 
    ~MetricErrorProcess() override = default;
 
 
    void operator()()
    {
        Execute();
    }
 
 
    void Execute() override;
 
    const Parameters GetDefaultParameters() const override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "MetricErrorProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "MetricErrorProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    ModelPart& mrThisModelPart; 
 
    double mMinSize;           
    double mMaxSize;           
 
    bool mSetElementNumber;    
    SizeType mElementNumber;   
    double mTargetError;       
    bool mAverageNodalH;       
 
    SizeType mEchoLevel;       
 
 
 
    void CalculateElementSize();
 
    void CalculateMetric();
 
 
 
 
 
    MetricErrorProcess& operator=(MetricErrorProcess const& rOther);
 
    //MetricErrorProcess(MetricErrorProcess const& rOther);
 
};// class MetricErrorProcess
 
 
 
 
 
template<unsigned int TDim, class TVarType>
inline std::istream& operator >> (std::istream& rIStream,
                                  MetricErrorProcess<TDim>& rThis);
 
template<unsigned int TDim, class TVarType>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const MetricErrorProcess<TDim>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
};// namespace Kratos.
#endif /* KRATOS_ERROR_METRICS_PROCESS defined */