有关使用OpenCL处理的大型数组的错误信息

时间:2012-05-27 16:46:33

标签: c opencl

我有一个内核执行问题,当我使用大型数组(1000 x 10000)时,它不会在正确的位置写入信息。但对于小阵列,没有问题,我检索正确的结果。对于内核执行,我使用的是ATI Mobility RADEON HD 4300系列的GPU。

C代码示例是:

#include <stdio.h>
#include <stdlib.h>

#ifdef __APPLE__
#include <OpenCL/opencl.h>
#else
#include <CL/cl.h>
#endif

#define MAX_SOURCE_SIZE (0x100000)
#define MAX_SIZE 108
#define NCOLS 1000
#define NROWS 10000

int main(void) {
    char* source_name = "mykernel.cl";
    char* source_code;
    size_t source_size;
    cl_platform_id platformId = NULL;
    cl_uint nbplatforms;
    cl_device_id deviceId = NULL;
    cl_uint nbdevices;
    cl_context context = NULL;
    cl_int errcode;
    cl_command_queue commandQueue = NULL;
    cl_program program;
    size_t global_work_size[2];
    size_t local_work_size[2];

    FILE* fh;

    //Retrieving platform information
    errcode = clGetPlatformIDs(1, &platformId, &nbplatforms);

    //Retrieving device (GPU) information
    errcode = clGetDeviceIDs(platformId, CL_DEVICE_TYPE_GPU, 1, &deviceId, &nbdevices);

    //Creation of a working context
    context = clCreateContext(NULL, 1, &deviceId, NULL, NULL, &errcode);

    commandQueue = clCreateCommandQueue(context, deviceId, 0, &errcode);

    //Opening and reading the kernel source file
    if((fh = fopen(source_name, "r")) == NULL){
        fprintf(stderr, "Failed to open the file containing the kernel source !\n");
        exit(EXIT_FAILURE);
    }
    source_code = (char*) malloc (MAX_SOURCE_SIZE * sizeof(char));
    source_size = fread(source_code, sizeof(char), MAX_SOURCE_SIZE, fh);
    fclose(fh);

    program = clCreateProgramWithSource(context, 1, (const char**) &source_code, (const size_t*) &source_size, &errcode);

    //Building kernel
    errcode = clBuildProgram(program, 1, &deviceId, NULL, NULL, NULL);

    //Creation of the kernel program
    cl_kernel kernel = clCreateKernel(program, "mykernel", &errcode);
    unsigned int *op1 = (unsigned int*) malloc (NCOLS * NROWS * sizeof(unsigned int));

    cl_mem op1buff = clCreateBuffer(context, CL_MEM_WRITE_ONLY, NCOLS * NROWS * sizeof(unsigned int), NULL, &errcode);

    clSetKernelArg(kernel, 0, sizeof(cl_mem), (void*) &op1buff);

    global_work_size[0] =  NCOLS;
    global_work_size[1] =  NROWS;
    local_work_size[0] = NCOLS;
    local_work_size[1] = 1;

    clEnqueueNDRangeKernel(commandQueue, kernel, 2, NULL, global_work_size, local_work_size, 0, NULL, NULL);

    errcode = clEnqueueReadBuffer(commandQueue, op1buff, CL_TRUE, 0, NCOLS * NROWS * sizeof(unsigned int), (void*)op1, 0, NULL, NULL);

    for(int i = 0; i < NROWS; i++){
        for(int j = 0; j < NCOLS; j++)
            printf("[index:%d - %u] ", i*NCOLS+j, op1[i*NCOLS+j]);
        printf("\n");
    }

    return EXIT_SUCCESS;
}

内核源代码放在名为mykernel.cl的文件中,如下所示:

__kernel void mykernel(__global unsigned int* op1buf){
    unsigned int index = get_group_id(1) * get_global_size(0) + get_local_id(0);
    op1buf[index] = index;
}

当我使用大型数组时,此程序的执行将返回从数组中读取的意外值。例如:

[index:0 - 16777215] [index:1 - 16777215] [index:2 - 16777215] [index:3 - 16777215] ...
[index:1000 - 3438339071] [index:1001 - 3941660159] [index:1002 - 1650092117] [index:1003 - 2529976771] ...
[index:1000 - 3438339071] [index:1001 - 3941660159] [index:1002 - 1650092117] [index:1003 - 2529976771] ...
[index:3000 - 16777215] [index:3001 - 16777215] [index:3002 - 16777215] [index:3003 - 16777215] ...
[index:4000 - 3438339071] [index:4001 - 3941660159] [index:4002 - 1650092117] [index:4003 - 2529976771] ...
....

我的代码可能是什么问题,或者是否有关于使用GPU的问题我没有考虑过?

提前致谢。

1 个答案:

答案 0 :(得分:1)

1000显然对您的设备来说太大了。使用clGetDeviceInfo和CL_DEVICE_MAX_WORK_GROUP_SIZE来确定可以使用的最大值。