point_3d.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Riccardo Rossi
//                   Janosch Stascheit
//                   Felix Nagel
//  contributors:    Hoang Giang Bui
//                   Josep Maria Carbonell
//
 
#if !defined(KRATOS_POINT_3D_H_INCLUDED )
#define  KRATOS_POINT_3D_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "geometries/geometry.h"
 
namespace Kratos
{
 
 
 
 
 
 
template<class TPointType>
class Point3D : public Geometry<TPointType>
{
public:
 
    typedef Geometry<TPointType> BaseType;
 
    KRATOS_CLASS_POINTER_DEFINITION(Point3D);
 
    typedef GeometryData::IntegrationMethod IntegrationMethod;
 
    typedef typename BaseType::GeometriesArrayType GeometriesArrayType;
 
    typedef TPointType PointType;
 
    typedef typename BaseType::IndexType IndexType;
 
 
    typedef typename BaseType::SizeType SizeType;
 
    typedef  typename BaseType::PointsArrayType PointsArrayType;
 
    typedef typename BaseType::IntegrationPointType IntegrationPointType;
 
    typedef typename BaseType::IntegrationPointsArrayType IntegrationPointsArrayType;
 
    typedef typename BaseType::IntegrationPointsContainerType IntegrationPointsContainerType;
 
    typedef typename BaseType::ShapeFunctionsValuesContainerType ShapeFunctionsValuesContainerType;
 
    typedef typename BaseType::ShapeFunctionsLocalGradientsContainerType ShapeFunctionsLocalGradientsContainerType;
 
    typedef typename BaseType::JacobiansType JacobiansType;
 
    typedef typename BaseType::ShapeFunctionsGradientsType ShapeFunctionsGradientsType;
 
    typedef typename BaseType::NormalType NormalType;
 
    typedef typename BaseType::CoordinatesArrayType CoordinatesArrayType;
 
 
//     Point3D(const PointType& FirstPoint)
//         : BaseType(PointsArrayType(), &msGeometryData)
//     {
//         BaseType::Points().push_back(typename PointType::Pointer(new PointType(FirstPoint)));
//     }
 
    Point3D(typename PointType::Pointer pFirstPoint)
        : BaseType(PointsArrayType(), &msGeometryData)
    {
        BaseType::Points().push_back(pFirstPoint);
    }
 
    Point3D(const PointsArrayType& ThisPoints)
        : BaseType(ThisPoints, &msGeometryData)
    {
        if( BaseType::PointsNumber() != 1)
            KRATOS_ERROR << "Invalid points number. Expected 2, given " << BaseType::PointsNumber() << std::endl;
    }
 
    explicit Point3D(
        const IndexType GeometryId,
        const PointsArrayType& rThisPoints
    ) : BaseType(GeometryId, rThisPoints, &msGeometryData)
    {
        KRATOS_ERROR_IF( this->PointsNumber() != 1 ) << "Invalid points number. Expected 1, given " << this->PointsNumber() << std::endl;
    }
 
    explicit Point3D(
        const std::string& rGeometryName,
        const PointsArrayType& rThisPoints
    ) : BaseType(rGeometryName, rThisPoints, &msGeometryData)
    {
        KRATOS_ERROR_IF(this->PointsNumber() != 1) << "Invalid points number. Expected 1, given " << this->PointsNumber() << std::endl;
    }
 
    Point3D(Point3D const& rOther)
        : BaseType(rOther)
    {
    }
 
 
    template<class TOtherPointType> Point3D(Point3D<TOtherPointType> const& rOther)
        : BaseType(rOther)
    {
    }
 
    ~Point3D() override {}
 
    GeometryData::KratosGeometryFamily GetGeometryFamily() const override
    {
        return GeometryData::KratosGeometryFamily::Kratos_Point;
    }
    GeometryData::KratosGeometryType GetGeometryType() const override
    {
        return GeometryData::KratosGeometryType::Kratos_Point3D;
    }
 
 
    Point3D& operator=(const Point3D& rOther)
    {
        BaseType::operator=(rOther);
 
        return *this;
    }
 
    template<class TOtherPointType>
    Point3D& operator=(Point3D<TOtherPointType> const & rOther)
    {
        BaseType::operator=(rOther);
 
        return *this;
    }
 
 
    typename BaseType::Pointer Create(
        const IndexType NewGeometryId,
        PointsArrayType const& rThisPoints
        ) const override
    {
        return typename BaseType::Pointer( new Point3D( NewGeometryId, rThisPoints ) );
    }
 
    typename BaseType::Pointer Create(
        const IndexType NewGeometryId,
        const BaseType& rGeometry
    ) const override
    {
        auto p_geometry = typename BaseType::Pointer( new Point3D( NewGeometryId, rGeometry.Points() ) );
        p_geometry->SetData(rGeometry.GetData());
        return p_geometry;
    }
 
//     /**
//      * @brief Lumping factors for the calculation of the lumped mass matrix
//      * @param rResult Vector containing the lumping factors
//      * @param LumpingMethod The lumping method considered. The three methods available are:
//      *      - The row sum method
//      *      - Diagonal scaling
//      *      - Evaluation of M using a quadrature involving only the nodal points and thus automatically yielding a diagonal matrix for standard element shape function
//      */
//     Vector& LumpingFactors(
//         Vector& rResult,
//         const typename BaseType::LumpingMethods LumpingMethod = BaseType::LumpingMethods::ROW_SUM
//         )  const override
//  {
//  }
 
 
    double Length() const override
    {
        return 0.00;
    }
 
    double Area() const override
    {
        return 0.00;
    }
 
 
    double DomainSize() const override
    {
        return 0.00;
    }
 
 
    //         virtual JacobiansType& Jacobian(JacobiansType& rResult, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
    //    virtual Matrix& Jacobian(Matrix& rResult, IndexType IntegrationPointIndex, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
//
    //     virtual Vector& DeterminantOfJacobian(Vector& rResult, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
    //     virtual double DeterminantOfJacobian(IndexType IntegrationPointIndex, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
//       virtual double DeterminantOfJacobian(const CoordinatesArrayType& rPoint) const
//  {
//  }
 
 
    //      virtual JacobiansType& InverseOfJacobian(JacobiansType& rResult, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
    //      virtual Matrix& InverseOfJacobian(Matrix& rResult, IndexType IntegrationPointIndex, IntegrationMethod ThisMethod) const
    //  {
    //  }
 
    //    virtual Matrix& InverseOfJacobian(Matrix& rResult, const CoordinatesArrayType& rPoint) const
    //  {
    //  }
 
 
    SizeType EdgesNumber() const override
    {
        return 1;
    }
 
    SizeType FacesNumber() const override
    {
        return 0;
    }
 
 
//      virtual double ShapeFunctionValue(IndexType ShapeFunctionIndex,
//                                            const CoordinatesArrayType& rPoint) const
//     {
//          KRATOS_ERROR << "This method is not implemented yet!" << *this << std::endl;
//       return 0;
//          }
//
    /*
 
 
         virtual void BoundingBox(Bounding_Box<TPointType>& rResult) const
                 {
                    this->BoundingBox(rResult.mLowPoint, rResult.mHighPoint);
                 }
 
 
          virtual void BoundingBox(TPointType& rLowPoint, TPointType& rHighPoint) const
                 {
 
                    array_1d<double,3 > point = ZeroVector(3);
 
            point[0] =  BaseType::GetPoint(0).X();
            point[1] =  BaseType::GetPoint(0).Y();
            point[2] =  BaseType::GetPoint(0).Z();
 
            rLowPoint.push_back(point);
            rHighPoint.push_back(point);
 
                }
 
    */
 
 
    std::string Info() const override
    {
        return "a point in 3D space";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "a point in 3D space";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
        rOStream << "a point in 3D space";
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
    static const GeometryData msGeometryData;
 
    static const GeometryDimension msGeometryDimension;
 
 
 
 
 
 
//      static Matrix CalculateShapeFunctionsIntegrationPointsValues(typename BaseType::IntegrationMethod ThisMethod)
//     {
//      }
 
    //      static ShapeFunctionsGradientsType CalculateShapeFunctionsIntegrationPointsLocalGradients(typename BaseType::IntegrationMethod ThisMethod)
    //  {
    //  }
 
    static const IntegrationPointsContainerType AllIntegrationPoints()
    {
        IntegrationPointsContainerType integration_points;
        return integration_points;
    }
 
    static const ShapeFunctionsValuesContainerType AllShapeFunctionsValues()
    {
        ShapeFunctionsValuesContainerType shape_functions_value;
        return shape_functions_value;
    }
    static const ShapeFunctionsLocalGradientsContainerType AllShapeFunctionsLocalGradients()
    {
        ShapeFunctionsLocalGradientsContainerType shape_functions_local_gradients;
        return shape_functions_local_gradients;
    }
 
 
 
 
 
 
    template<class TOtherPointType> friend class Point3D;
 
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, BaseType );
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, BaseType );
    }
 
    // Default constructor needed for serialization only
    Point3D():BaseType( PointsArrayType(), &msGeometryData ) {}
 
 
 
 
 
 
}; // Class Geometry
 
 
 
 
 
 
template<class TPointType>
inline std::istream& operator >> (std::istream& rIStream,
                                  Point3D<TPointType>& rThis);
 
template<class TPointType>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const Point3D<TPointType>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
template<class TPointType>
const GeometryData Point3D<TPointType>::msGeometryData(
        &msGeometryDimension,
        GeometryData::IntegrationMethod::GI_GAUSS_1,
        Point3D<TPointType>::AllIntegrationPoints(),
        Point3D<TPointType>::AllShapeFunctionsValues(),
        AllShapeFunctionsLocalGradients());
 
template<class TPointType>
const GeometryDimension Point3D<TPointType>::msGeometryDimension(3, 0);
 
}  // namespace Kratos.
 
#endif // KRATOS_LINE_2D_H_INCLUDED  defined