[Numpy-discussion] in place random generation

Robert Kern robert.kern at gmail.com
Thu Mar 8 12:48:06 EST 2007


Daniel Mahler wrote:
> On 3/8/07, Charles R Harris <charlesr.harris at gmail.com> wrote:

>> Robert thought this might relate to Travis' changes adding broadcasting to
>> the random number generator. It does seem certain that generating small
>> arrays of random numbers has a very high overhead.
> 
> Does that mean someone is working on fixing this?

It's not on the top of my list, no.

> Also what does 'adding broadcasting to the number generator' mean?

normal([[0.0], [0.5]], [1.0, 2.0, 3.0])

That gives you a (2, 3) array of random numbers drawn from 6 different normal
distributions: [[(mean=0, stdev=1), (mean=0, stdev=2), (mean=0, stdev=3)],
                [(mean=0.5, stdev=1), (mean=0.5, stdev=2), (mean=0.5, stdev=3)]]

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
 that is made terrible by our own mad attempt to interpret it as though it had
 an underlying truth."
  -- Umberto Eco



More information about the NumPy-Discussion mailing list