adjoint_finite_difference_potential_flow_element.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//
//  Main authors:    Marc Nunez, based on A. Geiser, M. Fusseder, I. Lopez and R. Rossi work
//
 
#if !defined(KRATOS_ADJOINT_FINITE_DIFFERENCE_POTENTIAL_FLOW_ELEMENT_H_INCLUDED )
#define KRATOS_ADJOINT_FINITE_DIFFERENCE_POTENTIAL_FLOW_ELEMENT_H_INCLUDED
 
 
// Project includes
#include "adjoint_base_potential_flow_element.h"
 
namespace Kratos
{
 
template <class TPrimalElement>
class AdjointFiniteDifferencePotentialFlowElement : public AdjointBasePotentialFlowElement<TPrimalElement>
{
public:
 
 
    typedef AdjointBasePotentialFlowElement<TPrimalElement> BaseType;
    using BaseType::NumNodes;
    using BaseType::Dim;
    typedef typename BaseType::GeometryType GeometryType;
    typedef typename BaseType::PropertiesType PropertiesType;
    typedef typename BaseType::NodesArrayType NodesArrayType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(AdjointFiniteDifferencePotentialFlowElement);
 
 
    explicit AdjointFiniteDifferencePotentialFlowElement(IndexType NewId = 0)
     : AdjointBasePotentialFlowElement<TPrimalElement>(NewId){};
 
    AdjointFiniteDifferencePotentialFlowElement(IndexType NewId,
                        typename GeometryType::Pointer pGeometry)
     : AdjointBasePotentialFlowElement<TPrimalElement>(NewId, pGeometry){};
 
    AdjointFiniteDifferencePotentialFlowElement(IndexType NewId,
                        typename GeometryType::Pointer pGeometry,
                        typename PropertiesType::Pointer pProperties)
     : AdjointBasePotentialFlowElement<TPrimalElement>(NewId, pGeometry, pProperties){};
    AdjointFiniteDifferencePotentialFlowElement(AdjointFiniteDifferencePotentialFlowElement const& rOther) {};
 
    ~AdjointFiniteDifferencePotentialFlowElement() override {};
 
 
    AdjointFiniteDifferencePotentialFlowElement & operator=(AdjointFiniteDifferencePotentialFlowElement const& rOther)
    {
        BaseType::operator=(rOther);
        Flags::operator =(rOther);
        return *this;
    }
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, typename PropertiesType::Pointer pProperties) const override;
 
    Element::Pointer Create(IndexType NewId, typename GeometryType::Pointer pGeom, typename PropertiesType::Pointer pProperties) const override;
 
    Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
 
    void CalculateSensitivityMatrix(const Variable<array_1d<double,3> >& rDesignVariable,
                                            Matrix& rOutput,
                                            const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateSensitivityMatrix(const Variable<double>& rDesignVariable,
                                            Matrix& rOutput,
                                            const ProcessInfo& rCurrentProcessInfo) override;
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
 
protected:
 
private:
 
    double GetPerturbationSize();
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
}; // Class AdjointFiniteDifferencePotentialFlowElement
 
 
} // namespace Kratos.
 
#endif