src/math/fdiml.c

   1 #include <math.h>
   2 #include <float.h>
   3
   4 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
   5 long double fdiml(long double x, long double y)
   6 {
   7         return fdim(x, y);
   8 }
   9 #else
  10 long double fdiml(long double x, long double y)
  11 {
  12         if (isnan(x))
  13                 return x;
  14         if (isnan(y))
  15                 return y;
  16         return x > y ? x - y : 0;
  17 }
  18 #endif