2

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

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 并行此代码？

（更多信息：）

在这篇文章之后，我首先尝试获取上三角距离矩阵值。这是内核代码：

__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( xx 等于 (1,2,3,4,5,6,7,8))。然后我应该有n(n-1)/2值（例如，总共 28 个）。我将全局和本地工作大小设置为n(8)。和are和(28)的cl_mem缓冲区对象的大小。但是当我读回结果时，我得到了 28 个值，但不是预期的值。我得到：（0,7,12,15,16,15,12,7,0,0,0...,0）。有什么帮助吗？xxautonn(n-1)/2

4

0 回答 0