24小时热门版块排行榜

返回列表

virtualzx

木虫 (著名写手)

应助: 263 (大学生)
金币: 7161.3
红花: 54
帖子: 1605
在线: 317.6小时
虫号: 2069080
注册: 2012-10-18
性别: GG
专业: 理论和计算化学

引用回帖:

8楼: Originally posted by hshx9291 at 2016-01-25 09:32:38
您好，我编的小程序如下：
program test
real*8 Dli
Dli=10**12
write(*,*)Dli
end

输出结果为-727379968.000000，仍不科学啊，有什么错误吗？
另外我原程序需要计算的是Dli(1)=(2.6 ...

看了你的程序，问题在于每一个常数都有类型的，10**12这里10是整型（缺省整数是integer*4)，这个类别只能达到10位。

取决于你需要的精度，如果你的程序需要完全精确的数值时才需要整数，整数没有任何误差，但能表示的大小范围很有限。你可以用10_8来指定10是8字节整数。

如果你不需要精确数值，而可以接受一定误差，那么就可以用浮点数。浮点数可以直接用科学计数法表示，比如10^15就是1D15

发自小木虫IOS客户端

赞一下

回复此楼

11楼2016-01-25 11:15:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hshx9291

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 264.5
帖子: 10
在线: 5.3小时
虫号: 4188722
注册: 2015-11-01
性别: GG
专业: 燃烧学

引用回帖:

11楼: Originally posted by virtualzx at 2016-01-25 11:15:31
看了你的程序，问题在于每一个常数都有类型的，10**12这里10是整型（缺省整数是integer*4)，这个类别只能达到10位。
取决于你需要的精度，如果你的程序需要完全精确的数值时才需要整数，整数没有任何误差，但能表示 ...

谢谢您的耐心解答，我原程序有大量的计算式，需要精确数值。如果要得到精确结果，每个常数后面都要写加上_8吗？例如程序中有一段计算为：
rholi(1)=410.70110+1.13071*(722-Tda)-9.5509*10**-4*(722-Tda)**2+6.7842*10**-7*(722-Tda)**3
      cpli(1)=0.00692*Tda**2-0.532*Tda+1720.3
      muli(1)=175.1*10**-3*exp(-10.31+1859/Tda)
      muli(2)=127.58*10**-3*exp(-9.47+1463/Tda)
      kli(1)=(0.4607*(1-0.00138*Tda)**0.38)/(Tda**(1/6))
      Dli(1)=(2.609*(10**-12)*Tda)/muli(2)
怎样才能让计算结果显示精确呢？我只在Dli(1)的赋值语句中加上_8显示结果仍为错误的。有没有全局设定什么的？谢谢您。

赞一下

回复此楼

12楼2016-01-25 14:41:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hshx9291

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 264.5
帖子: 10
在线: 5.3小时
虫号: 4188722
注册: 2015-11-01
性别: GG
专业: 燃烧学

引用回帖:

再比如，我又编了个小程序：
program test
real*8 Dli
real*8::Tda=300.00
real*8::muli=118.358
Dli=(2.609_8*(10_8**-12)*Tda)/muli
write(*,*)Dli
end

运算结果仍为0.00000000E+000

赞一下

回复此楼

13楼2016-01-25 15:00:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖