direction_damping_utilities.h

Go to the documentation of this file.

// ==============================================================================
//  KratosShapeOptimizationApplication
//
//  License:         BSD License
//                   license: ShapeOptimizationApplication/license.txt
//
//  Main authors:    Armin Geiser, https://github.com/armingeiser
//
// ==============================================================================
 
#ifndef DIRECTION_DAMPING_UTILITIES_H
#define DIRECTION_DAMPING_UTILITIES_H
 
// ------------------------------------------------------------------------------
// System includes
// ------------------------------------------------------------------------------
#include <iostream>
#include <string>
 
// ------------------------------------------------------------------------------
// Project includes
// ------------------------------------------------------------------------------
#include "includes/model_part.h"
#include "spatial_containers/spatial_containers.h"
#include "custom_utilities/filter_function.h"
 
// ==============================================================================
 
namespace Kratos
{
 
 
 
class KRATOS_API(SHAPE_OPTIMIZATION_APPLICATION) DirectionDampingUtilities
{
public:
 
    // Type definitions for better reading later
    typedef array_1d<double,3> array_3d;
    typedef ModelPart::NodeType NodeType;
    typedef ModelPart::NodeType::Pointer NodeTypePointer;
    typedef std::vector<NodeTypePointer> NodeVector;
    typedef std::vector<NodeTypePointer>::iterator NodeVectorIterator;
    typedef std::vector<double> DoubleVector;
    typedef std::vector<double>::iterator DoubleVectorIterator;
 
    // Type definitions for tree-search
    typedef Bucket< 3, NodeType, NodeVector, NodeTypePointer, NodeVectorIterator, DoubleVectorIterator > BucketType;
    typedef Tree< KDTreePartition<BucketType> > KDTree;
 
    KRATOS_CLASS_POINTER_DEFINITION(DirectionDampingUtilities);
 
 
    DirectionDampingUtilities( ModelPart& modelPartToDamp, Parameters DampingSettings );
 
    virtual ~DirectionDampingUtilities()
    {
    }
 
 
 
 
 
    void CreateListOfNodesOfModelPart();
 
    void CreateSearchTreeWithAllNodesOfModelPart();
 
    void InitalizeDampingFactorsToHaveNoInfluence();
 
    void SetDampingFactors();
 
    FilterFunction::Pointer CreateDampingFunction( std::string damping_type ) const;
 
    void ThrowWarningIfNodeNeighborsExceedLimit( const ModelPart::NodeType& given_node, const unsigned int number_of_neighbors ) const;
 
    void DampNodalVariable( const Variable<array_3d> &rNodalVariable );
 
 
 
    virtual std::string Info() const
    {
        return "DirectionDampingUtilities";
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const
    {
        rOStream << "DirectionDampingUtilities";
    }
 
    virtual void PrintData(std::ostream& rOStream) const
    {
    }
 
private:
 
    // Initialized by class constructor
    ModelPart& mrModelPartToDamp;
    Parameters mDampingSettings;
    array_3d mDirection;
    std::vector<double> mDampingFactors;
 
    // Variables for spatial search
    unsigned int mBucketSize = 100;
    unsigned int mMaxNeighborNodes = 10000;
    NodeVector mListOfNodesOfModelPart;
    KDTree::Pointer mpSearchTree;
 
 
}; // Class DirectionDampingUtilities
 
 
}  // namespace Kratos.
 
#endif // DIRECTION_DAMPING_UTILITIES_H