CGAL users discussion list

[Panagiotis Foteinos <>]

Thank you Sebastien.
I assumed that exact constructions were not necessary for the 4D kernels, simply because this is the case for the 3D kernels. I was wrong. The Cartesian_d<CGAL::Gmpq> kernel works fine now.

I encountered, however, another error which sadly does not trigger any assertion, but rather it throws a segmentation fault.

Consider the following insertion sequence of 4D points:
    1st point:   6.8,    6.8,    6.8,    6.8
    2nd point:  6.8,    6.8,    6.8,   43.2
    3rd point:   6.8,    6.8,   43.2,   6.8
    4th point:   6.8,    6.8,   43.2,  43.2

The function DT.lookup() on the 4th point before its insertion throws a segmentation fault in the file CGAL/Handle_for.h (Line 125) .

I looked into the code and I observed that CGAL checks the points of at least one simplex. However, no simplex has to exist in the triangulation so far: only 3 points have been inserted. Also, no extra (i.e., more than it is needed to describe the actual triangulation) simplices are created (for example infinite simplices as is the case in 3D), and therefore I do not understand why CGAL finds a simplex to check. Maybe this observation is related to the segmentation fault occurred.

Any input? I am really wondering if someone could reproduce the error I get, i.e. a segmentation fault during the insertion of the 4th point... I am using version 3.6.1 .

Thank you in advance,
Panagiotis Foteinos

On Tue, Aug 24, 2010 at 8:39 AM, Sebastien Loriot (GeometryFactory) <sloriot.ml@gmail.com> wrote:

Panagiotis Foteinos wrote:

Hello.

I am using the Delaunay_d class for the triangulation of 4D points. I am also using the Cartesian_d<double> kernel.

Looking at the documentation here:

http://www.cgal.org/Manual/latest/doc_html/cgal_manual/Kernel_d_ref/Class_Cartesian_d.html#Index_anchor_736

you'll see that you need to use an exact number type such as CGAL::Gmpq (in particular with collinear points) to compute safely the triangulation.

S.

I observed that the insertion of the following simple point sequence:

   1st point:  6.8,    6.8,    6.8,    6.8
   2nd point: 6.8,    6.8,    6.8,   43.2
   3rd point:  6.8,    6.8,   43.2,   6.8
   4th point:  6.8,    6.8,   43.2,  43.2

triggers, during the insertion of the 4th point,  the following assertion:

   terminate called after throwing an instance of 'CGAL::Assertion_exception'
   what():  CGAL ERROR: assertion violation!
   Expr: j != dim
   File: .../CGAL-3.6.1/include/CGAL/Kernel_d/HyperplaneCd.h
   Line: 103
   Explanation: HyperplaneCd::constructor: cannot use o to determine side.



Is there a workaround? Why is this behaviour happening?


Thank you,
Panagiotis Foteinos

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