关于浮点运算的两个函数：fetestexcept、feclearexcept

fetestexcept和feclearexcept是两个处理浮点运算异常的函数，在头文件是<fenv.h>中，该头文件中包含一系
列的处理浮点运算异常的函数，而我们今天介绍的是两个最简单，又是最常用的函数。

详细请看：fenv.h

1.fetestexcept介绍

1.1函数原型

int fetestexcept(int excepts);

1.2函数描述

用于查询指定异常excepts的状态值，其中：
excepts是异常的宏，它可以是：FE_INEXACT,FE_UNDERFLOW,FE_OVERFLOW,FE_DIVBYZERO,FE_INVALID的组合值，
也可以是FE_ALL_EXCEPT。
这些异常的具体含义如下：
FE_INEXACT：表示精度损失的异常
FE_UNDERFLOW：表示下溢
FE_OVERFLOW：表示上溢
FE_DIVBYZERO：表示被0除
FE_INVALID：表示不合法的浮点运算
FE_ALL_EXCEPT：表示所有异常

2.feclearexcept介绍

2.1函数原型

int feclearexcept(int excepts);

2.2函数描述

主要用以清除指定异常excepts的状态值，其中入参excepts同fetestexcept的入参。

3.一个简单的例子

源码如下：

#include <stdio.h>
#include <math.h>       //pow , sqrt
#include <fenv.h>       //fetestexcept , feclearexcept , FE_INEXACT...

using namespace std;

double divide(double x, double y)
{
	double res = x/y;
	return res;
}

void show_all_except(void)
{
	//设置该参数表示可以访问浮点运算异常的状态值
	#pragma STDC FENV_ACCESS ON
	
	//获取所有可能的异常
	int res = fetestexcept(FE_ALL_EXCEPT);
	
	printf( "Status: " );
	if(res & FE_INEXACT)//判断是否有精度损失的异常
	{
		printf( " inexact ");
	}
	if(res & FE_UNDERFLOW)//判断是否有下溢
	{
		printf( " underflow ");
	}
	if(res & FE_OVERFLOW)//判断是否有上溢
	{
		printf( " overflow ");
	}
	if(res & FE_DIVBYZERO)//判断是否有被0除
	{
		printf( " zero-divide ");
	}
	if(res & FE_INVALID)//判断是否有不合法的浮点运算
	{
		printf( " invalid ");
	}
	printf( " \n ");
	
	//清除所有异常
	feclearexcept(FE_ALL_EXCEPT);

}
int main( void )
{
	double a = 1e-40, b;
	float x, y;
	
	//没有异常
	show_all_except();
	
	//精度损失 and 下溢
	y = a;
	show_all_except();
	
	//精度损失 and 上溢
	double p = pow((double)1000,(int)1000);
	show_all_except();
	
	//被0除
	double d = divide(2.0, 0.0);
	show_all_except();
	
	//不合法的浮点运算
	double i = sqrt((double)-1);
	show_all_except();
	
	return 0;
}

输出如下：

 Status:  
 Status:  inexact  underflow  
 Status:  inexact  overflow  
 Status:  zero-divide  
 Status:  invalid