2

我有我这样实例化的二维数组 C：

const int wA = 16;
float * C[wA];
for(int i = 0; i < hA; i++)
{
C[i] = new float[hA];
for(int i2 = 0; i2 < hA; i2++)
    C[i][i2] = 0;
}

/* looks like this:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
*/

我创建了一个在 C 上运行的内核：

__kernel void simpleMultiply(__global float* outputC, 
                            int widthA, 
                            int heightA, 
                            int widthB, 
                            int heightB, 
                            __global float * inputA, 
                            __global float * inputB)
{
    int row = get_global_id(1);
    int col = get_global_id(0);
    float sum = 0.0f;
    for(int i = 0; i < widthA; i++)
    {
        sum += inputA[row*widthA+i] * inputB[i*widthB+col];
    }
    outputC[row*widthB+col] = sum;
}

一切顺利。从设置上下文到创建缓冲区、创建内核、程序、clEnqueueNDRangeKernel、clEnqueueReadBuffer 等，我一直获得状态的 CL_SUCCESS。

但是当我去阅读输出时它崩溃了。

status = clEnqueueNDRangeKernel(cmdQueue, kernel, 2, NULL, globalws, localws, 0, NULL, NULL);
cout << "\nclEnqueueNDRangeKernel: " << (status == CL_SUCCESS ? "SUCCESS" : "FAIL"); // prints SUCCESS
status = clEnqueueReadBuffer(cmdQueue, bufferC, CL_TRUE, 0, wC*hC*sizeof(float), (void*)C, 0, NULL, NULL);
cout << "\nclEnqueueReadBuffer: " << (status == CL_SUCCESS ? "SUCCESS" : "FAIL"); // prints SUCCESS
cout << "\nC[0][0]: " << C[0][0]; // <--crash

在此处输入图像描述

我对 C++ 和 OpenCL 一样陌生，所以这可能是由于对 C++ 中的数组和指针理解不足。

整个代码在这里

4

1 回答 1

4

数组

float * C[wA];

是一维float *指针数组。因此，您还没有在内存中创建具有连续行和列的二维数组。但是您已经创建了一个指向行的指针数组。

所以你应该在主机数组 C 上展平并像在内核上一样索引它。

float * C;

c = new float[ha * ha]; // Create a  contiguous memory area  to be addressed  in a 2D pattern

memset( C, 0, ha * ha  * sizeof(float) ); // Set all bytes to zero


...

现在您可以在运行内核后寻址 C 数组

cout << C[ icolumn  + irow * ha ];  // icolumn and irow are your row and columns indices
于 2013-07-01T14:36:18.840 回答