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

TetGenInterface()

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

inline

Default constructor.

Definition at line 65 of file TetGenInterface.h.

References Extract(), freetet(), GetNewPoints(), InitialMesh(), and RefineMesh().

{};

Member Function Documentation

Extract()

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

(

)

Get the list of connectivites of the nodes.

Definition at line 131 of file TetGenInterface.cpp.

Referenced by TetGenInterface().

{
    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;
}

freetet()

void Nektar::NekMeshUtils::TetGenInterface::freetet

(

)

Clear previous mesh.

Definition at line 147 of file TetGenInterface.cpp.

Referenced by TetGenInterface().

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

GetNewPoints()

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 98 of file TetGenInterface.cpp.

References CG_Iterations::loc.

Referenced by TetGenInterface().

{
    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);
    }
}

InitialMesh()

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

Assign parameters for meshing.

Definition at line 47 of file TetGenInterface.cpp.

References CG_Iterations::loc, CellMLToNektar.cellml_metadata::p, and REAL.

Referenced by TetGenInterface().

{
    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;
    }

    string cmd = "pYzqQ";
    char *cstr = new char[cmd.length() + 1];
    strcpy(cstr, cmd.c_str());

    tetrahedralize(cstr, &surface, &output);
}

RefineMesh()

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 111 of file TetGenInterface.cpp.

References REAL.

Referenced by TetGenInterface().

{
    input = output;

    input.numberofpointmtrs = 1;

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

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

    string cmd = "pYrmzq1.1/0QO2/7";
    char *cstr = new char[cmd.length() + 1];
    strcpy(cstr, cmd.c_str());

    tetrahedralize(cstr, &input, &output);
}

Friends And Related Function Documentation

MemoryManager< TetGenInterface >

friend class MemoryManager< TetGenInterface >

friend

Definition at line 60 of file TetGenInterface.h.

Member Data Documentation

additional

tetgenio Nektar::NekMeshUtils::TetGenInterface::additional

private

Definition at line 96 of file TetGenInterface.h.

input

tetgenio Nektar::NekMeshUtils::TetGenInterface::input

private

Definition at line 96 of file TetGenInterface.h.

output

tetgenio Nektar::NekMeshUtils::TetGenInterface::output

private

Definition at line 96 of file TetGenInterface.h.

surface

tetgenio Nektar::NekMeshUtils::TetGenInterface::surface

private

TetGen objects.

Definition at line 96 of file TetGenInterface.h.