dam_nodal_reference_temperature_process.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    D.J. Vicente
//
//
 
#if !defined(KRATOS_DAM_NODAL_REFERENCE_TEMPERATURE_PROCESS )
#define  KRATOS_DAM_NODAL_REFERENCE_TEMPERATURE_PROCESS
 
#include <cmath>
 
// Project includes
#include "includes/kratos_flags.h"
#include "includes/kratos_parameters.h"
#include "processes/process.h"
 
// Application include
#include "dam_application_variables.h"
 
namespace Kratos
{
 
class DamNodalReferenceTemperatureProcess : public Process
{
 
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(DamNodalReferenceTemperatureProcess);
 
    typedef Table<double,double> TableType;
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    DamNodalReferenceTemperatureProcess(ModelPart& rModelPart, TableType& Table,
                                Parameters& rParameters
                                ) : Process(Flags()) , mrModelPart(rModelPart) , mrTable(Table)
    {
        KRATOS_TRY
 
        //only include validation with c++11 since raw_literals do not exist in c++03
        Parameters default_parameters( R"(
            {
                "model_part_name":"PLEASE_CHOOSE_MODEL_PART_NAME",
                "variable_name"      : "PLEASE_PRESCRIBE_VARIABLE_NAME",
                "initial_value"      : 0.0,
                "input_file_name"    : "",
                "interval":[
                0.0,
                0.0
                ]
            }  )" );
 
        // Some values need to be mandatorily prescribed since no meaningful default value exist. For this reason try accessing to them
        // So that an error is thrown if they don't exist
        rParameters["variable_name"];
        rParameters["model_part_name"];
 
        // Now validate agains defaults -- this also ensures no type mismatch
        rParameters.ValidateAndAssignDefaults(default_parameters);
 
        mVariableName = rParameters["variable_name"].GetString();
        mInitialValue = rParameters["initial_value"].GetDouble();
        mInputFile = rParameters["input_file_name"].GetString();
 
        KRATOS_CATCH("");
    }
 
 
    virtual ~DamNodalReferenceTemperatureProcess() {}
 
 
    //----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void Execute() override
    {
    }
 
    //----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void ExecuteInitialize() override
    {
 
        KRATOS_TRY;
 
        const Variable<double>& var = KratosComponents< Variable<double> >::Get(mVariableName);
        const int nnodes = mrModelPart.GetMesh(0).Nodes().size();
 
        if(nnodes != 0)
        {
            ModelPart::NodesContainerType::iterator it_begin = mrModelPart.GetMesh(0).NodesBegin();
 
            std::string no_file_str = "- No file";
            std::string add_file_str = "- Add new file";
 
            if ((mInputFile == "") || strstr(mInputFile.c_str(),no_file_str.c_str()) || strstr(mInputFile.c_str(),add_file_str.c_str()))
            {
                #pragma omp parallel for
                for(int i = 0; i<nnodes; i++)
                {
                    ModelPart::NodesContainerType::iterator it = it_begin + i;
 
                    it->FastGetSolutionStepValue(var) = mInitialValue;
 
                }
            }
            else
            {
                #pragma omp parallel for
                for(int i = 0; i<nnodes; i++)
                {
                    ModelPart::NodesContainerType::iterator it = it_begin + i;
 
                    it->FastGetSolutionStepValue(var) = mrTable.GetValue(it->Id());
 
                }
            }
        }
 
        KRATOS_CATCH("");
    }
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void ExecuteInitializeSolutionStep() override
    {
 
        KRATOS_TRY;
 
        const Variable<double>& var = KratosComponents< Variable<double> >::Get(mVariableName);
        const int nnodes = mrModelPart.GetMesh(0).Nodes().size();
 
 
        if(nnodes != 0)
        {
            ModelPart::NodesContainerType::iterator it_begin = mrModelPart.GetMesh(0).NodesBegin();
 
            std::string no_file_str = "- No file";
            std::string add_file_str = "- Add new file";
 
            if ((mInputFile == "") || strstr(mInputFile.c_str(),no_file_str.c_str()) || strstr(mInputFile.c_str(),add_file_str.c_str()))
            {
                #pragma omp parallel for
                for(int i = 0; i<nnodes; i++)
                {
                    ModelPart::NodesContainerType::iterator it = it_begin + i;
 
                    it->FastGetSolutionStepValue(var) = mInitialValue;
 
                }
            }
            else
            {
                #pragma omp parallel for
                for(int i = 0; i<nnodes; i++)
                {
                    ModelPart::NodesContainerType::iterator it = it_begin + i;
 
                    it->FastGetSolutionStepValue(var) = mrTable.GetValue(it->Id());
 
                }
            }
        }
 
        KRATOS_CATCH("");
    }
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    std::string Info() const override
    {
        return "DamNodalReferenceTemperatureProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "DamNodalReferenceTemperatureProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
protected:
 
 
    ModelPart& mrModelPart;
    TableType& mrTable;
    std::string mVariableName;
    double mInitialValue;
    std::string mInputFile;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
private:
 
    DamNodalReferenceTemperatureProcess& operator=(DamNodalReferenceTemperatureProcess const& rOther);
 
};//Class
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                    DamNodalReferenceTemperatureProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const DamNodalReferenceTemperatureProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
} /* namespace Kratos.*/
 
#endif /* KRATOS_DAM_NODAL_REFERENCE_TEMPERATURE_PROCESS defined */