Random number generator: cpu/gpu

[navidcy]

If dev=GPU() I think what this does is that it creates random numbers on CPU and transfers then to GPU

GeophysicalFlows.jl/examples/barotropicqg_betaforced.jl

Line 72 in de360d7

ξ = ArrayType(dev)(exp.(2π*im*rand(eltype(grid), size(sol)))/sqrt(clock.dt))

Is this true? If so, we want to be creating the random numbers on the GPU itself, e.g.,

Random.rand(CURAND.default_rng(), Float32, 1)

?

[glwagner]

Yes definitely. When this example was written, random number generators on the GPU were not available. Let's update it.

[navidcy]

Closed by #225.