[Numpy-discussion] multivariate_normal issue with 'size' argument

[Peter Cock]

On Fri, May 24, 2013 at 2:47 PM, Warren Weckesser
<warren.weckesser at gmail.com> wrote:
>
>Peter wrote:
>> -------------------------------------------
>> Successes
>> -------------------------------------------
>>
>> 64 bit Linux:
>>
>> $ python2.6
>> Python 2.6.6 (r266:84292, Sep 11 2012, 08:34:23)
>> [GCC 4.4.6 20120305 (Red Hat 4.4.6-4)] on linux2
>> Type "help", "copyright", "credits" or "license" for more information.
>>>>> import sys;print("%x" % sys.maxsize, sys.maxsize > 2**32)
>> ('7fffffffffffffff', True)
>>>>> import numpy as np
>>>>> print(np.random.multivariate_normal(mean=np.zeros(2), cov=np.eye(2),
>>>>> size=1))
>> [[-0.27469218 -2.12911784]]
>>>>> print(np.random.multivariate_normal(mean=np.zeros(2), cov=np.eye(2),
>>>>> size=np.int64(1)))
>> [[ 0.02609307  0.32485211]]
>>>>> np.__version__
>> '1.6.2'
>>>>> quit()
>>
>
> Peter: wow, that's a lot of tests!

I try to keep a broad range on hand for testing my own code.

> Two more data points:
> On Ubuntu 12.04, using 64 bit builds of Python 2.7.4 (from Anaconda
> 1.5.0), and numpy built from source: numpy 1.6.1 gives the error, but
> 1.6.2 does not.
>
> Warren

That's interesting - and matches my only success being with NumPy 1.6.2

This suggests this was broken to up 1.6.1, but fixed in the 1.6.2 branch
and not the 1.7 branch. Have anyone tried the current master branch?

Peter