Question

我需要一种算法来计算数组的并行前缀和，而不使用共享内存。如果没有其他替代方法可以使用共享内存，那么解决冲突问题的最佳方法是什么？

Answer 1

此链接包含并行前缀和的顺序算法和并行算法的详细分析。

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

它还包含用于实现并行前缀算法的C代码片段以及避免共享内存冲突的详细说明。

您可以将代码移植到CUDAfy，也可以只定义C的区域，并将它们用作应用程序中的非托管代码。但是CUDA C代码中存在一些错误。我正在使用Cudafy.NET

编写更正版本的代码

[Cudafy]
public static void prescan(GThread thread, int[] g_odata, int[] g_idata, int[] n)
{
    int[] temp = thread.AllocateShared<int>("temp", threadsPerBlock);//threadsPerBlock is user defined
    int thid = thread.threadIdx.x;
    int offset = 1;
    if (thid < n[0]/2)
    {
        temp[2 * thid] = g_idata[2 * thid]; // load input into shared memory  
        temp[2 * thid + 1] = g_idata[2 * thid + 1];

        for (int d = n[0] >> 1; d > 0; d >>= 1)                    // build sum in place up the tree  
        {
            thread.SyncThreads();
            if (thid < d)
            {
                int ai = offset * (2 * thid + 1) - 1;
                int bi = offset * (2 * thid + 2) - 1;
                temp[bi] += temp[ai];
            }
            offset *= 2;
        }
        if (thid == 0)
        {
            temp[n[0] - 1] = 0;
        } // clear the last element  


        for (int d = 1; d < n[0]; d *= 2) // traverse down tree & build scan  
        {
            offset >>= 1;
            thread.SyncThreads();
            if (thid < d)
            {
                int ai = offset * (2 * thid + 1) - 1;
                int bi = offset * (2 * thid + 2) - 1;
                int t = temp[ai];
                temp[ai] = temp[bi];
                temp[bi] += t;
            }
        }
        thread.SyncThreads();
        g_odata[2 * thid] = temp[2 * thid]; // write results to device memory  
        g_odata[2 * thid + 1] = temp[2 * thid + 1];
    }
}

您可以使用上面修改的代码而不是链接中的代码。

并行前缀和CUDAfy

1 个答案: