[Octave-bug-tracker] [bug #60797] sqrtm: returns nan for matrix of ones with rows and columns >=4

[Markus Mützel]

Update of bug #60797 (project octave):

                  Status:          Ready For Test => Patch Submitted        
                 Release:                   6.2.0 => dev                    

    _______________________________________________________

Follow-up Comment #3:

The attached patch adds the shortcut for Hermitian matrices along the lines
described in comment #2.
For performance reasons, I opted to take the square root of the diagonal
elements in-place and don't touch the off-diagonal elements of the Schur
matrix. Those off-diagonal elements are checked to be below the tolerance. But
they might be non-zero.

Seems to be working correctly for me.


(file #51725)
    _______________________________________________________

Additional Item Attachment:

File name: bug60797-sqrtm-hermitian.patch Size:4 KB
   
<https://file.savannah.gnu.org/file/bug60797-sqrtm-hermitian.patch?file_id=51725>



    _______________________________________________________

Reply to this item at:

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

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