Verbessern Präzision auf Variablen durch ganzzahlige Quotienten definiert

Question

Say I haben folgendes Programm:

program derp 
    implicit none 
    integer, parameter :: ikind = selected_real_kind(18) 
    real (kind = ikind) :: a = 2.0/3.0 
    print*, a 
end program derp 

Das Programm derp0.6666666865348815917 Ausgänge, die eindeutig nicht mehr als 18 Ziffern der Präzision sind. Wenn ich jedoch a=2.0 und mit der gleichen Methode und definieren, dann definieren c=a/b ich bekomme eine Ausgabe von 0.666666666666666666685, was gut ist. Wie definiere ich einfach eine Variable als Quotient von ganzen Zahlen und habe sie alle Stellen der Genauigkeit, die ich will, von selected_real_kind?

Answer 1

Versuchen: real (kind = ikind) :: a = 2.0_ikind/3.0_ikind

Der Grund dafür ist, während die LHS hohe Präzision, die RHS in Ihrem Code Beispiel 2.0/3.0, nicht der Fall. Fortran führt diese Berechnung in einfacher Genauigkeit durch und weist das Ergebnis dann der LHS zu. Die RHS-Seite wird nicht mit höherer Genauigkeit berechnet, da die LHS eine hohe Genauigkeit aufweist. digits_kind ist die Art der Angabe des Typs einer Konstanten digits.