CGAL users discussion list

[Camille Wormser <>]

> I had a code (see below) that
> does that and it worked nicely using a homogeneous kernel with Gmpz.  
> Now I have
> changed the kernel to be Cartesian<double> and the function does not  
> work
> anymore.

I did not read your code in detail, but you say

> I tried different Kernels types (including
> Exact_predicates_inexact_constructions_kernel) but the problem  
> remains.

Then it probably means that you are performing some constructions.

I see that you test first "f1->plane()==f2->plane()", and then "f2- 
>halfedge()->next()->opposite()->facet()==f1".

So the supporting planes are constructed before being compared. You  
can avoid this construction by first testing that the two faces are  
adjacent, and then checking whether the four vertices are coplanar or  
not.

This may solve your problem, if you use  
Exact_predicates_inexact_constructions_kernel again.
--
Camille