neil@daikokuya.co.uk said:
> I'm afraid I don't see your point; your code seems to support
> everything I said. 

So I've appearently mistaken what you said:

neil@daikokuya.co.uk said:
> d) Despite many claims to the contrary that you will read on the
>    net, and on GCC lists, there is no double-rounding problem.
>    All results are correctly rounded to 53-bits given the value
>    they operated on

I had read this as "gcc will emit code which makes sure that
double results are correctly rounded even if the FPU is in
extended precision mode".
My code was just to show that this is not the case.

So we can't have good "long double" support and correct
"double" arithmetics at the same time. (at least as long
as we have to deal with pre-SSE2 CPUs)

best regards
Matthias




-----------------------------------------------------------------------------------------
-----------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich

Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzende des Aufsichtsrats: MinDirig'in Baerbel Brumme-Bothe
Vorstand: Prof. Dr. Achim Bachem (Vorsitzender), Dr. Ulrich Krafft (stellv. 
Vorsitzender)
-----------------------------------------------------------------------------------------
-----------------------------------------------------------------------------------------