feclearexcept
提供: cppreference.com
ヘッダ
<fenv.h> で定義
int feclearexcept( int excepts );
(C99以上)
ビットマスク引数 excepts で指定された浮動小数点例外のクリアを試みます。 excepts は浮動小数点例外マクロのビット単位の論理和です。
[編集] 引数
excepts
-
クリアする例外フラグを指定するビットマスク
[編集] 戻り値
指定された例外すべてのクリアに成功した場合、または excepts がゼロの場合は 0。 エラーが発生した場合は非ゼロの値を返します。
[編集] 例
Run this code
#include <fenv.h> #include <stdio.h> #include <math.h> #include <float.h> /* * A possible implementation of hypot which makes use of many advanced * floating point features. */ double hypot_demo(double a, double b) { const int range_problem = FE_OVERFLOW | FE_UNDERFLOW ; feclearexcept(range_problem); // try a fast algorithm double result = sqrt (a * a + b * b); if (!fetestexcept (range_problem)) // no overflow or underflow return result; // return the fast result // do a more complicated calculation to avoid overflow or underflow 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); // we can ignore the smaller value // scale so that fabs(a) is near 1 double a_scaled = scalbn (a, -a_exponent); double b_scaled = scalbn (b, -a_exponent); // overflow and underflow is now impossible result = sqrt (a_scaled * a_scaled + b_scaled * b_scaled); // undo scaling return scalbn (result, a_exponent); } int main(void) { // Normal case takes the fast route printf ("hypot(%f, %f) = %f\n", 3.0, 4.0, hypot_demo(3.0, 4.0)); // Extreme case takes the slow but more accurate route 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