ループを安全に並列化できることをコンパイラに伝える方法はありますか？

大きな配列のいくつかの要素を更新しています。ループを安全に並列化できることをコンパイラに伝える方法はありますか？

更新がからなる：（それはゼロではない場合）アレイ内の新しい位置に更新された値を移動する電流値

をクリア

10によって電流値を乗算

私は、移動が発生したときに衝突がないことを知っています。
ループを安全に並列化できることをコンパイラにどのように伝えることができますか？

do i = 1, 1e6 
    if (v[i] /= 0) then 
     temp = v[i] * 10 
     v[i] = 0 
     ndx = get_move_to_ndx(i) 
     v[ndx] = temp 
    end if 
end do

私はifortを利用していますが、これはコンパイラに依存しないと思います。

出典

2016-11-28 Anonymous

FortranとCの奇妙な混合物...削除された機能であるDOループの実際の最終値の使用については言及していません。 CONCURRENTは、ループが並列化に適しているが、PUREプロシージャにコールする必要があることをコンパイラに示す手がかりになります。 ifortに固有の何かを気にしなければ、！DIR $ IVDEPが役に立ちます。 –

上記は私にとって 'python'のようなものです。 – ja72

コンパイラが衝突がないと自動的に推測することはできません。 OpenMPを使用している場合、たとえば、共有配列を使用している場合は正常に動作します。 –

ここでは雑種のアプローチです。そのため、一時的なベクターを使用したアイデアがあります。どこが正しいかもしれない、あなたはそれを試さなければならないでしょう。 WHERE/ELSEWHEREの主な利点は、通常はループと同じくらい速くないため、読みやすさです。読みやすくなります。

!DIR$ SIMD 
FillTemp: Do I = 1, 1000000 
    Temps(I) = v(I)*10 
ENDDO FillTemp 

!$OMP PARALLEL DO 
FindIndex: Do I = 1, 1000000 
    ndx_vect(I) = get_move_to_ndx(i) 
ENDDO FindIndex 

WHERE(Temps /= 0) 
    V = 0 
ELSEWHERE 
    v(ndx_Vect) = tempz 
ENDWHERE

出典

2016-11-30 13:55:14 Holmz

ループを安全に並列化できることをコンパイラに伝える方法はありますか？

答えて

関連する問題