OpenMP：縮小と共有

Question

次のコードがあり、誤った出力になっています。

#pragma omp parallel private(i,piold) shared(pi,sign) 
{ 
#pragma omp for reduction(+:pi) schedule (static) 
    for (i = 0; i < 100000; i++){ 
     piold = pi; 
     pi += sign/(2*i+1); 
     sign=-sign; 
    } 
} 
    pi = 4*pi; 

私はOpenMPの新機能として少し失われています。私が混乱しているのは、スレッド間でsign,pi、pioldを渡す方法です。別の反復には異なる値が必要なので、それらを共有することはできません。したがって、1つの方法は、繰り返しを奇数と偶数に分けることですが、それは非常に非効率的です。

この場合のご提案はありますか？

Answer 1

円の近似にはLeiniz formulaを使用してみてください。

• C++でスコープの反復が宣言されているため、プライベートからiを削除します。用

  （unsigned int型i = 0; I < 10000; I ++）

• 符号はIの奇数/偶数プロパティ上のベースを知っているので、（iは& 1）の代わりにチェック使用してから符号を削除することができますプライベート

  int sign = 1 - (i & 1); 
      pi += (sign == 1 ? 1f : -1f)/(2 * i + 1); 
  

• pioldは、最終結果のパイに影響を与えるあなたはを読むことができるプライベート

から削除されることはありません詳しくはをご覧ください。私はコンパイラがOpenMPをサポートしていないので、私はテストすることができず、サンプルコードを与えることができません。しかし、私は2つのケースを提案します：

• openmpセクションを使用してiの奇数/偶数で2つのスレッドに式を分けます（+）と（ - ）。 2つのセクションがあります。

double plus = 0.0;

for (int i = 0; i < 5000; i++) 
     plus += (4 * i + 1); 

、次いで

double minus = 0.0; 
for (int i = 0; i < 5000; i++) 
    minus += (4 * i + 3); 

pi = plus - minus; 

• 使用i`は整数変数と `の両方である` `仮定sign`およびPI