|
From: | anonymous |
Subject: | [Octave-bug-tracker] [bug #60738] logm returning incorrect result with some real non-symmetric matrices |
Date: | Mon, 14 Jun 2021 21:58:47 -0400 (EDT) |
User-agent: | Mozilla/5.0 (Windows NT 6.3; Win64; x64; rv:89.0) Gecko/20100101 Firefox/89.0 |
Follow-up Comment #6, bug #60738 (project octave): Patch attached, incorporated tolerance check and 10x speedup for Hermitian matrices. For random matrices, I'm just using A=s*rand(n) and A=s*randn(n) or complex(rand(n),randn(n)) of different sizes (n) and scales (s). A*=A' for Hermitian matrices or A+=A' for symmetric not nessaraly positive definate along with A=[21 17 6;-5 -1 -6;4 4 16]; or A=[5 4 2 1;0 1 -1 -1;-1 -1 3 0;1 1 -1 2]; for repeated eigenvalues. I'm placing the matrix at the top of the file and checking with the max elementwise norm at the bottom max(abs(expm(logm(A))-A)(:)) (file #51565) _______________________________________________________ Additional Item Attachment: File name: bug60738logm.patch Size:3 KB <https://file.savannah.gnu.org/file/bug60738logm.patch?file_id=51565> _______________________________________________________ Reply to this item at: <https://savannah.gnu.org/bugs/?60738> _______________________________________________ Message sent via Savannah https://savannah.gnu.org/
[Prev in Thread] | Current Thread | [Next in Thread] |