在Fortran上使用openmp和fftw

Question

我目前正試圖在Fortran上運行帶有OpenMP的fftw，但是我在運行任何程序時遇到了一些問題。

我相信我已經安裝/配置FFTW正確：

./configure --enable-openmp --enable-threads 

，我似乎有所有正確的庫和文件，但我不能讓任何程序運行，我不斷收到錯誤

undefined reference to 'fftw_init_threads' 

我使用的代碼如下：

program trial 
    use omp_lib 
    implicit none 
    include "fftw3.f" 
    integer :: id, nthreads, void 
    integer :: error 

    call fftw_init_threads(void) 

    !$omp parallel private(id) 
    id = omp_get_thread_num() 
    write (*,*) 'Hello World from thread', id 
    !$omp barrier 

    if (id == 0) then 
     nthreads = omp_get_num_threads() 
     write (*,*) 'There are', nthreads, 'threads' 
    end if 

    !$omp end parallel 
    end program 

並運行它，我做

gfortran trial.f90 -I/home/files/include -L/home/files/lib -lfftw3_omp -lfftw3 -lm -fopenmp 

這將不勝感激，如果任何人都可以幫助我。

Answer 1

舊的FORTRAN界面似乎不支持OpenMP ...我建議你採用新的Fortran 2003界面。請注意，fftw_init_threads()是一個功能！

您還需要包括ISO_C_binding模塊：

program trial 
    use,intrinsic :: ISO_C_binding 
    use omp_lib 
    implicit none 
    include "fftw3.f03" 
    integer :: id, nthreads, void 
    integer :: error 

    void = fftw_init_threads() 

    !$omp parallel private(id) 
    id = omp_get_thread_num() 
    write (*,*) 'Hello World from thread', id 
    !$omp barrier 

    if (id == 0) then 
     nthreads = omp_get_num_threads() 
     write (*,*) 'There are', nthreads, 'threads' 
    end if 

    !$omp end parallel 
    end program