Jump to content

Guard digit

From Wikipedia, the free encyclopedia
(Redirected from Guard digits)

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

Example

[edit]

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.

In computing

[edit]

Guard digits are also used in floating point operations in most computer systems.

As an example, consider the subtraction . Here, the product notation indicates a binary floating point representation with the exponent of the representation given as a power of two and with the significand given with three bits after the binary point. To compute the subtraction it is necessary to change the forms of these numbers so that they have the same exponent, and so that when the product notation is expanded the resulting numbers have their binary points lined up with each other. Shifting the second operand into position, as , gives it a fourth digit after the binary point. This creates the need to add an extra digit to the first operand—a guard digit—putting the subtraction into the form .

Performing this operation gives as the result or .

Without using a guard digit the subtraction would be performed only to three bits of precision, as , yielding or , twice as large as the correct result. Thus, in this example, the use of a guard digit led to a more accurate result.

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.