Sun 08 Aug 2010 03:50:36 PM UTC, original submission:
Hi,
I am using octave 3.2 on ubuntu 10.04. I believe there is a bug in the interface of qhull. I ran the following code given as an example of the help manual.
x = rand (1, 10);
y = rand (size (x));
T = delaunay (x, y);
However this gives me following error:
qhull internal error (qh_qhull): temporary sets not empty (2)
While executing: | qhull d Qt Qbb Qc
Options selected for Qhull 2009.1 2009/06/14:
delaunay Qtriangulate Qbbound-last Qcoplanar-keep _pre-merge
_zero-centrum Pgood Qinterior-keep _max-width 0.93 Error-roundoff 1.3e-15
_one-merge 9.1e-15 Visible-distance 2.6e-15 U-coplanar-distance 2.6e-15
Width-outside 5.2e-15 _wide-facet 1.6e-14
Last point added to hull was p9.
At error exit:
Delaunay triangulation by the convex hull of 10 points in 3-d:
Number of input sites: 10
Number of Delaunay regions: 0
Statistics for: | qhull d Qt Qbb Qc
Number of points processed: 10
Number of hyperplanes created: 27
Number of facets in hull: 16
Number of distance tests for qhull: 48
error: _delaunayn_: qhull failed.
error: called from:
error: /usr/share/octave/3.2.3/m/geometry/delaunayn.m at line 55, column 5
error: /usr/share/octave/3.2.3/m/geometry/delaunay.m at line 55, column 11
Not sure about the possible cause or if wrong version of qhull is used in octave.
Thanks.
Sandeep
|