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

pasuka

依旧是则徐 发表于 2015-4-7 21:48
师兄，我现在是刚学用fortran，这个程序也是我自己写的，我也知道有很多错误问题。你说的matlab 验证程序 ...

个么改行来得及，matlab对照帮助文档的例子或者网上找些资料，三周就能入门，搞不定的问题就上网搜索，相关资料可谓铺天盖地，不少高校还购买正版了
Fortran的话，没有C、C++的基础，怎么也得花上三个月到半年的时间入门

依旧是则徐

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

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

pasuka

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

说实话，lz真不如直接matlab写个验证程序，调试起来也不会这么折磨人

依旧是则徐

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

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

楚香饭

依旧是则徐 发表于 2015-4-7 14:39
越界问题 是把variable模块中的sk3（3,0:2,0:22）改成了sk3（3,0：22,0:5） 改过之后，在运行就出现了sqr ...

sqrt domain 不一定就是负数，有时候本身就是 NaN 也会出错。
至少在我这里，发现很多容易出负数的地方。

例如 130.f90 文件
YITABOU3(i)=Amp*DCOS(WK*X3(i)-W1*T)&
                  -Amp**2*WK/(2*DSINH(2*WK*H)) &
这里的 WK 为 0，导致分母为0

wk 应该是在更早的代码计算：
wk=newton(a,func,dfunc)
这里计算结果为0

可能是由此导致的 NaN

正如 pasuka 所言，没有注释的代码比较难看懂。我想你学习一下如何 Debug 单步运行，或许会让你自己更容易检查错误。

本站的文章：http://debug.w.fcode.cn 会指导你如何学习 Debug

pasuka

依旧是则徐 发表于 2015-4-7 15:55
谢谢来帮我解决问题，http://pan.baidu.com/s/1o6tA04M这是链接。如果有空暇，请帮我看看问题在哪？ ...

粗看代码，没啥注释，也没发现啥问题，由于未装imsl所以没法编译测试
若是被导师强迫写Fortran程序的话，趁早和导师商量改matlab吧，相比Fortran往往事半功倍

依旧是则徐

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

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

依旧是则徐

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

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

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

楚香饭

楼主的完整代码在 7 楼。

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

sqrt domian 常见于负数开方。

pasuka

依旧是则徐 发表于 2015-4-7 12:54
是的，我改过来了！谢谢！然后我再请教一个问题，我在运行时出现error M6201： MATH     ——sqrt： DOMA ...

1、请贴出完整的代码，这样挤牙膏的话，就怕下次再发帖就风干喽
2、sqrt默认输入正实数，请检查以下输入值