- 2008-10-27 (Mon) 15:56
- 一般
OpenMPのreductionで||とか使ったら評価はどうなるんだ？ ||でのリダクションって並列化するとまずいんじゃないだろうか？
やってみた．関数fは引数がtrueならtrueを返して即座に抜け，それ以外は重たい処理をしてからfalseを返す．配列bsはbs[0]のみtrueで他はfalse．これで，逐次実行時にはi=0でrにtrueが入り以降はfの処理をしないので即座に計算が終わる．一方，並列実行時にi=0を含まない反復範囲の計算が行われるとしたら，その部分の計算では重たいfの処理が大量に行われ，計算時間がかさむはず．
#include<iostream> #include<vector> bool f(bool b) { unsigned int u = 0; if(!b) { for(int k = 0; k < 100000; k++) { u = (u + 2) % 7; } u = u ^ (((u/4)<<2) | (u&3)); } return b ^ (u!=0); } int k = 2; int main(int argc, char *argv[]) { int n = 10000000; std::vector<bool> bs(n); for(int i = 0; i < n; i++) { bs[i] = i == 0; } bool r = false; #pragma omp parallel for reduction(||:r) for(int i = 0; i < n; i++) { r = r || f(bs[i]); } int res = r; std::cout << res << std::endl; }
結果：逐次は一瞬，並列だと遅すぎて時間計るのも馬鹿らしい．
ということで，||でのリダクションを並列化するのはあまり意味がない気がする．さらに，本来評価されないものが評価される可能性が出てくるので，セグフォって落ちるプログラムが出来上がるかもしれない．なぜにOpenMPは||もリダクションに使えるようにしているのだろう？
- Newer: ことはじめ