While
it is possible to perform IEEE-754 double precision arithmetic
correctly on the IA32 architecture by writing every intermediate
result to memory,

This is not sufficient to simulate double x double -> double
operations---pre- and post-scalings, as well as setting the mantissa
precision to 53 bits, are necessary to eliminate double rounding and to
compute possible denormal results correctly.