使用 gfortran 编译 f77 内在函数时出错

Question

我正在开始使用 HPC 运行 Redhat 6.6。作为另一个程序的教程步骤，我被要求编译一些旧的 F77 代码。

我是运行行gfortran 3bar.f 我收到错误消息：

/tmp/ccgUvoGd.o: In function `MAIN__':
3bar.f:(.text+0x43): undefined reference to `mygetarg_'
3bar.f:(.text+0x61): undefined reference to `mygetarg_'
3bar.f:(.text+0x8a): undefined reference to `mygetarg_'
collect2: ld returned 1 exit status

我认为这与代码中 sqrt 内在函数的使用有关。我已经对这个错误和 gfortran 进行了一些搜索，但我对 Fortran 的特定版本为此和那个发生了一些争论。

如果有人能建议我如何让这段代码工作，我将不胜感激，因为这真的不是我正在做的事情的重点。

密码是：

      program Three_bar

      implicit none
      real*8  a1,a2, g1,g2, weight, sum,
     1        l0,d0,sqrt
      intrinsic sqrt
      character*80 inp_fname, out_fname
      character*4  card1
      integer i, ir,iw, numb_args, iargc, numb_iter

      numb_args = iargc()
      if ( numb_args .ge. 2 ) then
         call MYGETARG(1,inp_fname)
         call MYGETARG(2,out_fname)
      else if ( numb_args .eq. 1 ) then
         call MYGETARG(1,inp_fname)
         out_fname = 'rdcs_output'
      else
         inp_fname = 'rdcs_input'
         out_fname = 'rdcs_output'
      endif

      ir = 7
      open (ir, file=inp_fname,status='UNKNOWN',access='SEQUENTIAL')

      iw = 8
      open (iw, file=out_fname,status='UNKNOWN',access='SEQUENTIAL')

      l0 = 1.0
      d0 = 1.0

c loop to spend some cpu - to test recall feature

      sum = 1.0

      numb_iter = 20000
      do 1005 i=1,numb_iter
         sum = (sum*1.00001)**1.001
         if ( sum .gt. 1.0e6 ) sum = 1.0
 1005  continue

c read the data from param file

      weight = (2.d0*sqrt(2.d0)*a1 + a2)*l0*d0
      g1=(2.d0*a1+sqrt(2.d0)*a2)/(2.d0*a1*(a1+sqrt(2.d0)*a2))-1.0d0
      g2 = 1.0d0/(a1 + sqrt(2.d0)*a2) -1.0d0
      write (iw,601) g1,g2,weight

      close (unit=ir,status='KEEP')
      close (unit=iw,status='KEEP')

 501   format(a4,1x,f30.0)
 601   format("g1 = ",e25.16," ",/,"g2 = ",e25.16," ",/,
     1       "weight = ",e25.16," ")

      stop
      end

Answer 1

如评论中所述，编译器抱怨缺少 MYGETARG 的定义，而不是 SQRT。

您使用该函数的方式让我想起了 GETARG，它是一个 GNU 扩展。我的猜测是 MYGETARG 是一个包装函数，可以使代码与不同的编译器一起工作。

既然您使用的是 gfortran，我建议您直接使用 GETARG。然后，你的代码在我的机器上编译，尽管有一些警告。

在长运行上，您应该切换到符合标准的内容，例如 GET_COMMAND_ARGUMENT。

使用 gfortran 编译 f77 内在函数时出错

Error compiling f77 intrinsics with gfortran

fortran

gfortran