LibMeshPoly2TriTest Class Reference

Inheritance diagram for LibMeshPoly2TriTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (LibMeshPoly2TriTest)
	CPPUNIT_TEST (testLibMeshPoly2Tri)
	CPPUNIT_TEST (testLibMeshPoly2TriHole)
	CPPUNIT_TEST (testLibMeshPoly2TriSlivers)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	testLibMeshPoly2Tri ()
void	testLibMeshPoly2TriHole ()
void	testLibMeshPoly2TriSlivers ()

Detailed Description

Definition at line 12 of file libmesh_poly2tri.C.

Member Function Documentation

CPPUNIT_TEST() [1/3]

LibMeshPoly2TriTest::CPPUNIT_TEST

(

testLibMeshPoly2Tri

)

CPPUNIT_TEST() [2/3]

LibMeshPoly2TriTest::CPPUNIT_TEST

(

testLibMeshPoly2TriHole

)

CPPUNIT_TEST() [3/3]

LibMeshPoly2TriTest::CPPUNIT_TEST

(

testLibMeshPoly2TriSlivers

)

CPPUNIT_TEST_SUITE_END()

LibMeshPoly2TriTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

LibMeshPoly2TriTest::LIBMESH_CPPUNIT_TEST_SUITE

(

LibMeshPoly2TriTest

)

setUp()

void LibMeshPoly2TriTest::setUp ( )

inline

Definition at line 27 of file libmesh_poly2tri.C.

  {}

tearDown()

void LibMeshPoly2TriTest::tearDown ( )

inline

Definition at line 30 of file libmesh_poly2tri.C.

  {}

testLibMeshPoly2Tri()

void LibMeshPoly2TriTest::testLibMeshPoly2Tri ( )

inline

Definition at line 34 of file libmesh_poly2tri.C.

References libMesh::Utility::iota().

  {
    LOG_UNIT_TEST;

    std::vector<p2t::Point> pentagon {{0,0},{1,0},{1,1},{.5,1.5},{0,1}};

    std::vector<p2t::Point *> api_shim(pentagon.size());
    std::iota(api_shim.begin(), api_shim.end(), pentagon.data());

    p2t::CDT cdt(api_shim);

    cdt.Triangulate();

    auto tris = cdt.GetTriangles();

    // We mostly wanted to make sure this compiled, but might as well
    // make sure it gave us the expected triangle count while we're at
    // it.
    CPPUNIT_ASSERT_EQUAL(tris.size(), std::size_t(3));
  }

testLibMeshPoly2TriHole()

void LibMeshPoly2TriTest::testLibMeshPoly2TriHole ( )

inline

Definition at line 55 of file libmesh_poly2tri.C.

References libMesh::Utility::iota().

  {
//    const double eps = 0;  // This (or negative eps) succeeds
//    const double eps = 1e-12;  // This (or larger eps) succeeds
    const double eps = 1e-15; // This gave EdgeEvent - null triangle with older Poly2Tri

    std::vector<p2t::Point> outer_boundary
      {{0,0},{0.5,eps},{1,0},{1-eps,0.836541},
       {1,2},{.46,1.46+eps},{0,1},{eps,0.5}};

    std::vector<p2t::Point *> api_shim(outer_boundary.size());
    std::iota(api_shim.begin(), api_shim.end(), outer_boundary.data());

    const double r2o4 = std::sqrt(2.)/4;
    std::vector<p2t::Point> hole_boundary
      {{0.5+r2o4,0.5},{0.5,0.5+r2o4},{0.5-r2o4,0.5},{0.5-eps,0.5-r2o4}};

    std::vector<p2t::Point *> hole_shim(hole_boundary.size());
    std::iota(hole_shim.begin(), hole_shim.end(), hole_boundary.data());

    p2t::CDT cdt(api_shim);
    cdt.AddHole(hole_shim);

    std::vector<p2t::Point> interior_points
      {{0.21,0.79},{0.21,0.21},{0.79,0.21}};
    for (auto & p : interior_points)
      cdt.AddPoint(&p);

    cdt.Triangulate();

    auto tris = cdt.GetTriangles();

    std::cout << "With perturbation " << eps << " we got " << tris.size() << " triangles!" << std::endl;

    // We mostly wanted to make sure this compiled, but might as well
    // make sure it gave us the expected triangle count while we're at
    // it.
    CPPUNIT_ASSERT_EQUAL(tris.size(), std::size_t(18));
  }

testLibMeshPoly2TriSlivers()

void LibMeshPoly2TriTest::testLibMeshPoly2TriSlivers ( )

inline

Definition at line 95 of file libmesh_poly2tri.C.

References libMesh::Utility::iota().

  {
    // Points from a libMesh Poly2TriTriangulator failure case
    std::vector<p2t::Point> outer_boundary
    {{0, 0}, {1, 0}, {1, 2}, {0, 1}};

    std::vector<p2t::Point *> api_shim(outer_boundary.size());
    std::iota(api_shim.begin(), api_shim.end(), outer_boundary.data());

    const double r2o4 = std::sqrt(2.)/4;
    std::vector<p2t::Point> hole_boundary
      {{0.5,0.5+r2o4},{0.5-r2o4,0.5},{0.5,0.5}};

    std::vector<p2t::Point *> hole_shim(hole_boundary.size());
    std::iota(hole_shim.begin(), hole_shim.end(), hole_boundary.data());

    p2t::CDT cdt(api_shim);
    cdt.AddHole(hole_shim);

    // const double eps = 1e-12; // fails
    const double eps = 1e-11; // succeeds

    std::vector<p2t::Point> interior_points
    {{0.5-(r2o4/2)-eps, 0.5+(r2o4/2)+eps}};

    for (auto & p : interior_points)
      cdt.AddPoint(&p);

    cdt.Triangulate();
  }

---

The documentation for this class was generated from the following file:

• tests/mesh/libmesh_poly2tri.C