Float point issue

[Robert Dewar]

On 10/28/2011 12:35 AM, Zoltán Kócsi wrote:
> On Thu, 27 Oct 2011 23:31:14 -0400
> Robert Dewar<dewar@adacore.com>  wrote:
>
>>> - I am missing a gcc flag
>>
>> probably you should avoid extra precision and all the
>> issues it brings, as well as speed up your program, by
>> using SSE 64-bit arithmetic (using the appropriate gcc
>> flags)
>
> Indeed. -mpc64 fixes the issue and proper 53-bit rounding is applied.
> Thanks a lot.

If you are going to set 64-bit mode, you might as well specify
SSE artihmetic, it will be significantly faster for most applications.
>
> Zoltan