Fortran程序中无法分类的OpenMP指令

Question

我正在尝试使用openMP在Fortran中使用以下代码并行化代码：

program pigreco
    !----------------------------------------!
    use OMP_LIB
    implicit none
    !----------------------------------------!
    integer :: i
    integer, parameter :: N = 100000
    integer, parameter :: NCPU = 4 
    real*8 :: t0, t1 
    real :: h, totale, x, f
    !----------------------------------------!
    print '(a,2x,i15)', ' Number of intervals: ', N
    totale = 0.0
    h = 1. / N

    call OMP_SET_NUM_THREADS(NCPU)
    write(*, '(a,i10)') 'Numero di processori totali: ', NCPU

    t0 = OMP_GET_WTIME()
    !----------------------------------------!
#ifdef PARALLEL
    !
    print '(a)', "Scelta la versione parallela."
    !
    !$OMP PARALLEL DO PRIVATE(x, f) REDUCTION(+:totale)
    !
    do i = 1, N 
        x = (i - 0.5) * h
        f = (4 * h) / (1 + x**2)
        totale = totale + f
    enddo
    !$OMP END PARALLEL DO
    !
#endif
    !
    t1 = OMP_GET_WTIME()
    !
    PRINT '(a,2x,f30.25)', ' Computed PI =', totale
    PRINT '(a,2x,f30.25)', ' Total computational time =', t1 - t0
    !
end program pigreco

[然后我尝试使用以下代码行进行编译：gfortran prova.F90 -fopenmp -D PARALLEL，这给了我一个错误，提示“在（1）处无法分类的OpenMP指令”。

Answer 1

问题是您使用预处理器定义了PARALLEL，因此编译器将读取OMP PARALLEL DO，而不是读取OMP 1 DO，这当然没有意义。将#ifdef PARALLEL更改为#ifdef RUNPARALLEL，将-DPARALLEL更改为-DRUNPARALLEL，则编译器未给出错误。

或者，可以使用以下事实：在使用OpenMP进行编译时，将自动定义宏变量_OPENMP，因此可以使用#ifdef _OPENMP，并且没有-D标志。