3.1.6 MachinePrecision Numbers
Whenever machineprecision numbers appear in a calculation, the whole calculation is typically done in machine precision. Mathematica will then give machineprecision numbers as the result.
Whenever the input contains any machineprecision numbers, Mathematica does the computation to machine precision.
In[1]:= 1.4444444444444444444 ^ 5.7
Out[1]=
Zeta[5.6] yields a machineprecision result, so the N is irrelevant.
In[2]:= N[Zeta[5.6], 30]
Out[2]=
This gives a higherprecision result.
In[3]:= N[Zeta[56/10], 30]
Out[3]=
When you do calculations with arbitraryprecision numbers, as discussed in the previous section, Mathematica always keeps track of the precision of your results, and gives only those digits which are known to be correct, given the precision of your input. When you do calculations with machineprecision numbers, however, Mathematica always gives you a machineprecision result, whether or not all the digits in the result can, in fact, be determined to be correct on the basis of your input.
This subtracts two machineprecision numbers.
In[4]:= diff = 1.11111111  1.11111000
Out[4]=
The result is taken to have machine precision.
In[5]:= Precision[diff]
Out[5]=
Here are all the digits in the result.
In[6]:= InputForm[diff]
Out[6]//InputForm= 1.1099999999153454`*^6
The fact that you can get spurious digits in machineprecision numerical calculations with Mathematica is in many respects quite unsatisfactory. The ultimate reason, however, that Mathematica uses fixed precision for these calculations is a matter of computational efficiency.
Mathematica is usually set up to insulate you as much as possible from the details of the computer system you are using. In dealing with machineprecision numbers, you would lose too much, however, if Mathematica did not make use of some specific features of your computer.
The important point is that almost all computers have special hardware or microcode for doing floatingpoint calculations to a particular fixed precision. Mathematica makes use of these features when doing machineprecision numerical calculations.
The typical arrangement is that all machineprecision numbers in Mathematica are represented as "doubleprecision floatingpoint numbers" in the underlying computer system. On most current computers, such numbers contain a total of 64 binary bits, typically yielding 16 decimal digits of mantissa.
The main advantage of using the builtin floatingpoint capabilities of your computer is speed. Arbitraryprecision numerical calculations, which do not make such direct use of these capabilities, are usually many times slower than machineprecision calculations.
There are several disadvantages of using builtin floatingpoint capabilities. One already mentioned is that it forces all numbers to have a fixed precision, independent of what precision can be justified for them.
A second disadvantage is that the treatment of machineprecision numbers can vary slightly from one computer system to another. In working with machineprecision numbers, Mathematica is at the mercy of the floatingpoint arithmetic system of each particular computer. If floatingpoint arithmetic is done differently on two computers, you may get slightly different results for machineprecision Mathematica calculations on those computers.
Properties of numbers on a particular computer system.
Since machineprecision numbers on any particular computer system are represented by a definite number of binary bits, numbers which are too close together will have the same bit pattern, and so cannot be distinguished. The parameter $MachineEpsilon gives the distance between 1.0 and the closest number which has a distinct binary representation.
This gives the value of $MachineEpsilon for the computer system on which these examples are run.
In[7]:= $MachineEpsilon
Out[7]=
Although this prints as 1., Mathematica knows that the result is larger than 1.
In[8]:= 1. + $MachineEpsilon
Out[8]=
Subtracting 1 gives $MachineEpsilon.
In[9]:= %  1.
Out[9]=
This again prints as 1.
In[10]:= 1. + $MachineEpsilon/2
Out[10]=
In this case, however, subtracting 1 yields 0, since 1 + $MachineEpsilon/2 is not distinguished from 1. to machine precision.
In[11]:= %  1.
Out[11]=
Machine numbers have not only limited precision, but also limited magnitude. If you generate a number which lies outside the range specified by $MinMachineNumber and $MaxMachineNumber, Mathematica will automatically convert the number to arbitraryprecision form.
This is the maximum machineprecision number which can be handled on the computer system used for this example.
In[12]:= $MaxMachineNumber
Out[12]=
Mathematica automatically converts the result of this computation to arbitrary precision.
In[13]:= Exp[1000.]
Out[13]=
