I have problem with truly! random initialization of a sample code when I run it in parallel mode.
Please check the code in the folder ~/examples/codesByTopics/shanChenMultiPhase.
When I run it on a single core, everything is OK. But in parallel mode (like with 4 cores), the initial field is not truly random. Actually, the 4 blocks of the domain have the same random initial condition.
It seems the command line:
T rho = rho0 + ((T)random()/(T)RAND_MAX)*maxRho;
does not work properly in parallel mode.