fortran运行的时候出现sqrt:DOMAIN error问题

irides · 发表于 2016-4-30 18:17:03

您需要登录才可以下载或查看，没有帐号？立即注册

x

就是在t检验的这一部分出了问题

!t检验

do ix=1,x

do iy=1,y

t(ix,iy)=abs(r1(ix,iy))*sqrt(63/(1-r1(ix,iy)**2))

enddo

运行的结果是：

我尝试了定义a(ix,iy)=63/(1-r1(ix,iy)**2)，print出来之后显示全为正值，所以不应该是开方里面有负数的问题啊？

困扰我很久了，请各位大神帮帮忙啊！

freekiller · 发表于 2016-4-30 20:17:43

出现了负数，检查值

言深深 · 发表于 2016-5-1 08:29:49

输出，1-r1(ix,iy)**2，并看看结果，是否出现0，或者小于0的情况

irides · 发表于 2016-5-3 15:46:38

freekiller 发表于 2016-4-30 20:17
出现了负数，检查值

谢谢大神，已解决！

irides · 发表于 2016-5-3 15:47:58

言深深发表于 2016-5-1 08:29
输出，1-r1(ix,iy)**2，并看看结果，是否出现0，或者小于0的情况

谢谢版主，是NaN的问题，已解决

石榴之爱 · 发表于 2016-12-26 12:49:44

irides 发表于 2016-5-3 15:47
谢谢版主，是NaN的问题，已解决

你好，请问NaN的问题是怎么解决的？刚好也遇到相同的问题

irides · 发表于 2016-12-28 15:46:42

石榴之爱发表于 2016-12-26 12:49
你好，请问NaN的问题是怎么解决的？刚好也遇到相同的问题

将资料中NaN的部分赋值为0了

		自动登录	找回密码
密码			立即注册

[求助] fortran运行的时候出现sqrt:DOMAIN error问题