HLSL Tips 2：伪随机数生成器

Preface
PRNG on the GPU. Two pseudo random numbers per iteration (from the BA channels). Inspired by Pete Warden's fragment program for pseudo-random number generation. In turn based on an algorithm described by Francois Grieu, sci.crypt, 5th February 2004. Implementd in Cg by Alex Campbell. 7th August 2006. Thanks for their works. I implement it in HLSL now.

Function
#define cMult 0.0001002707309736288
#define aSubtract 0.2727272727272727
float4 randGrieu(float4 t)
{
float a=t.x+t.z*cMult+aSubtract-floor(t.x);
a*=a;
float b=t.y+a;
b-=floor(b);
float c=t.z+b;
c-=floor(c);
float d=c;
a+=c*cMult+aSubtract-floor(a);
a*=a;
b+=a;
b-=floor(b);
c+=b;
c-=floor(c);
return float4(a,b,c,d);
}

Application
Code snippet:
float4 color=tex2D(TextureSampler, In.TexCoord0);
color=randGrieu(color);
Out.color=color;

The original texture is as follows:

 

After called randGrieu() once, the result texture is: