フォートランの2.0d-5を出力したい。

Question

フォートランについてです。
real(8) dx
dx = 2.0d-5
write(,) dx
と入力したときに出力結果が 2.0000000000000002E-005とちょうど0.00002にならなくて困っています。
このdxに正しい値を入力する方法を教えていただきたいです。
もちろん
real(8) dx
dx = 2.0d0
write(,) dx
と入力した場合はちょうど2.0000000000000000と出力されました。
よろしくお願いします。

Answer 1

real*8の精度は十進表示の場合約15桁です。
writeのデフォルトが何故か17桁のようですが、下2－3桁はゴミです。

write(*,'(E20.14)') dx

Comments

[okkun-gg]

ありがとうございました。m(__)m

[otn]

> real*8の精度は十進表示の場合約15桁です。
正確には、「今のメジャーなCPUで採用されている浮動小数点形式では」という前提で。

[ikadzuchi]

> 何故か
17桁はdoubleの値を10進数で正確に記述するのに必要な桁数ですので、それででしょうね。
ここで「正確に」とは、任意のdoubleを10進数で表記した文字列を再度doubleに変換して元の値に戻るということを指します。

Answer 2

double precision dx
dx = 2.0d-5
write(,) dx

ではどうなりますか。

Comments

[okkun-gg]

otnさんの回答のように精度の桁数が問題であったようです。回答ありがとうございました。