Sourceware Bugzilla – Attachment 5991 Details for
Bug 13304
fma, fmaf, fmal produce wrong results
Home
|
New
|
Browse
|
Search
|
[?]
|
Reports
|
Requests
|
Help
|
New Account
|
Log In
[x]
|
Forgot Password
Login:
[x]
test case #1 for fma()
fma-bug1.c (text/x-csrc), 680 bytes, created by
Bruno Haible
on 2011-10-17 02:04:29 UTC
(
hide
)
Description:
test case #1 for fma()
Filename:
MIME Type:
Creator:
Bruno Haible
Created:
2011-10-17 02:04:29 UTC
Size:
680 bytes
patch
obsolete
>/* Bug in fma() of glibc 2.11.3 on x86_64. */ > >#include <math.h> >#include <stdio.h> > >int >main () >{ > volatile double x = 1.5; /* 3 * 2^-1 */ > volatile double y = x; > volatile double z = ldexp (1.0, 54); /* 2^54 */ > /* x * y + z with infinite precision: 2^54 + 9 * 2^-2. > Lies between (2^52 + 0) * 2^2 and (2^52 + 1) * 2^2 > and is closer to (2^52 + 1) * 2^2, therefore the rounding > must round up and produce (2^52 + 1) * 2^2. */ > volatile double expected = z + 4.0; > volatile double result = fma (x, y, z); > > if (result != expected) > { > printf ("Expected: %A\n", expected); > printf ("Actual: %A\n", result); > return 1; > } > > return 0; >}
<b>You cannot view the attachment while viewing its details because your browser does not support IFRAMEs. <a href="attachment.cgi?id=5991">View the attachment on a separate page</a>.</b>
View Attachment As Raw
Actions:
View
Attachments on
bug 13304
:
5990
| 5991 |
5992
|
5993
|
5994
|
5995
|
5996
|
5997
|
5998
|
5999
|
6000
|
6001
|
6002
|
6003
|
6004
|
6005
|
6006
|
6008
|
6014
|
6015
|
6017
|
6657