我正在尝试使用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
我一直试图解决这个问题,包括搜索互联网,但无济于事。有人有主意吗?这可能是一个非常简单的问题。谢谢你的帮助。
编辑:如果我将模块复制并粘贴到与子例程相同的文件中,我就可以使用它,但是将它作为单独的文件使用它会很好。
答案 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子程序也属于模块。