diff options
author | Christopher Faylor <cgf@redhat.com> | 2000-02-17 19:39:46 +0000 |
---|---|---|
committer | Christopher Faylor <cgf@redhat.com> | 2000-02-17 19:39:46 +0000 |
commit | 6d491d8ad9122891d1a48447e8c49702eb678921 (patch) | |
tree | 45619a131b9f211ceca9fce7a866328344bf00ec /newlib/libm/mathfp/s_fabs.c | |
parent | bd592806f45307757c24ec3f75bdefe13abfd51f (diff) | |
download | gdb-cvs/SNAPSHOT.tar.gz |
import newlib-2000-02-17 snapshotcvs/SNAPSHOT
Diffstat (limited to 'newlib/libm/mathfp/s_fabs.c')
-rw-r--r-- | newlib/libm/mathfp/s_fabs.c | 80 |
1 files changed, 80 insertions, 0 deletions
diff --git a/newlib/libm/mathfp/s_fabs.c b/newlib/libm/mathfp/s_fabs.c new file mode 100644 index 00000000000..9e1d75eb336 --- /dev/null +++ b/newlib/libm/mathfp/s_fabs.c @@ -0,0 +1,80 @@ + +/* @(#)z_fabs.c 1.0 98/08/13 */ + +/* +FUNCTION + <<fabs>>, <<fabsf>>---absolute value (magnitude) +INDEX + fabs +INDEX + fabsf + +ANSI_SYNOPSIS + #include <math.h> + double fabs(double <[x]>); + float fabsf(float <[x]>); + +TRAD_SYNOPSIS + #include <math.h> + double fabs(<[x]>) + double <[x]>; + + float fabsf(<[x]>) + float <[x]>; + +DESCRIPTION +<<fabs>> and <<fabsf>> calculate +@tex +$|x|$, +@end tex +the absolute value (magnitude) of the argument <[x]>, by direct +manipulation of the bit representation of <[x]>. + +RETURNS +The calculated value is returned. + +PORTABILITY +<<fabs>> is ANSI. +<<fabsf>> is an extension. + +*/ + +/****************************************************************** + * Floating-Point Absolute Value + * + * Input: + * x - floating-point number + * + * Output: + * absolute value of x + * + * Description: + * fabs computes the absolute value of a floating point number. + * + *****************************************************************/ + +#include "fdlibm.h" +#include "zmath.h" + +#ifndef _DOUBLE_IS_32BITS + +double +_DEFUN (fabs, (double), + double x) +{ + switch (numtest (x)) + { + case NAN: + errno = EDOM; + return (x); + case INF: + errno = ERANGE; + return (x); + case 0: + return (0.0); + default: + return (x < 0.0 ? -x : x); + } +} + +#endif /* _DOUBLE_IS_32BITS */ |