Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Squared Radius

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Squared Radius


Chronological Thread 
  • From: Panagiotis Foteinos <>
  • To:
  • Subject: Re: [cgal-discuss] Squared Radius
  • Date: Fri, 28 Aug 2009 17:55:38 -0400
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; b=Puol38FGUPyqZV5+pNcvhlTdsFVtYvwqaqjwhJ0rAQ/RTAKvdhLRQCDwnaIIEF5Y4Z 2nuZpziO9/pwysVdekEqwCVAu1koJEP3amzGJ0fi1FhfjGvj3OvM3IdhrsyHMB/qOxAl p9zj0taLFnVS5Ql5wg329+llV/Xhp/wcmWxxw=

Yes, it is an exact predicate, inexact construction kernel.

This is exactly what I am thinking... I am running out of precision. The exact construction kernel, however, is not an option for my code. :(

Thank you,
Panagiotis

On Fri, Aug 28, 2009 at 5:49 PM, Ben Supnik <> wrote:
The experts will correct me I am sure, but ...

This is an exact predicate, inexact construction kernel.

If I understand correctly, this kernel can tell you IF the triangle is degenerate, but if it is not degenerate it might still NOT be able to calculate the area properly, since the underlying numeric type (double) can round out.

Try using CGAL::Quotient<CGAL::MP_Float> instead of double...if things get better maybe it's a construction issue?

cheers
ben

Panagiotis Foteinos wrote:
I am using the  CGAL::Filtered_kernel<CGAL::Cartesian<double> > kernel.

On Fri, Aug 28, 2009 at 5:27 PM, Ben Supnik < <mailto:>> wrote:

   What Kernel are you using?


   Panagiotis Foteinos wrote:

       Hello everybody.

       I observed something strange.

       If a triangle t(p, q, r) is close to a degenerate one with the
       value, however, of the predicate t.is_degenerate() always being
       FALSE, then during the calculation of CGAL::squared_radius(p, q,
       r), the assertion:

       ********************************************************************************************
         terminate called after throwing an instance of
       'CGAL::Assertion_exception'
        what():  CGAL ERROR: assertion violation!
       Expr: ! CGAL_NTS is_zero(den)
       File:
       /scratch/local/CGAL-3.4/include/CGAL/constructions/kernel_ftC2.h
       Line: 60
       ********************************************************************************************

       is trigerred, and the calculation fails.

       Unfortunately, this error is difficult to be reproduced, but
       sometimes it happens.


       Is this behavior expected? Is it relevant to the fact that the
       triangle is almost degenerate?


       Thank you,
       Panagiotis


   --    Scenery Home Page: http://scenery.x-plane.com/
   Scenery blog: http://xplanescenery.blogspot.com/
   Plugin SDK: http://www.xsquawkbox.net/xpsdk/
   X-Plane Wiki: http://wiki.x-plane.com/
   Scenery mailing list:
   <mailto:>

   Developer mailing list:
   <mailto:>

   --    You are currently subscribed to cgal-discuss.
   To unsubscribe or access the archives, go to
   https://lists-sop.inria.fr/wws/info/cgal-discuss



--
Scenery Home Page: http://scenery.x-plane.com/
Scenery blog: http://xplanescenery.blogspot.com/
Plugin SDK: http://www.xsquawkbox.net/xpsdk/
X-Plane Wiki: http://wiki.x-plane.com/
Scenery mailing list:
Developer mailing list:
--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://lists-sop.inria.fr/wws/info/cgal-discuss




Archive powered by MHonArc 2.6.16.

Top of Page