Parallelising私はサブルーチンで次の操作を実行しようとしたとき、私は任意のスピードアップを取得しないループのFortran

Question

を行うネストされた：OpenMPの一部に

!$ call omp_set_num_threads(threadno) 
call system_clock(x1) 

!$OMP PARALLEL do private(i), reduction(+:total)  
    do i = 1,m 
     total = 0.d0 
     call result(a,l,b,qm,q,en) !here l is input for subroutine and en is output 
     qm(:,i) = q 
     qtv(i) = qt 
     mean = sum(q)/size(q) 
     do i2 = 1,k 
      total = total + ((mean-q(i2))**2)/(a+b) 
     end do 
     qvv(i1) = total 
    end do 
    call system_clock(x2) 
    print *, x2-x1 
!$OMP END PARALLEL do 

Answer 1

コメント：

1. totalがリセットされるべきでありませんループの中で、!$OMP節の前。
2. i2およびmeanはプライベートである必要があります。
3. qがループの繰り返しの間で変更されない場合は、sum(q)/size(q)を外側に配置する必要があります。

プライベート設定がないと、メモリアクセスの競合が発生する可能性があります。

あなたが表示するコードは、実際にコンパイルしたものに近いが実際にはそうではないと思う。より良いヘルプを提供するためにコンパイルされたコードを持つことは有益でしょう。

乾杯、

ピエール

EDIT：OpenMPのコードのタイミングのために、あなたは（https://gcc.gnu.org/onlinedocs/libgomp/omp_005fget_005fwtime.htmlを参照）walltime https://en.wikipedia.org/wiki/Wall-clock_timeを与えることomp_get_wtimeを使用する必要があります。 openmpルーチンのモジュールには、use omp_lib