我正在尝试使用 C++ 中的 OpenACC 并行化 Karatsuba 算法的迭代版本。我想问一下如何矢量化内部 for 循环
。我的编译器显示了关于该循环的消息:
526, Complex loop carried dependence of result-> prevents parallelization
Loop carried dependence of result-> prevents parallelization
Loop carried backward dependence of result-> prevents vectorization
这里是两个嵌套循环的代码:
#pragma acc kernels num_gangs(1024) num_workers(32) copy (result[0:2*size-1]) copyin(A[0:size],$
{
#pragma acc loop gang
for (TYPE position = 1; position < 2 * (size - 1); position++) {
// for even coefficient add Di/2
if (position % 2 == 0)
result[position] += D[position / 2];
TYPE start = (position >= size) ? (position % size ) + 1 : 0;
TYPE end = (position + 1) / 2;
// inner loop: sum (Dst) - sum (Ds + Dt) where s+t=i
#pragma acc loop worker
for(TYPE inner = start; inner < end; inner++){
result[position] += (A[inner] + A[position - inner]) * (B[inner] + B[position - inn$
result[position] -= (D[inner] + D[position - inner]);
}
}
}
实际上,我不确定是否可以对其进行矢量化。但如果是,我无法意识到我做错了什么。谢谢
最佳答案
“复杂循环携带结果依赖性”问题是由指针别名引起的。编译器无法判断“result”指向的对象是否与另一个指针的对象重叠。
作为 C++ 扩展,您可以将 C99“restrict”关键字添加到数组声明中。这将向编译器断言指针不会别名。
或者,您可以在循环指令中添加 OpenACC“独立”子句,以告诉编译器循环没有任何依赖性。
请注意,OpenACC 不支持数组缩减,因此您将无法并行化“内部”循环,除非您修改代码以使用标量。像这样的东西:
rtmp = result[position];
#pragma acc loop vector reduction(+:rtmp)
for(TYPE inner = start; inner < end; inner++){
rtmp += (A[inner] + A[position - inner]) * (B[inner] + B[position - inn$
rtmp -= (D[inner] + D[position - inner]);
}
result[position] = rtmp;
关于c++ - 在 C++ 中使用 OpenACC 并行化和矢量化迭代 Karatsuba 算法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/49817402/