CUDA プログラミング - asyncAPI 学習記録

CUDAプログラミング - asyncAPI学習記録

2024-07-12

非同期API

1. 完全なコード
2. 重要部分の分析

1. 完全なコード

#include <stdio.h>

// includes CUDA Runtime
#include <cuda_runtime.h>
#include <cuda_profiler_api.h>

// includes, project
#include <helper_cuda.h>
#include <helper_functions.h>  // helper utility functions

__global__ void increment_kernel(int *g_data, int inc_value) {
  int idx = blockIdx.x * blockDim.x + threadIdx.x;
  g_data[idx] = g_data[idx] + inc_value;
}

bool correct_output(int *data, const int n, const int x) {
  for (int i = 0; i < n; i++)
    if (data[i] != x) {
      printf("Error! data[%d] = %d, ref = %dn", i, data[i], x);
      return false;
    }

  return true;
}

int main(int argc, char *argv[]) {
  int devID;
  cudaDeviceProp deviceProps;

  printf("[%s] - Starting...n", argv[0]);

  // This will pick the best possible CUDA capable device
  devID = findCudaDevice(argc, (const char **)argv);

  // get device name
  checkCudaErrors(cudaGetDeviceProperties(&deviceProps, devID));
  printf("CUDA device [%s]n", deviceProps.name);

  int n = 16 * 1024 * 1024;
  int nbytes = n * sizeof(int);
  int value = 26;

  // allocate host memory
  int *a = 0;
  checkCudaErrors(cudaMallocHost((void **)&a, nbytes));
  memset(a, 0, nbytes);

  // allocate device memory
  int *d_a = 0;
  checkCudaErrors(cudaMalloc((void **)&d_a, nbytes));
  checkCudaErrors(cudaMemset(d_a, 255, nbytes));

  // set kernel launch configuration
  dim3 threads = dim3(512, 1);
  dim3 blocks = dim3(n / threads.x, 1);

  // create cuda event handles
  cudaEvent_t start, stop;
  checkCudaErrors(cudaEventCreate(&start));
  checkCudaErrors(cudaEventCreate(&stop));

  StopWatchInterface *timer = NULL;
  sdkCreateTimer(&timer);
  sdkResetTimer(&timer);

  checkCudaErrors(cudaDeviceSynchronize());
  float gpu_time = 0.0f;

  // asynchronously issue work to the GPU (all to stream 0)
  checkCudaErrors(cudaProfilerStart());
  sdkStartTimer(&timer);
  cudaEventRecord(start, 0);
  cudaMemcpyAsync(d_a, a, nbytes, cudaMemcpyHostToDevice, 0);
  increment_kernel<<<blocks, threads, 0, 0>>>(d_a, value);
  cudaMemcpyAsync(a, d_a, nbytes, cudaMemcpyDeviceToHost, 0);
  cudaEventRecord(stop, 0);
  sdkStopTimer(&timer);
  checkCudaErrors(cudaProfilerStop());

  // have CPU do some work while waiting for stage 1 to finish
  unsigned long int counter = 0;

  while (cudaEventQuery(stop) == cudaErrorNotReady) {
    counter++;
  }

  checkCudaErrors(cudaEventElapsedTime(&gpu_time, start, stop));

  // print the cpu and gpu times
  printf("time spent executing by the GPU: %.2fn", gpu_time);
  printf("time spent by CPU in CUDA calls: %.2fn", sdkGetTimerValue(&timer));
  printf("CPU executed %lu iterations while waiting for GPU to finishn",
         counter);

  // check the output for correctness
  bool bFinalResults = correct_output(a, n, value);

  // release resources
  checkCudaErrors(cudaEventDestroy(start));
  checkCudaErrors(cudaEventDestroy(stop));
  checkCudaErrors(cudaFreeHost(a));
  checkCudaErrors(cudaFree(d_a));

  exit(bFinalResults ? EXIT_SUCCESS : EXIT_FAILURE);
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105

2. 重要部分の分析

装置初期化:findCudaDevice 関数は、最適な CUDA デバイスを選択し、デバイス ID を返すために使用されます。

devID = findCudaDevice(argc, (const char **)argv);

1
2

デバイスプロパティの取得: cudaGetDeviceProperties 関数は、デバイス名やその他の情報を含む、指定されたデバイスのプロパティを取得します。

checkCudaErrors(cudaGetDeviceProperties(&deviceProps, devID));

1
2

メモリ割り当て: CPU 上でアクセス可能なページロックメモリを割り当てるには cudaMallocHost を使用し、デバイス上にメモリを割り当てるには cudaMalloc を使用します。

int *a = 0;
checkCudaErrors(cudaMallocHost((void **)&a, nbytes));

1
2
3

スレッドブロックとグリッドを設定します。ここでは、スレッドブロックサイズが 512 スレッドに設定され、グリッドサイズはデータサイズに基づいて動的に計算されます。

dim3 threads = dim3(512, 1);
dim3 blocks = dim3(n / threads.x, 1);

1
2
3

CUDA イベントとタイマーの作成: CUDA イベントは時間の記録に使用され、タイマーは CPU 実行時間の測定に使用されます。

cudaEvent_t start, stop;
checkCudaErrors(cudaEventCreate(&start));
checkCudaErrors(cudaEventCreate(&stop));

1
2
3
4

CUDA ストリーム処理: 非同期メモリコピーに cudaMemcpyAsync を使用、<<<blocks, threads> >>構文は、同時に実行されるCUDAカーネル関数increment_kernelを開始します。

cudaMemcpyAsync(d_a, a, nbytes, cudaMemcpyHostToDevice, 0);
increment_kernel<<<blocks, threads, 0, 0>>>(d_a, value);
cudaMemcpyAsync(a, d_a, nbytes, cudaMemcpyDeviceToHost, 0);

1
2
3
4

タイミングと待機: cudaEventRecord はイベントを記録し、GPU 実行時間を計算するために使用されます。 cudaEventQuery(stop) を通じて GPU 操作が完了するまで待ちます。

cudaEventRecord(start, 0);
// ...
cudaEventRecord(stop, 0);

1
2
3
4

結果の検証: GPU 計算結果の正確性を検証するには、correct_output 関数を使用します。

bool bFinalResults = correct_output(a, n, value);

1
2

リソース解放: 割り当てられたメモリと CUDA イベントを解放します。

checkCudaErrors(cudaEventDestroy(start));
checkCudaErrors(cudaEventDestroy(stop));
checkCudaErrors(cudaFreeHost(a));
checkCudaErrors(cudaFree(d_a));

1
2
3
4
5

CUDA カーネル関数 increment_kernel:

この単純な CUDA カーネル関数は、配列内の各要素を指定された値 inc_value だけインクリメントするために使用されます。 blockIdx.x と threadIdx.x は、各スレッドのグローバルインデックス idx を計算し、加算演算を実行するために使用されます。

__global__ void increment_kernel(int *g_data, int inc_value) {
  int idx = blockIdx.x * blockDim.x + threadIdx.x;
  g_data[idx] = g_data[idx] + inc_value;
}

1
2
3
4
5

その他の補助機能
CudaErrorsをチェック：CUDAをチェックする関数呼び出しにエラーがあるかどうか。
sdkCreateTimer および sdkResetTimer: タイマーの作成とリセットに使用されます。
sdkStartTimer および sdkStopTimer: タイマーの開始と停止、および CPU 実行時間を記録するために使用されます。

技術共有