Hum, slightly amended patch, after checking 32-bit results on another linux machine.
The test for FMA has been made a bit less strict, because otherwise we have surprised on 387 arithmetic due to excess precision.

Final patch is attached. Regression-tested on x86_64-pc-linux-gnu, both 32- and 64-bit.
OK to commit?

FX