arrays - 英特尔 Cilk Plus 数组符号的奇怪复制行为

Question

我正在使用 Intel Cilk Plus 数组表示法来练习矢量编程。但是，我遇到了数组赋值的奇怪复制行为。

要解决的问题是并行前缀。D是输入，P是输出。

//错误的结果代码。在英特尔 icc 14.0.2 上测试。

void vec_prefix(int n, int D[n], int P[n]) {
  P[0:n] = D[0:n]; //initial copy
  int bound=1;
  int * Buf1 = _mm_malloc(n*sizeof(int), 16);
  int * Buf2 = _mm_malloc(n*sizeof(int), 16);

  while(bound<n){
    Buf1[0:n-bound] = P[0:n-bound];
    Buf2[0:n-bound] = P[bound:n-bound];
    //printf("WHY??\n");  //Add this fence, the result will be correct.
    P[bound:n-bound] = Buf1[0:n-bound] + Buf2[0:n-bound];
    bound<<=1;
  }
  _mm_free(Buf1); _mm_free(Buf2);
}

如果我删除“printf”行的注释，结果是正确的。否则就错了。但如果所有副本都遵循英特尔文档中的描述，则代码应该是正确的。

好像没有内存栅栏之类的，前两行的Buf1/Buf2拷贝没有完成，后面的add操作使用了一些不稳定的值。或者，编译器优化器只是使用复制传播来删除副本，并创建“P[bound:n-bound] = P[0:n-bound] + P[bound:n-bound]”。这在英特尔的文档中未定义。

//正确的结果代码

void vec_prefix(int n, int D[n], int P[n]) {
  P[0:n] = D[0:n]; //initial copy
  int bound=1;
  int * Buf1 = _mm_malloc(n*sizeof(int), 16);

  while(bound<n){
    //direct copy part
    Buf1[0:bound] = P[0:bound];
    //add part
    Buf1[bound:n-bound] = P[bound:n-bound] + P[0:n-bound];
    //copy back
    P[0:n] = Buf1[0:n];
    bound<<=1;
  }
  _mm_free(Buf1);
}

参考：英特尔的文件

• http://software.intel.com/sites/products/documentation/studio/composer/en-us/2011Update/compiler_c/optaps/common/optaps_par_cean_prog.htm
• https://www.cilkplus.org/tutorial-array-notation

Answer 1

你的代码对我来说看起来是正确的。不需要“围栏”。我向英特尔编译器组提交了一份内部错误报告，其中包含所有示例作为输入。谢谢你的例子。

如果编译器工作正常，您甚至可以将 while 循环缩短为：

while(bound<n){
    Buf1[0:n-bound] = P[0:n-bound];
    P[bound:n-bound] = Buf1[0:n-bound] + P[bound:n-bound];
    bound<<=1;
}

第二个数组部分分配是可以的，因为 P 的重叠是完全正确的。唉 icc 在这个例子中也展示了它的错误。