F2PY에서 모듈 범위 변수를 볼 수 없습니다

Question

Fortran 90과 f2py를 처음 사용하게되어 죄송합니다.

저는 Windows 64 비트, Python 3.4 64 비트, gfortran을 사용하고 있습니다. NumPy와 버전은 1.9.1이며,이 링크의 지시에 따라 나는 "NotImplementedError 인상의 gnu.py에서 ("오직 MS 컴파일러의 Win64에 gfortran와 지원 ")"의 주석 부분 : http://scientificcomputingco.blogspot.com.au/2013/02/f2py-on-64bit-windows-python27.html

나는이 로 작성 포트란의 모듈은 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 

을 그럼, 그것은 그 DP 선언되지 않은 진술 오류를보고 할 f2py -c testf2py.f90 -m testf2py

를 실행합니다.

module-scope를 function-scope에 복사하면 작동합니다. 꽤 "젖은"입니다 로

! 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?