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

[Rik]

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

The diagonal matrix class is supposed to be a memory saving optimization.  All
calculations on diagonal matrices should result in exactly the same behavior
as when a full matrix is used.  My suggestion would be to change the diagonal
class to return an Inf array for a singular matrix the same as is done fore a
full matrix.

    _______________________________________________________

Reply to this item at:

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

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