Warum wird die gleiche Gleitkommakonstante in Fortran und in C anders gedruckt?

Question

Druck konstant 0.8804418 in Fortran

program hello 
     print *, 0.8804418 
    end program hello 

gibt 0.880441785 während die C-Version

#include <stdio.h> 

int main() { 
    printf("%.10g", 0.8804418); 
} 

druckt 0.8804418 auf demselben System (Linux x86-64 mit gfortran und gcc). Warum ist die Ausgabe anders? Beachten Sie, dass die Erhöhung der Genauigkeit in printf die Ausgabe nicht ändert.

Dies ist kein Duplikat von Is floating point math broken? oder ähnlichem. Diese Frage bezieht sich speziell auf Unterschiede in Fortran und C-Darstellung (oder Formatierung).

Answer 1

Standardmäßig sind Fortran Anzahl Konstanten einfachpräzise; In C haben jedoch Fließkomma-Literale doppelte Präzision.

Wenn Sie 0.8804418 einfacher Genauigkeit übersetzen und dann als double in C drucken, erhalten Sie 0.8804417849 gedruckt (demo)

float x = 0.8804418f; 
printf("%.10g\n", x); 

Fortran der Ausdruck die gleiche Zahl aufzurunden zu sein scheint.

Fortran Syntax für doppelte Genauigkeit REAL Zahlen verwendet Suffix d:

print *, 0.8804418d+0 

Dieser druckt 0.88044180000000005 (demo).