[Octave-bug-tracker] [bug #64801] besselj(1, z) is not pure imaginary for pure imaginary inputs

[Lapresté Jean-Thierry]

Follow-up Comment #4, bug #64801 (project octave):

[(Erreur - Introuvable)]
I don't know either. But you are right a floating point problem.

In fact I am writing the kyosu library based on eve which implements in scalar
and simd complex and caley-dickson (quaternions octonions x-nions) and
classical functions on them.

I used octave to verify some of my ouputs and I spotted this problem.

What I think is that a complex function is meant to return pure imaginary or
real values in quite obvious an general circumstances it is better that it
does. 

However the problem can be complicated to fix.

I just wanted to signal the issue. Sorry for the noise.



    _______________________________________________________

Reply to this item at:

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

_______________________________________________
Message posté via Savannah
https://savannah.gnu.org/