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 58 of file TetGenInterface.h.

Constructor & Destructor Documentation

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

inline

Default constructor.

Definition at line 66 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 132 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 148 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 99 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 CellMLToNektar.cellml_metadata::p.

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

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

{
    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

friend class MemoryManager< TetGenInterface >

friend

Definition at line 61 of file TetGenInterface.h.

Member Data Documentation

tetgenio Nektar::NekMeshUtils::TetGenInterface::additional

private

Definition at line 97 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::input

private

Definition at line 97 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::output

private

Definition at line 97 of file TetGenInterface.h.

tetgenio Nektar::NekMeshUtils::TetGenInterface::surface

private

TetGen objects.

Definition at line 97 of file TetGenInterface.h.