newlib/libm/common/sf_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 #ifdef __STDC__
  10         float fdimf(float x, float y)
  11 #else
  12         float fdimf(x,y)
  13         float x;
  14         float y;
  15 #endif
  16 {
  17   if (__fpclassifyf(x) == FP_NAN)  return(x);
  18   if (__fpclassifyf(y) == FP_NAN)  return(y);
  19
  20   return x > y ? x - y : 0.0;
  21 }
  22
  23 #ifdef _DOUBLE_IS_32BITS
  24
  25 #ifdef __STDC__
  26         double fdim(double x, double y)
  27 #else
  28         double fdim(x,y)
  29         double x;
  30         double y;
  31 #endif
  32 {
  33   return (double) fdimf((float) x, (float) y);
  34 }
  35
  36 #endif /* defined(_DOUBLE_IS_32BITS) */