Nektar::NekMeshUtils::TriangleInterface Class Reference

class for interfacing with external library triangle More...

#include <TriangleInterface.h>

Collaboration diagram for Nektar::NekMeshUtils::TriangleInterface:

Public Member Functions
	TriangleInterface ()
	default constructor More...
void	Assign (std::vector< std::vector< NodeSharedPtr > > &boundingloops, std::vector< Array< OneD, NekDouble > > &centers, int i, NekDouble str=1.0)
	assign meshing paramters More...
void	AssignStiener (std::vector< NodeSharedPtr > stiner)
void	Mesh (bool Quality=false)
	Execute meshing. More...
void	Extract (std::vector< std::vector< NodeSharedPtr > > &Connec)
	Extract mesh. More...

Private Member Functions
void	SetUp ()
	Clear memory. More...

Private Attributes
std::vector< std::vector < NodeSharedPtr > >	m_boundingloops
	List of bounding nodes to the surface. More...
std::vector< NodeSharedPtr >	m_stienerpoints
	List of additional nodes. More...
std::vector< Array< OneD, NekDouble > >	m_centers
	Coordinates of the centers of the loops. More...
std::map< int, NodeSharedPtr >	nodemap
	Map from NekMesh id to triangle id. More...
int	sid
	ID of the surface. More...
NekDouble	m_str
	Stretching factor of parameter plane. More...
DelaunayTriangle	dt
	Triangle data strucutres. More...

Friends
class	MemoryManager< TriangleInterface >

Detailed Description

class for interfacing with external library triangle

Definition at line 55 of file TriangleInterface.h.

Constructor & Destructor Documentation

Nektar::NekMeshUtils::TriangleInterface::TriangleInterface ( )

inline

default constructor

Definition at line 63 of file TriangleInterface.h.

{};

Member Function Documentation

void Nektar::NekMeshUtils::TriangleInterface::Assign ( std::vector< std::vector< NodeSharedPtr > > &  boundingloops, std::vector< Array< OneD, NekDouble > > &  centers, int  i, NekDouble  str = 1.0  )

inline

assign meshing paramters

Definition at line 68 of file TriangleInterface.h.

References m_boundingloops, m_centers, m_str, and sid.

    {
        m_boundingloops = boundingloops;
        m_centers       = centers;
        m_str           = str;
        sid             = i;
    }

void Nektar::NekMeshUtils::TriangleInterface::AssignStiener ( std::vector< NodeSharedPtr >  stiner )

inline

Definition at line 79 of file TriangleInterface.h.

References m_stienerpoints.

    {
        m_stienerpoints = stiner;
    }

void Nektar::NekMeshUtils::TriangleInterface::Extract

(

std::vector< std::vector< NodeSharedPtr > > &

Connec

)

Extract mesh.

Definition at line 193 of file TriangleInterface.cpp.

References ASSERTL0, and Nektar::iterator.

{
    Connec.clear();
    for (int i = 0; i < dt.out.numberoftriangles; i++)
    {
        map<int, NodeSharedPtr>::iterator n1, n2, n3;
        n1 = nodemap.find(dt.out.trianglelist[i * 3 + 0]);
        n2 = nodemap.find(dt.out.trianglelist[i * 3 + 1]);
        n3 = nodemap.find(dt.out.trianglelist[i * 3 + 2]);

        ASSERTL0(n1 != nodemap.end() && n2 != nodemap.end() &&
                     n3 != nodemap.end(),
                 "node index error");

        vector<NodeSharedPtr> tri(3);
        tri[0] = n1->second;
        tri[1] = n2->second;
        tri[2] = n3->second;
        Connec.push_back(tri);
    }
}

void Nektar::NekMeshUtils::TriangleInterface::Mesh

(

bool

Quality = false

)

Execute meshing.

Definition at line 46 of file TriangleInterface.cpp.

References ASSERTL0.

{
    SetUp();

    int numPoints = 0;
    int numSeg = 0;
    for (int i = 0; i < m_boundingloops.size(); i++)
    {
        numSeg += m_boundingloops[i].size();
    }
    numPoints = numSeg + m_stienerpoints.size();

    stringstream ss;
    ss << "3 points required for triangulation, " << numPoints << " provided";

    ASSERTL0(numPoints > 2, ss.str());

    dt.in.numberofpoints          = numPoints;
    dt.in.numberofpointattributes = 0;
    dt.in.pointlist               = new double[dt.in.numberofpoints * 2];
    dt.in.pointmarkerlist         = new int [dt.in.numberofpoints];

    int pointc = 0;

    for (int i = 0; i < m_boundingloops.size(); i++)
    {
        for (int j = 0; j < m_boundingloops[i].size(); j++, pointc++)
        {
            nodemap[pointc] = m_boundingloops[i][j];

            Array<OneD, NekDouble> uv =
                m_boundingloops[i][j]->GetCADSurfInfo(sid);
            dt.in.pointlist[pointc * 2 + 0] = uv[0] * m_str;
            dt.in.pointlist[pointc * 2 + 1] = uv[1];
        }
    }

    for (int i = 0; i < m_stienerpoints.size(); i++, pointc++)
    {
        nodemap[pointc] = m_stienerpoints[i];

        Array<OneD, NekDouble> uv    = m_stienerpoints[i]->GetCADSurfInfo(sid);
        dt.in.pointlist[pointc * 2 + 0] = uv[0] * m_str;
        dt.in.pointlist[pointc * 2 + 1] = uv[1];
    }

    dt.in.numberofsegments = numSeg;
    dt.in.segmentlist      = new int[dt.in.numberofsegments * 2];
    pointc = 0;
    for (int i = 0; i < m_boundingloops.size(); i++, pointc++)
    {
        int first = pointc;
        for (int j = 0; j < m_boundingloops[i].size() - 1; j++, pointc++)
        {
            dt.in.segmentlist[pointc * 2 + 0] = pointc;
            dt.in.segmentlist[pointc * 2 + 1] = pointc + 1;
        }
        dt.in.segmentlist[pointc * 2 + 0] = pointc;
        dt.in.segmentlist[pointc * 2 + 1] = first;
    }

    dt.in.numberofregions = 0;
    dt.in.numberofholes   = m_centers.size() - 1;
    dt.in.holelist        = new double[dt.in.numberofholes * 2];

    for (int i = 1; i < m_centers.size(); i++)
    {
        dt.in.holelist[(i - 1) * 2 + 0] = m_centers[i][0] * m_str;
        dt.in.holelist[(i - 1) * 2 + 1] = m_centers[i][1];
    }

    string cmd;
    if (Quality)
    {
        cmd = "pqY";
    }
    else if (!Quality)
    {
        cmd = "pY";
    }
    char *cstr = new char[cmd.length() + 1];
    strcpy(cstr, cmd.c_str());

    dt.triangulate(cstr);
}

void Nektar::NekMeshUtils::TriangleInterface::SetUp ( )

private

Clear memory.

Definition at line 132 of file TriangleInterface.cpp.

{
    dt.in.pointlist                   = (double *)NULL;
    dt.in.pointattributelist          = (double *)NULL;
    dt.in.pointmarkerlist             = (int *)NULL;
    dt.in.numberofpoints              = 0;
    dt.in.numberofpointattributes     = 0;
    //
    dt.in.trianglelist                = (int *)NULL;
    dt.in.triangleattributelist       = (double *)NULL;
    dt.in.trianglearealist            = (double *)NULL;
    dt.in.neighborlist                = (int *)NULL;
    dt.in.numberoftriangles           = 0;
    dt.in.numberofcorners             = 0;
    dt.in.numberoftriangleattributes  = 0;
    //
    dt.in.segmentlist                 = (int *)NULL;
    dt.in.segmentmarkerlist           = (int *)NULL;
    dt.in.numberofsegments            = 0;
    //
    dt.in.holelist                    = (double *)NULL;
    dt.in.numberofholes               = 0;
    //
    dt.in.regionlist                  = (double *)NULL;
    dt.in.numberofregions             = 0;
    //
    dt.in.edgelist                    = (int *)NULL;
    dt.in.edgemarkerlist              = (int *)NULL;
    dt.in.normlist                    = (double *)NULL;
    dt.in.numberofedges               = 0;
    //
    dt.out.pointlist                  = (double *)NULL;
    dt.out.pointattributelist         = (double *)NULL;
    dt.out.pointmarkerlist            = (int *)NULL;
    dt.out.numberofpoints             = 0;
    dt.out.numberofpointattributes    = 0;
    //
    dt.out.trianglelist               = (int *)NULL;
    dt.out.triangleattributelist      = (double *)NULL;
    dt.out.trianglearealist           = (double *)NULL;
    dt.out.neighborlist               = (int *)NULL;
    dt.out.numberoftriangles          = 0;
    dt.out.numberofcorners            = 0;
    dt.out.numberoftriangleattributes = 0;
    //
    dt.out.segmentlist                = (int *)NULL;
    dt.out.segmentmarkerlist          = (int *)NULL;
    dt.out.numberofsegments           = 0;
    //
    dt.out.holelist                   = (double *)NULL;
    dt.out.numberofholes              = 0;
    //
    dt.out.regionlist                 = (double *)NULL;
    dt.out.numberofregions            = 0;
    //
    dt.out.edgelist                   = (int *)NULL;
    dt.out.edgemarkerlist             = (int *)NULL;
    dt.out.normlist                   = (double *)NULL;
    dt.out.numberofedges              = 0;
}

Friends And Related Function Documentation

friend class MemoryManager< TriangleInterface >

friend

Definition at line 58 of file TriangleInterface.h.

Member Data Documentation

DelaunayTriangle Nektar::NekMeshUtils::TriangleInterface::dt

private

Triangle data strucutres.

Definition at line 113 of file TriangleInterface.h.

std::vector<std::vector<NodeSharedPtr> > Nektar::NekMeshUtils::TriangleInterface::m_boundingloops

private

List of bounding nodes to the surface.

Definition at line 101 of file TriangleInterface.h.

Referenced by Assign().

std::vector<Array<OneD, NekDouble> > Nektar::NekMeshUtils::TriangleInterface::m_centers

private

Coordinates of the centers of the loops.

Definition at line 105 of file TriangleInterface.h.

Referenced by Assign().

std::vector<NodeSharedPtr> Nektar::NekMeshUtils::TriangleInterface::m_stienerpoints

private

List of additional nodes.

Definition at line 103 of file TriangleInterface.h.

Referenced by AssignStiener().

NekDouble Nektar::NekMeshUtils::TriangleInterface::m_str

private

Stretching factor of parameter plane.

Definition at line 111 of file TriangleInterface.h.

Referenced by Assign().

std::map<int, NodeSharedPtr> Nektar::NekMeshUtils::TriangleInterface::nodemap

private

Map from NekMesh id to triangle id.

Definition at line 107 of file TriangleInterface.h.

int Nektar::NekMeshUtils::TriangleInterface::sid

private

ID of the surface.

Definition at line 109 of file TriangleInterface.h.

Referenced by Assign().