support_lagrange_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
 
#if !defined(KRATOS_SUPPORT_LAGRANGE_CONDITION_H_INCLUDED )
#define  KRATOS_SUPPORT_LAGRANGE_CONDITION_H_INCLUDED
 
// System includes
#include "includes/define.h"
#include "includes/condition.h"
 
// External includes
 
// Project includes
#include "iga_application_variables.h"
 
 
namespace Kratos
{
 
 
    class SupportLagrangeCondition
        : public Condition
    {
    public:
 
        KRATOS_CLASS_POINTER_DEFINITION(SupportLagrangeCondition);
 
        typedef std::size_t SizeType;
        typedef std::size_t IndexType;
 
 
        SupportLagrangeCondition(
            IndexType NewId,
            GeometryType::Pointer pGeometry)
            : Condition(NewId, pGeometry)
        {};
 
        SupportLagrangeCondition(
            IndexType NewId,
            GeometryType::Pointer pGeometry,
            PropertiesType::Pointer pProperties)
            : Condition(NewId, pGeometry, pProperties)
        {};
 
        SupportLagrangeCondition()
            : Condition()
        {};
 
        virtual ~SupportLagrangeCondition() = default;
 
 
        Condition::Pointer Create(
            IndexType NewId,
            GeometryType::Pointer pGeom,
            PropertiesType::Pointer pProperties
            ) const override
        {
            return Kratos::make_intrusive<SupportLagrangeCondition>(
                NewId, pGeom, pProperties);
        };
 
        Condition::Pointer Create(
            IndexType NewId,
            NodesArrayType const& ThisNodes,
            PropertiesType::Pointer pProperties
            ) const override
        {
            return Kratos::make_intrusive< SupportLagrangeCondition >(
                NewId, GetGeometry().Create(ThisNodes), pProperties);
        };
 
 
        void CalculateRightHandSide(
            VectorType& rRightHandSideVector,
            const ProcessInfo& rCurrentProcessInfo) override
        {
            const SizeType mat_size = GetNumberOfNonZeroNodes() * 6;
 
            if (rRightHandSideVector.size() != mat_size)
                rRightHandSideVector.resize(mat_size);
            noalias(rRightHandSideVector) = ZeroVector(mat_size);
 
            MatrixType left_hand_side_matrix;
 
            CalculateAll(left_hand_side_matrix, rRightHandSideVector,
                rCurrentProcessInfo, false, true);
        }
 
        void CalculateLeftHandSide(
            MatrixType& rLeftHandSideMatrix,
            const ProcessInfo& rCurrentProcessInfo) override
        {
            const SizeType mat_size = GetNumberOfNonZeroNodes() * 6;
 
            VectorType right_hand_side_vector;
 
            if (rLeftHandSideMatrix.size1() != mat_size && rLeftHandSideMatrix.size2())
                rLeftHandSideMatrix.resize(mat_size, mat_size);
            noalias(rLeftHandSideMatrix) = ZeroMatrix(mat_size, mat_size);
 
            CalculateAll(rLeftHandSideMatrix, right_hand_side_vector,
                rCurrentProcessInfo, true, false);
        }
 
        void CalculateLocalSystem(
            MatrixType& rLeftHandSideMatrix,
            VectorType& rRightHandSideVector,
            const ProcessInfo& rCurrentProcessInfo) override
        {
            const SizeType mat_size = GetNumberOfNonZeroNodes() * 6;
 
            if (rRightHandSideVector.size() != mat_size)
                rRightHandSideVector.resize(mat_size);
            noalias(rRightHandSideVector) = ZeroVector(mat_size);
 
            if (rLeftHandSideMatrix.size1() != mat_size)
                rLeftHandSideMatrix.resize(mat_size, mat_size);
            noalias(rLeftHandSideMatrix) = ZeroMatrix(mat_size, mat_size);
 
            CalculateAll(rLeftHandSideMatrix, rRightHandSideVector,
                rCurrentProcessInfo, true, true);
        }
 
        void EquationIdVector(
            EquationIdVectorType& rResult,
            const ProcessInfo& rCurrentProcessInfo
            ) const override;
 
        void GetDofList(
            DofsVectorType& rElementalDofList,
            const ProcessInfo& rCurrentProcessInfo
            ) const override;
 
        void CalculateAll(
            MatrixType& rLeftHandSideMatrix,
            VectorType& rRightHandSideVector,
            const ProcessInfo& rCurrentProcessInfo,
            const bool CalculateStiffnessMatrixFlag,
            const bool CalculateResidualVectorFlag
            );
 
        /* @brief This function computes the determinants of jacobian for
         *         the initial configuration. Required for conservative integration.
         * @param rGeometry: corresponding geometry
         * @param rDeterminantOfJacobian: output determinants
         */
        void DeterminantOfJacobianInitial(
            const GeometryType& rGeometry,
            Vector& rDeterminantOfJacobian);
 
 
        std::string Info() const override
        {
            std::stringstream buffer;
            buffer << "\"SupportLagrangeCondition\" #" << Id();
            return buffer.str();
        }
 
        void PrintInfo(std::ostream& rOStream) const override
        {
            rOStream << "\"SupportLagrangeCondition\" #" << Id();
        }
 
        void PrintData(std::ostream& rOStream) const override {
            pGetGeometry()->PrintData(rOStream);
        }
 
 
    private:
 
        const double shape_function_tolerance = 1e-6;
 
 
        /* @brief checks all shape functions and
         *        returns the number of non zero nodes.
         */
        SizeType GetNumberOfNonZeroNodes() const;
 
 
        friend class Serializer;
 
        virtual void save(Serializer& rSerializer) const override
        {
            KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Condition);
        }
 
        virtual void load(Serializer& rSerializer) override
        {
            KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Condition);
        }
 
 
    }; // Class SupportLagrangeCondition
 
}  // namespace Kratos.
 
#endif // KRATOS_SUPPORT_LAGRANGE_CONDITION_H_INCLUDED  defined