cad_tessellation_modeler.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Dagmawi Bekel
//                   Ruben Zorrilla
//
 
#if USE_TRIANGLE_NONFREE_TPL
 
#if !defined(KRATOS_CAD_TESSELLATION_MODELER_INCLUDED)
#define KRATOS_CAD_TESSELLATION_MODELER_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "modeler.h"
#include "utilities/tessellation_utilities/curve_tessellation.h"
#include "geometries/brep_surface.h"
#include "geometries/brep_curve_on_surface.h"
#include "integration/triangle_gauss_legendre_integration_points.h"
 
namespace Kratos {
 
 
 
class KRATOS_API(KRATOS_CORE) CadTessellationModeler : public Modeler {
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(CadTessellationModeler);
 
    typedef std::size_t IndexType;
 
    typedef std::size_t SizeType;
 
    typedef Node NodeType;
 
    typedef Point EmbeddedNodeType;
 
    typedef PointerVector<NodeType> ContainerNodeType;
 
    typedef PointerVector<EmbeddedNodeType> ContainerEmbeddedNodeType;
 
    typedef BrepSurface<ContainerNodeType, ContainerEmbeddedNodeType> BrepSurfaceType;
 
    typedef BrepCurveOnSurface<ContainerNodeType, ContainerEmbeddedNodeType> BrepCurveOnSurfaceType;
 
 
    CadTessellationModeler() : Modeler()
    {
    }
 
    CadTessellationModeler(
        Model& rModel,
        Parameters ModelerParameters = Parameters())
        : Modeler(rModel, ModelerParameters)
        , mpModel(&rModel)
    {
    }
 
    virtual ~CadTessellationModeler() = default;
 
 
 
 
    Modeler::Pointer Create(
        Model& rModel,
        const Parameters ModelParameters) const override;
 
 
    void SetupModelPart() override;
 
 
    std::string Info() const override
    {
        return "CadTessellationModeler";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << Info();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
private:
 
 
 
 
    Model* mpModel = nullptr;
 
 
    friend class Serializer;
 
 
 
 
    std::vector<array_1d<double, 2>> ComputeBoundaryTessellation(const BrepCurveOnSurfaceType& rBoundarySegment);
 
    std::vector<BoundedMatrix<double,3,3>> ComputeSurfaceTriangulation(
        const BrepSurfaceType& rSurfaceGeometry,
        const std::vector<array_1d<double, 2>>& rBoundaryLoop);
 
    std::vector<BoundedMatrix<double,3,3>> InsertGaussPointsExactSurface(
        const BrepSurfaceType& rSurfaceGeometry,
        const std::vector<double>& rPointsCoordinates,
        const std::vector<IndexType>& rTriangleConnectivities);
 
    std::vector<BoundedMatrix<double,3,3>> InsertGaussPointsApproxSurface(
        const BrepSurfaceType& rSurfaceGeometry,
        const std::vector<double>& rPointsCoordinates,
        const std::vector<IndexType>& rTriangleConnectivities);
 
    double ComputeDiscretizationError(
        const std::vector<BoundedMatrix<double,3,3>>& rGaussPointsExact,
        const std::vector<BoundedMatrix<double,3,3>>& rGaussPointsApprox);
 
}; // Class CadTessellationModeler
 
 
inline std::istream& operator>>(std::istream& rIStream,
    CadTessellationModeler& rThis)
{
    return rIStream;
}
 
inline std::ostream& operator<<(std::ostream& rOStream, const CadTessellationModeler& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
} // namespace Kratos.
 
#endif // KRATOS_CAD_TESSELLATION_MODELER_INCLUDED  defined
#endif