toolchain/gcc/newlib/libm/common/s_fdim.c

   1 /* Copyright (C) 2002 by  Red Hat, Incorporated. All rights reserved.
   2  *
   3  * Permission to use, copy, modify, and distribute this software
   4  * is freely granted, provided that this notice is preserved.
   5  */
   6
   7 #include "fdlibm.h"
   8
   9 #ifndef _DOUBLE_IS_32BITS
  10
  11 #ifdef __STDC__
  12         double fdim(double x, double y)
  13 #else
  14         double fdim(x,y)
  15         double x;
  16         double y;
  17 #endif
  18 {
  19   int c = __fpclassifyd(x);
  20   if (c == FP_NAN || c == FP_INFINITE)
  21     return HUGE_VAL;
  22
  23   return x > y ? x - y : 0.0;
  24 }
  25
  26 #endif /* _DOUBLE_IS_32BITS */