[Octave-bug-tracker] [bug #52285] inv(0) gives 0 (instead of inf)

[Rik]

Follow-up Comment #4, bug #52285 (project octave):

I stepped through the code in dMatrix.cc and in fact the MatrixType is set to
Diagonal for the case of a scalar.  It looks like the general question is why
there isn't an early decision to take care of diagonal matrices which would be
very easy


  if (typ == MatrixType::Diagonal)
    ret = 1 / (*this);  
  else if (...)


Find attached a quick test (inv.diff).  It seems to pass "make check" with no
problems.  To be complete there should also be some BIST tests for the correct
behavior of scalars of all types (double, single, double complex, single
complex) with respect to the various operations '1 / x', 'inv (x)', x^-1,
'pinv (x)'. 


(file #42320)
    _______________________________________________________

Additional Item Attachment:

File name: inv.diff                       Size:4 KB


    _______________________________________________________

Reply to this item at:

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

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