Re: INTERP2 / GRIDDATA : problems using those functions

[Pierre Baldensperger]

Thank you Jonathan,

I tried your fix and Octave no longer crashes, but the problems just moved
a few lines further : now the "cross" function complains that its argument
vectors should have 3 components, even though it is obviously the case, 
since
it is programmed that way... I don't get it : my Octave installation must be
totally corrupt or something...

I haven't been able to solve my "interp2" problem either.

Thanks again anyway,
-Pierre.



Jonathan Stickel wrote:

> More than a year ago I had some trouble with the griddata function, 
> program crashes like yours if I recall.  It turned out to be the call 
> for delaunay triangulation.  Adding the "QJ" option fixed the problem 
> (diff attached).  I tried to get this solution checked into CVS, but I 
> met a lot of resistance and I don't know why.
>
> Jonathan
>
> --- griddata_orig.m     2004-09-01 12:04:42.000000000 -0700
> +++ griddata.m  2004-09-01 12:06:10.228458056 -0700
> @@ -49,7 +49,7 @@
>   [nr,nc]=size(xi);
>   
>   ## triangulate data
> -  tri=delaunay(x,y);
> +  tri=delaunay(x,y,"QJ");
>   zi=zeros(size(xi));
>   
>   if strcmp(method,'cubic')




-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------