浮点数的比较

简介

运算的修正

修正的方法是选取一个极小的区间，如果数a,落在[ b - eps, b + eps] 上则可以认为是数a和数b是相等的。这里选择eps = 1e-8(10的-8次方)，因为这个区间可以兼容大多数数的比较

• 等于
  
  当 b - eps < a < b + eps 说明 a == b 成立

bool equal(double a, double b)
{
	return fabs(a - b) < eps; 
}

写成宏定义：

#define equal(a,b) ( (abs((a) - (b))) < (eps) )

• 小于
  
  当 a < b - eps 说明 a < b 成立

bool less(double a, double b)
{
	return a - b < -eps; 
}

写成宏定义：

#define less(a,b) ( (a) - (b) < (-eps) )

• 大于
  
  当 a > b + eps 说明 a > b 成立

bool more(double a, double b)
{
	return a - b > eps; 
}

写成宏定义：

#define more(a,b) ( (a) - (b) > (eps) )

• 大于等于
  
  当 a > b - eps 说明 a ≥ b 成立

bool more_equ(double a, double b)
{
	return a - b > -eps; 
}

写成宏定义：

#define more_equ(a,b) ( (a) - (b) > (-eps) )

• 小于等于
  
  当 a < b + eps 说明 a ≤ b 成立

bool less_equ(double a, double b)
{
	return a - b < eps; 
}

写成宏定义：

#define less_que(a,b) ( (a) - (b) < (eps) )

关于Π的计算

因为 cos( pai ) = -1 所以有 pai = arccos(-1)
直接得到pai的大小为： double pai = acos( -1 )