triangle_2d_3_ausas_incised_shape_functions.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Franziska Wahl
//
 
#if !defined(KRATOS_TRIANGLE_2D_3_AUSAS_INCISED_SHAPE_FUNCTIONS)
#define KRATOS_TRIANGLE_2D_3_AUSAS_INCISED_SHAPE_FUNCTIONS
 
// System includes
 
// External includes
 
// Project includes
#include "modified_shape_functions/triangle_2d_3_ausas_modified_shape_functions.h"
 
namespace Kratos
{
 
 
 
 
class KRATOS_API(KRATOS_CORE) Triangle2D3AusasIncisedShapeFunctions : public Triangle2D3AusasModifiedShapeFunctions
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(Triangle2D3AusasIncisedShapeFunctions);
 
    // General type definitions
    typedef AusasModifiedShapeFunctions                        BaseType;
    typedef BaseType::GeometryPointerType                      GeometryPointerType;
 
 
    Triangle2D3AusasIncisedShapeFunctions(const GeometryPointerType rpInputGeometry,
        const Vector& rNodalDistancesWithExtrapolated, const Vector& rExtrapolatedEdgeRatios);
 
    ~Triangle2D3AusasIncisedShapeFunctions();
 
 
 
 
    std::string Info() const override;
 
    void PrintInfo(std::ostream& rOStream) const override;
 
    void PrintData(std::ostream& rOStream) const override;
 
 
 
    const Vector& GetExtrapolatedEdgeRatios() const ;
 
 
protected:
 
 
    // Arrays to get edge and node IDs of geometry from edge ID of splitting utility
    const std::array<size_t, 3> edge_id_for_geometry {{2, 0, 1}};
    const std::array<std::array<size_t,2>, 3> node_ids_for_geometry {{{{0,1}}, {{1,2}}, {{2,0}}}};
 
 
 
    void SetPositiveSideCondensationMatrix(Matrix& rPosSideCondMatrix) override;
 
    void SetNegativeSideCondensationMatrix(Matrix& rNegSideCondMatrix) override;
 
    void SetPositiveSideCondensationMatrix(
        Matrix& rPosSideCondMatrix,
        const std::vector<int>& rEdgeNodeI,
        const std::vector<int>& rEdgeNodeJ,
        const std::vector<int>& rSplitEdges);
 
    void SetNegativeSideCondensationMatrix(
        Matrix& rNegSideCondMatrix,
        const std::vector<int>& rEdgeNodeI,
        const std::vector<int>& rEdgeNodeJ,
        const std::vector<int>& rSplitEdges);
 
 
 
 
 
private:
 
 
    const Vector mExtraEdgeRatios;
 
 
 
 
 
 
 
    Triangle2D3AusasIncisedShapeFunctions& operator=(Triangle2D3AusasModifiedShapeFunctions const& rOther);
 
    Triangle2D3AusasIncisedShapeFunctions(Triangle2D3AusasIncisedShapeFunctions const& rOther) :
        Triangle2D3AusasModifiedShapeFunctions(rOther.GetInputGeometry(), rOther.GetNodalDistances()),
        mExtraEdgeRatios(rOther.mExtraEdgeRatios) {
    };
 
 
};// class Triangle2D3AusasIncisedShapeFunctions
 
}
#endif /* KRATOS_TRIANGLE_2D_3_AUSAS_INCISED_SHAPE_FUNCTIONS defined */