Fortran оператор с самой высокой появляются две точности в системе

Question

Кто-то хочет менее точно бы написать

999 format ('The answer is x = ', F8.3) 

Другие желающих более высокую выходную точность может написать

999 format ('The answer is x = ', F18.12) 

Таким образом, это полностью зависит от того, какие желания пользователя. Каков формат , который точно соответствует точности, используемой при расчете? (Примечание может отличаться от системы к системе)

Answer 1

Это сложный вопрос, потому что вы запрашиваете «точность расчета», которая зависит от многих факторов. Например: если я разрешаю f (x) = 0 методом Ньютона до допуска 1E-6, вы хотите формат с семью цифрами?

С другой стороны, если вы имеете в виду «наивысшую точность, достигаемую по типу» (например, двойную или одинаковую точность), вы можете просто найти соответствующий epsilon (машинный eps или точность) и использовать это как формат флаг. Если epsilon равен 1E-15, вы можете использовать флаг формата, который не должен содержать более 16 цифр.

В Fortran вы можете использовать функцию EPSILON (X) для получения этого номера (ответ будет зависеть от типа X), вы можете взять слово абсолютного значения логарифма (база 10) epsilon , и сделайте это число десятичных знаков в вашем поплавком представлении.

Например, если epsilon равен 1E-12, то регистр равен -12, абс составляет 12, а пол 12, поэтому вам нужен формат, такой как 15.12F (12 знаков после запятой + 1 точка + нуль + знак = 15 мест)

Answer 2

Проблема с числами с плавающей запятой заключается в том, что нет такой точности: только значимые цифры.

Например, если вы вычисляете долготы в реальном * 1, недалеко от Великобритании, вы были бы точны до шести знаков после запятой, но если бы вы были в Колорадо-Спрингсе, это было бы точно только до 4 знаков после запятой. Было бы бессмысленно печатать номер в формате F, это просто мусор после четвертого знака после запятой.

Если вы хотите печатать до максимальной точности, напечатайте в формате E. Так как это всегда n.nn..nEnn, вы получаете все значимые цифры.

Edit - запрос user4050 в Попробуйте следующий пример

program main 
    real intpart, multiplier 
    integer ii 

    multiplier = 1 
    do ii = 1, 6 
     intpart = 9.87654321 
     intpart = intpart * multiplier 
     print '(F15.7 E15.7 G15.8)', intpart, intpart, intpart 
     multiplier = multiplier * 10 
    end do 
    stop 
end program 

Что вы получите что-то вроде

 9.8765430 0.9876543E+01 9.8765430 
    98.7654266 0.9876543E+02 98.765427 
    987.6542969 0.9876543E+03 987.65430 
    9876.5429688 0.9876543E+04 9876.5430 
    98765.4296875 0.9876543E+05 98765.430 
987654.3125000 0.9876543E+06 987654.31 

Обратите внимание, что точность изменения, как число становится больше, потому что только поплавок имеет 7 значительных цифр.