diff options
author | Jakub Jelinek <jj@ultra.linux.cz> | 2011-01-16 16:44:35 +0000 |
---|---|---|
committer | Jakub Jelinek <jj@ultra.linux.cz> | 2011-01-16 16:44:35 +0000 |
commit | 90df67fd97d65ca9739a483499c8979d79f0dc99 (patch) | |
tree | cbbd98da06ef06f334e26d4cd08eea03e60befdb /libquadmath/math/fdimq.c | |
parent | faf99708104ef936eecaaa20e221eb54b20ce5b5 (diff) |
PR fortran/46416
* quadmath.h (cbrtq, finiteq, isnanq, signbitq, sqrtq): Remove
const from prototype argument.
(cimagq, conjq, cprojq, crealq, fdimq, fmaxq, fminq, ilogbq,
llrintq, log2q, lrintq, nearbyintq, remquoq): New prototypes.
(__quadmath_extern_inline): Define.
(cimagq, conjq, crealq): New inlines.
* Makefile.am (libquadmath_la_SOURCES): Add math/cimagq.c,
math/conjq.c, math/cprojq.c, math/crealq.c, math/fdimq.c,
math/fmaxq.c, math/fminq.c, math/ilogbq.c, math/llrintq.c,
math/log2q.c, math/lrintq.c, math/nearbyintq.c and math/remquoq.c.
* Makefile.in: Regenerated.
* quadmath_weak.h (cimagq, conjq, cprojq, crealq, fdimq, fmaxq,
fminq, ilogbq, llrintq, log2q, lrintq, nearbyintq, remquoq): Add.
* quadmath-imp.h (__LITTLE_ENDIAN__): Don't define.
(ieee854_float128): Use __BYTE_ORDER == __ORDER_BIG_ENDIAN__
tests instead of __BIG_ENDIAN__.
* quadmath.map (QUADMATH_1.0): Add cimagq, conjq, cprojq, crealq,
fdimq, fmaxq, fminq, ilogbq, llrintq, log2q, lrintq, nearbyintq
and remquoq.
* libquadmath.texi (cimagq, conjq, cprojq, crealq, fdimq, fmaxq,
fminq, ilogbq, llrintq, log2q, lrintq, nearbyintq, remquoq): Add.
* math/cprojq.c: New file.
* math/ilogbq.c: New file.
* math/fminq.c: New file.
* math/llrintq.c: New file.
* math/log2q.c: New file.
* math/lrintq.c: New file.
* math/crealq.c: New file.
* math/nearbyintq.c: New file.
* math/fmaxq.c: New file.
* math/conjq.c: New file.
* math/remquoq.c: New file.
* math/cimagq.c: New file.
* math/fdimq.c: New file.
* math/ldexpq.c: Include errno.h. Set errno to ERANGE if needed.
git-svn-id: https://gcc.gnu.org/svn/gcc/trunk@168854 138bc75d-0d04-0410-961f-82ee72b054a4
Diffstat (limited to 'libquadmath/math/fdimq.c')
-rw-r--r-- | libquadmath/math/fdimq.c | 43 |
1 files changed, 43 insertions, 0 deletions
diff --git a/libquadmath/math/fdimq.c b/libquadmath/math/fdimq.c new file mode 100644 index 00000000000..539fb08c641 --- /dev/null +++ b/libquadmath/math/fdimq.c @@ -0,0 +1,43 @@ +/* Return positive difference between arguments. + Copyright (C) 1997, 2004, 2009 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, write to the Free + Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA + 02111-1307 USA. */ + +#include <errno.h> +#include "quadmath-imp.h" + +__float128 +fdimq (__float128 x, __float128 y) +{ + int clsx = fpclassifyq (x); + int clsy = fpclassifyq (y); + + if (clsx == QUADFP_NAN || clsy == QUADFP_NAN + || (y < 0 && clsx == QUADFP_INFINITE && clsy == QUADFP_INFINITE)) + /* Raise invalid flag. */ + return x - y; + + if (x <= y) + return 0.0Q; + + __float128 r = x - y; + if (isinfq (r)) + errno = ERANGE; + + return r; +} |