bug in __lin_interpn__.cc for decreasing coordinates

[David Bateman]

Alexander Barth wrote:
> Hi all again,
> this code,
>
> [x,y,z] = ndgrid(0:2);
> f = x.^2+y.^2+z.^2;
> fi = interpn(x,y,-z,f,1.5,1.5,-1.5)
>
> returns NA (octave 3.0.0) while the correct answer is 7.5. The
> attached patch, handles the case for decreasing coordinate values in
> __lin_interpn__.cc
>
> Proposed changelog:
> __lin_interpn__.cc: handle decreasing coordinate values
>
> Thanks,
> Alex
>
Its a bit hard to read your patch due to changes in the tabulation.
Though this is your code, and it addresses a real issue, so I'm all for
applying it..

D.

-- 
David Bateman                                David.Bateman at motorola.com
Motorola Labs - Paris                        +33 1 69 35 48 04 (Ph) 
Parc Les Algorithmes, Commune de St Aubin    +33 6 72 01 06 33 (Mob) 
91193 Gif-Sur-Yvette FRANCE                  +33 1 69 35 77 01 (Fax) 

The information contained in this communication has been classified as: 

[x] General Business Information 
[ ] Motorola Internal Use Only 
[ ] Motorola Confidential Proprietary