helmholtz_surface_shape_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   license: OptimizationApplication/license.txt
//
//  Main authors:    Reza Najian Asl
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "includes/condition.h"
#include "includes/model_part.h"
#include "utilities/integration_utilities.h"
#include "utilities/geometry_utilities.h"
#include "optimization_application_variables.h"
 
 
namespace Kratos
{
 
 
 
class KRATOS_API(OPTIMIZATION_APPLICATION) HelmholtzSurfaceShapeCondition
    : public Condition
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(HelmholtzSurfaceShapeCondition);
 
 
    HelmholtzSurfaceShapeCondition(
    IndexType NewId,
    GeometryType::Pointer pGeometry);
 
    HelmholtzSurfaceShapeCondition(
    IndexType NewId,
    GeometryType::Pointer pGeometry,
    PropertiesType::Pointer pProperties);
 
    virtual ~HelmholtzSurfaceShapeCondition();
 
 
    Condition::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties
        ) const override;
 
    Condition::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties
        ) const override;
 
    Condition::Pointer Clone (
        IndexType NewId,
        NodesArrayType const& rThisNodes
        ) const override;
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    void GetDofList(
        DofsVectorType& rConditionalDofList,
        const ProcessInfo& rCurrentProcessInfo
        ) const override;
 
    void GetValuesVector(
        Vector& rValues,
        int Step = 0
        ) const override;
 
    void CalculateLocalSystem(
        MatrixType& rLeftHandSideMatrix,
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void CalculateLeftHandSide(
        MatrixType& rLeftHandSideMatrix,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void CalculateRightHandSide(
        VectorType& rRightHandSideVector,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
 
    void CalculateMassMatrix(
        MatrixType& rMassMatrix,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    void Calculate(
        const Variable<double>& rVariable,
        double& rOutput,
        const ProcessInfo& rCurrentProcessInfo
        ) override;
 
    int Check( const ProcessInfo& rCurrentProcessInfo ) const override;
 
 
protected:
    // Protected default constructor necessary for serialization
    HelmholtzSurfaceShapeCondition() : Condition()
    {
    }
 
private:
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Condition);
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Condition);
    }
 
 
 
    void CalculateStiffnessMatrix(
        MatrixType& rStiffnessMatrix,
        const ProcessInfo& rCurrentProcessInfo) const;
 
    void CalculateNormal(
        VectorType & rN) const;
 
    void GetParentElementShapeFunctionsValues(
        MatrixType& rNMatrix,
        const IntegrationMethod& rIntegrationMethod,
        const ProcessInfo& rCurrentProcessInfo) const;
 
    void GetParentElementShapeFunctionsGlobalGradients(
        MatrixType& rDN_DX,
        const IndexType PointNumber,
        const IntegrationMethod& rIntegrationMethod,
        const ProcessInfo& rCurrentProcessInfo) const;
 
}; // Class HelmholtzSurfaceShapeCondition
 
 
}  // namespace Kratos.