[Octave-bug-tracker] [bug #48307] sinc loses precision for large arguments

[Colin Macdonald]

Follow-up Comment #23, bug #48307 (project octave):

Here's one concrete approach: we could find a library for approximating
spherical bessel functions jn(z).  And j0(z) is sinc(z)...

(We could relate jn(x) to besselj(1/2,x)*1/sqrt(x)*sqrt(pi/2) but I don't
think this would be much or any improvement over sin(x)/x.  Also in my very
limited tests, Octave's non-integer-order bessel functions seem to have poorer
relative error compared to integers.)

    _______________________________________________________

Reply to this item at:

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

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