1

以下に示すように、次のスニペットを並列化するために openmp を適用しようとしています。ただし、並列化されたコードは、openmp の派生物を使用しないシリアル化されたコードよりも実行速度が遅くなることが判明しました。私は 64 Linux プラットフォームで実行しており、gfortran でコンパイルしています。適切に行う方法についてのコメントや提案をお待ちしております。

        call omp_set_num_threads(4)

    do i = 2, natoms - 1

          rti(1:3) = R_for(i,1:3)
          fti(1:3) = ftmp(i,1:3)

    !$OMP PARALLEL DO DEFAULT(SHARED)&
    !$OMP& private(rtij,rsqij,rsqijinv,sr2,sr6,sr12,vij,wij,fij,ftij,ncut)&
    !$OMP& REDUCTION(+:vtmp,wtmp,ftmp,fti) &
    !$OMP& firstprivate(i,rti,R_for)

    do j = i + 1, natoms
        rtij = rti - R_for(j,1:3)
        rtij = rtij - boxl*idnint( rtij*boxlinv )
        rsqij = sum(rtij**2)

                if(rsqij.lt.rcutsq) then
                   rsqijinv = 1d0/rsqij
                   sr2 = sigsq*rsqijinv
                   sr6 = sr2*sr2*sr2
                   sr12 = sr6*sr6
                   vij = sr12 - sr6
                   vtmp = vtmp + vij
                   wij = vij + sr12
                   wtmp = wtmp + wij
                   fij = wij*rsqijinv

                   ftij = fij*rtij
                   fti = fti + ftij
                   ftmp(j,1:3) = ftmp(j,1:3) - ftij(1:3)
                   ncut = ncut + 1
                endif
enddo
    !$OMP END PARALLEL DO

        ftmp(i,1:3) = fti(1:3) 
enddo
4

1 に答える 1