gpurng

Control random number generation on the GPU

collapse all in page

Syntax

gpurng("default")

gpurng(seed)

gpurng(seed,generator)

gpurng(generator)

gpurng(S)

S = gpurng

S = gpurng(___)

Description

gpurng("default") initializes the GPU random number generator using the default algorithm and seed. The default algorithm is the Threefry generator with seed 0. The random numbers produced are the same as if you had restarted MATLAB^®.

The gpurng function controls the global GPU stream, which determines how the rand, randi, randn, and randperm functions produce a sequence of random numbers on the GPU. To create one or more independent streams separate from the global GPU stream, see parallel.gpu.RandStream.

gpurng(seed) specifies the seed for the GPU random number generator using the current generator algorithm.

Specify seed as a nonnegative integer, such as gpurng(1), to initialize the GPU random number generator with that seed.
Specify seed as "shuffle" to initialize the generator seed based on the current time so that rand, randi, randn, and randperm produce different sequences of numbers after each time you call gpurng.

gpurng(seed,generator) also specifies the algorithm for the GPU random number generator to use. For example, gpurng(2,"philox") initializes the Philox 4x32 generator with a seed of 2.

example

gpurng(generator) specifies the algorithm for the GPU random number generator to use with a seed of 0. This syntax is equivalent to gpurng(0,generator). (since R2023b)

gpurng(S) initializes the state of the random number generator based on settings contained in a structure S with fields Type, Seed, and State. The structure S must be a structure that is returned by a previous call to S = gpurng or S = gpurng(__).

S = gpurng returns the current state of the random number generator as a structure S with fields Type, Seed, and State.

S = gpurng(___) returns the current state of the random number generator in a structure S before changing the settings using the specified arguments.

Note

The default algorithm and seed for the random number generator in the MATLAB Preferences window affect random numbers generated on the CPU only and do not affect calls to gpurng("default").

Examples

collapse all

Create Predictable Arrays of Random Numbers on the GPU and CPU

Open Live Script

Capture the GPU generator settings, and set the state of the CPU random number generator to match the GPU generator settings. Create predictable arrays of random numbers on the CPU and GPU.

Restore the generator type and seed to their default values on both the CPU and the GPU.

gpurng("default") 
rng("default")

Save the default seed and generator type of the GPU random number generator.

GPUdef = gpurng

GPUdef = struct with fields:
     Type: 'threefry'
     Seed: 0
    State: [17×1 uint32]

Set the CPU random number generator to match the default GPU settings.

rng(GPUdef)

Create an array of uniformly distributed random numbers on the GPU.

rGPU = rand(1,10,"gpuArray")

rGPU =

    0.3640    0.5421    0.6543    0.7436    0.0342    0.8311    0.7040    0.2817    0.1163    0.5671

Create an array of random numbers on the CPU.

rCPU = rand(1,10)

rCPU = 1×10

    0.3640    0.5421    0.6543    0.7436    0.0342    0.8311    0.7040    0.2817    0.1163    0.5671

The seed and generator type are the same for both the GPU and the CPU, so the arrays are the same.

isequal(rGPU,rCPU)

ans = logical
   1

The gpurng state does not save the settings for the transformation applied to generate a normally distributed set of random numbers. Even though the seed and the generator type are the same on the GPU and the CPU, the set of normally distributed random numbers is different.

nGPU = randn(1,1000,"gpuArray");
nCPU = randn(1,1000);

figure
hold on
histogram(nGPU)
histogram(nCPU)
legend("GPU","CPU")
title("Normally Distributed Random Numbers")
xlabel("Value")
ylabel("Count")
hold off

Figure contains an axes object. The axes object with title Normally Distributed Random Numbers, xlabel Value, ylabel Count contains 2 objects of type histogram. These objects represent GPU, CPU.

The statistics of the normal distribution of random numbers are the same on the GPU and the CPU.

By default, the CPU uses the Ziggurat transformation, while the GPU uses the BoxMuller algorithm for the Threefry generator. The only transformation method supported on both the CPU and GPU is the Inversion transform.

You can change the transformation method on the GPU using parallel.gpu.RandStream.

Input Arguments

collapse all

`seed` — Random number seed
0 (default) | nonnegative integer | "shuffle"

Random number seed, specified as a nonnegative integer or "shuffle". The seed specifies the starting point for the algorithm to generate random numbers. Specify the seed as a nonnegative integer when you want reproducible results. The default seed is 0.

When you specify the seed as "shuffle", the software initializes the generator seed based on the current time, resulting in a different sequence of random numbers after each call to gpurng.

Example: gpurng(7)

`generator` — Random number generator
`"Threefry"` (default) | `"Philox` | `"CombRecursive"`

Random number generator, specified as a character vector or string for any valid random number generator that supports multiple streams and substreams. Three random number generator algorithms are supported on the GPU.

Keyword	Generator	Multiple Stream and Substream Support	Approximate Period in Full Precision
`"Threefry"` or `"Threefry4x64_20"`	Threefry 4x64 generator with 20 rounds	Yes	2⁵¹⁴ (2²⁵⁶ streams of length 2²⁵⁸)
`"Philox"` or `"Philox4x32_10"`	Philox 4x32 generator with 10 rounds	Yes	2¹⁹³ (2⁶⁴ streams of length 2¹²⁹)
`"CombRecursive"` or `"mrg32k3a"`	Combined multiple recursive generator	Yes	2¹⁹¹ (2⁶³ streams of length 2¹²⁷)

The default generator is Threefry.

For more information on the differences between generating random numbers on the GPU and CPU, see Control Random Number Streams on Workers.

Example: gpurng("Philox")

`S` — Previous random number generator state
structure

Previous random number generator state, specified as a structure previously created using S = gpurng.

Example: S = gpurng captures the current state of the random number generator, and gpurng(S) restores the generator to those settings.

Data Types: struct

Output Arguments

collapse all

`S` — Random number generator state
structure

Random number generator state, returned as a structure with fields Type, Seed, and State.

Example: S = gpurng captures the current state of the random number generator, and gpurng(S) restores the generator to those settings.

Data Types: struct

Extended Capabilities

Thread-Based Environment
Run code in the background using MATLAB® `backgroundPool` or accelerate code with Parallel Computing Toolbox™ `ThreadPool`.

This function fully supports thread-based environments. For more information, see Run MATLAB Functions in Thread-Based Environment.

Version History

Introduced in R2011b

expand all

R2023b: Specify random number algorithm without specifying seed

Use the new syntax gpurng(generator) to specify the algorithm for the random number generator to use. This syntax allows you to set the random number algorithm without specifying the seed, where gpurng uses a seed of 0. This syntax is equivalent to gpurng(0,generator). For example, gpurng("philox") initializes the Philox 4x32 generator with a seed of 0.

R2019a: Default random number generator change for `gpurng`

Starting in R2019a, the default random number generator for parallel computations is changed to Threefry. This generator offers performance enhancements for parallel calculations over the previous default. In releases up to R2018b, the default random number generator for parallel computations is CombRecursive.

With a different default generator, MATLAB generates different random numbers sequences by default in the context of parallel computations. However, statistics of these calculations remain unaffected. Therefore, you might want to update any code that relies on the specific random numbers being generated, but most calculations on the random numbers are unaffected.

To set the generator to the settings used by default in R2018b and earlier on GPU arrays, use the following command.

gpurng(0,"CombRecursive")

R2018a: `parallel.gpu.rng` is renamed to `gpurng`

Starting in R2018a, the function parallel.gpu.rng is renamed to gpurng. Replace all instances of parallel.gpu.rng with gpurng.

parallel.gpu.rng will continue to work but is not recommended.

gpurng

Syntax

Description

Examples

Create Predictable Arrays of Random Numbers on the GPU and CPU

Input Arguments

`seed` — Random number seed
0 (default) | nonnegative integer | "shuffle"

`generator` — Random number generator
`"Threefry"` (default) | `"Philox` | `"CombRecursive"`

`S` — Previous random number generator state
structure

Output Arguments

`S` — Random number generator state
structure

Extended Capabilities

Thread-Based Environment
Run code in the background using MATLAB® `backgroundPool` or accelerate code with Parallel Computing Toolbox™ `ThreadPool`.

Version History

R2023b: Specify random number algorithm without specifying seed

R2019a: Default random number generator change for `gpurng`

R2018a: `parallel.gpu.rng` is renamed to `gpurng`

See Also

Topics

gpurng

Syntax

Description

Examples

Create Predictable Arrays of Random Numbers on the GPU and CPU

Input Arguments

seed — Random number seed 0 (default) | nonnegative integer | "shuffle"

generator — Random number generator "Threefry" (default) | "Philox | "CombRecursive"

S — Previous random number generator state structure

Output Arguments

S — Random number generator state structure

Extended Capabilities

Thread-Based Environment Run code in the background using MATLAB® backgroundPool or accelerate code with Parallel Computing Toolbox™ ThreadPool.

Version History

R2023b: Specify random number algorithm without specifying seed

R2019a: Default random number generator change for gpurng

R2018a: parallel.gpu.rng is renamed to gpurng

See Also

Topics

`seed` — Random number seed
0 (default) | nonnegative integer | "shuffle"

`generator` — Random number generator
`"Threefry"` (default) | `"Philox` | `"CombRecursive"`

`S` — Previous random number generator state
structure

`S` — Random number generator state
structure

Thread-Based Environment
Run code in the background using MATLAB® `backgroundPool` or accelerate code with Parallel Computing Toolbox™ `ThreadPool`.

R2019a: Default random number generator change for `gpurng`

R2018a: `parallel.gpu.rng` is renamed to `gpurng`