[Octave-bug-tracker] [bug #46103] inversion of a singular diagonal matrix A produces inv(A)==A

[Rik]

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

I'll ask on the Maintainer's list for someone to run the Matlab code.  I knew
the mutation would be hard so I think we should avoid it in dDiagMatrix.cc. 
Maybe, as I said, we can figure it out in inv.cc and do the mutation there.

But, I don't think we should always promote to full.  Lots, even most, of the
diagonal matrices will not be singular and so we don't want to go from linear
memory consumption (N) to quadratic consumption (N^2) every time.

    _______________________________________________________

Reply to this item at:

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

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