Re: griddata

[Paul Kienzle]

There is a bug in octave introduced sometime between octave
2.1.40 and 2.1.48 where the wrong thing happens if the shape
of a boolean indexing matrix does not match the shape of the
matrix being indexed.

E.g.,

octave-2.1.40> x=[1,2;3,4]; x(:)'(x!=4)
ans =

  1  3  2

octave-2.1.48> x=[1,2;3,4]; x(:)'(x!=4)
ans =

  1  2
  3  5

It is still present in octave-2.1.50.

The following patch to griddata works around this problem.

I will post a new version shortly.

Paul Kienzle
address@hidden

diff -c -p -r1.5 griddata.m
*** griddata.m  27 Dec 2002 19:19:04 -0000      1.5
--- griddata.m  10 Nov 2003 16:46:32 -0000
*************** function [rx, ry, rz] = griddata (x,y,z,
*** 67,73 ****
  
      ## keep non zero values before overwriting zeros with 1
      valid=!isnan(reshape(tri_list,size(xi)));
!     tri_list = tri_list(valid);
      nr_t=rows(tri_list);
  
      ## assign x,y,z for each point of triangle
--- 67,73 ----
  
      ## keep non zero values before overwriting zeros with 1
      valid=!isnan(reshape(tri_list,size(xi)));
!     tri_list = tri_list(!isnan(tri_list));
      nr_t=rows(tri_list);
  
      ## assign x,y,z for each point of triangle

On Mon, Nov 10, 2003 at 12:36:17PM +0100, Claudio Belotti wrote:
> Hi all,
>    has anyone succeeded in using griddata from octave.sourceforge.net?
> 
> Claudio
> 
> 
> 
> -------------------------------------------------------------
> 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
> -------------------------------------------------------------
> 
> 
> 



-------------------------------------------------------------
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
-------------------------------------------------------------