The quadruple precision available in ifort-9.0 is supposed to
be a IEEE-style floating point format with 113 bits of accuracy, right?
It seems that sometimes it doesn't seem to have the full 113 bits
of accuracy; the following program finds the smallest i such that
fl(2^i + 1) - 2^i != 1, but if fails around i = 81 already.
Is there a but in the quad precision emulation somewhere?