[Octave-bug-tracker] [bug #62321] Wrong result for Airy funtion

[Arun Giridhar]

Follow-up Comment #7, bug #62321 (project octave):

I found some strange behavior. Consider the input matrix V from comment #2:

V =
   -0.4649 -       0i   -0.3135 -       0i   -0.2155 -       0i
   -4.5181 -  0.0036i   -9.2814 -  0.0080i  -19.2341 -  0.0171i

where the first row is real and the second row is complex. In that case,
calling "P = airy (0, V)" gives different results for the first row than going
element by element in a loop with "Q(i) = airy (0, V(i))". The second row
agrees between P and Q, only the first row is incorrect for P.

The bizarre part is this: trying to isolate the bug, I used only the first row
of V where the discrepancy was, and magically P and Q both agreed! The
presence of the completely unrelated second row of V causes the *first* row to
give wrong results.

When real and complex values are mixed in a matrix, then along the line
something breaks down in applying the wrong function (real for complex or vice
versa) and it returns wrong results.


    _______________________________________________________

Reply to this item at:

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

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