cuda のバックアップ(No.2)

バックアップ一覧
差分を表示
現在との差分を表示
ソースを表示
cuda へ行く。
- 1 (2009-07-07 (火) 12:04:04)
- 2 (2009-07-07 (火) 12:05:37)
- 3 (2010-05-21 (金) 08:26:12)

CUDA プログラミング †

ランダム関数(GAP のルーチンを移植) †

#define r228 (1<<28)

__device__ int myrandom(int max, int rn, long rx[])
{
    rn = (rn%55) + 1;
    rx[rn-1] = (rx[rn-1] + rx[(rn+30)%55])%r228;
    return rx[rn-1]*max/r228;
}

__device__ void myseed(int n, int rn, long rx[])
{
  int i;
  rn = 1;
  rx[0] = n%r228;
  for(i = 2 ; i <= 55 ; i++)rx[i-1]=(1664525*rx[i-2]+1)%r228;
  for(i = 1 ; i < 100 ; i++){
    rn = (rn%55) + 1;
    rx[rn-1] = (rx[rn-1]+rx[(rn+30)%55])%r228;
  }

}