|
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] |