tetrahedra_3d_4_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_TETRAHEDRA_3D_4_AUSAS_INCISED_SHAPE_FUNCTIONS)
#define KRATOS_TETRAHEDRA_3D_4_AUSAS_INCISED_SHAPE_FUNCTIONS
 
// System includes
 
// External includes
 
// Project includes
#include "modified_shape_functions/tetrahedra_3d_4_ausas_modified_shape_functions.h"
 
namespace Kratos
{
 
 
 
 
class KRATOS_API(KRATOS_CORE) Tetrahedra3D4AusasIncisedShapeFunctions : public Tetrahedra3D4AusasModifiedShapeFunctions
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(Tetrahedra3D4AusasIncisedShapeFunctions);
 
    // General type definitions
    typedef AusasModifiedShapeFunctions                        BaseType;
    typedef BaseType::GeometryPointerType                      GeometryPointerType;
 
 
    Tetrahedra3D4AusasIncisedShapeFunctions(
        const GeometryPointerType rpInputGeometry,
        const Vector& rNodalDistancesWithExtrapolated,
        const Vector& rExtrapolatedEdgeRatios);
 
    ~Tetrahedra3D4AusasIncisedShapeFunctions();
 
 
 
 
    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, 6> edge_id_for_geometry {{0, 2, 3, 1, 4, 5}};
    const std::array<std::array<size_t,2>, 6> node_ids_for_geometry {{{{0,1}}, {{2,0}}, {{0,3}}, {{1,2}}, {{1,3}}, {{2,3}}}};
 
 
 
    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;
 
 
 
 
 
 
 
    Tetrahedra3D4AusasIncisedShapeFunctions& operator=(Tetrahedra3D4AusasIncisedShapeFunctions const& rOther);
 
    Tetrahedra3D4AusasIncisedShapeFunctions(Tetrahedra3D4AusasIncisedShapeFunctions const& rOther) :
        Tetrahedra3D4AusasModifiedShapeFunctions(rOther.GetInputGeometry(), rOther.GetNodalDistances()),
        mExtraEdgeRatios(rOther.mExtraEdgeRatios) {
    };
 
 
};// class Tetrahedra3D4AusasIncisedShapeFunctions
 
}
#endif /* KRATOS_TETRAHEDRA_3D_4_AUSAS_INCISED_SHAPE_FUNCTIONS defined */