[SciPy-dev] Instability in numpy.random.logseries

[Alan Jackson]

For parameter values close to one, the returned valued can be bad.

In [242]: np.random.logseries(.999999999,10)
Out[242]: 
array([        910, -2147483648,       13814,       17661,          14,
              3783,         335,   180167317,          38,   256949708])

For that matter, the function accepts an input of 1.0, but the output 
is not useful. It should probably throw an error.

In [225]: np.random.logseries(1.,10)
Out[225]: 
array([-2147483648, -2147483648, -2147483648, -2147483648, -2147483648,
       -2147483648, -2147483648, -2147483648, -2147483648, -2147483648])


-- 
-----------------------------------------------------------------------
| Alan K. Jackson            | To see a World in a Grain of Sand      |
| alan at ajackson.org          | And a Heaven in a Wild Flower,         |
| www.ajackson.org           | Hold Infinity in the palm of your hand |
| Houston, Texas             | And Eternity in an hour. - Blake       |
-----------------------------------------------------------------------