bugGNU Octave - Bugs: bug #53897, poly of conjugate complex pairs...


bug #53897: poly of conjugate complex pairs should return a real polynomial

Submitted by:  Marco Caliari <caliari>
Submitted on:  Mon 14 May 2018 09:11:56 AM UTC  
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Inaccurate Result
Status:  None Assigned to:  None
Originator Name:  Open/Closed:  Open
Release:  dev Operating System:  GNU/Linux

Add a New Comment(Rich Markup)

You are not logged in

Please log in, so followups can be emailed to you.


Mon 14 May 2018 09:11:56 AM UTC, original submission:

If you try the following

you see that N has spurious imaginary parts (quite large, in absolute value). This should not happen, since the roots x are complex conjugate pairs. We should first detect this fact (I thought to imag (sum (sort (x))) == 0, it is probably possible to do without computations, but only with sorting and comparison). Then, either we use the formula

(y-x(1))*(y-conj(x(1)) = y^2-2*real(x(1))*y+abs(x(1))^2

for each complex conjugate pair or we simply remove the spurious imaginary parts after the computation of the coefficients.

Marco Caliari <caliari>
Project Member


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by caliari (Submitted the item)

    Do you think this task is very important?
    If so, you can add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.




    No changes have been made to this item

    Back to the top

    Powered by Savane 3.3