diff options
Diffstat (limited to 'newlib/libm/common/sf_fdim.c')
-rw-r--r-- | newlib/libm/common/sf_fdim.c | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/newlib/libm/common/sf_fdim.c b/newlib/libm/common/sf_fdim.c new file mode 100644 index 000000000..80221ace1 --- /dev/null +++ b/newlib/libm/common/sf_fdim.c @@ -0,0 +1,37 @@ +/* Copyright (C) 2002 by Red Hat, Incorporated. All rights reserved. + * + * Permission to use, copy, modify, and distribute this software + * is freely granted, provided that this notice is preserved. + */ + +#include "fdlibm.h" + +#ifdef __STDC__ + float fdimf(float x, float y) +#else + float fdimf(x,y) + float x; + float y; +#endif +{ + int c = __fpclassifyf(x); + if (c == FP_NAN || c == FP_INFINITE) + return HUGE_VAL; + + return x > y ? x - y : 0.0; +} + +#ifdef _DOUBLE_IS_32BITS + +#ifdef __STDC__ + double fdim(double x, double y) +#else + double fdim(x,y) + double x; + double y; +#endif +{ + return (double) fdimf((float) x, (float) y); +} + +#endif /* defined(_DOUBLE_IS_32BITS) */ |