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

楚香饭 · 发表于 2015-4-7 13:36:05

楼主的完整代码在 7 楼。

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

sqrt domian 常见于负数开方。

依旧是则徐 · 发表于 2015-4-7 14:39:10

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

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

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

依旧是则徐 · 发表于 2015-4-7 15:55:07

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

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

pasuka · 发表于 2015-4-7 19:15:49

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

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

楚香饭 · 发表于 2015-4-7 20:05:05

依旧是则徐发表于 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

依旧是则徐 · 发表于 2015-4-7 20:15:56

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

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

pasuka · 发表于 2015-4-7 20:58:31

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

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

依旧是则徐 · 发表于 2015-4-7 21:48:10

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

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

pasuka · 发表于 2015-4-7 22:14:55

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

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

[数学库] IMSL编译错误，This name has already been used as an external function name.

新人勋章

爱心勋章

水王勋章

元老勋章

热心勋章