[Nek5000-users] Initialise with random noise

[nek5000-users at lists.mcs.anl.gov]

Hi,
We had the same issue a while back; I guess the simple solution is to 
use different seeds for the pseudo random number generator depending on 
the processor you are on; in this way each processor will generate an 
independent series of random numbers. Depending on which generator you 
use, and what requirements you have (always the same sequence of numbers 
or completely random, etc.), maybe it is enough to just call 
random_seed() on each process before you pull out the numbers.

Philipp

On 2014-07-17 12:40, nek5000-users at lists.mcs.anl.gov wrote:
> Hi again Nek's,
>
> I am working on the simple pipe flow as well, and I have a question
> regarding the initialisation of the perturbation mode. I initialise my
> perturbation with call random_number(vxp) call random_number(vyp) and
> call random_number(vzp) and run the calculation on 4 processors for
> validation and tests purposes. The domain is divided in four partitions
> in the streamwise direction. Unfortunately, the call random_number
> repoduces four times the same random field and as a consequence
> prescribes a streamwise wavenumber which I'd like to get rid of. Any
> ideas on how to do that?
>
> Regards,
> JC
>
> --
> Jean-Christophe Loiseau
> Homepage <https://sites.google.com/site/loiseaujc/>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>