F2PY не может видеть переменные области видимости

Question

Извините, что вы новичок как в Fortran 90, так и в f2py.

Я использую Windows 64 бит, Python 3.4 64 бит, gfortran. Numpy версия 1.9.1, и я откомментировал «Поднятие NotImplementedError (» Только MS компилятор поддерживал с gfortran на Win64 «)» в gnu.py, как указано на этой ссылке: http://scientificcomputingco.blogspot.com.au/2013/02/f2py-on-64bit-windows-python27.html

У меня есть модуль Fortran, записывается следующим образом, с модулем-Scope переменной dp:

! testf2py.f90 
module testf2py 
    implicit none 
    private 
    public dp, i1 
    integer, parameter :: dp=kind(0.d0) 
contains 
    real(dp) function i1(m) 
     real(dp), intent(in) :: m(3, 3) 
     i1 = m(1, 1) + m(2, 2) + m(3, 3) 
     return 
    end function i1 
end module testf2py 

Тогда, если я бегу f2py -c testf2py.f90 -m testf2py

Он сообщит об ошибке, о том, что дп не была объявлена.

Если я копирую область модуля в область функций, она будет работать.

! testf2py.f90 
module testf2py 
    implicit none 
    private 
    public i1 
    integer, parameter :: dp=kind(0.d0) 
contains 
    real(dp) function i1(m) 
     integer, parameter :: dp=kind(0.d0) 
     real(dp), intent(in) :: m(3, 3) 
     i1 = m(1, 1) + m(2, 2) + m(3, 3) 
     return 
    end function i1 
end module testf2py 

Однако, это не выглядит как лучшая практика кодирования, хотя, как это довольно «мокрая».

Любые идеи?

Answer 1

Вот обходные, в котором dp перемещается в types модуля, а use types заявление добавляется к функции i1.

! testf2py.f90 

module types 
    implicit none 
    integer, parameter :: dp=kind(0.d0) 
end module types 

module testf2py 
    implicit none 
    private 
    public i1 
contains 
    real(dp) function i1(m) 
     use types 
     real(dp), intent(in) :: m(3, 3) 
     i1 = m(1, 1) + m(2, 2) + m(3, 3) 
     return 
    end function i1 
end module testf2py 

В действии:

In [6]: import numpy as np 

In [7]: m = np.array([[10, 20, 30], [40, 50, 60], [70, 80, 90]]) 

In [8]: import testf2py 

In [9]: testf2py.testf2py.i1(m) 
Out[9]: 150.0 

Изменение похож на третий вариант, который я описал в этом ответе: f2py: Specifying real precision in fortran when interfacing with python?