[Octave-bug-tracker] [bug #51246] pinv(0) different from Matlab

[Rik]

Follow-up Comment #8, bug #51246 (project octave):

Any fix also needs to respect bug #49690 where the inverse of Inf should be 0
for Matlab compatibility.  My guess is that is where the problem in inverse
for the input 0 has crept in.  See cset 22154.


changeset:   22154:1d242ae72240
user:        Rik <address@hidden>
date:        Tue Nov 29 16:46:35 2016 -0800
files:       libinterp/corefcn/inv.cc liboctave/array/CMatrix.cc
liboctave/array/dMatrix.cc liboctave/array/fCMatrix.cc
liboctave/array/fMatrix.cc
description:
Return 0 for special case of scalar Inf input to inverse (bug #49690).

* inv.cc (Finv): Add BIST tests for new behavior.

* CMatrix.cc, fCMatrix.cc, dMatrix.cc, fMatrix.cc (::inverse):
If function has had to calculate the reciprocal condition number, and
the rcond value is zero, then return 0 if the input was a scalar, or
a matrix of all Inf values.




    _______________________________________________________

Reply to this item at:

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

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