C++ OpenCL内核没有';我不能被处决
有一组规则,其中每个规则对应于ceratain细胞自动机。我需要检查每个规则的双射性属性。由于它们太多了(准确地说是2^32),我决定使用我的GPU来实现这个目的。但大约一周后,我仍在与一个bug作斗争。C++ OpenCL内核没有';我不能被处决,c++,c,opencl,C++,C,Opencl,有一组规则,其中每个规则对应于ceratain细胞自动机。我需要检查每个规则的双射性属性。由于它们太多了(准确地说是2^32),我决定使用我的GPU来实现这个目的。但大约一周后,我仍在与一个bug作斗争。 简单地说,当内核排队时,它的执行应该在GPU上执行,GPU的使用就好像它是空闲的一样。此外,在我向内核代码中添加了几个语句以查看内核是否正在执行之后,我发现没有任何迹象表明这些语句以及内核本身已经执行。此外,所有错误代码都等于CL_成功。我可能会出错,因为我是OpenCL编程新手,非常感谢您的
简单地说,当内核排队时,它的执行应该在GPU上执行,GPU的使用就好像它是空闲的一样。此外,在我向内核代码中添加了几个语句以查看内核是否正在执行之后,我发现没有任何迹象表明这些语句以及内核本身已经执行。此外,所有错误代码都等于CL_成功。我可能会出错,因为我是OpenCL编程新手,非常感谢您的帮助。
这是主机端代码,带有一些缩写:
#define CL_USE_DEPRECATED_OPENCL_2_0_APIS
//some includes here
#define GLOBAL_SIZE 4096
#define LOCAL_SIZE 256
#define GLOBAL_SCOPE 0xffffffff
int main()
{
//we assume that global_scope divides evenly into global_size
//and therefore there is no need in processing remainder
long rules_per_thread = GLOBAL_SCOPE / GLOBAL_SIZE;
int * starts = new int[GLOBAL_SIZE];
int * stops = new int[GLOBAL_SIZE];
int count = 0;
for (int i = 0; i < GLOBAL_SIZE; i++) {
starts[i] = count;
count += rules_per_thread;
stops[i] = count;
count++;
}
...
/*obtainig platform, device, building program*/
...
/*====CREATING BUFFERS====*/
//buffer for storing approved automata
const int bufSize = 10000; //size of buffer picked at random guess; might need to add some space later
uint32_t* bijective_aut = new uint32_t[bufSize];
std::fill(&bijective_aut[0], &bijective_aut[bufSize - 1], 0);
//first value in array serves as global iterator over array
//and initially is set to base offset
bijective_aut[0] = 3;
//second value serves as indicator of array length
bijective_aut[1] = bufSize;
cl::Buffer buf(context, CL_MEM_READ_WRITE, sizeof(uint32_t) * bufSize);
cl::Buffer starts_buf(context, CL_MEM_READ_ONLY, sizeof(int) * GLOBAL_SIZE);
cl::Buffer stops_buf(context, CL_MEM_READ_ONLY, sizeof(int) * GLOBAL_SIZE);
/*====SETTING UP COMMAND QUEUE====*/
cl::CommandQueue queue(context, device);
err = queue.enqueueWriteBuffer(buf, CL_FALSE, 0, sizeof(uint32_t) * bufSize, bijective_aut);
err = queue.enqueueWriteBuffer(starts_buf, CL_FALSE, 0, sizeof(int) * GLOBAL_SIZE, starts);
err = queue.enqueueWriteBuffer(stops_buf, CL_FALSE, 0, sizeof(int) * GLOBAL_SIZE, stops);
/*====CREATING KERNEL, SETTING ITS VARIABLES====*/
cl::Kernel bc_kernel(program, "bijection_check", &err);
err = bc_kernel.setArg(0, buf);
err = bc_kernel.setArg(1, starts_buf);
err = bc_kernel.setArg(2, stops_buf);
/*====EXECUTING KERNEL====*/
cl::Event event;
err = queue.enqueueNDRangeKernel(bc_kernel, cl::NullRange, cl::NDRange(GLOBAL_SIZE), cl::NDRange(LOCAL_SIZE), nullptr, &event);
event.wait();
err = queue.enqueueReadBuffer(buf, CL_FALSE, 0, sizeof(uint32_t) * bufSize, bijective_aut);
cl::finish();
}
\define CL\u USE\u DEPRECATED\u OPENCL\u 2\u 0\u API
//有些包括在这里
#定义全局_大小4096
#定义本地_大小256
#定义全局范围0xffffffff
int main()
{
//我们假设全局_范围平均划分为全局_大小
//因此,不需要处理剩余部分
每个线程的长规则=全局范围/全局大小;
int*开始=新的int[全局大小];
int*stops=newint[GLOBAL_SIZE];
整数计数=0;
对于(int i=0;i
然后是内核代码:
__kernel void bijection_check (
__global uint * bijective_rules, //stores approved bijective rules
__global const uint * starts,
__global const uint * stops
)
{
__private int idx = get_global_id(0);
int iterator = bijective_rules[0]; //inditates next free cell to write in
int start = starts[idx];
int stop = stops[idx];
bool check = true;
//there is some variables required for test
//iterating over rules between *start* and *stop*
for (uint rule = start; rule < stop; rule++)
{
...
/*then there goes test of the rule for bijectivity*/
...
//writing current rule to general list if it turned to be bijective
if ((check == true) && (iterator < 10000))
{
bijective_rules[iterator] = rule;
bijective_rules[0]++;
}
else
{
bijective_rules[2]++;
}
}
bijective_rules[3]++;
}
\u内核无效双射\u检查(
__全局uint*双射_规则,//存储批准的双射规则
__全局consuint*开始,
__全局consuint*停止
)
{
__private int idx=get\u global\u id(0);
int迭代器=双射_规则[0];//指示下一个要写入的空闲单元格
int start=start[idx];
int-stop=stops[idx];
布尔检查=真;
//测试需要一些变量
//在*开始*和*停止*之间迭代规则*
对于(uint规则=开始;规则<停止;规则++)
{
...
/*然后对双射性规则进行了检验*/
...
//如果当前规则变为双射,则将其写入常规列表
if((check==true)&&(迭代器<10000))
{
双射_规则[迭代器]=规则;
双射_规则[0]+;
}
其他的
{
双射_规则[2]+;
}
}
双射_规则[3]+;
}
根据执行后从缓冲区读取的数组判断,最后的两条语句都没有一次执行。也就是说,在内核执行之后,双射_规则数组的状态与以前在主机端定义的完全相同。您有一个竞争条件:您读取
双射_规则[0]
,但其他线程可能同时执行双射_规则[0]++代码>,从而读取和写入该内存位置。如果两个线程将不同的数据写入同一个内存地址,则存在竞争条件,这两个线程中的哪一个将随机决定结果。因此,您的结果将是随机和不可复制的
如果多个线程需要在同一内存位置增加一个值,请使用atoimic函数。原子函数在一个线程运行时阻塞内存位置,而所有其他线程都必须等待
要消除竞争条件,请从缓冲区的一个副本(或一个特定内存地址)读取,然后写入第二个副本(或地址)。这样,您就不会写入其他并发线程正在读取的内存。谢谢您的回复!这绝对是我自己不会发现的问题。但我不确定这是否是导致程序奇怪行为的唯一原因。