facet - CGAL: Access to Results of find_conflicts() -

facet - CGAL: Access to Results of find_conflicts() -

so confused find_conflicts function in cgal. thought knew std::pair, , thought knew going on in find_conflicts(), life of me, not sure how access results. thought iterator passed find_conflicts plenty access values directly. (i.e., want @ facets set in "vector facets,") , appears if i'm doing because can

typedef std::pair<std::vector<facet>, std::vector<cell> > facetandcell; * * * cell_handle cell = t.locate(curr_point); std::vector<facet> facets; t.find_conflicts(curr_point, cell, std::back_inserter(facets), cgal::emptyset_iterator()); cgal::first_of_pair_property_map<facetandcell> my_p_map(); delaunay::finite_facets_iterator ff_iter; std::vector<facet>::iterator facet_iter; // here, i'm trying accomplish figuring out of facets // in conflict finite. had wanted utilize kind of test // "is_infinite()" on facet @ hand, isn't available // facet out of context of triangulation it's part of. for(facet_iter = facets.begin(); facet_iter != facets.end(); facet_iter++){ for(ff_iter = t.finite_facets_begin(); ff_iter != t.finite_facets_end(); ff_iter++){ // since error facet_iter of type pair, thought work, doesn't. // error! cout << facet_iter->first << endl; // works, led me believe comparing 1 facet facet. /* if(*facet_iter == *ff_iter){ cout << "finite facet!" << endl; break; }*/ } }

in summary:

1) overall, want know facets result of find_conflicts() finite. if there easier way this, sense free allow me know.

2) otherwise, more specific problem here need @ vector of facets results find_conflicts() , @ vertices of each facet. supposed working returned "pair" of cells , facets or can access vector directly, i'm trying do?

help, please, thanks.

for finiteness testing, utilize is_infinite(). see http://doc.cgal.org/latest/triangulation_3/classcgal_1_1triangulation__3.html#af024721d3ae4b9327ffe442fb828935c

otherwise, maybe confused because facet type typedef pair (unfortunately).

ps: alternatively, can utilize marc glisse's suggestion in other question (using locate(midpoint(p, sphere center))). might easier. cgal using locate() find cell on triangulation surface

facet triangulation cgal vertices