CUDA：アクセス違反の読み取り場所

Question

CUDAを初めて使用していて、場所のランタイム例外を読み取ってアクセス違反が発生したときに、.tgaファイルをぼかすための簡単なプログラムを作成しようとしていました。私はCUDAには新しかったので、私はそれを修正する方法を見つけ出すことができず、Googleは役に立たなかった。 は、ここでは、コードです：

#include "cuda_runtime.h" 
#include "device_launch_parameters.h" 
#include <iostream> 
#include <math.h> 
#include <time.h> 

#define println(...) std::cout << __VA_ARGS__ << "\n"; 
#define WIDTH 1920 
#define HEIGHT 1080 
#define BLUR 5 

unsigned char in[HEIGHT][WIDTH][3]; 
unsigned char out[HEIGHT][WIDTH][3]; 
unsigned char header[18]; 

void read(char input[256]) { 
    FILE* f; 
    f = fopen(input, "rb"); 
    if (!f) { 
     printf("File Reading Failed\n"); 
    } 
    fread(&header, 1, 18, f); 
    fread(&in, 1, HEIGHT*WIDTH * 3, f); 
    fclose(f); 
} 

void write(char output[256]) { 
    FILE* fw; 
    fw = fopen(output, "wb+"); 
    if (!fw) { 
     printf("File Writing Failed\n"); 
    } 
    header[16] = 24; 
    header[13] = WIDTH/256; 
    header[12] = WIDTH % 256; 
    header[15] = HEIGHT/256; 
    header[14] = HEIGHT % 256; 
    fwrite(&header, 1, 18, fw); 
    fwrite(&out, 1, HEIGHT*WIDTH * 3, fw); 
    fclose(fw); 
} 

__device__ 
int toIndex(int x, int y) { 
    return x + y/WIDTH; 
} 

__device__ 
void doPixel(int x, int y, char *red, char *green, char *blue) { 
    int avgRed = 0; 
    int avgGreen = 0; 
    int avgBlue = 0; 
    for (int i = -BLUR; i <= BLUR; i++) { 
     for (int j = -BLUR; j <= BLUR; j++) { 
      avgRed += red[toIndex(i, j)]; 
      avgBlue += blue[toIndex(i, j)]; 
      avgGreen += green[toIndex(i, j)]; 
     } 
    } 
    red[toIndex(x, y)] = avgRed/(BLUR*BLUR); 
    green[toIndex(x, y)] = avgGreen/(BLUR*BLUR); 
    blue[toIndex(x, y)] = avgBlue/(BLUR*BLUR); 
} 

__global__ 
void setValue(char *red, char *green, char *blue) { 
    int x; 
    int y; 
    for (int i = threadIdx.x; i < WIDTH * HEIGHT; i += 1024) { 
     x = i % WIDTH; 
     y = i/WIDTH; 
     doPixel(x, y, red, green, blue); 
    } 
} 

int main(void) { 
    char *red, *green, *blue; 
    double time; 
    read("test.tga"); 
    cudaMallocManaged(&red, WIDTH * HEIGHT); 
    cudaMallocManaged(&green, WIDTH * HEIGHT); 
    cudaMallocManaged(&blue, WIDTH * HEIGHT); 
    for (int x = 0; x < WIDTH; x++) { 
     for (int y = 0; y < HEIGHT; y++) { 
      red[x + y*WIDTH] = in[y][x][2]; 
      green[x + y*WIDTH] = in[y][x][1]; 
      blue[x + y*WIDTH] = in[y][x][0]; 
     } 
    } 
    time = clock(); 
    setValue<<<1, 1024>>>(red, green, blue); 
    cudaDeviceSynchronize(); 
    println((clock() - time)/CLOCKS_PER_SEC); 
    int x; 
    int y; 
    for (int i = 0; i < WIDTH * HEIGHT; i++) { 
     x = i % WIDTH; 
     y = i/WIDTH; 
     out[y][x][0] = blue[i];  //Program gives error here 
     out[y][x][1] = green[i]; 
     out[y][x][2] = red[i]; 
    } 
    write("test.tga"); 
    cudaFree(red); 
    cudaFree(green); 
    cudaFree(blue); 
    getchar(); 
} 

私はcudaDeviceSynchronize（）がこの問題を解決する方法だったことを読んで、それが動作していないようです。 cudaThreadSynchronize（）も問題を解決しません。

Answer 1

不正なメモリアクセスエラーを見つける最も簡単な方法は、バイナリをcuda-gdbで実行することです。あなたは、このコードセグメントでいくつかのエラーを見つけるかもしれないあなたのケースでは

をコンパイルするとき、あなたが-g -G -O0フラグを与えていることを確認してください

for (int i = -BLUR; i <= BLUR; i++) { 
    for (int j = -BLUR; j <= BLUR; j++) { 
     avgRed += red[toIndex(i, j)]; 
     avgBlue += blue[toIndex(i, j)]; 
     avgGreen += green[toIndex(i, j)]; 
    } 
}