CGAL users discussion list

[Nicholas Mario Wardhana <>]

On 7 March 2012 00:33, Sebastien Loriot (GeometryFactory)
<>
 wrote:
> Thanks for the update.
>
> Can you try before the assertion
>
> CGAL_assertion(new_rd >= copy_rd);
>
> to change
> new_rd=Orthogonal_distance_instance.new_distance(copy_rd,old_off,new_off,new_cut_dim);
> into
> new_rd=old_off==new_off?copy_rd:Orthogonal_distance_instance.new_distance(copy_rd,old_off,new_off,new_cut_dim);
>
> I still suspect an optimization problem.
>
> Sebastien.
>
>

Hi Dr Loriot,

After reverting to the previous assertions and changing the lines as
you suggested, I still got the same error. For example, in one
occasion, I got the following values:

new_rd = 3240.3969698751616
copy_rd = 3240.3969698751621 ,

so I suppose it was still the same floating point problem. Anyway, I
am wondering if you mind telling me a bit about what that new line
does?

Thank you.

Best regards,
Nicholas