Nektar::NekMeshUtils::TetGenInterface Class Reference

Class for interacting with the external library TetGen. More...

#include <TetGenInterface.h>

Public Member Functions
	TetGenInterface ()
	Default constructor. More...
void	InitialMesh (std::map< int, NodeSharedPtr > tgidton, std::vector< Array< OneD, int > > tri)
	Assign parameters for meshing. More...
void	GetNewPoints (int num, std::vector< Array< OneD, NekDouble > > &newp)
	Gets the locations of the Stiener points added by TetGen. More...
void	RefineMesh (std::map< int, NekDouble > delta)
	Refines a previously made tetmesh with node delta information from the Octree. More...
std::vector< Array< OneD, int > >	Extract ()
	Get the list of connectivites of the nodes. More...
void	freetet ()
	Clear previous mesh. More...

Private Attributes
tetgenio	surface
	TetGen objects. More...
tetgenio	output
tetgenio	input
tetgenio	additional

Friends
class	MemoryManager< TetGenInterface >

Detailed Description

Class for interacting with the external library TetGen.

Definition at line 57 of file TetGenInterface.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::TetGenInterface::TetGenInterface ( )

inline

Default constructor.

Definition at line 65 of file TetGenInterface.h.

{};

Member Function Documentation

vector< Array< OneD, int > > Nektar::NekMeshUtils::TetGenInterface::Extract

(

)

Get the list of connectivites of the nodes.

Definition at line 124 of file TetGenInterface.cpp.

{
    vector<Array<OneD, int> > tets;
    for (int i = 0; i < output.numberoftetrahedra; i++)
    {
        Array<OneD, int> tet(4);
        tet[0] = output.tetrahedronlist[i * 4 + 0];
        tet[1] = output.tetrahedronlist[i * 4 + 1];
        tet[2] = output.tetrahedronlist[i * 4 + 2];
        tet[3] = output.tetrahedronlist[i * 4 + 3];
        tets.push_back(tet);
    }
    return tets;
}

void Nektar::NekMeshUtils::TetGenInterface::freetet

(

)

Clear previous mesh.

Definition at line 139 of file TetGenInterface.cpp.

{
    surface.deinitialize();
    input.deinitialize();
    output.deinitialize();
}

void Nektar::NekMeshUtils::TetGenInterface::GetNewPoints	(	int	num,
		std::vector< Array< OneD, NekDouble > > &	newp
	)

Gets the locations of the Stiener points added by TetGen.

Definition at line 95 of file TetGenInterface.cpp.

{
    for (int i = num; i < output.numberofpoints; i++)
    {
        Array<OneD, NekDouble> loc(3);
        loc[0] = output.pointlist[i * 3 + 0];
        loc[1] = output.pointlist[i * 3 + 1];
        loc[2] = output.pointlist[i * 3 + 2];
        newp.push_back(loc);
    }
}

void Nektar::NekMeshUtils::TetGenInterface::InitialMesh	(	std::map< int, NodeSharedPtr >	tgidton,
		std::vector< Array< OneD, int > >	tri
	)

Assign parameters for meshing.

Definition at line 48 of file TetGenInterface.cpp.

References Nektar::iterator, and REAL.

{
    surface.initialize();
    output.initialize();

    // build surface input
    tetgenio::facet *f;
    tetgenio::polygon *p;

    surface.firstnumber    = 0;
    surface.numberofpoints = tgidton.size();
    surface.pointlist      = new REAL[surface.numberofpoints * 3];

    map<int, NodeSharedPtr>::iterator it;
    for (it = tgidton.begin(); it != tgidton.end(); it++)
    {
        Array<OneD, NekDouble> loc = it->second->GetLoc();

        surface.pointlist[it->first * 3 + 0] = loc[0];
        surface.pointlist[it->first * 3 + 1] = loc[1];
        surface.pointlist[it->first * 3 + 2] = loc[2];
    }

    surface.numberoffacets  = tri.size();
    surface.facetlist       = new tetgenio::facet[surface.numberoffacets];
    surface.facetmarkerlist = new int[surface.numberoffacets];

    for (int i = 0; i < tri.size(); i++)
    {
        f                          = &surface.facetlist[i];
        f->numberofpolygons        = 1;
        f->polygonlist             = new tetgenio::polygon[f->numberofpolygons];
        f->numberofholes           = 0;
        f->holelist                = NULL;
        p                          = &f->polygonlist[0];
        p->numberofvertices        = 3;
        p->vertexlist              = new int[p->numberofvertices];
        p->vertexlist[0]           = tri[i][0];
        p->vertexlist[1]           = tri[i][1];
        p->vertexlist[2]           = tri[i][2];
        surface.facetmarkerlist[i] = 0;
    }

    tetrahedralize("pYzqQ", &surface, &output);
}

void Nektar::NekMeshUtils::TetGenInterface::RefineMesh

(

std::map< int, NekDouble >

delta

)

Refines a previously made tetmesh with node delta information from the Octree.

Definition at line 108 of file TetGenInterface.cpp.

References REAL.

{
    input = output;

    input.numberofpointmtrs = 1;

    input.pointmtrlist = new REAL[input.numberofpoints];

    for (int i = 0; i < input.numberofpoints; i++)
    {
        input.pointmtrlist[i] = delta[i];
    }

    tetrahedralize("pYrmzqQO2/7o/120", &input, &output);
}

Friends And Related Function Documentation

friend class MemoryManager< TetGenInterface >

friend

Definition at line 60 of file TetGenInterface.h.

Member Data Documentation

tetgenio Nektar::NekMeshUtils::TetGenInterface::additional

private

Definition at line 96 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::input

private

Definition at line 96 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::output

private

Definition at line 96 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::surface

private

TetGen objects.

Definition at line 96 of file TetGenInterface.h.