GeomElements.cpp File Reference

#include <LibUtilities/Python/NekPyConfig.hpp>
#include <SpatialDomains/Python/SpatialDomains.h>
#include <SpatialDomains/HexGeom.h>
#include <SpatialDomains/MeshGraph.h>
#include <SpatialDomains/PointGeom.h>
#include <SpatialDomains/PrismGeom.h>
#include <SpatialDomains/PyrGeom.h>
#include <SpatialDomains/QuadGeom.h>
#include <SpatialDomains/SegGeom.h>
#include <SpatialDomains/TetGeom.h>
#include <SpatialDomains/TriGeom.h>

Go to the source code of this file.

Typedefs
using	NekError = Nektar::ErrorUtil::NekError

Functions
template<class T , class S >
unique_ptr_objpool< T >	Geometry_Init (int id, py::list &facets)
template<class T , class S >
unique_ptr_objpool< T >	Geometry_Init_Curved (int id, py::list &facets, Curve *curve)
template<class T , class S >
void	export_Geom_2d (py::module &m, const char *name)
template<class T , class S >
void	export_Geom_3d (py::module &m, const char *name)
unique_ptr_objpool< SegGeom >	SegGeom_Init (int id, int coordim, py::list &points, Curve *curve)
unique_ptr_objpool< PointGeom >	PointGeom_Init (int coordim, int id, NekDouble x, NekDouble y, NekDouble z)
py::tuple	PointGeom_GetCoordinates (const PointGeom &self)
void	export_GeomElements (py::module &m)

Typedef Documentation

NekError

using NekError = Nektar::ErrorUtil::NekError

Definition at line 51 of file GeomElements.cpp.

Function Documentation

export_Geom_2d()

template<class T , class S >

void export_Geom_2d	(	py::module &	m,
		const char *	name
	)

Definition at line 90 of file GeomElements.cpp.

{
    py::classh<T, Geometry2D>(m, name)
        .def(py::init<>(&Geometry_Init<T, S>), py::arg("id"),
             py::arg("segments") = py::list())
        .def(py::init<>(&Geometry_Init_Curved<T, S>), py::arg("id"),
             py::arg("segments"), py::arg("curve"));
}

export_Geom_3d()

template<class T , class S >

void export_Geom_3d	(	py::module &	m,
		const char *	name
	)

Definition at line 99 of file GeomElements.cpp.

{
    py::classh<T, Geometry3D>(m, name).def(py::init<>(&Geometry_Init<T, S>),
                                           py::arg("id"),
                                           py::arg("segments") = py::list());
}

export_GeomElements()

void export_GeomElements

(

py::module &

m

)

Definition at line 143 of file GeomElements.cpp.

{
    // Geometry dimensioned base classes
    py::classh<Geometry1D, Geometry>(m, "Geometry1D");
    py::classh<Geometry2D, Geometry>(m, "Geometry2D")
        .def("GetCurve", &Geometry2D::GetCurve,
             py::return_value_policy::reference);
    py::classh<Geometry3D, Geometry>(m, "Geometry3D");
 
    // Point geometries
    py::classh<PointGeom, Geometry>(m, "PointGeom")
        .def(py::init<>(&PointGeom_Init), py::arg("coordim"), py::arg("id"),
             py::arg("x"), py::arg("y"), py::arg("z"))
        .def("GetCoordinates", &PointGeom_GetCoordinates);
 
    // Segment geometries
    py::classh<SegGeom, Geometry>(m, "SegGeom")
        .def(py::init<>(&SegGeom_Init), py::arg("id"), py::arg("coordim"),
             py::arg("points") = py::list(), py::arg("curve") = nullptr)
        .def("GetCurve", &SegGeom::GetCurve,
             py::return_value_policy::reference);
 
    export_Geom_2d<TriGeom, SegGeom>(m, "TriGeom");
    export_Geom_2d<QuadGeom, SegGeom>(m, "QuadGeom");
    export_Geom_3d<TetGeom, TriGeom>(m, "TetGeom");
    export_Geom_3d<PrismGeom, Geometry2D>(m, "PrismGeom");
    export_Geom_3d<PyrGeom, Geometry2D>(m, "PyrGeom");
    export_Geom_3d<HexGeom, QuadGeom>(m, "HexGeom");
}

References Nektar::SpatialDomains::Geometry2D::GetCurve(), Nektar::SpatialDomains::SegGeom::GetCurve(), PointGeom_GetCoordinates(), PointGeom_Init(), and SegGeom_Init().

Referenced by PYBIND11_MODULE().

Geometry_Init()

template<class T , class S >

unique_ptr_objpool< T > Geometry_Init	(	int	id,
		py::list &	facets
	)

Definition at line 54 of file GeomElements.cpp.

{
    std::array<S *, T::kNfacets> geomArr;
 
    if (py::len(facets) != T::kNfacets)
    {
        throw new NekError("Incorrect number of facets for this geometry");
    }
 
    for (int i = 0; i < T::kNfacets; ++i)
    {
        geomArr[i] = py::cast<S *>(facets[i]);
    }
 
    return ObjPoolManager<T>::AllocateUniquePtr(id, geomArr);
}

References Nektar::ObjPoolManager< DataType >::AllocateUniquePtr().

Geometry_Init_Curved()

template<class T , class S >

unique_ptr_objpool< T > Geometry_Init_Curved	(	int	id,
		py::list &	facets,
		Curve *	curve
	)

Definition at line 72 of file GeomElements.cpp.

{
    std::array<S *, T::kNfacets> geomArr;
 
    if (py::len(facets) != T::kNfacets)
    {
        throw new NekError("Incorrect number of facets for this geometry");
    }
 
    for (int i = 0; i < T::kNfacets; ++i)
    {
        geomArr[i] = py::cast<S *>(facets[i]);
    }
 
    return ObjPoolManager<T>::AllocateUniquePtr(id, geomArr, curve);
}

References Nektar::ObjPoolManager< DataType >::AllocateUniquePtr().

PointGeom_GetCoordinates()

py::tuple PointGeom_GetCoordinates

(

const PointGeom &

self

)

Definition at line 138 of file GeomElements.cpp.

{
    return py::make_tuple(self.x(), self.y(), self.z());
}

References Nektar::NekPoint< data_type >::x(), Nektar::NekPoint< data_type >::y(), and Nektar::NekPoint< data_type >::z().

Referenced by export_GeomElements().

PointGeom_Init()

unique_ptr_objpool< PointGeom > PointGeom_Init	(	int	coordim,
		int	id,
		NekDouble	x,
		NekDouble	y,
		NekDouble	z
	)

Definition at line 132 of file GeomElements.cpp.

{
    return ObjPoolManager<PointGeom>::AllocateUniquePtr(coordim, id, x, y, z);
}

References Nektar::ObjPoolManager< DataType >::AllocateUniquePtr().

Referenced by export_GeomElements().

SegGeom_Init()

unique_ptr_objpool< SegGeom > SegGeom_Init	(	int	id,
		int	coordim,
		py::list &	points,
		Curve *	curve
	)

Definition at line 106 of file GeomElements.cpp.

{
    std::array<PointGeom *, 2> geomArr;
 
    if (py::len(points) != 2)
    {
        throw new NekError("Incorrect number of facets for this geometry");
    }
 
    for (int i = 0; i < 2; ++i)
    {
        geomArr[i] = py::cast<PointGeom *>(points[i]);
    }
 
    if (!curve)
    {
        return ObjPoolManager<SegGeom>::AllocateUniquePtr(id, coordim, geomArr);
    }
    else
    {
        return ObjPoolManager<SegGeom>::AllocateUniquePtr(id, coordim, geomArr,
                                                          curve);
    }
}

References Nektar::ObjPoolManager< DataType >::AllocateUniquePtr().

Referenced by export_GeomElements().