[Octave-bug-tracker] [bug #63384] Error for "inv(NaN)" with LAPACK 3.11.0

[Rik]

Follow-up Comment #35, bug #63384 (project octave):

Should we bother with Matlab compatibility here?  Or maybe this is a different
issue.  See test code below:


octave:1> x = nan (2,2)
x =

   NaN   NaN
   NaN   NaN

octave:2> [y,rc] = inv (x)
y =

   NaN   NaN
   NaN   NaN

rc = NaN
octave:3> [y,rc] = inv (single (x))
y =

   NaN     0
   NaN   NaN

rc = NaN
octave:4> [y,rc] = inv (complex (single (x), ones (2,2)))
y =

   Inf   Inf
   Inf   Inf

rc = 0
octave:5> [y,rc] = inv (complex (x, ones (2,2)))
y =

   Inf   Inf
   Inf   Inf

rc = 0


In Matlab, the return value is always all NaN for every entry.  Complex
matrices are returned with "NaN + NaNi" for each element.


    _______________________________________________________

Reply to this item at:

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

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