CGAL users discussion list

[Olivier Devillers <>]

pam a écrit :
> Hi all!
> I'm working with delaunay 3D and I was able to view cells ans faces
> but when I try to go trougth the edge I get some rare points that are
> not part of my points.
>
> Here is my code. Any idea of what is wrong?
>
> Delaunay::Finite_edges_iterator edge;
>   Point pt0;
>   Point pt1;
>   for(edge=T.finite_edges_begin(); edge != T.finite_edges_end();edge++)
>   {
>      std::cout<<" Edges "<<j<<std::endl;
>      pt0=edge->first->vertex(Delaunay::cw(edge->second))->point();
>      pt1=edge->first->vertex(Delaunay::ccw(edge->second))->point();
>          std::cout<<pt0<<std::endl;
>          std::cout<<pt1<<std::endl;
>           j++;  
>         
>   }
>
> Thanks in advance
>
> Pamela
>   
from
http://www.cgal.org/Manual/3.3/doc_html/cgal_manual/TriangulationDS_3_ref/Concept_TriangulationDataStructure_3.html#Cross_link_anchor_986
typedef Triple<Cell_handle, int, int>
   Edge;     (c,i,j) is the edge of cell c whose vertices indices are i 
and j. (See Section 28.1.)

thus it seems that you should use

pt0=edge->first->vertex(edge->second)->point();
pt1=edge->first->vertex(edge->third)->point();