feclearexcept
来自cppreference.com
在标头
<fenv.h> 定义
int feclearexcept( int excepts );
(C99 起)
尝试清除列于位掩码实参 excepts 的浮点数异常,该实参值为浮点数异常宏的逐位或。
[编辑] 参数
excepts
-
列出欲清除的异常标志的位掩码
[编辑] 返回值
若成功清除所有指明的异常,或若 excepts 为零则返回 0。错误时返回非零。
[编辑] 示例
运行此代码
#include <fenv.h> #include <stdio.h> #include <math.h> #include <float.h> /* * 可能的 hypot 实现,它会活用许多高级浮点数特性。 */ double hypot_demo(double a, double b) { const int range_problem = FE_OVERFLOW | FE_UNDERFLOW ; feclearexcept(range_problem); // 尝试快速算法 double result = sqrt (a * a + b * b); if (!fetestexcept (range_problem)) // 未上溢或下溢 return result; // 返回结果 // 做更多复杂计算以避免上溢或下溢 int a_exponent,b_exponent; frexp (a, &a_exponent); frexp (b, &b_exponent); if (a_exponent - b_exponent > DBL_MAX_EXP ) return fabs (a) + fabs (b); // 我们可以忽略小值 // 缩放以令 fabs(a) 接近1 double a_scaled = scalbn (a, -a_exponent); double b_scaled = scalbn (b, -a_exponent); // 现在上溢和下溢是不可能的 result = sqrt (a_scaled * a_scaled + b_scaled * b_scaled); // 撤销缩放 return scalbn (result, a_exponent); } int main(void) { // 通常情况选择较快的路线 printf ("hypot(%f, %f) = %f\n", 3.0, 4.0, hypot_demo(3.0, 4.0)); // 极限情形会选择较慢但更准确的路线 printf ("hypot(%e, %e) = %e\n", DBL_MAX / 2.0, DBL_MAX / 2.0, hypot_demo(DBL_MAX / 2.0, DBL_MAX / 2.0)); return 0; }
输出:
hypot(3.000000, 4.000000) = 5.000000 hypot(8.988466e+307, 8.988466e+307) = 1.271161e+308