CGAL users discussion list

[Adam Getchell <>]

I'm looking to use the flip function in a 3D Delaunay triangulation in order to convert 2 simplices into 3, starting from a seed with 2 cells:

S3->insert(typename T::Point(0,0,1));
  S3->insert(typename T::Point(1, 0, 1));
  S3->insert(typename T::Point(0, 1, 1));
  S3->insert(typename T::Point(0, 0, 2));
  /// This point gives us two cells, 5 vertices, 9 edges, and 7 faces
  /// These two cells are the basis for the 2-6 3D ergodic move
  S3->insert(typename T::Point(0, 0, 0));

It looks like I need to obtain an edge to do so. I have this code so far:

/// Get a cell
  Delaunay::Finite_cells_iterator cit;
  cit = S3->finite_cells_begin();
  assert(S3->is_cell(cit));
  /// Get an edge, hopefully in the same cell
  Delaunay::Finite_facets_iterator fit;
  fit = S3->finite_facets_begin();
  /// Now flip them
  //Delaunay::flip(fit);

What I'd like to do is get the internal facet of the two cells and then do a

Delaunay::flip(fit);

After determining that fit is indeed the internal facet to the cells initialized above.

Pointers?

I'll keep reading the manuals.

The source code is here:

https://github.com/acgetchell/CDT-plusplus/blob/master/spherical_3_complex.h

--

Adam Getchell

about.me/adamgetchell

"Invincibility is in oneself, vulnerability in the opponent." -- Sun Tzu