Question

我有CERN程序来计算复杂参数的gamma函数，但由于缺乏Fortran知识，我无法运行该程序。

我有以下程序：

PROGRAM Console1
  IMPLICIT REAL *8 (A-H,O-Z) 
  COMPLEX *16 gama,z,w 
  z=cmplx(0,0.707106781186548d0) 
  gama=cgamma(0,z,w) 
END

SUBROUTINE cgamma(mo, z, w) 
  INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(14, 60)
  INTEGER, INTENT(IN)       :: mo
  COMPLEX (dp), INTENT(IN)  :: z
  COMPLEX (dp), INTENT(OUT) :: w
... the subroutine code here 
END SUBROUTINE cgamma

然后出现错误

错误＃6553：函数引用正在调用子例程子程序。 [CGAMMA]

Answer 1

从技术上讲，它不是一个程序。我将在下面的链接中重新发布部分文档，请注意SUBROUTINE cgamma和END SUBROUTINE cgamma作为第一个和非常最后行。从技术上讲，这个文件不是fortran程序，它只是一个子程序，它是更大程序的一部分。

SUBROUTINE cgamma(mo, z, w)

   IMPLICIT NONE

!  variables and code not worth posting to illustrate what is happening

   60 w = CMPLX(w1, w2, KIND=dp)
   RETURN

   70 w = (0.0_dp, 0.0_dp)
   RETURN

   CONTAINS

   FUNCTION rexp(x) RESULT(fn_val)
!     extra code here, do not need to post it
      fn_val = -1.0_dp
      RETURN
   END FUNCTION rexp

END SUBROUTINE cgamma

在fortran中，文件的第一行和最后一行是PROGRAM whatever和END。

因此您需要获取完整的程序并使用调用此cgamma子例程的程序，或编写您自己的fortran PROGRAM，然后调用此cgamma子例程。

Answer 2

查看错误消息：

函数引用正在调用子例程子程序。

SUBROUTINE cgamma(mo, z, w)
END SUBROUTINE

将cgamma定义为（在本例中为外部）子例程。

在程序中

gama=cgamma(0,z,w)

引用cgamma作为函数（返回要分配给gama的结果）。这是错误消息所指的不兼容性。

要解决此问题，您需要使用call语句引用子例程，或将cgamma更改为函数，并返回适当的结果。

鉴于程序的设计（三个参数），我希望该功能是有意的。

函数引用正在调用子例程

2 个答案: