[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: trivial delaunay example fails
From: |
Ben Abbott |
Subject: |
Re: trivial delaunay example fails |
Date: |
Thu, 20 Oct 2011 16:09:26 -0400 |
On Oct 20, 2011, at 11:16 AM, Ben Abbott wrote:
> I've filed a bug report.
>
> https://savannah.gnu.org/bugs/index.php?34604
>
> I'm not sure if this is a problem local to my OS, my package manager, or
> something I am doing. Can someone try the example below?
>
> x = [0 0 1 1];
> y = [0 1 1 0];
> tri = delaunay (x, y)
>
> The answer should be ...
>
> tri =
>
> 3 2 4
> 2 1 4
>
> However, I'm getting a qhull error (see the bug report for more detail).
>
> Ben
Adding the "QJ" option resolves the error.
tri = delaunay (x, y, {"QJ"})
qhull precision warning:
The initial hull is narrow (cosine of min. angle is 1.0000000000000000).
A coplanar point may lead to a wide facet. Options 'QbB' (scale to unit box)
or 'Qbb' (scale last coordinate) may remove this warning. Use 'Pp' to skip
this warning. See 'Limitations' in qh-impre.htm.
tri =
2 4 1
2 4 3
Does anyone know if there is a problem with turning these options on by default?
Ben