我目前正在使用一本名为 Cuda by Example 的书参加并行计算课程。在本书的第 4 章中,我使用了一些包含“GL/glut.h”和“GL/gext.h”的 .h 文件,我有在线安装 GLUT 的步骤,并遵循了这些步骤。我认为这行得通,但我不确定。然后我试图找到 glext 的方向,但我似乎找不到那么多。我确实找到了一个 .h 文件,并尝试通过将其包含在 GL 文件夹中来使用它。这似乎不起作用,因为我在编译类似的东西时收到错误:
错误 1 错误:不允许从设备调用主机函数(“cuComplex::cuComplex”)/_全局_ 函数(“julia”) C:\Users\Laptop\Documents\Visual Studio 2010\Projects\Lab1\Lab1 \lab1.cu 29 1 实验室1
我认为这是因为我需要更多的 glext.h,比如 .dll 和类似于 glut 的东西,但我不确定。对此的任何帮助将不胜感激。谢谢你。
编辑: - 这是我正在使用的代码,除了前两个包含语句和 .h 文件来自谷歌代码之外,我没有改变我在书中看到的内容:谢谢你的帮助
#include "cuda_runtime.h"
#include "device_launch_parameters.h"
#include "book.h"
#include "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, .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 threadIdx/BlockIdx to pixel position
int x = blockIdx.x;
int y = blockIdx.y;
int offset = x + y * gridDim.x;
//now claculate 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;
}
int main( void ) {
CPUBitmap bitmap(DIM, DIM);
unsigned char *dev_bitmap;
HANDLE_ERROR(cudaMalloc((void**)&dev_bitmap, bitmap.image_size()));
dim3 grid(DIM,DIM);
kernel<<<grid,1>>>( dev_bitmap );
HANDLE_ERROR( cudaMemcpy( bitmap.get_ptr(), dev_bitmap, bitmap.image_size(), cudaMemcpyDeviceToHost));
bitmap.display_and_exit();
HANDLE_ERROR( cudaFree( dev_bitmap ));
}