Question

我无法为此序列代码创建内核：

void distlag(double *x, double *y,double *data, double *ncoord, double *res)
{
    // x: x-coordinate
    // y: y-coordinate
    // data: value associated with x-y coordinate
    // ncoord: number of coordinates
    int i=0, j=0;
    double u=0.0, v=0.0, lags=0.0;
    for(i=0;i<(ncoord[0]-1);i++)
        for(j=(i+1); j<ncoord[0];j++)
        {
            // Pairwise distances
            lags=hypot(x[i]-x[j],y[i]-y[j]);
            u=data[i]; 
            v=data[j]; 
            *res+= u+v+lags;
        } 
    }

我认为主要问题是在hypot中，距离取决于同一向量中的值。或者，j值可能取决于i。如何使用OpenCL并行此代码？

（更多信息：）

在this post之后，我首先尝试获得上三角距离矩阵值。这是内核代码：

__kernel void totexp(__global float *x, __global float *xauto, int n)
{
int num_wrk_items  = get_local_size(0);
int local_id       = get_local_id(0);
int group_id       = get_group_id(0);
int j, gid = get_global_id(0);
float sum = 0.0;


for (j = 0; j < n; j++) 
{
if ( ((gid+j)!= j) && ((gid+j) < n)) 
       {
   sum = fabs(x[j]-x[gid+j]);
        }
        else
            continue;
  }
  xauto[gid] = sum;
}

int n是length的{{1}}（x等于（1,2,3,4,5,6,7,8））。然后我应该有x个值（示例中总共28个）。我将全局和本地工作大小设置为n(n-1)/2（8）。 n和cl_mem的{{1}}缓冲区对象的大小分别为x和xauto（28）。但当我读回结果时，我得到28个值但不是预期值。我得到：（0,7,12,15,16,15,12,7,0,0,0 ......，0）。有什么帮助吗？

OpenCL：并行化滞后向量运算（三角距离矩阵）

0 个答案: