Andrea Riciputi wrote:
> Why not arctan2(ff.imag, ff.real)?
Yes, for the love of all that is holy, please use arctan2()!
For one thing, arctan2() solves other problems besides the divide-by-zero.
In [1]: from numarray import *
In [2]: arctan(-1/-1) == arctan(1/1)
Out[2]: True
In [3]: arctan2(-1,-1) == arctan2(1,1)
Out[3]: False
--
Robert Kern
rkern at ucsd.edu
"In the fields of hell where the grass grows high
Are the graves of dreams allowed to die."
-- Richard Harter