explicit_filter.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   license: OptimizationApplication/license.txt
//
//  Main authors:    Baumgaertner Daniel
//                   Aditya Ghantasala
//                   Suneth Warnakulasuriya
//
 
#pragma once
 
// System includes
#include <string>
#include <variant>
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
#include "spatial_containers/spatial_containers.h"
#include "expression/container_expression.h"
 
// Application includes
#include "entity_point.h"
#include "filter_function.h"
#include "damping_function.h"
 
namespace Kratos {
 
 
template<class TContainerType>
class KRATOS_API(OPTIMIZATION_APPLICATION) ExplicitFilter
{
public:
 
    using EntityType = typename TContainerType::value_type;
 
    using EntityPointVector = std::vector<typename EntityPoint<EntityType>::Pointer>;
 
    // Type definitions for tree-search
    using BucketType = Bucket<3, EntityPoint<EntityType>, EntityPointVector>;
 
    using KDTree = Tree<KDTreePartition<BucketType>>;
 
    KRATOS_CLASS_POINTER_DEFINITION(ExplicitFilter);
 
 
    ExplicitFilter(
        const ModelPart& rModelPart,
        const std::string& rKernelFunctionType,
        const IndexType MaxNumberOfNeighbours);
 
    ExplicitFilter(
        const ModelPart& rModelPart,
        const ModelPart& rFixedModelPart,
        const std::string& rKernelFunctionType,
        const std::string& rDampingFunctionType,
        const IndexType MaxNumberOfNeighbours);
 
 
    void SetFilterRadius(const ContainerExpression<TContainerType>& rContainerExpression);
 
    void Update();
 
    ContainerExpression<TContainerType> FilterField(const ContainerExpression<TContainerType>& rContainerExpression) const;
 
    ContainerExpression<TContainerType> FilterIntegratedField(const ContainerExpression<TContainerType>& rContainerExpression) const;
 
    void GetIntegrationWeights(ContainerExpression<TContainerType>& rContainerExpression) const;
 
    std::string Info() const;
 
private:
 
    const ModelPart& mrModelPart;
 
    const ModelPart* mpFixedModelPart = nullptr;
 
    FilterFunction::UniquePointer mpKernelFunction;
 
    DampingFunction::UniquePointer mpDampingFunction;
 
    typename ContainerExpression<TContainerType>::Pointer mpFilterRadiusContainer;
 
    Expression::ConstPointer mpNodalDomainSizeExpression;
 
    EntityPointVector mEntityPointVector;
 
    EntityPointVector mFixedModelPartEntityPointVector;
 
    IndexType mBucketSize = 100;
 
    IndexType mMaxNumberOfNeighbors;
 
    typename KDTree::Pointer mpSearchTree;
 
    typename KDTree::Pointer mpFixedModelPartSearchTree;
 
 
    template<class TWeightIntegrationType>
    ContainerExpression<TContainerType> GenericFilterField(const ContainerExpression<TContainerType>& rContainerExpression) const;
 
};
 
 
template<class TContainerType>
inline std::ostream& operator<<(
    std::ostream& rOStream,
    const ExplicitFilter<TContainerType>& rThis)
{
    return rOStream << rThis.Info();
}
 
 
} // namespace Kratos