[Numpy-discussion] numarray: how to make sinc(x) function - "divide by zero" !

[Sebastian H. <ha...@ms...>]

Hi,
The best I could come up with was:

def sinc(r):
    na.Error.pushMode(all="ignore")
    a = na.where(r, na.divide(na.sin(r),r), 1)
    na.Error.popMode()
    return a

but I still seem to get  a warning...
>>> F.sinc(0)
Warning: Encountered invalid numeric result(s)  in divide
1.0

Thanks,
Sebastian Haase