Re: griddata segmentation fault

[Chenjie Gu]

Thanks, David.

If x,y,z are column vectors, it works on one of my intel machines. (which encountered an error without segmentation fault before.)

However, still I got a segmentation fault on one of my amd64 machine (running 32-bit linux).

Is there anything to do with the processor?

Chenjie

On Mon, Dec 29, 2008 at 2:55 PM, David Bateman <address@hidden> wrote:

Chenjie Gu wrote:

Hi there,

The error after I run the griddata is
====================================
panic: Segmentation fault -- stopping myself...
attempting to save variables to `octave-core'...
save to `octave-core' complete
Segmentation fault
====================================

To reproduce the error, run the following script:
====================================
x = rand(1,100);
y = rand(1,100);
z = rand(1,100);
[xx yy] = meshgrid(0.2:0.1:0.8,0.2:0.1:0.8);
zz = griddata(x,y,z,xx,yy);
surf(xx,yy,zz)
====================================

I am using octave 3.0.3, I have installed qhull package (2003.1), and the above script runs correctly under matlab.
(I have tried to change "tri = delaunay (x, y);" to "tri = delaunay (x, y, "QJ");" in "octave/3.0.3/m/geometry/griddata.m", as mentioned in some posts, but it did not work.)

Any ideas to fix this problem?

Chenjie

Works here without a seg-fault if x,y and z are column rather than row vectors. That is

x = rand(100,1);
y = rand(100,1);
z = rand(100,1);

[xx yy] = meshgrid(0.2:0.1:0.8,0.2:0.1:0.8);
zz = griddata(x,y,z,xx,yy);
surf(xx,yy,zz)

However, with your code even though there is an error there is no seg-fault..

D.

--
David Bateman                                address@hidden
35 rue Gambetta                              +33 1 46 04 02 18 (Home)
92100 Boulogne-Billancourt FRANCE            +33 6 72 01 06 33 (Mob)