|
| From: | Marco Caliari |
| Subject: | [Octave-bug-tracker] [bug #53897] poly of conjugate complex pairs should return a real polynomial |
| Date: | Wed, 17 Oct 2018 14:37:57 -0400 (EDT) |
| User-agent: | Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:62.0) Gecko/20100101 Firefox/62.0 |
Follow-up Comment #6, bug #53897 (project octave):
As far as I know, eig (that is, the undelying LAPACK routines) gives exact
complex conjugates for real matrices. This is confirmed by experiments like
octave:1> A = randn (100);
octave:2> x = eig (A);
octave:3> all(sort(x(imag(x)>0))==sort(conj(x(imag(x)<0))))
ans = 1
Therefore, I would not introduce the tolerance.
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?53897>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
| [Prev in Thread] | Current Thread | [Next in Thread] |