[英]undefined symbol compiling fortran 66/77 in c++ to octave?
當我在ising3d_fort2oct.cc內部執行ising3d.f時,我不斷收到“未定義符號:ising3d_”錯誤。 誰能提供關於我可以更改ising3d.f或ising3d_fort2oct.cc來解決此問題的建議?
八度終端命令:
octave:12> system("gfortran ising3d.f -o ising3d")
ans = 0
octave:13> mkoctfile ising3d_fort2oct.cc
octave:14> ising3d_fort2oct
error: /home/b/Desktop/Umaine_Classes/CHY_573/3disingmodelf/ising3d_fort2oct.oct: failed to load: /home/b/Desktop/Umaine_Classes/CHY_573/3disingmodelf/ising3d_fort2oct.oct: undefined symbol: ising3d_
octave:14>
ising3d_fort2oct.cc:動態鏈接的C ++函數,調用fortran子例程ising3d.f:
#include <octave/oct.h>
#include <octave/f77-fcn.h>
extern "C" void F77_FUNC (ising3d,ISING3D)();
DEFUN_DLD (ising3d_fort2oct, args , ,"work in progress")
{
octave_value_list retval;
F77_FUNC(ising3d,ISING3D)();
return retval;
}
ising3d.f:確實很老的fortran子例程:
SUBROUTINE ising3d()
C 3D ISING MODEL
C Critical temperature TC = 4.5116=1/.22165
DIMENSION IS(10,10,10), EX(13)
REAL *8 R(1)
DATA IS/1000*1/
ITMAX=5000
ISTART=4000
L=10
NR=1
ISEED=768521034
M=L*L*L
DO 1000 K=1,24
TR=0.05+(K-1)*.05
T=TR/.221655
MR=0.
DO 3 I=1, 13, 2
3 EX(I)=EXP(-2*(I-7.)/T)
DO 2 ITIME=1,ITMAX
DO 1 K1=1,L
K1P1=K1+1
K1M1=K1-1
IF(K1.EQ.1) K1M1=L
IF(K1.EQ.L) K1P1=1
DO 1 K2=1,L
K2P1=K2+1
K2M1=K2-1
IF(K2.EQ.1) K2M1=L
IF(K2.EQ.L) K2P1=1
DO 1 K3=1,L
K3P1=K3+1
K3M1=K3-1
IF(K3.EQ.1) K3M1=L
IF(K3.EQ.L) K3P1=1
IEN=7+IS(K3,K2,K1)*(IS(K3M1,K2,K1)+IS(K3P1,K2,K1)+IS(K3,K2M1,K1)
& +IS(K3,K2P1,K1)+IS(K3,K2,K1M1)+IS(K3,K2,K1P1))
CALL GGUBS(ISEED,NR,R)
IF(EX(IEN).LT.R(1)) GOTO 1
IS(K3,K2,K1)=-IS(K3,K2,K1)
M=M+2*IS(K3,K2,K1)
C WRITE(*,*) M,ITIME
1 CONTINUE
IF (ITIME.GT.ISTART) MR=MR+M
2 CONTINUE
WRITE(*,*) FLOAT(MR)/1000./FLOAT((ITMAX-ISTART)),TR
WRITE(1,*) FLOAT(MR)/1000./FLOAT((ITMAX-ISTART)),TR
1000 CONTINUE
STOP
END
RETURN
END
SUBROUTINE GGUBS(ISEED,NR,R)
IMPLICIT REAL *8(A-H,O-Z)
DIMENSION R(NR)
DATA D2P31M/2147483647.D0/
DATA D2P31/2147483648.D0/
DO 7 I=1,NR
ISEED=MOD(16807.D0*ISEED, D2P31M)
R(I)=ISEED/D2P31
7 CONTINUE
RETURN
END
問題不在於代碼,而在於您編譯和鏈接要從Octave調用的oct文件的方式。 構建ising3d_fort2oct
,沒有指示它鏈接或如何找到ising3d
Fortran函數。
您可以通過以下兩種方式之一將這兩個源文件鏈接在一起。 您可以單獨編譯每個文件,然后將生成的目標文件鏈接在一起,或者可以對兩個源文件調用一次mkoctfile
,然后讓它為您做所有事情。
要使用一個命令從兩個源構建oct文件,請執行
octave:1> mkoctfile ising3d_fort2oct.cc ising3d.f
在Octave Shell中,告訴mkoctfile
第二個要編譯的源文件並與第一個源文件鏈接在一起。
如要嘗試分別構建Fortran子例程,然后再將其與oct-file包裝器鏈接,則可以執行以下操作:
octave:1> system ("gfortran -c ising3d.f -o ising3d.o");
octave:2> mkoctfile ising3d_fort2oct.cc ising3d.o
請注意,您仍然需要將oct文件與目標文件明確鏈接。 否則,假定函數由Octave本身或系統庫提供,則將編譯oct文件,這就是為什么會出現未定義的符號錯誤的原因。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.