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

[Marco Caliari]

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

Hi,

I teach my students that the inverse of a diagonal matrix is simply a diagonal
matrix with the inverse diagonal. So, I prefer solution a) in #6, I consider
it the most natural and the least surprising. I agree with #5 that Octave
should give the same result as with a regular matrix, but the challenge should
be have Octave to recognize a regular diagonal matrix A as a diagonal matrix
and compute the inverse by diag(1./diag(A)). My 2 cents...

Marco

    _______________________________________________________

Reply to this item at:

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

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