亚洲av成人无遮挡网站在线观看,少妇性bbb搡bbb爽爽爽,亚洲av日韩精品久久久久久,兔费看少妇性l交大片免费,无码少妇一区二区三区

Chinaunix

標題: 【原創(chuàng)】跟劉峰六學C語言(5) 線程柵欄 [打印本頁]

作者: sohu2000000    時間: 2010-10-02 16:56
標題: 【原創(chuàng)】跟劉峰六學C語言(5) 線程柵欄
某些并行計算需要面臨某些在計算進行前的某些單通瓶頸點,這種情況下,當然可以使用信號量的方式來進行處理,但是還存在著另外的一種處理方式是更加方便的,它就是:柵欄(在pthread庫里面被定義成為類型 pthread_barrier_t),下面我們來看一段程序作為示例
  1. #define _XOPEN_SOURCE 600

  2. #include <pthread.h>
  3. #include <stdlib,h>
  4. #include <stdio.h>

  5. #define ROWS        10000
  6. #define COLS        10000
  7. #define THREADS        10

  8. double initial_matrix[ROWS][COLS];
  9. double final_matrix[ROWS][COLS];

  10. //Barrier variable

  11. pthread_barrier_t barr

  12. extern void DotProduct(int row, int col, double source[ROWS][COLS], double destination[ROWS][COLS]);
  13. extern void determinant(double matrix[ROWS][COLS]);

  14. void * entry_point(void * arg)
  15. {
  16.     int rank = (int)arg;
  17.     int row;
  18.     for(row=rank*ROWS/THREADS; row < (rank+1)*THREADS;++row)
  19.         for(int col=0;col<COLS;++col)
  20.             DotProduct(row,col,initial_matrix,final_matrix);

  21.     //synchronization pointer

  22.     int rc = pthread_barrier_wait(&barr);
  23.     if(rc!=0 && rc != PTHREAD_BARRIER_SERIAL_THREAD)
  24.     {        
  25.         printf("Could not wait on barrier\n");
  26.         exit(-1);
  27.     }

  28.     for(row=rank*ROWS/THREADS; row < (rank+1)*THREADS;++row)
  29.         for(int col=0;col<COLS;++col)
  30.             DotProduct(row,col,final_matrix,initial_matrix);
  31. }

  32. int main(int argc, char* argv[])
  33. {
  34.     int i;
  35.     pthread_t thr[THREADS];

  36.     //Barrier initialization

  37.     if(pthread_barrierattr_init(&barr,NULL,THREADS))
  38.     {
  39.         
  40.         printf("Could not create a barrier\n");
  41.         return -1;
  42.     }

  43.     for(i=0;i<THREADS;++i)
  44.     {
  45.         if(pthread_create(&thr[i],NULL,&entry_point, (void*)i))
  46.         {
  47.             
  48.             printf("Could not create thread %d\n", i);
  49.             return -1;
  50.         }
  51.     }

  52.     for(i=0;i<THREADS;++i)
  53.     {
  54.         if(pthread_join(thr[i],NULL))
  55.         {
  56.             
  57.             printf("Could not join thread %d\n", i);
  58.             return -1;
  59.         }
  60.     }

  61.     double det = determinant(initial_matrix);
  62.     printf("The determinant of M^4 = %f\n", det);

  63.     return 0;
  64.         
  65. }
復制代碼
這段程序產生出許多個線程,并且分配給每個線程計算矩陣乘法的一部分,然后每個線程使用這次計算的結果,繼續(xù)進行下一步的計算:另一個矩陣的乘法

幾點關于API的說明:

barrier 變量必須在最開始聲名為全局變量
barrier 變量的初始化必須在main函數(shù)里進行初始化
在點上每一個線程都會等待它的對端完成工作
注意
在程序頂部的宏定義 _XOPEN_SOURCE 是非常重要的;如果沒有這個變量,那么barrier類型就會在pthread.h中被屏蔽掉,這個定義必須在所有的頭文件引用之前被定義出來
作者: sohu2000000    時間: 2010-10-02 16:58





歡迎光臨 Chinaunix (http://72891.cn/) Powered by Discuz! X3.2