3.1.4 Numerical Precision
As discussed in Section 1.1.2, Mathematica can handle approximate real numbers with any number of digits. In general, the precision of an approximate real number is the effective number of decimal digits in it which are treated as significant for computations. The accuracy is the effective number of these digits which appear to the right of the decimal point. Note that to achieve full consistency in the treatment of numbers, precision and accuracy often have values that do not correspond to integer numbers of digits.
Precision and accuracy of real numbers.
This generates a number with 30digit precision.
In[1]:= x = N[Pi^10, 30]
Out[1]=
This gives the precision of the number.
In[2]:= Precision[x]
Out[2]=
The accuracy is lower since only some of the digits are to the right of the decimal point.
In[3]:= Accuracy[x]
Out[3]=
This number has all its digits to the right of the decimal point.
In[4]:= x / 10^6
Out[4]=
Now the accuracy is larger than the precision.
In[5]:= {Precision[%], Accuracy[%]}
Out[5]=
An approximate real number always has some uncertainty in its value, associated with digits beyond those known. One can think of precision as providing a measure of the relative size of this uncertainty. Accuracy gives a measure of the absolute size of the uncertainty.
Mathematica is set up so that if a number has uncertainty , then its true value can lie anywhere in an interval of size from to . An approximate number with accuracy is defined to have uncertainty , while a nonzero approximate number with precision is defined to have uncertainty .
Definitions of precision and accuracy in terms of uncertainty.
Adding or subtracting a quantity smaller than the uncertainty has no visible effect.
In[6]:= {x  10^26, x, x + 10^26}
Out[6]=
Numerical evaluation with arbitraryprecision and machineprecision numbers.
Mathematica distinguishes two kinds of approximate real numbers: arbitraryprecision numbers, and machineprecision numbers or machine numbers. Arbitraryprecision numbers can contain any number of digits, and maintain information on their precision. Machine numbers, on the other hand, always contain the same number of digits, and maintain no information on their precision.
Here is a machinenumber approximation to .
In[7]:= N[Pi]
Out[7]=
These are both arbitraryprecision numbers.
In[8]:= {N[Pi, 4], N[Pi, 20]}
Out[8]=
As discussed in more detail below, machine numbers work by making direct use of the numerical capabilities of your underlying computer system. As a result, computations with them can often be done more quickly. They are however much less flexible than arbitraryprecision numbers, and difficult numerical analysis can be needed to determine whether results obtained with them are correct.
Machine numbers.
This returns the symbol MachinePrecision to indicate a machine number.
In[9]:= Precision[ N[Pi] ]
Out[9]=
On this computer, machine numbers have slightly less than 16 decimal digits.
In[10]:= $MachinePrecision
Out[10]=
When you enter an approximate real number, Mathematica has to decide whether to treat it as a machine number or an arbitraryprecision number. Unless you specify otherwise, then if you give less than $MachinePrecision digits, Mathematica will treat the number as machine precision, and if you give more digits, it will treat the number as arbitrary precision.
Input forms for numbers.
When Mathematica prints out numbers, it usually tries to give them in a form that will be as easy as possible to read. But if you want to take numbers that are printed out by Mathematica, and then later use them as input to Mathematica, you need to make sure that no information gets lost.
In standard output form, Mathematica prints a number like this to six digits.
In[11]:= N[Pi]
Out[11]=
In input form, Mathematica prints all the digits it knows.
In[12]:= InputForm[%]
Out[12]//InputForm= 3.141592653589793
Here is an arbitraryprecision number in standard output form.
In[13]:= N[Pi, 20]
Out[13]=
In input form, Mathematica explicitly indicates the precision of the number, and gives extra digits to make sure the number can be reconstructed correctly.
In[14]:= InputForm[%]
Out[14]//InputForm= 3.1415926535897932384626433832795028842`20.
This makes Mathematica not explicitly indicate precision.
In[15]:= InputForm[%, NumberMarks>False]
Out[15]//InputForm= 3.1415926535897932385
Controlling printing of numbers.
The default setting for the NumberMarks option, both in InputForm and in functions such as ToString and OpenWrite is given by the value of $NumberMarks. By resetting $NumberMarks, therefore, you can globally change the way that numbers are printed in InputForm.
This makes Mathematica by default always include number marks in input form.
In[16]:= $NumberMarks = True
Out[16]=
Even a machineprecision number is now printed with an explicit number mark.
In[17]:= InputForm[N[Pi]]
Out[17]//InputForm= 3.141592653589793`
Even with no number marks, InputForm still uses *^ for scientific notation.
In[18]:= InputForm[N[Exp[600], 20], NumberMarks>False]
Out[18]//InputForm= 3.7730203009299398234*^260
In doing numerical computations, it is inevitable that you will sometimes end up with results that are less precise than you want. Particularly when you get numerical results that are very close to zero, you may well want to assume that the results should be exactly zero. The function Chop allows you to replace approximate real numbers that are close to zero by the exact integer 0.
Removing numbers close to zero.
This computation gives a small imaginary part.
In[19]:= Exp[ N[2 Pi I] ]
Out[19]=
You can get rid of the imaginary part using Chop.
In[20]:= Chop[%]
Out[20]=
