CUDAでのシーケンシャルプログラミング

Question

私はこの単純なループをCUDAで実現したいと思います。

for (int i = 1; i <= N; i++) 
    { 
     for (int j = 1; j <= N; j++) 
     { 
     x[i, j] = (x0[i, j] + a*(x[i - 1, j] + x[i + 1, j] + x[i, j - 1] + x[i, j + 1]))/c; 
     } 
    } 

問題は次のとおりです。 X [i、j] X（i-1、j）とX [i、j-1]の新しい値を知る必要がある。しかし、GPUはすべてを並列に計算しているので、CPUとGPUから受け取った結果は異なります。私はCUDAとcudaDeviceSynchronize（）について動的パラレル化についての情報を見つけましたが、それは有用かもしれませんが、とにかくカーネルにこのループを実装する方法はまだ分かりません。私はあなたの助けに感謝します。

Answer 1

上記のコメントは、純粋な順次実装であっても、データをコピーする必要があります。あなたのhardareに合ったどんな値で打ち上げ境界を置き換える

update<<<dim3(16, 16), dim3(64, 64)>>>(input, output, width, height); 

使用して、

__global__ void update(...) 
    for(int i = threadIdx.x + blockDim.x * blockIdx.x; i <= N; i += blockDim.x * gridDim.x) 
{ 
    for(int j = threadIdx.y + blockDim.y * blockIdx.y; j <= N; j += blockDim.y * gridDim.y) 
    { 
     output[i,j] = update_func(input, i, j); 
    } 
} 

あなたは（ホストから）呼び出すことができます。ここで（メモリ管理コードまたは任意のさらなる詳細なし）カーネルかもしれません