Skip to Content.
Sympa Menu

cgal-discuss - [cgal-discuss] Floating Point Exception Handling

Subject: CGAL users discussion list

List archive

[cgal-discuss] Floating Point Exception Handling


Chronological Thread 
  • From: Panagiotis Foteinos <>
  • To:
  • Subject: [cgal-discuss] Floating Point Exception Handling
  • Date: Mon, 1 Mar 2010 14:31:59 -0500
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:date:message-id:subject:from:to:content-type; b=lsatJpOwaM38aN/kk+iQGMy1FPIeSuN/4MWE2ePgNzX7c1j+K+e4vpPYLpTykIKABU niqOnEfegB4K6zpW+A+uIhNz8vIDfXcOBna7fWPE4IVhlSdu3ozuicBUBa58YYlPEKuy C/f2oCmftkHDy7biapoxYYHtj1Wr6yYQzpHos=

Hello all.

I am using the 3D Delaunay Triangulation package with a custom refinement algorithm. In my program, there are almost flat elements which I want to identify as flat (by checking the volume for example), and leave them alone without performing any calculations on them any more. The problem is that the calculation of the circumcenter or the radius (and the volume as a consequence) of an almost flat element gives a "floating point" exception (in release mode) and terminations follows.

A workaround is to hard-code a couple of lines inside CGAL package's code that throw an exception if there is division with zero , and then catch it from my code. If the catch statement is activated that means that I am possibly dealing with an almost element and I should just overpass it.

Is there any better way to do this?


Regards,
Panagiotis



Archive powered by MHonArc 2.6.16.

Top of Page