edge_based_gradient_recovery_element.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//
//
 
#if !defined(KRATOS_EDGE_BASED_GRADIENT_RECOVIERY_ELEMENT_H_INCLUDED )
#define  KRATOS_EDGE_BASED_GRADIENT_RECOVIERY_ELEMENT_H_INCLUDED
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "includes/element.h"
 
namespace Kratos
{
 
 
 
 
 
 
 
 
template<std::size_t TDim>
class KRATOS_API(KRATOS_CORE) EdgeBasedGradientRecoveryElement : public Element
{
public:
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(EdgeBasedGradientRecoveryElement);
 
    constexpr static std::size_t NumNodes = 2;
 
    constexpr static std::size_t BlockSize = TDim;
 
    constexpr static std::size_t LocalSize = NumNodes * BlockSize;
 
    using NodeType = Node;
 
    using GeometryType = Geometry<NodeType>;
 
    using NodesArrayType = Geometry<NodeType>::PointsArrayType;
 
    using VectorType = Vector;
 
    using MatrixType = Matrix;
 
    using IndexType = Element::IndexType;
    using DofsArrayType = Element::DofsArrayType;
    using DofsVectorType = Element::DofsVectorType;
    using EquationIdVectorType = Element::EquationIdVectorType;
 
 
    //Constructors.
 
 
    EdgeBasedGradientRecoveryElement(IndexType NewId = 0)
        : Element(NewId)
    {}
 
 
    EdgeBasedGradientRecoveryElement(
        IndexType NewId,
        const NodesArrayType& ThisNodes)
        : Element(NewId, ThisNodes)
    {}
 
 
    EdgeBasedGradientRecoveryElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry)
        : Element(NewId, pGeometry)
    {}
 
 
    EdgeBasedGradientRecoveryElement(
        IndexType NewId,
        GeometryType::Pointer pGeometry,
        PropertiesType::Pointer pProperties)
        : Element(NewId, pGeometry, pProperties)
    {}
 
    ~EdgeBasedGradientRecoveryElement() override = default;
 
 
 
 
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties) const override
    {
        return Kratos::make_intrusive<EdgeBasedGradientRecoveryElement<TDim>>(NewId, GetGeometry().Create(ThisNodes), pProperties);
    }
 
    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;
 
    int Check(const ProcessInfo &rCurrentProcessInfo) const override;
 
    void EquationIdVector(
        EquationIdVectorType& rResult,
        const ProcessInfo& rCurrentProcessInfo) const override;
 
    void GetDofList(
        DofsVectorType& rElementalDofList,
        const ProcessInfo& rCurrentProcessInfo) const override;
 
 
 
 
 
 
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "EdgeBasedGradientRecoveryElement #" << Id();
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "EdgeBasedGradientRecoveryElement";
    }
 
    virtual void PrintData(std::ostream& rOStream) const override
    {}
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
 
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, Element );
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, Element);
    }
 
 
 
 
 
 
 
 
 
 
    EdgeBasedGradientRecoveryElement & operator=(EdgeBasedGradientRecoveryElement const& rOther) = delete;
 
    EdgeBasedGradientRecoveryElement(EdgeBasedGradientRecoveryElement const& rOther) = delete;
 
 
}; // Class EdgeBasedGradientRecoveryElement
 
 
 
 
 
template<std::size_t TDim>
inline std::istream &operator>>(
    std::istream &rIStream,
    EdgeBasedGradientRecoveryElement<TDim> &rThis)
{
    return rIStream;
}
 
template<std::size_t TDim>
inline std::ostream &operator<<(
    std::ostream &rOStream,
    const EdgeBasedGradientRecoveryElement<TDim> &rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
} // namespace Kratos.
 
#endif // KRATOS_EDGE_BASED_GRADIENT_RECOVIERY_ELEMENT_H_INCLUDED  defined