'delauney3' fails

[Jaroslav Hajek]

On Wed, May 27, 2009 at 9:00 AM, Søren Hauberg <soren at hauberg.org> wrote:
> Hi All
>
> I've tried the following example from the Mathworks webpage [1]
>
>        d = [-1 1];
>        [x,y,z] = meshgrid(d,d,d);  % A cube
>        x = [x(:);0];
>        y = [y(:);0];
>        z = [z(:);0];
>        % [x,y,z] are corners of a cube plus the center.
>        Tes = delaunay3(x,y,z)
>

Works for me, apparently:
Tes =

   4   9   3   1
   4   9   2   1
   7   9   3   1
   7   9   5   1
   7   4   9   3
   7   4   8   9
   6   9   2   1
   6   9   5   1
   6   4   9   2
   6   4   8   9
   6   7   9   5
   6   7   8   9

 IIRC, miscompiled qhull is not uncommon (it's sensitive to flags), so
maybe that's the source of your problem.

cheers

-- 
RNDr. Jaroslav Hajek
computing expert & GNU Octave developer
Aeronautical Research and Test Institute (VZLU)
Prague, Czech Republic
url: www.highegg.matfyz.cz