我正在使用Fortran 90编写代码，现在需要使用amos Fortran 77库中的特殊函数(http://www.netlib.org/amos/)。我找到了这些例程的模块接口(https://github.com/certik/fortran-utils/blob/master/src/amos.f90)。
我的问题是：如何将它们结合起来并在我的Fortran 90程序中使用，并如何正确编译它们？
我已经为此苦苦挣扎了一整天，仍然无法弄清楚。
以下是我的测试代码：
我得到的结果如下：
我的问题是：如何将它们结合起来并在我的Fortran 90程序中使用，并如何正确编译它们？
我已经为此苦苦挣扎了一整天，仍然无法弄清楚。
以下是我的测试代码：
PROGRAM TEST_ZBESI
USE set_precisions
USE amos
IMPLICIT NONE
INTEGER :: n, i, nz, ierr
!double precision :: zr,zi, cyr(5), cyi(5)
REAL(kind=DBL) :: zr, zi, cyr(5), cyi(5)
n=5
zr=1.0_DBL
zi=2.0_DBL
call ZBESI(zr,zi,0.0_DBL,1,n,cyr,cyi,nz,ierr)
print *,' '
do i=1, n
write(*,10) i-1, cyr(i)
write(*,11) i-1, cyi(i)
end do
print *,' NZ=', NZ
print *,' Error code:', ierr
print *,' '
10 format(' zr(',I1,') = ',F10.6)
11 format(' zi(',I1,') = ',F10.6)
END PROGRAM TEST_ZBESI
我得到的结果如下：
zr(0) = 0.000000
zi(0) = 0.000000
zr(1) = 0.000000
zi(1) = 0.000000
zr(2) = 0.000000
zi(2) = 0.000000
zr(3) = 0.000000
zi(3) = 0.000000
zr(4) = 0.000000
zi(4) = 0.000000
NZ= 0
Error code: 4
无论如何我都似乎无法得到正确的答案。
我尝试手动将ZBESI.f Fortran 77代码转换为Fortran 90代码。但是代码太长了，这是一场灾难。