[Octave-bug-tracker] [bug #53897] poly of conjugate complex pairs should return a real polynomial

[Robert Jenssen]

Follow-up Comment #13, bug #53897 (project octave):

This patch fails for zp2tf from the signal package with values that are
"nearly" real. Copying from line 84 of
share/octave/5.1.0/m/polynomial/poly.m:

octave:14> [z1,p1,k1]=cheby2(5,40,0.1);
octave:15> p1
p1 =
   9.345539872652292e-01 - 1.826618879469199e-01i
   8.372613608092301e-01 - 1.304057652412007e-01i
   7.781404349401341e-01 - 1.537554717468438e-17i
   8.372613608092301e-01 + 1.304057652412007e-01i
   9.345539872652292e-01 + 1.826618879469199e-01i

octave:16> tmp = sort (p1(imag (p1) > 0)) == sort (conj (p1(imag (p1) < 0)))
error: mx_el_eq: nonconformant arguments (op1 is 2x1, op2 is 3x1)
octave:16> 


    _______________________________________________________

Reply to this item at:

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

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