关于新旧版本的GCC对浮点数运算的一些行为（读《CS:APP》）

《CS:APP》第二章中说到GCC对浮点运算的一些处理方式，在Intel IA32的处理器中浮点数寄存器使用的是80位的扩展精度格式，而在GCC中float和double类型的浮点数长度为32和64。这就导致在浮点数寄存器向内存中写入数据的时候会产生精度丢失。书中的例子如下:

#include <stdio.h>

double recip(int denom)
{
return 1.0/(double) denom;
}

void do_nothing(){}

void test1(int denom)
{
double r1,r2;
int t1,t2;

r1=recip(denom);
r2=recip(denom);
t1=r1==r2;
do_nothing();
t2=r1==r2;
printf("test1 t1:r1 %f %c=r2 %f/n",r1,t1?'=':'!',r2);
printf("test1 t2:r1 %f %c=r2 %f/n",r1,t2?'=':'!',r2);
}

main()
{
test1(10);
}

 

按照书中描述的GCC的版本得出的结果应该是(没有试过)

test1 t1: r1 0.100000!=r2 0.100000

test1 t2: r1 0.100000==r2 0.100000

 

但是在我的系统下，gcc 版本 4.2.3 (Ubuntu 4.2.3