Cet*_*ert 2 floating-point fortran notation double-precision
有没有人知道后缀" d + 0 "在下面的M1,M2和M4的作业中是什么意思,或者网上是否有任何资源或一本很可能找到这些信息的书?
subroutine plot( t, x, p, q, nga, nt, wron,
& ngq, gq, ngaq1, ngaq2, gaq, rwh, iwh )
implicit none
integer*4 nga, nt, ngq, ngaq1, ngaq2, iwh(*)
real*8 t, x(*), p(*), q(*), wron(nga,*),
& gq(ngq,*), gaq(ngaq1,ngaq2,*), rwh(*)
real*8 M1, M2, M3, M4, mr, mst, h3, Tc
integer*8 iflag
c DISCRETIZE1( Tc, rwh, iwh )
M1 = 0.1362d+0
M2 = 0.09806d+0
M3 = M1 + M2
M4 = 0.236d+0
mr = M1*x(1) + M2*x(2) + M3*x(3) + M4*x(4)
mst = M1*q(1) + M2*q(2) + M4*q(3)
call mess3( t, x, h3, p, q, rwh, iwh, iflag )
write (10,100) t, x(1), x(2), x(3), x(4), Tc, h3, mr-mst
100 format( e20.10, 7(1x,e20.10) )
end
Run Code Online (Sandbox Code Playgroud)
它表示"双精度",+ 0表示第10次幂(零是指数).
这是我能找到的第一个参考文献的引用:
存储在计算机中的数字在数量和精度上受到限制.限制取决于特定的计算机.因此,REAL数字只有一定数量的有效数字.如果计算需要更多有效数字,则必须使用DOUBLE PRECISION编号.DOUBLE PRECISION常量以与单精度REAL常数相同的指数形式写入,除了D而不是E将尾数与指数分开.
这里引用了更好的参考资料,www.fortran.com,特别是F77标准.
4.5.2双精度常数.
双精度常数的形式是:
- 基本实常数后跟双精度指数
- 整数常量后跟双精度指数
双精度常数的值是先于d和10的通过遵循D.形式(2)的整数恒定部分可以以多位数被写入比一个处理器将使用整数指示的功率的常数的乘积近似常数的值.