c++ - undefined symbol 将 c++ 中的 fortran 66/77 编译为 Octave ？

Question

当我在 ising3d_fort2oct.cc 中执行 ising3d.f 时，我一直收到“undefined symbol: ising3d_”错误。任何人都可以提供有关我可以更改 ising3d.f 或 ising3d_fort2oct.cc 以解决此问题的建议吗？

Octave 终端命令:

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

Answer 1

问题不在于代码，而在于您编译和链接要从 Octave 调用的八进制文件的方式。当您构建 ising3d_fort2oct 时，您没有指示它链接或如何找到 ising3d Fortran 函数。

您可以通过以下两种方式之一将这两个源文件链接在一起。您可以单独编译每个文件，然后将生成的目标文件链接在一起，或者您可以对两个源文件调用一次 mkoctfile，让它为您做所有事情。

要使用一个命令从两个源构建 oct 文件，请执行

octave:1> mkoctfile ising3d_fort2oct.cc ising3d.f

在 Octave shell 中，mkoctfile 被告知要编译并与第一个源文件链接在一起的第二个源文件。

要像您尝试做的那样单独构建 Fortran 子例程，然后稍后将其与 oct 文件包装器链接，您可以改为做

octave:1> system ("gfortran -c ising3d.f -o ising3d.o");
octave:2> mkoctfile ising3d_fort2oct.cc ising3d.o

请注意，您仍然需要明确地将 oct 文件与目标文件链接起来。否则，将假设该函数是由 Octave 本身或系统库提供的，编译 oct 文件，这就是为什么会出现 undefined symbol 错误的原因。