f2py: using openMP parallel in fortran fails

血红的双手。 提交于 2019-12-08 06:52:51

问题


I am trying to compile a fortran routine that uses openMP for python using f2py.

This is the file bsp.f90:

module OTmod
  !$ use omp_lib
  implicit none

  public :: get_threads

contains

  function get_threads() result(nt)
    integer :: nt

    nt = 0
    !$ nt = omp_get_max_threads()

    !$omp parallel num_threads(nt)
        write( *, * ) 'hello world!'
    !$omp end parallel

  end function get_threads

end module OTmod

If I compile it with

f2py -m testmod --fcompiler=gfortran --f90flags='-fopenmp' -lgomp -c bsp.f90

compilation works, but importing it to python fails with the error

ImportError: dlopen(/Users/USER/omp_py/testmod/testmod.cpython-36m-darwin.so, 2): Symbol not found: _GOMP_parallel
  Referenced from: /Users/USER/omp_py/testmod/testmod.cpython-36m-darwin.so
  Expected in: flat namespace
 in /Users/USER/omp_py/testmod/testmod.cpython-36m-darwin.so

However if I remove only the two lines that make the parallel loop

!$omp parallel num_threads(nt)

and

!$omp end parallel

then importing to python and executing it works, so omp_get_max_threads() works, but the parallel loop causes it to fail. Is this some mismatch between Anaconda and my brew install gcc? How can I fix this?


System setup information below


OS

macOS 10.13

python

Python 3.6.2 |Anaconda custom (x86_64)| (default, Sep 21 2017, 18:29:43) 
[GCC 4.2.1 Compatible Clang 4.0.1 (tags/RELEASE_401/final)] on darwin

gfortran -v

Using built-in specs.
COLLECT_GCC=gfortran
COLLECT_LTO_WRAPPER=/usr/local/Cellar/gcc/7.2.0/libexec/gcc/x86_64-apple-darwin17.0.0/7.2.0/lto-wrapper
Target: x86_64-apple-darwin17.0.0
Configured with: ../configure --build=x86_64-apple-darwin17.0.0 --prefix=/usr/local/Cellar/gcc/7.2.0 --libdir=/usr/local/Cellar/gcc/7.2.0/lib/gcc/7 --enable-languages=c,c++,objc,obj-c++,fortran --program-suffix=-7 --with-gmp=/usr/local/opt/gmp --with-mpfr=/usr/local/opt/mpfr --with-mpc=/usr/local/opt/libmpc --with-isl=/usr/local/opt/isl --with-system-zlib --enable-checking=release --with-pkgversion='Homebrew GCC 7.2.0' --with-bugurl=https://github.com/Homebrew/homebrew-core/issues --disable-nls
Thread model: posix
gcc version 7.2.0 (Homebrew GCC 7.2.0) 

f2py -c --help-fcompiler

Gnu95FCompiler instance properties:
  archiver        = ['/usr/local/bin/gfortran', '-cr']
  compile_switch  = '-c'
  compiler_f77    = ['/usr/local/bin/gfortran', '-Wall', '-g', '-ffixed-
                    form', '-fno-second-underscore', '-m64', '-fPIC', '-O3', '
                    -funroll-loops']
  compiler_f90    = ['/usr/local/bin/gfortran', '-Wall', '-g', '-fno-second-
                    underscore', '-m64', '-fPIC', '-O3', '-funroll-loops']
  compiler_fix    = ['/usr/local/bin/gfortran', '-Wall', '-g', '-ffixed-
                    form', '-fno-second-underscore', '-Wall', '-g', '-fno-
                    second-underscore', '-m64', '-fPIC', '-O3', '-funroll-
                    loops']
  libraries       = ['gfortran']
  library_dirs    = ['/usr/local/Cellar/gcc/7.2.0/lib/gcc/7/gcc/x86_64-apple
                    -darwin17.0.0/7.2.0']
  linker_exe      = ['/usr/local/bin/gfortran', '-Wall', '-Wall']
  linker_so       = ['/usr/local/bin/gfortran', '-Wall', '-g', '-m64', '-
                    Wall', '-g', '-undefined', 'dynamic_lookup', '-bundle']
  object_switch   = '-o '
  ranlib          = ['/usr/local/bin/gfortran']
  version         = LooseVersion ('7.2.0')
  version_cmd     = ['/usr/local/bin/gfortran', '-dumpversion']
Fortran compilers found:
  --fcompiler=gnu95  GNU Fortran 95 compiler (7.2.0)

otool -L testmod.cpython-36m-darwin.so

/Users/USER/omp_py/testmod/testmod.cpython-36m-darwin.so:
    @rpath/libgomp.1.dylib (compatibility version 2.0.0, current version 2.0.0)
    @rpath/libgfortran.3.dylib (compatibility version 4.0.0, current version 4.0.0)
    /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 1252.0.0)
    /usr/local/lib/gcc/7/libgcc_s.1.dylib (compatibility version 1.0.0, current version 1.0.0)
    @rpath/libquadmath.0.dylib (compatibility version 1.0.0, current version 1.0.0)

回答1:


I am not sure how to solve this issue for dynamic linking. It mostly has to do with something related to LD_LIBRARY_PATH or something.

But as a workaround you can do static linking of the library

Add the -static option to the linker options in setup.py

extra_link_args=['-lgomp -static']

This will remove the need for locating the dependency at run-time



来源:https://stackoverflow.com/questions/46726276/f2py-using-openmp-parallel-in-fortran-fails

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!