apply_perturbation_function_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Miguel Maso Sotomayor
//
 
#ifndef KRATOS_APPLY_PERTURBATION_FUNCTION_PROCESS_H_INCLUDED
#define KRATOS_APPLY_PERTURBATION_FUNCTION_PROCESS_H_INCLUDED
 
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "includes/model_part.h"
#include "processes/process.h"
#include "includes/kratos_parameters.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
template<class TVarType>
class KRATOS_API(SHALLOW_WATER_APPLICATION) ApplyPerturbationFunctionProcess : public Process
{
public:
 
    typedef std::size_t                     IndexType;
    typedef Node                         NodeType;
    typedef Node::Pointer                NodePointerType;
    typedef ModelPart::NodesContainerType   NodesArrayType;
    typedef NodesArrayType::iterator        NodeIteratorType;
 
    KRATOS_CLASS_POINTER_DEFINITION(ApplyPerturbationFunctionProcess);
 
 
    ApplyPerturbationFunctionProcess(
        ModelPart& rThisModelPart,
        NodePointerType pNode,
        TVarType& rThisVariable,
        Parameters& rThisParameters);
 
    ApplyPerturbationFunctionProcess(
        ModelPart& rThisModelPart,
        NodesArrayType& rSourcePoints,
        TVarType& rThisVariable,
        Parameters& rThisParameters);
 
    ~ApplyPerturbationFunctionProcess() override {}
 
 
    /*
     * @brief This operator is provided to call the process as a function and simply calls the Execute method.
     */
    void operator()()
    {
        Execute();
    }
 
 
    void Execute() override;
 
    int Check() override;
 
    void ExecuteBeforeSolutionLoop() override;
 
 
 
 
 
 
    virtual std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "ApplyPerturbationFunctionProcess" ;
        return buffer.str();
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const override {rOStream << "ApplyPerturbationFunctionProcess";}
 
    virtual void PrintData(std::ostream& rOStream) const override {}
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
    ModelPart& mrModelPart;
    NodesArrayType mSourcePoints;
    TVarType& mrVariable;
    double mDefaultValue;
    double mInfluenceDistance;
    double mPerturbation;
    double mHalfWaveNumber;
 
 
 
 
    void ValidateParameters(Parameters& rParameters);
 
    double ComputeDistance(NodeType& rNode);
 
    double PointPointSquareDistance(array_1d<double, 3>& rCoordA, array_1d<double, 3>& rCoordB);
 
    double ComputeInitialValue(double& rDistance);
 
 
 
 
 
 
    // ApplyPerturbationFunctionProcess& operator=(ApplyPerturbationFunctionProcess const& rOther);
 
    // ApplyPerturbationFunctionProcess(ApplyPerturbationFunctionProcess const& rOther);
 
 
 
}; // Class ApplyPerturbationFunctionProcess
 
 
 
 
 
// /// input stream function
// inline std::istream& operator >> (std::istream& rIStream,
//                 ApplyPerturbationFunctionProcess& rThis);
 
// /// output stream function
// inline std::ostream& operator << (std::ostream& rOStream,
//                 const ApplyPerturbationFunctionProcess& rThis)
// {
//     rThis.PrintInfo(rOStream);
//     rOStream << std::endl;
//     rThis.PrintData(rOStream);
 
//     return rOStream;
// }
 
 
 
}  // namespace Kratos.
 
#endif // KRATOS_APPLY_PERTURBATION_FUNCTION_PROCESS_H_INCLUDED  defined