问使用OpenCL的Random123随机数

EN

是。该示例代码在每次调用时都生成相同的随机数序列。

要获得不同的随机数流，只需以不同的方式初始化任意值K1.3和/或c1..3。您可以通过命令行参数、环境变量、每日时间、保存状态、/dev/urandom或任何其他源来初始化它们。只是要意识到：

( a)如果在两次不同的运行中以完全相同的方式初始化它们，那么这两次运行将得到相同的随机数流。

如果在两次不同的运行中对它们进行不同的初始化，那么这两次运行将得到不同的随机数流。

有时你想要财产a)。有时你想要财产b)。花点时间想想你想要什么，并确保你正在做你想做的事情。

更普遍地说，库中的函数(例如threefry4x32 )具有无状态。如果您更改输入中的任何位(即c或k的任何元素中的任何位)，则会得到完全不同的随机、统计独立、均匀分布的输出。

P.S.我是该图书馆的作者之一，论文“并行数字:和1，2，3一样容易”：http://dl.acm.org/citation.cfm?id=2063405

如果您不是ACM数字图书馆的订阅者，上面的链接可能会碰到付费墙。或者，您可以通过以下链接免费获得纸张：

http://www.thesalmons.org/john/random123/index.html

我无法帮助您处理库本身，但我可以告诉您，在OpenCL中生成随机数的最常见方法是在调用内核之间保存一些状态。

随机数发生器通常使用一种状态，由此生成一个新的状态和一个随机数。实际上，这一点也不复杂:您只需传递一个包含状态的额外数组。在我的代码中，我实现了以下随机数：

uint rand_uint(uint2* rvec) {  //Adapted from http://cas.ee.ic.ac.uk/people/dt10/research/rngs-gpu-mwc64x.html
    #define A 4294883355U
    uint x=rvec->x, c=rvec->y; //Unpack the state
    uint res = x ^ c;          //Calculate the result
    uint hi = mul_hi(x,A);     //Step the RNG
    x = x*A + c;
    c = hi + (x<c);
    *rvec = (uint2)(x,c);      //Pack the state back up
    return res;                //Return the next result
    #undef A
}
inline float rand_float(uint2* rvec) {
    return (float)(rand_uint(rvec)) / (float)(0xFFFFFFFF);
}
__kernel void my_kernel(/*more arguments*/ __global uint2* randoms) {
    int index = get_global_id(0);
    uint2 rvec = randoms[index];

    //Call rand_uint or rand_float a number of times with "rvec" as argument.
    //These calls update "rvec" with new state, and return a random number

    randoms[index] = rvec;
}

。。。然后，您所要做的就是传递一个额外的数组，该数组将RNG的状态随机保存起来。实际上，您需要为每个工作项设置不同的数组种子。

0xdecafbad，0xfacebead，0x12345678和0xf00dcafe，0xdeadbeef，0xbeeff00d只是任意选择的数字，它们不是特殊的。任何其他数字(甚至0)都可以用来代替它们--我将在示例代码中添加一个注释。

您可以用传入的变量替换它们中的任何一个；在输出随机“流”中避免不必要的重复的唯一要求是避免重复(c，k)输入元组。示例代码使用线程id和循环索引来确保唯一性，但您可以轻松地添加更多变量以确保唯一性--例如，在主机代码中计算内核调用并传入计数器，使用它代替k或c的元素之一。

顺便说一句，尽管名字是“基于计数器的随机数生成器”，但没有要求输入(c，k)是“计数器”，只是计数器恰好是确保输入不重复的最方便的成语。