Как использовать и проверять четырехкратную точность в gfortran?

Я пытаюсь использовать четырехкратную точность в gfortran, но похоже, что real*16не работает. После некоторого поиска я обнаружил, что это может быть реализовано как real*10. Естьreal*10на самом деле четырехкратная точность?

Как я могу проверить точность моего кода? Существует ли стандартный простой алгоритм проверки точности? Например, когда я хочу выяснить, что такое компьютерный ноль, я продолжаю делить на 2,0, пока не достигну 0,0. Отслеживание значений позволяет мне узнать, когда компьютер «думает», что мое не-нулевое число равно нулю -, что дает мне компьютерный ноль.

Есть ли хороший способ определить точность с помощью алгоритма, который я описал?