Оператор Fortran fomat с высочайшей точностью в системе

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

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

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

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

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

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

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

Изменить – запрос пользователя 4050 Попробуйте следующий пример.

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 значащих цифр.