[Octave-bug-tracker] [bug #49506] interp2(..., 'linear') gives incorrect result

[Nir Krakauer]

Follow-up Comment #6, bug #49506 (project octave):

Senseless undocumented Matlab behavior.

As a work-around,  we can replace line 228 in scripts/general/interp2.m


ZI = a(idx) + b(idx).*Xsc + c(idx).*Ysc + d(idx).*Xsc.*Ysc;


with

fis = @(x, i) reshape(x(i), size(i));
ZI = fis(a,idx) + fis(b,idx).*Xsc + fis(c,idx).*Ysc + fis(d,idx).*Xsc.*Ysc;


Then we can check that

%!assert (interp2 ([2 3], [2 3 4], [1 2; 3 4; 5 6], [2 3], 3, "linear"), [3
4])


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?49506>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/