max_overhang_response_utils.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   license: OptimizationApplication/license.txt
//
//  Main author:     Reza Najian Asl
//
 
#pragma once
 
// System includes
#include <unordered_map>
#include <variant>
#include <vector>
 
// Project includes
#include "includes/define.h"
#include "includes/model_part.h"
 
// Application includes
 
namespace Kratos
{
 
 
class KRATOS_API(OPTIMIZATION_APPLICATION) MaxOverhangAngleResponseUtils
{
public:
 
    using IndexType = std::size_t;
 
    using GeometryType = ModelPart::ElementType::GeometryType;
 
    using array_3d = array_1d<double, 3>;
 
    using SensitivityFieldVariableTypes = std::variant<const Variable<double>*, const Variable<array_3d>*>;
 
    using SensitivityVariableModelPartsListMap = std::unordered_map<SensitivityFieldVariableTypes, std::vector<ModelPart*>>;
 
 
    static double CalculateValue(const std::vector<ModelPart*>& rModelParts, 
                                 const Parameters ResponseSettings);
 
    static void CalculateSensitivity(
        const std::vector<ModelPart*>& rEvaluatedModelParts,
        const SensitivityVariableModelPartsListMap& rSensitivityVariableModelPartInfo,
        const Parameters ResponseSettings);
 
private:
 
    static double CalculateConditionValue(const Condition& rCondition, const Parameters ResponseSettings);
 
    static void CalculateConditionFiniteDifferenceShapeSensitivity(
        Condition& rCondition,
        Condition::Pointer& pThreadLocalCondition,
        ModelPart& rModelPart,
        std::vector<std::string>& rModelPartNames,
        const Parameters ResponseSettings,
        const IndexType MaxNodeId,
        const Variable<array_3d>& rOutputSensitivityVariable);
 
    static void CalculateFiniteDifferenceShapeSensitivity(
       ModelPart& rModelPart,
       const Parameters ResponseSettings,
       const Variable<array_3d>& rOutputSensitivityVariable);
 
};
 
}