Re: bug in float on Pentium

From: Matt Billenstein (mbillens@mbillens.dyndns.org)
Date: Sun Apr 15 2001 - 01:02:04 EST


It's not that you found a new bug or that floats are inaccurate (they are
just less exact than doubles)... For example, if you make your program
print some more digits, you'll get:

5483.98999999999978172127157449722290039062500000000000
5483.99023437500000000000000000000000000000000000000000

#include <stdio.h>

int main() {

    unsigned int *X;
    unsigned int *Y;

 double x = 5483.99;
 float y = 5483.99;

    X = (unsigned int *)&x;
    Y = (unsigned int *)&y;

 printf ("%60.50lf\n%60.50f\n", x, y);

 printf("%lf %x%x %f %x\n", x, X[1], X[0], y, *Y);
 return 0;
}

which you can verify by hand:

5483.990000 40b56bfd70a3d70a
0 10000001011 1.0101 0110 1011 1111 1101 0111 0000 1010 0011 1101 0111 0000
1010
+ 2^12 * 1.3388647460937499467092948179925 =
5483.989999999999781721271574497222900390625

5483.990234 45ab5fec
0 10001011 1.010 1011 0101 1111 1110 1100
+ 2^12 * 1.338864803314208984375 = 5489.990234375

check out:
http://www.psc.edu/general/software/packages/ieee/ieee.html

l8r

m

Matt Billenstein
mbillens (at) one (dot) net
http://w3.one.net/~mbillens/

----- Original Message -----
From: "Joe" <joeja@mindspring.com>
To: <linux-kernel@vger.kernel.org>
Sent: Friday, April 13, 2001 7:23 PM
Subject: bug in float on Pentium

| Not sure but I think I found a NEW bug.
|
| I know that there have been some issues with pentiums and floating point
| arrithmatic, but this takes the cake...
|
| Linux Lserver.org 2.2.18 #43 SMP Fri Mar 9 14:19:41 EST 2001 i586
| unknown
|
| >kgcc --version
| egcs-2.91.66
|
| RH 6.2.x / 7.0
|
| try this program
|
| #include <stdio.h>
|
| int main() {
|
| char tmpx[100];
| char tmpy[100];
|
| double x = 5483.99;
| float y = 5483.99;
|
| sprintf (tmpx, "%f",x );
| sprintf (tmpy, "%f",y );
|
| printf ("%s\n%s\n", tmpx, tmpy);
| return 0;
| }
|
|
| I am getting the following as output
|
| joeja@Lserver$ ./testf
| 5483.990000
| 5483.990234
|
|
| what is with the .990234?? it should be .990000
|
| any ideas on this??
|
| --
| Joe Acosta ........
| home: joeja@mindspring.com
|
|
|
| -
| To unsubscribe from this list: send the line "unsubscribe linux-kernel" in
| the body of a message to majordomo@vger.kernel.org
| More majordomo info at http://vger.kernel.org/majordomo-info.html
| Please read the FAQ at http://www.tux.org/lkml/

-
To unsubscribe from this list: send the line "unsubscribe linux-kernel" in
the body of a message to majordomo@vger.kernel.org
More majordomo info at http://vger.kernel.org/majordomo-info.html
Please read the FAQ at http://www.tux.org/lkml/



This archive was generated by hypermail 2b29 : Sun Apr 15 2001 - 21:00:23 EST