Форум программистов, компьютерный форум CyberForum.ru

Разбор программы - C++

Восстановить пароль Регистрация
 
Danzel
6 / 6 / 3
Регистрация: 21.08.2012
Сообщений: 45
30.08.2012, 21:20     Разбор программы #1
После недели освоения cuda и разбора программ,наткнулся на урок по куда,прилагается программа перемножения матриц,собственно код:

C++
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
// Perform "naive" square matrix multiplication
//
#include <stdio.h>
 
#define BLOCK_SIZE  16          // submatrix size
#define N           1024        // matrix size is N*N
 
 
// KERNEL //
__global__ void matMult ( float * a, float * b, int n, float * c ){
  int   bx  = blockIdx.x;       // block index
  int   by  = blockIdx.y;
  int   tx  = threadIdx.x;      // thread index
  int   ty  = threadIdx.y;
  float sum = 0.0f;         // подэлемент
  int     ia  = n * BLOCK_SIZE * by + n * ty;   // a [i][0]
  int   ib  = BLOCK_SIZE * bx + tx;
    
  // умножение матриц
  for ( int k = 0; k < n; k++ )
    sum += a [ia + k] * b [ib + k*n];
            
  // блок глобальной памяти для матрицы
  // каждый поток в один элемент
  int ic = n * BLOCK_SIZE * by + BLOCK_SIZE * bx;
    
  c [ic + n * ty + tx] = sum;
}
 
 
// HOST CODE //
int main ( int argc, char *  argv [] ){
  int   numBytes = N * N * sizeof ( float );
  
  // выделение памяти хоста
  float * a = new float [N*N];
  float * b = new float [N*N];
  float * c = new float [N*N];
 
  for ( int i = 0; i < N; i++ )
    for ( int j = 0; j < N; j++ ){
      a [i] = 0.0f;
      b [i] = 1.0f;
    }
        
  // выделение памяти устройства
  float * adev = NULL;
  float * bdev = NULL;
  float * cdev = NULL;
    
  cudaMalloc ( (void**)&adev, numBytes );
  cudaMalloc ( (void**)&bdev, numBytes );
  cudaMalloc ( (void**)&cdev, numBytes );
 
  // конфигурация запуска ядра
  dim3 threads ( BLOCK_SIZE, BLOCK_SIZE );
  dim3 blocks  ( N / threads.x, N / threads.y);
 
  // создание куда таймера
  cudaEvent_t start, stop;
  float gpuTime = 0.0f;
 
  cudaEventCreate ( &start );
  cudaEventCreate ( &stop );
    
  // asynchronously issue work to the GPU (all to stream 0)
  cudaEventRecord ( start, 0 );
  cudaMemcpy      ( adev, a, numBytes, cudaMemcpyHostToDevice );
  cudaMemcpy      ( bdev, b, numBytes, cudaMemcpyHostToDevice );
  
  matMult<<<blocks, threads>>> ( adev, bdev, N, cdev );
    
  cudaMemcpy      ( c, cdev, numBytes, cudaMemcpyDeviceToHost );
  cudaEventRecord ( stop, 0 );
 
  cudaEventSynchronize ( stop );
  cudaEventElapsedTime ( &gpuTime, start, stop );
 
  
  
 
  // время работы
  printf("time spent executing by the GPU: %.2f millseconds\n", gpuTime );
 
  // освобождение ресурсов
  cudaEventDestroy ( start );
  cudaEventDestroy ( stop  );
  cudaFree         ( adev  );
  cudaFree         ( bdev  );
  cudaFree         ( cdev  );
 
  delete a;
  delete b;
  delete c;
  system("pause");
 
  return 0;
}

Возникли два вопроса:
1)Как заполняется матрица?
2)Как вывести на экран эту матрицу?(Пробовал вывести с и sum выдает ошибку,с остальной частью вроде все понятно)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
30.08.2012, 21:20     Разбор программы
Посмотрите здесь:

Разбор кода программы C++
C++ Классы, конструктор копирования (разбор куска программы)
Разбор кода простой программы новичка C++
разбор программы C++
C++ разбор программы
C++ Разбор алгоритма программы
C++ Разбор алгоритма программы
C++ Разбор кода программы

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
AEXks
24 / 3 / 1
Регистрация: 28.10.2012
Сообщений: 35
02.11.2012, 23:34     Разбор программы #2
Danzel, еще актуально, или уже решил вопрос?
Yandex
Объявления
02.11.2012, 23:34     Разбор программы
Ответ Создать тему
Опции темы

Текущее время: 09:32. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru