pseudo-random number generator rand() in .mex
2 views (last 30 days)
Show older comments
In standard C, the following code generate 5 random numbers.
#include <stdio.h>
#include <stdlib.h>
int main (int argc, char *argv[])
{
/* Now generate 5 pseudo-random numbers */
int i;
for (i=0; i<5; i++)
{
printf ("rand[%d]= %u\n",
i, rand ());
}
return 0;
}
Each time you run the program, same serial of 5 numbers will be generated since the seed is retested. However, the same code as a mex function does not behave like this. Here is the mex code:
#include <stdio.h>
#include <stdlib.h>
#include "mex.h"
void mexFunction(int nlhs, mxArray *plhs[],
int nrhs, const mxArray *prhs[]) {
/* Now generate 5 pseudo-random numbers */
int i;
for (i=0; i<5; i++)
{
mexPrintf ("rand[%d]= %u\n",
i, rand ());
}
return;
}
Every time this mex function produces different serial of 5 random numbers. Does anybody know why? How the mex function resets the seed number each time? Thanks in advance.
2 Comments
Elliott Rachlin
on 12 Oct 2018
I notice that in a mex file, the code "rand()" seems to return an integer between 0 and 32767, not the value between 0.0 and 1.0 that I was expecting. Does anybody know why this is (and why the documentation does not seem to mention this)?
Walter Roberson
on 12 Oct 2018
That is, you are invoking C's rand(), not invoking MATLAB's rand()
Accepted Answer
James Tursa
on 3 May 2011
When you first run a mex function it loads into memory and stays in memory until it gets cleared. So any variables that are global or persistent in the mex function, like a random number generator seed that the rng internally maintains, will retain the last value between calls. So you get a new set of numbers with each call. If you want to start fresh each time you can either clear the mex function from memory or in the mex code itself reinitialized the rng each time you enter the mex function.
0 Comments
More Answers (1)
James
on 3 May 2011
1 Comment
James Tursa
on 3 May 2011
Yes and no. Using srand (or rand) in one mex routine will affect other mex routines that use rand. However, there is apparently another piece of the library hanging around in memory even after you clear the mex routines, so at least in my testing you will not get a fresh start this way by just doing clear. You must use srand.
See Also
Categories
Find more on Write C Functions Callable from MATLAB (MEX Files) in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!