为什么 ifort 不并行化这段代码?它一直说“循环未并行化:存在并行依赖”。我不明白依赖在哪里。gfortran 会生成并行代码,但速度提升不是很高。
PROGRAM erat
IMPLICIT NONE
INTEGER*8 :: i, rm, sn=1000000000
LOGICAL*1 , ALLOCATABLE, DIMENSION(:) :: nums
rm = INT( DBLE(sn)**0.5) + 1
ALLOCATE(nums(sn))
nums = .TRUE. !This line not parallelized
PRINT *, 'Doing initial sieve...'
nums(1) = .FALSE.
DO i = 2,rm
nums(i**2:sn:i) = .FALSE. !This line not parallelized
END DO
END PROGRAM erat