CGAL users discussion list

[pam <>]

Thanks a lot. It works.


2008/10/8 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();
>
>
>
>
>
>
>
> --
> You are currently subscribed to cgal-discuss.
> To unsubscribe or access the archives, go to
> https://lists-sop.inria.fr/wws/info/cgal-discuss
>