fdiml.c raw

   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
  19