Some functions don't handle inifinities and nans

[GoogleCodeExporter]

Powers, in particular, lack correct handling of special numbers.

Original issue reported on code.google.com by fredrik....@gmail.com on 16 Feb 2008 at 12:20

[GoogleCodeExporter]

>>> for f in [lambda x: nthroot(x, 3), sqrt, cos, log, atan, cosh]:
...   print f(mpf('inf'))
...   print f(mpf('nan'))
...
+inf
nan
+inf
nan
nan
nan
+inf
nan
1.5707963267949
nan
+inf
nan
>>> mpf('nan')**2
mpf('nan')
>>> mpf('inf')**2
mpf('+inf')

Are there still some functions left?

Original comment by Vinzent.Steinberg@gmail.com on 9 Oct 2008 at 2:33

[GoogleCodeExporter]

Yes, plenty.

>>> acosh(nan)
mpc(real='0.0', imag='nan')

>>> hyp2f1(nan,nan,nan,nan)
(infinite recursion)

Original comment by fredrik....@gmail.com on 9 Oct 2008 at 3:26

[GoogleCodeExporter]

There are still some functions that will crash or return strange things at 
infs/nans, but the situation is much better, and there are lots of tests.

Original comment by fredrik....@gmail.com on 14 Jun 2010 at 7:51

• Changed state: Fixed