mapper_vertex_morphing_matrix_free.h

Go to the documentation of this file.

// ==============================================================================
//  KratosShapeOptimizationApplication
//
//  License:         BSD License
//                   license: ShapeOptimizationApplication/license.txt
//
//  Main authors:    Baumgaertner Daniel, https://github.com/dbaumgaertner
//
// ==============================================================================
 
#ifndef MAPPER_VERTEX_MORPHING_MATRIX_FREE_H
#define MAPPER_VERTEX_MORPHING_MATRIX_FREE_H
 
// ------------------------------------------------------------------------------
// System includes
// ------------------------------------------------------------------------------
#include <iostream>
#include <string>
 
// ------------------------------------------------------------------------------
// Project includes
// ------------------------------------------------------------------------------
#include "includes/model_part.h"
#include "spatial_containers/spatial_containers.h"
#include "mapper_base.h"
#include "custom_utilities/filter_function.h"
 
// ==============================================================================
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(SHAPE_OPTIMIZATION_APPLICATION) MapperVertexMorphingMatrixFree : public Mapper
{
public:
 
    // Type definitions for better reading later
    typedef Node NodeType;
    typedef Node ::Pointer NodeTypePointer;
    typedef std::vector<NodeType::Pointer> NodeVector;
    typedef std::vector<NodeType::Pointer>::iterator NodeIterator;
    typedef std::vector<double>::iterator DoubleVectorIterator;
    typedef ModelPart::ConditionsContainerType ConditionsArrayType;
 
    // Type definitions for tree-search
    typedef Bucket< 3, NodeType, NodeVector, NodeTypePointer, NodeIterator, DoubleVectorIterator > BucketType;
    typedef Tree< KDTreePartition<BucketType> > KDTree;
 
    KRATOS_CLASS_POINTER_DEFINITION(MapperVertexMorphingMatrixFree);
 
 
    MapperVertexMorphingMatrixFree( ModelPart& rOriginModelPart, ModelPart& rDestinationModelPart, Parameters MapperSettings )
        : mrOriginModelPart( rOriginModelPart ),
          mrDestinationModelPart( rDestinationModelPart ),
          mMapperSettings( MapperSettings ),
          mFilterRadius( MapperSettings["filter_radius"].GetDouble() ),
          mMaxNumberOfNeighbors( MapperSettings["max_nodes_in_filter_radius"].GetInt())
    {
    }
 
    virtual ~MapperVertexMorphingMatrixFree()
    {
    }
 
 
 
 
 
    // --------------------------------------------------------------------------
    void Initialize() override;
 
    // --------------------------------------------------------------------------
    void Map( const Variable<array_3d> &rOriginVariable, const Variable<array_3d> &rDestinationVariable ) override;
 
    // --------------------------------------------------------------------------
    void Map( const Variable<double> &rOriginVariable, const Variable<double> &rDestinationVariable ) override;
 
    // --------------------------------------------------------------------------
    void InverseMap( const Variable<array_3d> &rDestinationVariable, const Variable<array_3d> &rOriginVariable ) override;
 
    // --------------------------------------------------------------------------
    void InverseMap( const Variable<double> &rDestinationVariable, const Variable<double> &rOriginVariable ) override;
 
    // --------------------------------------------------------------------------
    void Update() override;
 
    // --------------------------------------------------------------------------
 
 
    FilterFunction::UniquePointer mpFilterFunction;
 
 
 
 
    virtual std::string Info() const override
    {
        return "MapperVertexMorphingMatrixFree";
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "MapperVertexMorphingMatrixFree";
    }
 
    virtual void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
    // Initialized by class constructor
    ModelPart& mrOriginModelPart;
    ModelPart& mrDestinationModelPart;
    Parameters mMapperSettings;
    double mFilterRadius;
    unsigned int mMaxNumberOfNeighbors;
 
    // Variables for spatial search
    unsigned int mBucketSize = 100;
    NodeVector mListOfNodesInOriginModelPart;
    KDTree::Pointer mpSearchTree;
 
    // Variables for mapping
    std::vector<Vector> mValuesOrigin;
    std::vector<Vector> mValuesDestination;
    bool mIsMappingInitialized = false;
 
 
 
 
    // --------------------------------------------------------------------------
    void CreateListOfNodesInOriginModelPart();
 
    // --------------------------------------------------------------------------
    void CreateFilterFunction();
 
    // --------------------------------------------------------------------------
    void InitializeMappingVariables();
 
    // --------------------------------------------------------------------------
    void AssignMappingIds();
    // --------------------------------------------------------------------------
    void CreateSearchTreeWithAllNodesInOriginModelPart();
 
    // --------------------------------------------------------------------------
    void ThrowWarningIfNumberOfNeighborsExceedsLimit(ModelPart::NodeType& given_node, unsigned int number_of_neighbors);
 
    // --------------------------------------------------------------------------
    virtual void ComputeWeightForAllNeighbors(  const ModelPart::NodeType& design_node,
                                        const NodeVector& neighbor_nodes,
                                        const unsigned int number_of_neighbors,
                                        std::vector<double>& list_of_weights,
                                        double& sum_of_weights );
 
    // --------------------------------------------------------------------------
    double GetVertexMorphingRadius(const NodeType& rNode) const override
    {
        return mFilterRadius;
    }
 
 
 
 
 
 
//      MapperVertexMorphingMatrixFree& operator=(MapperVertexMorphingMatrixFree const& rOther);
 
//      MapperVertexMorphingMatrixFree(MapperVertexMorphingMatrixFree const& rOther);
 
 
 
}; // Class MapperVertexMorphingMatrixFree
 
 
 
 
 
 
 
}  // namespace Kratos.
 
#endif // MAPPER_VERTEX_MORPHING_MATRIX_FREE_H