feclearexcept

#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;
}