[Octave-bug-tracker] [bug #64361] sinc() return NaN on sinc(-1:0.20:1-0.20) expression

[Dmitri A. Sergatskov]

Follow-up Comment #6, bug #64361 (project octave):

I can reproduce this problem with octave compiled with AVX/AVX2:


octave:1> sinc(-1:0.2:1-0.2)
ans =

 Columns 1 through 6:

   3.8982e-17   2.3387e-01   5.0455e-01   7.5683e-01   9.3549e-01         
NaN

 Columns 7 through 10:

   9.3549e-01   7.5683e-01   5.0455e-01   2.3387e-01



octave:2> format bit
octave:3> x = -1:0.20:1-0.20;
octave:4> x(6)
ans = 0000000000000000000000000000000000000000000000000000000000000000
octave:5> x(6) == 0
ans = 00000001
octave:6> x(6) != 0
ans = 00000000
octave:7> xpi = pi * x(6)
xpi = 0000000000000000000000000000000000000000000000000000000000000000
octave:8> xs = sin (xpi)
xs = 0000000000000000000000000000000000000000000000000000000000000000
octave:9> y = xs / xpi
y = 1111111111111000000000000000000000000000000000000000000000000000


Markus patch helps.

Dmitri.
-- 





    _______________________________________________________

Reply to this item at:

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

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