Random number for GPU

最新推荐文章于 2024-10-20 20:55:43 发布

原创最新推荐文章于 2024-10-20 20:55:43 发布 · 863 阅读

0 ·

CC 4.0 BY-SA版权

计算机图形技术基础专栏收录该内容

41 篇文章

订阅专栏

本文探讨了使用CUDA进行高效伪随机数生成的方法和技术，包括cuRAND库、快速GPU随机数生成、高质量伪随机数生成器以及GPU伪随机数生成库GPUTEA。此外，介绍了基于“TEA”加密算法的伪随机数生成函数，并提供了Morton编码和无符号整数位反转等技巧。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Efficient random number generation and application using CUDA

http://http.developer.nvidia.com/GPUGems3/gpugems3_ch37.html

cuRAND: Nvidia random number generation library

https://developer.nvidia.com/cuRAND

Quick and easy gpu random numbers in D3D11:

http://www.reedbeta.com/blog/2013/01/12/quick-and-easy-gpu-random-numbers-in-d3d11/

A Fast high quality pseudo random number generator for nVidia CUDA

http://www0.cs.ucl.ac.uk/staff/ucacbbl/ftp/papers/langdon_2009_CIGPU.pdf

GPU pseudo random number: GPU TEA

http://www.csee.umbc.edu/~olano/papers/GPUTEA.pdf

/*
* Pseudo random number generator, based on "TEA, a tiny Encrytion Algorithm"
* http://citeseer.ist.psu.edu/viewdoc/download?doi=10.1.1.45.281&rep=rep1&type=pdf
* @param v - old seed (full 32bit range)
* @param IterationCount - >=1, bigger numbers cost more performance but improve quality
* @return new seed
*/
uint2 ScrambleTEA(uint2 v, uint IterationCount = 3)
{
// Start with some random data (numbers can be arbitrary but those have been used by others and seem to work well)
uint k[4] = { 0xA341316Cu, 0xC8013EA4u, 0xAD90777Du, 0x7E95761Eu };

uint y = v[0];
uint z = v[1];
uint sum = 0;

for (uint i = 0; i < IterationCount; ++i)
{
sum += 0x9e3779b9;
y += (z << 4u) + k[0] ^ z + sum ^ (z >> 5u) + k[1];
z += (y << 4u) + k[2] ^ y + sum ^ (y >> 5u) + k[3];
}

return uint2(y, z);
}

uint MortonCode(uint x)
{
x = (x ^ (x << 2)) & 0x33333333;
x = (x ^ (x << 1)) & 0x55555555;
return x;
}

uint ReverseUIntBits(uint bits)
{
bits = ((bits & 0x33333333) << 2) | ((bits & 0xcccccccc) >> 2);
bits = ((bits & 0x55555555) << 1) | ((bits & 0xaaaaaaaa) >> 1);
return bits;
}