Fortran'da çift kesinlik değerleri nasıl ayarlanır?

Son zamanlarda, FORTRAN95 ile garip bir sorunla karşılaştım. X ve Y değişkenlerini aşağıdaki gibi başlattım:

X=1.0
Y=0.1

Daha sonra bunları bir araya getirip sonucu yazdırıyorum:

1.10000000149012

Değişkenler incelendikten sonra, 0.1 tam hassasiyetle çift hassasiyetle temsil edilmemiş gibi görünüyor. Bundan kaçınmamın bir yolu var mı?

Bunu yapmanın başka bir yolu da, önce SELECTED_REAL_KIND intrinsic'i kullanarak değişkente istediğiniz hassasiyeti açıkça belirtmek ve daha sonra bunu değişkenleri tanımlamak ve başlatmak için kullanmaktır. Gibi bir şey:

INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(15)
REAL(dp) :: x
x = 1.0_dp

Bunu bu şekilde yapmanın güzel bir avantajı, dpmodülün tanımını , sonra USEda gerekli olduğunda bu modülde depolayabilmenizdir . Şimdi programınızın hassasiyetini değiştirmek isterseniz , değişken başlatmalarınızın sonundaki dptüm D0leri aramak ve değiştirmek yerine yalnızca o yerde tanımını değiştirmeniz gerekir. (Bu yüzden 1.0D-1Y'yi önerilen şekilde tanımlamak için sözdizimini kullanmamanızı tavsiye ederim . Çalışır, ancak gelecekteki tüm örnekleri bulmayı ve değiştirmeyi zorlaştırır.)

Fortran Wiki'deki bu sayfa hakkında iyi ek bilgiler vermektedir SELECTED_REAL_KIND.

Değişkenleri çift kesinlik olarak beyan ettiniz, ancak tek kesinlik değerleriyle başlattınız.

Yazmış olabilirsiniz:

X=1.0d0
Y=1.0d-1

Barron'un aşağıdaki yanıtı , değişkenlerinizin hassasiyetini daha sonra değiştirmenize izin vermesi avantajı ile gerçek bir çift hassasiyet yapmanın başka bir yoludur.