diff options
Diffstat (limited to 'newlib/libm/common/s_fdim.c')
-rw-r--r-- | newlib/libm/common/s_fdim.c | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/newlib/libm/common/s_fdim.c b/newlib/libm/common/s_fdim.c new file mode 100644 index 000000000..0010bf59d --- /dev/null +++ b/newlib/libm/common/s_fdim.c @@ -0,0 +1,26 @@ +/* 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" + +#ifndef _DOUBLE_IS_32BITS + +#ifdef __STDC__ + double fdim(double x, double y) +#else + double fdim(x,y) + double x; + double y; +#endif +{ + int c = __fpclassifyd(x); + if (c == FP_NAN || c == FP_INFINITE) + return HUGE_VAL; + + return x > y ? x - y : 0.0; +} + +#endif /* _DOUBLE_IS_32BITS */ |