Speedily create a random array

Question

I have a for-loop:

for (i=0; i<N; i++) a[i] = rand()

But that's too slow, and BLAS/vDSP/vecLib/etc are at my disposal. I'm looking for a function to create an array of random numbers; it seems like it should be included.

Besides the name for this function, I would also appreciate a some good vDSP/vecLib/etc documentation. I've found some documentation, but it's not good -- just general descriptions of the functions.

Yes, that's okay. I'll take any random array -- I can easily transform that with BLAS/etc. — Scott, Nov 07 '13 at 17:41

score 1 · Accepted Answer · answered Nov 07 '13 at 18:43

1

dlarnv does what you describe. It populates a double array with random (uniform or normal) values, and it is included in vecLib.

answered Nov 07 '13 at 18:43

finnw

47,861
24
143
221

1

And there's no mentions of the functions included anywhere in Apple's docs. – Scott Nov 08 '13 at 00:44

score 0 · Answer 2 · answered Nov 07 '13 at 04:56

Fast and well-distributed random number generators are useful for many purposes, so I'm sure you can find some great ones out there. I would start by looking at the Mersenne Twister, and in particular implementations that are optimized for SIMD, like this.

The next step in speed would come from GPU or DSP based algorithms, but I'm not aware of any in the public domain at the moment.

Speedily create a random array

2 Answers2