assign_scalar_input_to_entities_process.h

Go to the documentation of this file.

//
//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Vicente Mataix Ferrandiz
//
 
#if !defined(KRATOS_ASSIGN_SCALAR_INPUT_TO_ENTITIES_PROCESS_H_INCLUDED )
#define  KRATOS_ASSIGN_SCALAR_INPUT_TO_ENTITIES_PROCESS_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/model_part.h"
#include "includes/kratos_parameters.h"
#include "processes/process.h"
#include "utilities/result_dabatase.h"
 
namespace Kratos
{
 
 
struct AssignScalarInputToEntitiesProcessSettings
{
    // Defining clearer options
    constexpr static bool SaveAsHistoricalVariable = true;
    constexpr static bool SaveAsNonHistoricalVariable = false;
};
 
template<class TEntity, bool THistorical = false>
class KRATOS_API(KRATOS_CORE) AssignScalarInputToEntitiesProcess
    : public Process
{
public:
 
    typedef Node NodeType;
 
    typedef std::size_t IndexType;
 
    typedef std::size_t SizeType;
 
    typedef PointerVectorSet<TEntity, IndexedObject> EntityContainerType;
 
    KRATOS_CLASS_POINTER_DEFINITION(AssignScalarInputToEntitiesProcess);
 
    KRATOS_DEFINE_LOCAL_FLAG( GEOMETRIC_DEFINITION ); 
 
 
    enum class Algorithm {
        NEAREST_NEIGHBOUR = 0
    };
 
 
    AssignScalarInputToEntitiesProcess(
        ModelPart& rModelPart,
        Parameters rParameters
        );
 
    ~AssignScalarInputToEntitiesProcess() override {}
 
 
 
    void ExecuteInitializeSolutionStep() override;
 
    const Parameters GetDefaultParameters() const override;
 
 
 
 
 
 
    std::string Info() const override
    {
        return "AssignScalarInputToEntitiesProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "AssignScalarInputToEntitiesProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
protected:
 
 
private:
 
 
    ModelPart& mrModelPart;                                                  
 
    const Variable<double>* mpVariable = nullptr;                            
 
    ResultDatabase mDatabase;                                                
 
    std::vector<std::unordered_map<IndexType, double>> mWeightExtrapolation; 
 
    std::vector<array_1d<double, 3>> mCoordinates;                           
 
    Algorithm mAlgorithm = Algorithm::NEAREST_NEIGHBOUR;                     
 
 
    void IdentifyDataTXT(const std::string& rFileName);
 
    void IdentifyDataJSON(const std::string& rFileName);
 
    void ReadDataTXT(const std::string& rFileName);
 
    void ReadDataJSON(const std::string& rFileName);
 
    void ComputeExtrapolationWeight();
 
    void InternalAssignValue(
        const Variable<double>& rVariable,
        const double Value
        );
 
    Algorithm ConvertAlgorithmString(const std::string& Str)
    {
        if(Str == "NEAREST_NEIGHBOUR" || Str == "nearest_neighbour")
            return Algorithm::NEAREST_NEIGHBOUR;
        else
            return Algorithm::NEAREST_NEIGHBOUR;
    }
 
    array_1d<double, 3> GetCoordinatesEntity(const IndexType Id);
 
    EntityContainerType& GetEntitiesContainer();
 
    void ResetValues();
 
    void SetValue(
        TEntity& rEntity,
        const Variable<double>& rVariable,
        const double Value
        );
 
    double& GetValue(
        TEntity& rEntity,
        const Variable<double>& rVariable
        );
 
 
}; // Class AssignScalarInputToEntitiesProcess
 
 
 
 
 
 
 
template<class TEntity>
inline std::istream& operator >> (std::istream& rIStream,
                                  AssignScalarInputToEntitiesProcess<TEntity>& rThis);
 
template<class TEntity>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const AssignScalarInputToEntitiesProcess<TEntity>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_ASSIGN_SCALAR_INPUT_TO_ENTITIES_PROCESS_H_INCLUDED  defined