错误:遇到非法内存访问

Error: an illegal memory access was encountered

本文关键字:内存 访问 非法 遇到 错误      更新时间:2023-10-16

我问过疑问错误:不允许从__global__函数调用__host__函数,我得到了 ans . 因此,我修改了我的代码 BT 我无法访问d_point[i]。 我怎么能访问它。

  __global__ void densefun(int *d_counters,float2 *d_points,int d_x_max,int d_y_max,int width,int height, int min_distance,int size)
{
  int i = blockDim.x * blockIdx.x + threadIdx.x;
   if(i <= size)
   {
        float2 point = (d_points)[i];
        int x = floorf(point.x);
        int y = floorf(point.y);
        printf(" ( %d  %d )",x,y);
        if(x < d_x_max && y < d_y_max)
        {  
            x /= min_distance;
            y /= min_distance;
            (d_counters)[y*width+x]++;
            __syncthreads();
        }
    }
}

void DenseSample(const Mat& grey, std::vector<Point2f>& points, const double quality, const int min_distance)
{
    int width = grey.cols/min_distance;
    int height = grey.rows/min_distance;
    Mat eig;
    cornerMinEigenVal(grey, eig, 3, 3);
    double maxVal = 0;
    minMaxLoc(eig, 0, &maxVal);
    const double threshold = maxVal*quality;
    std::vector<int> counters(width*height);
    int x_max = min_distance*width;
    int y_max = min_distance*height; 
    printf("in descriptor size:%ld ",points.size());   
    int *d_counters;
    float2 *d_points;
    cudaMalloc(&d_counters,counters.size()*width*height*sizeof(int));
    printf("in cuda point size:%d       ",points.size());
    cudaMalloc(&d_points,points.size()*sizeof(float2));
    cout<<"points.size() : "<<points.size()<<endl;
    cudaMemcpy(d_points, &points, points.size()*sizeof(float2), cudaMemcpyHostToDevice);
    int blk=cvFloor(points.size()/1024)+1;
    cout<<"blk : "<<blk<<endl;
    if(points.size()>0)
    {
        densefun<<<blk,1024>>>(d_counters,d_points,x_max,y_max,width,height,min_distance, points.size());
        cudaError_t err = cudaGetLastError();
        if (err != cudaSuccess) 
        printf("Error: %sn", cudaGetErrorString(err));
        cudaMemcpy(&counters, d_counters, counters.size()* width*height*sizeof(int), cudaMemcpyDeviceToHost);
    }
    cudaFree(d_counters);
    cudaFree(d_points);
    points.clear();
    int index = 0;
    int offset = min_distance/2;
    for(int i = 0; i < height; i++)
    for(int j = 0; j < width; j++, index++) 
    {
        if(counters[index] <= 0)
        {
            int x = j*min_distance+offset;
            int y = i*min_distance+offset;
            if(eig.at<float>(y, x) > threshold)
            points.push_back(Point2f(float(x), float(y)));
        }
    }
}

输出为:在描述符大小:1605 在 CUDA 点大小:1605 点.大小() : 1605

大牌 : 2错误:遇到非法内存访问在描述符大小:918 在 CUDA 点大小:918 点.大小() : 918

大牌 : 1错误:遇到非法内存访问

创建块长度1024且网格长度等于

int blk=cvFloor(points.size()/1024)+1;

这基本上意味着线程数将是 1024 的倍数,大于 points.size() . 在这种情况下,请使用:

int i = blockDim.x * blockIdx.x + threadIdx.x;
float2 point = (d_points)[i];

不能成功,因为您几乎可以肯定您将获得越界内存访问。添加一些条件以确保它不会发生。

__global__ void densefun(int *d_counters,float2 *d_points,int d_x_max,int d_y_max,int width, int height, int min_distance)
{
  int i = blockDim.x * blockIdx.x + threadIdx.x;
  if(i < width * height)
  {
      //rest of the code
  }
}

此外,您不会为d_points分配 enugh 内存:

float2 *d_points;
cudaMalloc(&d_points,points.size()*sizeof(float));

如果要分配float2数组(或复制到其中),则需要使用 sizeof(float2) .