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

[Markus Mützel]

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

Anyway, it might be a good idea to not compare floating point values exactly
here:
https://hg.savannah.gnu.org/hgweb/octave/file/4555f9918009/scripts/signal/sinc.m#l47
Maybe, something like this would be better?

 
   y = ones (size (x));
 
-  i = (x != 0);
+  i = (abs (x) > eps);
 
   if (any (i(:)))
     t = pi * x(i);




    _______________________________________________________

Reply to this item at:

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

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