IMSL编译错误，This name has already been used as an external function name.

依旧是则徐

出现这样的错误，修改思路是啥
Error: The same named entity from different modules and/or program units cannot be referenced.   [SXYZ]
   SXYZ(1,0,j)=-15+1*J
---^
Error: This name has already been used as an external function name.   [SXYZ]
   SXYZ(1,0,j)=-15+1*J
---^
Error: This function, which is specified as the left side of an assignment statement, is invalid.   [SXYZ]
   SXYZ(1,0,j)=-15+1*J
这是程序的一部分。

[Fortran] 纯文本查看 复制代码

!主程序
program MAIN
use     NUMERICAL
use     CONSTANT
USE     VARIABLE
USE     IMSL

IMPLICIT NONE
external FCN
                                      !!!主程序出现的变量        
real :: a=20.0                          !牛顿迭代求波数时的初始值        
INTEGER IDO ,ISTEP,NOUT,PSTEP2
real*8  TOL,PARAM(MXPARM),Y(N),TEND,T
INTEGER  I,J, K
  !浮子浮标坐标划分
    zxyz(1,0)=-15
    zxyz(2,0)=0.0
    zxyz(3,0)=0.0

    zxyz(1,22)=15
    zxyz(2,22)=0.0
    zxyz(3,22)=0.0

    do i=1,21
       zxyz(1,i)=-10+1*(i-1)
       zxyz(2,i)=0
       zxyz(3,i)=0
    end do

     ! 梗绳单元划分
    
   do J=1,4
   SXYZ(1,0,j)=-15+1*J
   SXYZ(2,0,j)=0
   SXYZ(3,0,j)=0
   SXYZ(1,0,5)=zxyz(1,1)
   SXYZ(1,0,0)=zxyz(1,0)
   SXYZ(2,0,5)=zxyz(2,1)
   SXYZ(2,0,0)=zxyz(2,0)
   SXYZ(3,0,5)=zxyz(3,1)
   SXYZ(3,0,0)=zxyz(3,0)
end do
   
   do I=1,20 
   do J=1,4
      sxyz(1,i,j)=-10+1*(I-1)+0.2*J
      sxyz(2,i,j)=0
      sxyz(3,i,j)=0
      sxyz(1,i,5)=zxyz(1,i+1)
      sxyz(1,i,0)=zxyz(1,i)
      sxyz(2,i,5)=zxyz(2,i+1)
      sxyz(2,i,0)=zxyz(2,i)
      sxyz(3,i,5)=zxyz(3,i+1)
      sxyz(3,i,0)=zxyz(3,i)
   end do
   end do

I=21
do J=1,4
   sxyz(1,i,j)=-10+1*(I-1)+1*J
   sxyz(2,i,j)=0
   sxyz(3,i,j)=0
   sxyz(1,i,5)=zxyz(1,i+1)
   sxyz(1,i,0)=zxyz(1,i)
   sxyz(2,i,5)=zxyz(2,i+1)
   sxyz(2,i,0)=zxyz(2,i)
   sxyz(3,i,5)=zxyz(3,i+1)
   sxyz(3,i,0)=zxyz(3,i)
end do

依旧是则徐

楚香饭 发表于 2015-4-7 09:26
use     NUMERICAL
use     CONSTANT
USE     VARIABLE

对啊，我就是在variable 里面定义了sxyz  其他的都没有定义。而且对于imsl 我不太懂，当我把IMSL去掉之后时，上述error就没有了 ，反倒出现了这个错误error LNK2001: unresolved external symbol _IVPRK@32
Debug/130.exe : fatal error LNK1120: 1 unresolved externals

依旧是则徐

楚香饭 发表于 2015-4-7 10:28
试试
use IMSL
改为

谢谢，确实有帮助，这个问题解决了！然后又出现了这个问题，请帮忙看下d:\files\ipv\130.f90(634) : Error: The type of the actual argument differs from the type of the dummy argument.   [T]
call IVPRK(IDO, N, FCN, T, TEND, TOL, PARAM, Y)
------------------------^
d:\files\ipv\130.f90(634) : Error: The type of the actual argument differs from the type of the dummy argument.   [TEND]
call IVPRK(IDO, N, FCN, T, TEND, TOL, PARAM, Y)
我看了下我的主程序和子程序里面都定义这个量了（下面贴出来了），为啥还出现这样的问题？

!主程序
program MAIN
use     NUMERICAL
use     CONSTANT
USE     VARIABLE
use     IMSL,sxyz_1=>sxyz

IMPLICIT NONE
external FCN
                                      !!!主程序出现的变量       
real :: a=20.0                          !牛顿迭代求波数时的初始值       
INTEGER IDO ,ISTEP,NOUT,PSTEP2
real*8  TOL,PARAM(MXPARM),Y(N),TEND,T
INTEGER  I,J, K

！子程序
SUBROUTINE FCN(NN,T1,Y,YPRIME)
use CONSTANT
USE VARIABLE
IMPLICIT NONE
INTEGER NN
REAL*8 T1,Y(N),YPRIME(N)

依旧是则徐

楚香饭 发表于 2015-4-7 11:09
参数类型不一致。ivprk 用于处理单精度，而 T 是双精度。
尝试使用 DIVPRK 函数 ...

非常感谢！编译没问题了.我在运行时又出错了, 提示我的问题是数组fz超出。可是我不太理解，我设定的FUFSUM(3，0：22)没有问题吧？请楼主如果有空暇时，帮我看看下面的文件（链接是文件所在位置），http://pan.baidu.com/s/1o6tA04M

依旧是则徐

楚香饭 发表于 2015-4-7 12:32
VARIABLE 模块中，
real*8  sk3(3,0:2,0:22)
而在  shouli.f90 第 50 行：

是的，我改过来了！谢谢！然后我再请教一个问题，我在运行时出现error M6201： MATH     ——sqrt： DOMAIN error。我看了下 shouli.f90中的sqrt公式，没发现错误，帮我检查一下吧！

依旧是则徐

楚香饭 发表于 2015-4-7 13:36
楼主的完整代码在 7 楼。

但是我不知道他之前的越界怎么改的，所以也不知道问题在哪儿。

越界问题 是把variable模块中的sk3（3,0:2,0:22）改成了sk3（3,0：22,0:5） 改过之后，在运行就出现了sqrt domain，我检查了一下，sqrt里面的值都是平方项的和，不应该出现负值啊，请楚香饭师兄帮忙看看。看能否发现错误，

依旧是则徐

pasuka 发表于 2015-4-7 13:22
1、请贴出完整的代码，这样挤牙膏的话，就怕下次再发帖就风干喽
2、sqrt默认输入正实数，请检查以下输入 ...

谢谢来帮我解决问题，http://pan.baidu.com/s/1o6tA04M这是链接。如果有空暇，请帮我看看问题在哪？

依旧是则徐

楚香饭 发表于 2015-4-7 20:05
sqrt domain 不一定就是负数，有时候本身就是 NaN 也会出错。
至少在我这里，发现很多容易出负数的地方。 ...

确实是有很多NaN错误，谢谢师兄！我这就看看Debug单步运行的文章。有问题我再来请教！

依旧是则徐

pasuka 发表于 2015-4-7 20:58
说实话，lz真不如直接matlab写个验证程序，调试起来也不会这么折磨人

师兄，我现在是刚学用fortran，这个程序也是我自己写的，我也知道有很多错误问题。你说的matlab 验证程序我还不会呢。正想办法学，