CGAL Using Locate() to Find Cell on Triangulation Surface -

CGAL Using Locate() to Find Cell on Triangulation Surface -

using cgal, have 3d delaunay triangulation of set of random points on unit sphere, obtained via:

delaunay t(points.begin(), points.end());

now, able query t (using locate() or that) find 3 vertices of surface facet arbitrary point (also on unit sphere) contained inside.

when utilize locate(), interior cells results sometimes, include infinite vertex. don't want of these. want surface facets , able arbitrary point seek find on unit sphere. trying figure out has taken lot longer thought would.

any help much obliged. thanks.

so utilize find_conflit(), cgal::emptyset_iterator cit because don't need these. in bfit, facets of boundary of "hole", , hole tetrahedra in conflict point (whose circumscribing sphere contains point, natural extension infinite vertex). so, bfit, set them in standard container using std::back_inserter example. then, iterate on these, tests if finite facets. finite facets separate point rest of triangulation, so, can orientation() tests center of sphere 1 interested in.

triangulation cgal delaunay