I was being stupid. Sorry. My reply was correct but unhelpful. You were using -mfpmath=387 (possible as the default), and that uses Intel native arithmetic, which is NOT standard IEEE 754 double, but an extended format. If you use -mfpmath=sse, you will get what you expect. Regards, Nick Maclaren, University of Cambridge Computing Service, New Museums Site, Pembroke Street, Cambridge CB2 3QH, England. Email: nmm1@xxxxxxxxx Tel.: +44 1223 334761 Fax: +44 1223 334679
