f2py的问题:对子程序中模块的未定义引用

时间:2012-07-30 10:02:53

标签: python fortran f2py

我正在尝试使用f2py编译Fortran f90文件,以便在Python中使用。该文件是从另一个文件调用模块的子例程。该模块基本上用于分配。我可以在命令窗口中使用'gfortran my_dec.f90'编译模块,但在尝试编译子例程文件时出错。这特别困难,因为我几乎没有使用Fortran,这是别人的代码。

这是模块和子程序的一部分,因为它很长,包括它的开始和结束:

 module my_dec

    integer ndir, nfreq
    integer ihmax,ier
    integer nk,nth,nspec


    real hspmin
    real wsmult
    real wscut

    logical flcomb, flc 

    parameter(ndir=24)
    parameter(nfreq=23)

    parameter(nk=nfreq)
    parameter(nth=ndir)
    parameter(nspec=nk*nth)

    REAL DTH, SIG(0:nk+1), DSII(0:nk+1), DSIP(0:nk+1)
    REAL ECOS(nspec+nth), ESIN(nspec+nth), XFR 
    REAL FACHFE, TH(nth), FTE

    REAL ES2(nspec+NTH),EC2(nspec+NTH),ESC(nspec+NTH)
    REAL DDEN(NK),DDEN2(nspec)
    REAL SIG2(nspec)

    INTEGER IAPROC, NAPERR, NDSE, NDST

    INTEGER year, TIME 

    real pcg ! percentage either side of peakfor gamma estimate

    data pcg/0.3/ 

end module my_dec

子程序:

子程序my_init

use my_dec
use constants

iaproc=1
naperr=1
ndset=1
ndst=1

IHM    = 100
HSPM   = 0.05
WSM    = 1.7
WSC    = 0.333
FLC    = .true.


IHMAX  = MAX ( 50, IHM )
HSPMIN = MAX ( 0.0001 , HSPM )
WSMULT = MAX ( 1. , WSM )
WSCUT  = MIN ( 1.0001 , MAX ( 0. , WSC ) )
FLCOMB = FLC

...      

    return

    end

当我尝试使用'f2py -c my_init.f90 -m my_init_m'编译子程序文件'my_init.f90'时,我在子程序中获得了大量关于模块中变量引用的消息:

c:\users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o:my_init.f90:(.text+0
xb): undefined reference to `__my_dec_MOD_iaproc'
c:\users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o:my_init.f90:(.text+0
x15): undefined reference to `__my_dec_MOD_naperr'
c:\users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o:my_init.f90:(.text+0
x26): undefined reference to `__my_dec_MOD_ndst'
c:\users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o:my_init.f90:(.text+0
x4f): undefined reference to `__my_dec_MOD_flc'

然后是错误,这对我来说并不多见:

collect2: ld returned 1 exit status
error: Command "C:\Python27\Scripts\gfortran.exe -Wall -Wall -shared c:\users\lw
l\appdata\local\temp\tmptlve6z\Release\users\lwl\appdata\local\temp\tmptlve6z\sr
c.win-amd64-2.7\my_init_mmodule.o c:\users\lwl\appdata\local\temp\tmptlve6z\Rele
ase\users\lwl\appdata\local\temp\tmptlve6z\src.win-amd64-2.7\fortranobject.o c:\
users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o -Lc:\python27\egg-info\
mingw\usr\lib\gcc\x86_64-w64-mingw32\4.5.2 -LC:\Python27\libs -LC:\Python27\PCbu
ild\amd64 -lpython27 -lgfortran -o .\my_init_m.pyd" failed with exit status 1

我一直试图解决这个问题,包括搜索互联网,但无济于事。有人有主意吗?这可能是一个非常简单的问题。谢谢你的帮助。

编辑:如果我将模块复制并粘贴到与子例程相同的文件中,我就可以使用它,但是将它作为单独的文件使用它会很好。

1 个答案:

答案 0 :(得分:2)

对不起,如果这个解释涵盖了你已经知道的理由,但是你写的是你几乎不了解Fortran。

您的例程my_init使用名为my_dec的模块(以及一个名为constants的模块)。这就是use陈述的含义。错误消息,例如

c:\users\lwl\appdata\local\temp\tmptlve6z\Release\my_init.o:my_init.f90:(.text+0
xb): undefined reference to `__my_dec_MOD_iaproc'
如果你试图编译my_init而没有提供链接到my_dec的编译版本,那么我期望看到

。诸如__my_dec_MOD_iaproc之类的名称由编译器生成,您可以将该名称读取为在iaproc MODule中标识名为my_dec的实体。如果my_dec根本没有定义iaproc,您会收到类似的消息,但事实并非如此。

暂且离开f2py你(在大多数情况下)只需确保首先编译由另一个模块或子程序或程序used的任何模块,链接器将发挥其魔力(提供路径)设置等等。

我不知道你如何告诉f2py在哪里寻找my_dec的编译版本。

我看到你现在已经在你的编辑中找到了解决方案。我很困惑为什么你认为将子程序的源文件放在一个单独的文件中会很好。如果你真的想写Fortran 90子程序也属于模块。