This paper discusses a FORTRAN subprogram library developed primarily to support extended-precision floating-point arithmetic. The general strategy, which makes limited use of guard digits, is developed to achieve high accuracy with reasonable execution time and storage space. In addition to describing some previously unpublished algorithms, the authors present subprograms for simulating extended-precision arithmetic and for input and output conversion.
PDF