OpenMP* のリダクション操作

OpenMP* のリダクションは、共有数値変数のインクリメントや配列の共有数値変数への累積など、単純な操作に使用できます。リダクション操作を行うには、reduction 節を並列領域内に追加して、指定した操作と変数で並列に合計操作を行うようにコンパイラーに指示します。

次のようなアノテート付きの C/C++ シリアルコードについて考えてみます。

  int  i, n=500000; 
  float *array, total=0.0;
     ...
     for (i=0; i <n ; ++i { 
        ANNOTATE_LOCK_ACQUIRE(0);
        total+ = array[i];
        ANNOTATE_LOCK_RELEASE(0);
      }
...

次に、#include <omp.h> と #pragma omp parallel for simd を追加した C/C++ 並列コードを示します。

#include <omp.h> // .dll を検出できないロード時の問題を解決します 
     int i, n=500000; 
     float *array, total=0.0; 
     ...
     #pragma omp parallel for reduction(+:total)
     for (i=0; i < 100; i++) { 
          total+ = array[i];
        }
. . .

次のアノテーション付きの Fortran シリアルコードについて考えてみます。

  integer(4) n
  real(4) array(50000), total = 0.0
  n = 500000
  ...
     do i=1, n
     call annotate_lock_acquire(0)
        total = total + array(i)
     call annotate_lock_release(0)
 .. .
     end do

次に、use omp_lib、!$omp parallel do reduction(+:ret) および !$omp end parallel do を追加した C/C++ 並列コードを示します。

  use omp_lib
  integer(4) n
  real(4) array(50000), total = 0.0
  n = 500000
  ...
  !$omp parallel do reduction(+:total)
     do i=1, n
        total = total + array(i)
  !$omp end parallel do
 .. .
     end do