[Octave-bug-tracker] [bug #49900] cplxpair return error message where data seems good

[Henk Borsje]

Follow-up Comment #2, bug #49900 (project octave):

It looks like there still is a problem with cplxpair().  It seems unrelated to
units of tolerance, as someone suggested, it is a different bug.

When the function is entered with:

z = -1.0000 + 1.2246e-16

tol = 2.2204e-10

the same error still appears in the following section of code:

  ## For each remaining z, place the value and its conjugate at the start of
  ## the returned list, and remove them from further consideration.
  
  for j = 1:m
    p = n - nr(j);
    for i = 1:2:p
      if (i+1 > p)
        error ("cplxpair: could not pair all complex numbers");
      endif



'nr' seems to be an empty sparse matrix and 'p' gets the value of 1 and always
has a value less than 'i+1'

The input arguments come from a call to 'butter' and are not made up.

Thank you for your attention.

Henk Borsje


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?49900>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/