Guard digit

From Wikipedia the free encyclopedia

In numerical analysis, one or more guard digits can be used to reduce the amount of roundoff error.

For example, suppose that the final result of a long, multi-step calculation can be safely rounded off to N decimal places. That is to say, the roundoff error introduced by this final roundoff makes a negligible contribution to the overall uncertainty.

However, it is quite likely that it is not safe to round off the intermediate steps in the calculation to the same number of digits. Be aware that roundoff errors can accumulate. If M decimal places are used in the intermediate calculation, we say there are M−N guard digits.

Guard digits are also used in floating point operations in most computer systems. Given we have to line up the binary points. This means we must add an extra digit to the first operand—a guard digit. This gives us . Performing this operation gives us or . Without using a guard digit we have , yielding or . This gives us a relative error of 1. Therefore, we can see how important guard digits can be.

An example of the error caused by floating point roundoff is illustrated in the following C code.

int main(){    double a;    int i;     a = 0.2;     a += 0.1;     a -= 0.3;     for (i = 0; a < 1.0; i++)         a += a;     printf("i=%d, a=%f\n", i, a);     return 0; } 

It appears that the program should not terminate. Yet the output is:

i=54, a=1.000000

Another example is:

Take two numbers:

and

We bring the first number to the same power of as the second one:

The addition of the two numbers is:

0.0256*10^2  2.3400*10^2 +   ____________  2.3656*10^2  

After padding the second number (i.e., ) with two s, the bit after is the guard digit, and the bit after is the round digit. The result after rounding is as opposed to , without the extra bits (guard and round bits), i.e., by considering only . The error therefore is .

See also[edit]

References[edit]

  • Forman S. Acton. Numerical Methods that Work, The Mathematical Association of America (August 1997).
  • Higham, Nicholas J. Accuracy and Stability of Numerical Algorithms, Washington D.C.: Society for Industrial & Applied Mathematics, 2002.