1

大きな配列 (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の使用に関して考慮していないことがありますか?

前もって感謝します。

4

1 に答える 1

1

1000 は明らかにデバイスには大きすぎます。使用できる最大値を決定するには、CL_DEVICE_MAX_WORK_GROUP_SIZE で clGetDeviceInfo を使用します。

于 2012-05-30T20:32:51.480 に答える