[Octave-bug-tracker] [bug #61863] interp2 gives wrong results for complex matrices

[Christof Kaufmann]

Follow-up Comment #11, bug #61863 (project octave):

Regarding the performance of `complex()`: In my patch it was `NA + 1i*NA`, but
I think it doesn't matter, since it is just a single scalar.

Yes, I agree, that it is reasonable to fix it in the C++ code. I just have not
set up the tools for compilation and have not looked at the C++ API. There is
even one more file, that has to be considered. Here is an overview:
* `interp1()` uses `pchip()` and pchip()` uses `__pchip_deriv__()`
* `interp2()` uses `__pchip_deriv__()` directly
* `interpn()` might use `__pchip_deriv__()` in future, but "pchip" is not yet
implemented.
* `interpn()` uses `__lin_interpn__()`, which also ignores the imaginary
part.
So libinterp/corefcn/__lin_interpn__.cc and
libinterp/corefcn/__pchip_deriv__.cc would have to be fixed.

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?61863>

_______________________________________________
  Message sent via Savannah
  https://savannah.gnu.org/