[Octave-bug-tracker] [bug #55862] error when using polynomial function poly.m

[anonymous]

URL:
  <https://savannah.gnu.org/bugs/?55862>

                 Summary: error when using polynomial function poly.m
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Fri 08 Mar 2019 08:13:22 AM UTC
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Unexpected Error
                  Status: None
             Assigned to: None
         Originator Name: Gerhard Doblinger
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 5.1.0
        Operating System: GNU/Linux

    _______________________________________________________

Details:

Checking of conjugate complex roots (in order to return real-valued results)
fails in certain cases due to a logical comparison with == instead of
isequal(). 
Reason: comparison with == requires (matrix) arguments of same dimension
Example:
[b,a] = cheby2(5,70,1/4)
works with version 4.x, fails with version 5.1.0

Solution:
replace 
tmp = sort (v(imag (v) > 0)) == sort (conj (v(imag (v) < 0)));
with
tmp = isequal(sort (v(imag (v) > 0)), sort (conj (v(imag (v) < 0))));
in function poly.m






    _______________________________________________________

Reply to this item at:

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

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