トップQs
タイムライン
チャット
視点

OpenMP

ウィキペディアから

OpenMP
Remove ads

OpenMP（オープンエムピー）は、並列計算機環境において共有メモリ・マルチスレッド型の並列アプリケーションソフトウェア開発をサポートするために標準化されたAPIである[3]。「OpenMP」は「open multiprocessing」の略である[3]。

概要 開発元, 初版 ...

同様に並列コンピューティングに利用されるMPIでは、メッセージの交換をプログラム中に明示的に記述しなければならないが、OpenMPではディレクティブ（指令）を挿入することによって並列化を行う。OpenMPが使用できない環境では、このディレクティブは無視されるため、並列環境と非並列環境でほぼ同一のソースコードを使用できるという利点がある。また、プラットフォーム固有のスレッドAPIを使わず、コンパイラによって暗黙的に生成されたスレッド[注釈 1]を利用してタスクを振り分けることになるため、並列プログラムを簡潔に記述できるだけでなく、複数の環境に移植しやすくなる。

MPIとの比較では、OpenMPは異なるスレッドが同一のデータを同じアドレスで参照できるのに対して、MPIでは明示的にメッセージ交換を行わなければならない。そのため、OpenMPは、SMP環境においては大きなデータの移動を行なわずにすむので高い効率が期待できる。ただし並列化の効率はコンパイラに依存するので、チューニングによる性能改善がMPIほど高くならないという問題がある。また、OpenMPはMPIに比べてメモリアクセスのローカリティが低くなる傾向があるので、頻繁なメモリアクセスがあるプログラムでは、MPIの方が高速な場合が多い[要出典]。

OpenMPは、並列プログラミングにおいて最も広く利用されているAPIであるが、共有メモリに対してUniform Memory Access（英語版）に近いアクセスができるハードウェアシステムアーキテクチャでは、スケーラビリティに限界がある[3]。そのため、現在のほとんどのスーパーコンピューターでは、OpenMP単独ではなく、分散メモリ環境で高いスケーラビリティを発揮するMPIと組み合わせた、ハイブリッドMPI+OpenMPが利用されている[3][6]。

2013年にリリースされたOpenMP 4.0では、多数の先進的な機能が追加された[7]。SIMD命令を使った自動ベクトル化機能（omp simd）や、GPUなどのアクセラレータに並列処理を委譲する分散メモリ型のオフロード機能などがサポートされている[8][9]。

2023年現在、FORTRANとC/C++について標準化が行われている。

Remove ads

OpenMPを用いたコード例

要約
視点

以下はC言語における for ループを並列処理させる例である。

int main(int argc, char *argv[])
{
    int i;
#pragma omp parallel for
    for (i = 0; i < 10000; ++i)
    {
        /* (並列処理させたいプログラム) */
    }
    return 0;
}

OpenMPはループの反復処理を自動的に複数のスレッドに分割して並行処理できるようにする。例えば4つのスレッドを用いて処理を分割する場合、上記例ではインデックス[0, 2499], [2500, 4999], [5000, 7499], [7500, 9999]の各範囲をそれぞれのスレッドに分担させる、といった具合である。実際にいくつのスレッドを起動するのか、また各スレッドに対してどのように処理を振り分けるのかはOpenMP処理系（コンパイラ）およびプログラム実行環境などの条件に依存する[10]。

以下は区分求積法を用いた円周率πの数値計算を、OpenMP並列リダクションを用いて行うC++のコード例である。一部にC++11の機能が使われているが、OpenMPのディレクティブ自体は言語バージョンとは無関係であり、C++98/C++03でも利用できる。

#include <iostream>
#include <chrono>
#include <cmath>
#include <iomanip>
#include <omp.h>

const double D_PI = 3.1415926535897932384626433832795;

// 区分求積法で π の近似値を求める。
// 1 / (x^2 + 1) を区間 [0, 1] で積分すると π/4 になるという定積分を利用する。

int main()
{
  const int DivNum = 1000 * 1000 * 1000;
  const double delta = 1.0 / DivNum;

  std::cout << "OpenMP max threads count = " << omp_get_max_threads() << std::endl;

  const auto startTime = std::chrono::system_clock::now();
  double sum = 0;
#pragma omp parallel for reduction(+ : sum)
  for (int i = 0; i < DivNum; ++i)
  {
    const double x = (delta * i);
    const double area = delta *  1.0 / (x * x + 1.0);
    sum += area;
  }
  const double pi = sum * 4.0;
  const auto endTime = std::chrono::system_clock::now();
  std::cout << std::setprecision(15) << "PI ~= " << pi << std::endl;
  std::cout << "Error [%] = " << (100.0 * std::fabs(D_PI - pi) / D_PI) << std::endl;
  std::cout << "Elapsed time [ms] = " << std::chrono::duration_cast<std::chrono::milliseconds>(endTime - startTime).count() << std::endl;
  return 0;
}

#pragma omp parallel forは並列ループのディレクティブである。直後に続くreductionはclause[注釈 2]と呼ばれ、並列処理の動作を調整することができる。ここでは総和を格納するスレッド共有変数sumに対する更新操作の演算子（加算）を指定している。異なるスレッドから共有変数にアクセスするときは排他制御やアトミック操作が必要となるが、OpenMPのclauseを使用することでそのような煩雑なコードを記述する必要がなくなり、詳細を処理系に任せて隠蔽することができる。

OpenMPコンパイルオプションの有無を切り替えるか、OpenMPディレクティブをコメントアウト／コメント解除してからコンパイル・実行することで、マルチスレッド版およびシングルスレッド版の速度性能比較を簡単に行なうことができるのがOpenMPプログラムの特徴である。

Remove ads

対応コンパイラ

  • GCC：バージョン4.1から対応[15]。バージョン4.9でC/C++用のOpenMP 4.0を、バージョン4.9.1でFortran用のOpenMP 4.0をサポートした[16]。GCC 5ではオフロード機能のサポートが追加された。GCC 6以降でC/C++用のOpenMP 4.5を、GCC 7以降でFortran用のOpenMP 4.5をサポートしている。GCC 9以降でOpenMP 5.0の初期サポートが始まっている。
  • Clang：バージョン3.7でOpenMP 3.1に対応した[17]。Clang 3.7以前は派生プロジェクトが存在した[18]。Clang 3.9でオフロード以外のOpenMP 4.5機能をすべてサポートした[19]。
  • Microsoft Visual C++：Visual C++ 2017時点でOpenMP 2.0をサポートしている[20]。Visual C++ 2019ではSIMDベクトル化機能を実験的にサポートする[21][22]。
  • Intel C++ Compiler：バージョン12.1においてOpenMP 3.1をサポートしている。また、バージョン14.0においてOpenMP 4.0の機能を一部サポートしている[23]。
  • Intel Fortran Compiler（英語版）: バージョン18.0以降でOpenMP 5.0の機能の大部分をサポートしている[24]。
Remove ads

歴史

  • 1997年10月 - OpenMP Fortran API 1.0[25][1]
  • 1998年10月 - OpenMP C/C++ API 1.0[25][1]
  • 1999年11月 - OpenMP Fortran API 1.1[1]
  • 2000年11月 - OpenMP Fortran API 2.0[1]
  • 2002年3月 - OpenMP C/C++ API 2.0[1]
  • 2005年5月 - OpenMP Fortran C/C++ API Version 2.5[1]
  • 2007年10月 - OpenMP Fortran C/C++ API Version 3.0 Draft[1]
  • 2008年 - OpenMP 3.0[26]
  • 2011年7月21日 - OpenMP 3.1[26]
  • 2012年11月13日 - OpenMP 4.0 Draft[27]
  • 2013年7月26日 - OpenMP 4.0[28]
  • 2015年11月16日 - OpenMP 4.5[29]
  • 2016年11月10日 - OpenMP 5.0 Preview 1[30]
  • 2018年11月8日 - OpenMP 5.0[31]
  • 2020年11月13日 - OpenMP 5.1[32]
  • 2021年11月9日 - OpenMP 5.2[33]
  • 2024年8月1日 - OpenMP 6.0 Draft[34]
  • 2024年11月14日 - OpenMP 6.0[2]

関連書籍

要約
視点

洋書：

  • Rohit Chandra, Ramesh Menon, Leo Dagum, David Kohr, Dror Maydan and Jeff McDonald: Parallel Programming in OpenMP, Morgan Kaufmann, ISBN 978-1-55860-671-5 (2000年10月).
  • Barbara Chapman, Gabriele Jost and Ruud Van der Pas: Using OpenMP: Portable Shared Memory Parallel Programming, MIT Press, ISBN 978-0-262-53302-7 (2007年10月).
  • Ruud van der Pas, Eric Stotzer and Christian Terboven: Using OpenMP -- The Next Step: Affinity, Accelerators, Tasking, and SIMD, The MIT Press, ISBN 978-0-262-53478-9 (2017年10月27日). ※ OpenMP 4.5 の仕様を記述。
  • Timothy G. Mattson, Yun (Helen) He, and Alice Evelyn Konigs: The OpenMP Common Core: Making OpenMP Simple Again, The MIT Press, ISBN 978-0-262-53886-2 (2019年11月19日).
  • Tom Deakin and Timothy G. Mattson: Programming Your GPU with OpenMP: Performance Portability for GPUs, The MIT Press, ISBN 978-0-262-54753-6 (2023年11月7日).

和書等：

Remove ads

脚注

関連項目

外部リンク

Loading related searches...

Wikiwand - on

Seamless Wikipedia browsing. On steroids.

Remove ads
Remove ads