Fortran代码循环依赖分析

1）代码

    !$omp parallel private(n, j, dx, dy, dz, r, a)
    do n = 1, nsteps
        !$omp do
        do i = 0, nparticles - 1
            x_tmp(i) = x(i) + vx(i) * dt
            y_tmp(i) = y(i) + vy(i) * dt
            z_tmp(i) = z(i) + vz(i) * dt

            do j = 0, nparticles - 1
                dx = x(j) - x(i)
                dy = y(j) - y(i)
                dz = z(j) - z(i)

                R = sqrt(dx ** 2 + dy ** 2 + dz ** 2 + softening)
                a = G * m(j) / R ** 3 * dt

                vx(i) = vx(i) + a * dx
                vy(i) = vy(i) + a * dy
                vz(i) = vz(i) + a * dz
            end do
        end do
        !$omp do
        do i = 0, nparticles - 1
            x(i) = x_tmp(i)
            y(i) = y_tmp(i)
            z(i) = z_tmp(i)
        end do
        !$omp master
        print *, "n =", n
        do i = 0, nparticles - 1
            print *, i, ":", x(i), y(i), z(i)
        end do
        !$omp end master
    end do
    !$omp end parallel

2）结果：

 [8, 19, 1, None]

 [3, 20, 4, True]
OMP程序分析结果：
私有变量： ['dy', 'i', 'a', 'dz', 'n', 'dx', 'r', 'j']
归约变量： []

 [22, 26, 4, True]
OMP程序分析结果：
私有变量： ['n', 'i', 'j']
归约变量： []

 [29, 31, 4, True]
OMP程序分析结果：
私有变量： ['n', 'i', 'j']
归约变量： []

 [1, 33, -1, False]

3)分析

对于循环，变量都多了n，但是应该不影响结果的正确性。