我是OpenCL的新手,可以从Matthew Scarpino's “OpenCL in Action”中学习。我研究了一个矩阵与矢量相乘的示例(第11-13页)。由于某种原因,该示例在我的计算机上不起作用。内核未返回值。我开始寻找从内核输出数据的简单示例。
我在Wesley Shillingford的youtube频道中发现,输出字符串“ Hello world!”的example。从内核。在我的家用计算机上,该示例有效。但是,自提供的示例用C ++编写以来,OpenCL“厨房”一直关闭。代码的简短性掩盖了正在发生的事情的概念。因此,我开始进一步寻找C代码中的示例。
在Stackoverflow的答案中,我发现了一个example of a minimal OpenCL program,它会增加内核中的值。我将这段代码作为编写程序的基础,因为它对于初学者来说既简单又方便。后来我发现,该示例包含一个错误。
另一个great example说服我使用指针从内核返回数据。使用数组存储内核的输出值会导致以下事实:目标数组的值不会更改,并且内核的数据在输出期间会消失。我意识到我们需要使用指针从内核输出数据。但是,这没有帮助我。从内核向主机程序传输数据的问题仍然存在。如果我在某些方面弄错了,请纠正我。主题的实质:内核不会将字符数据返回到主机程序。可能是什么问题?
#include <CL/cl.h>
#include <stdio.h>
#include <stdlib.h>
int main(){
cl_platform_id *platforms =NULL;
cl_device_id *devices=NULL;
cl_context context;
cl_command_queue cmdQueue;
cl_program program;
cl_kernel kernel = NULL;
char *cpOutputData;
int output_size = 8;
cl_mem output_buff;
cl_int status; // to check the output of each API call
const char *source =
"__kernel void Hello( __global char* ch) {\n"
" ch[0]='P';"
" ch[1]='r';"
" ch[2]='i';"
" ch[3]='v';"
" ch[4]='e';"
" ch[5]='t';"
" ch[6]='!';"
" ch[7]='\0';"
"}\0";
printf("GetPlatformIDs... ");
cl_uint numPlatforms = 0;
//STEP 1: Discover and initialize platforms
// Use clGetPlatformIDs to retreive the number of platforms
status = clGetPlatformIDs(0,
NULL,
&numPlatforms);
// Allocate enough space for each platform
platforms = (cl_platform_id*)malloc(numPlatforms*sizeof(cl_platform_id));
// Fill in platforms with clGetPlatformIDs()
status=clGetPlatformIDs(numPlatforms,
platforms,
NULL);
printf("\nNumber of discovered platforms is %d. ", numPlatforms);
// STEP 2: Discover and initialize devices
printf("OK.\nGetDeviceIDs... ");
cl_uint numDevices = 0;
// Use clGetDeviceIDs() to retrieve the number of devices present
status = clGetDeviceIDs(platforms[0],
CL_DEVICE_TYPE_CPU,
0,
NULL,
&numDevices);
// Allocate enough space for each device
devices = (cl_device_id*)malloc(numDevices*sizeof(cl_device_id));
// Fill in devices with clGetDeviceIDs()
clGetDeviceIDs(platforms[0],
CL_DEVICE_TYPE_CPU,
numDevices,
devices,
NULL);
printf("\nNumber of discovered devices is %d. ", numDevices);
// STEP 3: Create a context
printf("OK.\nCreating context... ");
// Create context using clCreateContext() and associate it with the devices
context = clCreateContext(NULL,
numDevices,
devices,
NULL,
NULL,
&status);
// STEP 4: Create a command queue
printf("OK.\nQueue creating... ");
cmdQueue = clCreateCommandQueue(context,
devices[0],
CL_QUEUE_PROFILING_ENABLE,
&status);
// STEP 5: Create device buffers
printf("OK.\nOutput buffer creating... ");
output_buff = clCreateBuffer(context,
CL_MEM_WRITE_ONLY,
sizeof(char)*output_size,
NULL,
&status);
// STEP 6: Create and compile program
printf("OK.\nBuilding program... ");
// Create a program using clCreateProgramWithSource()
program = clCreateProgramWithSource(context,
1,
(const char**)&source,
NULL,
&status);
// Build (compile) the program for the devices with clBuildProgram()
status=clBuildProgram(program,
numDevices,
devices,
NULL,
NULL,
NULL);
// STEP 7: Create a kernel
printf("OK.\nCreating kernel... ");
kernel = clCreateKernel(program,
"Hello",
&status);
// STEP 8: Set kernel arguments
// Associate ouput buffer with the kernel
printf("OK.\nSetting kernel arguments... ");
status = clSetKernelArg(kernel,
0,
sizeof(cl_mem),
&output_buff);
// STEP 9: Configure the work-item structure
// Define an index space (global work size) of work itmes for execution.
// A workgroup size (local work size) is not required, but can be used.
size_t globalWorkSize[1];
// There are 'elements' work-items
globalWorkSize[0] = output_size;
// STEP 10: Enqueue the kernel for execution
printf("OK.\nExecuting kernel... ");
//Execute the kernel by using clEnqueueNDRangeKernel().
// 'globalWorkSize' is the 1D dimension of the work-items
clEnqueueNDRangeKernel(cmdQueue,
kernel,
1,
NULL,
globalWorkSize,
NULL,
0,
NULL,
NULL);
clFinish(cmdQueue);
// STEP 11: Read the ouput buffer back to the host
printf("OK.\nReading buffer... ");
// Allocate space for the data to be read
cpOutputData = (char*)malloc(output_size*sizeof(char));
// Use clEnqueueReadBuffer() to read the OpenCL ouput buffer to the host ouput array
clEnqueueReadBuffer(cmdQueue,
output_buff,
CL_TRUE,
0,
output_size,
cpOutputData,
0,
NULL,
NULL);
printf("\nPrinting output data: \n");
printf(cpOutputData);
// STEP 12: Releasing resources
printf("\n...Releasing OpenCL resources... ");
clReleaseKernel(kernel);
clReleaseProgram(program);
clReleaseCommandQueue(cmdQueue);
clReleaseMemObject(output_buff);
clReleaseContext(context);
printf("OK.\n...Releasing host resources... ");
free(cpOutputData);
free(platforms);
free(devices);
printf("OK.\nEnd of program. Bey!\n");
system("PAUSE");
return 0;
}
我的程序的执行输出is here。
您遇到的问题非常细微,不幸的是,您没有在一个可以抓住它的地方进行错误检查。具体来说,使用clBuildProgram
编译内核的源代码将失败,不幸的是,未选中status
。我不确定为什么程序的其余部分不会在您的实现上产生错误,而肯定会在我的系统上产生错误。
您的内核源无效的原因是此行:
" ch[7]='\0';"
// ^^---- This terminates the string early!
基本上,您的内核源代码在OpenCL编译器中看起来像这样:
__kernel void Hello( __global char* ch) {
ch[0]='P';
ch[1]='r';
ch[2]='i';
ch[3]='v';
ch[4]='e';
ch[5]='t';
ch[6]='!';
ch[7]='
因为字符串文字中的转义码\0
在source
变量最终指向的内存中插入了一个实际的nul字符,导致该字符被视为内核源代码的结尾。
您真正想要的是让转义序列出现在您的OpenCL内核代码中,因此您需要对它进行两次转义:一次是对主机程序的C编译器,而第二次是对OpenCL编译器。那将是:
" ch[7]='\\0';"
// ^--- note second backslash
在您的source
字符串中,双反斜杠被转换为单个反斜杠,在其中OpenCL编译器将其与随后的零组合以将字符文字转换为nul字符。
有了该修复程序,一切正常!
我建议在单独的文件中编写内核源代码。您可以使用程序中的文件I / O加载该文件,也可以自动生成用于将数据嵌入到源代码中的文字。 unix tool xxd
可以使用xxd
标志来执行此操作,您可能可以找到Windows等效版本,甚至可以找到该工具本身的Windows版本。