Я недавно начал изучать CUDA и интегрировал CUDA в MS Visual Studio 2010 с помощью Nsight. Я также приобрел книгу «CUDA by Example», и я изучаю все примеры и собираю их. Однако я столкнулся с ошибкой, которую не понимаю.
Программа взята из главы 4 и является примером julia_gpu. Оригинальный код:
#include "../common/book.h"
#include "../common/cpu_bitmap.h"
#define DIM 1000
struct cuComplex {
float r;
float i;
cuComplex( float a, float b ) : r(a), i(b) {}
__device__ float magnitude2( void ) {
return r * r + i * i;
}
__device__ cuComplex operator*(const cuComplex& a) {
return cuComplex(r*a.r - i*a.i, i*a.r + r*a.i);
}
__device__ cuComplex operator+(const cuComplex& a) {
return cuComplex(r+a.r, i+a.i);
}
};
__device__ int julia( int x, int y ) {
const float scale = 1.5;
float jx = scale * (float)(DIM/2 - x)/(DIM/2);
float jy = scale * (float)(DIM/2 - y)/(DIM/2);
cuComplex c(-0.8, 0.156);
cuComplex a(jx, jy);
int i = 0;
for (i=0; i<200; i++) {
a = a * a + c;
if (a.magnitude2() > 1000)
return 0;
}
return 1;
}
__global__ void kernel( unsigned char *ptr ) {
// map from blockIdx to pixel position
int x = blockIdx.x;
int y = blockIdx.y;
int offset = x + y * gridDim.x;
// now calculate the value at that position
int juliaValue = julia( x, y );
ptr[offset*4 + 0] = 255 * juliaValue;
ptr[offset*4 + 1] = 0;
ptr[offset*4 + 2] = 0;
ptr[offset*4 + 3] = 255;
}
// globals needed by the update routine
struct DataBlock {
unsigned char *dev_bitmap;
};
int main( void ) {
DataBlock data;
CPUBitmap bitmap( DIM, DIM, &data );
unsigned char *dev_bitmap;
HANDLE_ERROR( cudaMalloc( (void**)&dev_bitmap, bitmap.image_size() ) );
data.dev_bitmap = dev_bitmap;
dim3 grid(DIM,DIM);
kernel<<<grid,1>>>( dev_bitmap );
HANDLE_ERROR( cudaMemcpy( bitmap.get_ptr(), dev_bitmap,
bitmap.image_size(),
cudaMemcpyDeviceToHost ) );
HANDLE_ERROR( cudaFree( dev_bitmap ) );
bitmap.display_and_exit();
}
Моя Visual Studio, однако, вынуждает меня приукрашивать конструктор cuComplex на устройство , иначе он не скомпилируется (он говорит мне, что я не могу использовать его позже в функции julia), что, я думаю, достаточно справедливо , Итак, у меня есть:
__device__ cuComplex( float a, float b ) : r(a), i(b) {}
Но когда я запускаю пример (добавив необходимые включения для запуска через VS, то есть cuda_runtime.h и device_launch_parameters.h, а также скопировав glut32.dll в ту же папку, что и exe), он быстро выходит из строя, убивая мой драйвер устройства и говоря, что это происходит из-за неизвестной ошибки в строке 94, которая является вызовом cudaMemcpy в main. Чтобы быть точным, это фактическая строка, содержащая вызов "cudaDeviceToHost". Честно говоря, я попытался создать несколько точек останова строка за строкой, и драйвер умирает при вызове ядра.
Может кто-нибудь сказать мне, что может быть не так? Я новичок в CUDA и понятия не имею, почему тривиальный пример так себя убивает. Что я могу делать не так? Потому что, честно говоря, я даже не знаю, что расследовать.
У меня есть набор инструментов CUDA 4.1, NSight 2.1 и GeForce GT445M с вычислительными возможностями, равными 2.1, и версия драйверов 295.