feclearexcept
来自cppreference.com
| 在标头 <fenv.h> 定义
|
||
| |
(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
引用
- C11 标准(ISO/IEC 9899:2011):
- 7.6.2.1 The feclearexcept function (第 209 页)
- C99 标准(ISO/IEC 9899:1999):
- 7.6.2.1 The feclearexcept function (第 190 页)
参阅
(C99) |
确认设置了哪些浮点数异常状态标志 (函数) |
feclearexcept 的 C++ 文档
| |