Question

几乎就是在这个问题上撕掉我的头发。

我有一个CUDA内核，可以对存储在3D数组中的数据进行一些数学计算。在测试时，我曾经为数组分配了一些值（非零）并观察结果。我注释掉那些行，但结果仍然相同。就好像它完全忽略了我正在做一个memset到0的事实。

当我在Debug中逐步执行代码时代码正常工作...但不是在Release中！ 我的猜测是这个矩阵有内存泄漏。

我将此数组分配为：

cudaExtent m_extent = make_cudaExtent(sizeof(float)*matdim.x, matdim.y, matdim.z); // width, height, depth
cudaPitchedPtr m_device;
cudaMalloc3D(&m_device, m_extent);
cudaMemset3D(m_device, 0, m_extent);

我在这样的循环中调用内核：

for (int iter = 0; iter < gpu_iterations; iter++)
    {
        PF_iteration_kernel<<<grids,threads>>>(m_device, m_extent, matdim);
        cudaDeviceSynchronize(); 
    }

之后我释放了m_device投球指针：

cudaFree(m_device.ptr);

matdim只是由dim3持有的矩阵维度。

在内核中我执行以下操作（好吧，我评论了一切功能......）：

__global__ void PF_iteration_kernel(cudaPitchedPtr mPtr, cudaExtent mExt, dim3 matrix_dimensions)
{
int x = threadIdx.x + blockIdx.x * blockDim.x;
int y = threadIdx.y + blockIdx.y * blockDim.y;

// Find location within the pitched memory
char *m = (char*)mPtr.ptr;

int sof = sizeof(float);
size_t pitch = mPtr.pitch;
size_t slice_pitch = pitch*mExt.height;
char* m_addroff = m + y * pitch + x * sof;
printf("m(%d,%d) is %f \n", x, y, *m_addroff); // display the slice

*m_addroff = 0; // WILL THIS RESET IT?!

__syncthreads();
}

那应该只是显示0，但它会显示我的旧值（25,26,27,28等）。

我已经清理并重新清洁并重建了几次。我重新启动了IDE。

我的IDE是使用NSight 4.6的Visual Studio 2010（CUDA 7.0）。我在Windows 7 x64上

Answer 1

考虑一下

char* m_addroff = m + y * pitch + x * sof;
printf("m(%d,%d) is %f \n", x, y, *m_addroff);

编译器会看到char并将其推广到int推送到堆叠 - 而不是float提升为格式所需的double。

编译器不提供适合格式规范的参数，但是一些编译器会检查格式规范并警告问题。

我建议你提出这个论点。我冒险猜测和失败，但这样的事情

printf("m(%d,%d) is %f \n", x, y, *(float*)m_addroff);

Herer是一个简单的例子。

#include <stdio.h>
int main()
{
    char car [4] = {0};
    char *cptr = car;
    printf ("Hello %f\n", *(float*)cptr);
    return 0;
}

为什么我的CUDA内核返回旧值？

1 个答案: