Computational Domain
Range of numbers (from /usr/include/limits.h):
FLT_DIG 6 /* decimal digits of precision of a float */
DBL_DIG 15 /* decimal digits of precision of a double */
DBL_MAX 1.7976931348623157E+308
DBL_MIN 2.2250738585072014E-308
LONGLONG_MIN -9223372036854775807LL-1LL
LONGLONG_MAX 9223372036854775807LL
ULONGLONG_MAX 18446744073709551615LLU
The extended precision (REAL*16) is available and supported by the compiler. But this mode of calculation is slow (by factor ~40)