Fortran Coder

标题: sqrt 定义域出错的问题 [打印本页]

作者: wansq 时间: 2014-3-28 12:42
标题: sqrt 定义域出错的问题
还是这段老程序，出了好几次错了，出一次改一次，又才出一次。。。运行出来，是开方的问题，但是不知道具体是哪里，怎么改啊，求助求助！

QQ截图20140328124341.png (11.04 KB, 下载次数: 329)

11.7 KB, 下载次数: 1

70 Bytes, 下载次数: 2

28 Bytes, 下载次数: 1

166 Bytes, 下载次数: 1

作者: 楚香饭 时间: 2014-3-28 13:15
检查一下 144，145 行。

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

PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*RO1(J)))

      CSI=2*ASIN((SQRT(XXXI-XXX(J))**2+(YYYI-YYY(J))**2)/(2*RO2(J)))

你确定是要对 XXI-XX(j) 先开根号，再平方？？
而不是平方和，然后开根号？

作者: wansq 时间: 2014-3-28 13:26

chuxf 发表于 2014-3-28 13:15
检查一下 144，145 行。
[mw_shl_code=fortran,true]PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*R ...

再次谢谢你，我改过来了

作者: wansq 时间: 2014-3-28 13:29

chuxf 发表于 2014-3-28 13:15
检查一下 144，145 行。
[mw_shl_code=fortran,true]PSI=2*ASIN((SQRT(XXI-XX(J))**2+(YYI-YY(J))**2)/(2*R ...

请问一下，您是经过调试，查找到这一行的是吧

作者: 楚香饭 时间: 2014-3-28 13:33

wansq 发表于 2014-3-28 13:29
请问一下，您是经过调试，查找到这一行的是吧

是的。直接调试，出错后就自动停在这一行

作者: wansq 时间: 2014-3-28 13:56

chuxf 发表于 2014-3-28 13:33
是的。直接调试，出错后就自动停在这一行

哦，明白了，如果下次出问题，我试着调试一下，就知道问题出在哪里了

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)