1 /* Compute remainder and a congruent to the quotient.
2 Copyright (C) 1997 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4 Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
6 The GNU C Library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2.1 of the License, or (at your option) any later version.
11 The GNU C Library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details.
16 You should have received a copy of the GNU Lesser General Public
17 License along with the GNU C Library; if not, write to the Free
18 Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
23 #include "math_private.h"
26 static const long double zero = 0.0;
30 __remquol (long double x, long double p, int *quo)
36 GET_LDOUBLE_WORDS (ex, hx, lx, x);
37 GET_LDOUBLE_WORDS (ep, hp, lp, p);
39 qs = (sx ^ (ep & 0x8000)) >> 15;
43 /* Purge off exception values. */
44 if ((ep | hp | lp) == 0)
45 return (x * p) / (x * p); /* p = 0 */
46 if ((ex == 0x7fff) /* x not finite */
47 || ((ep == 0x7fff) /* p is NaN */
49 return (x * p) / (x * p);
52 x = __ieee754_fmodl (x, 8 * p); /* now x < 8p */
54 if (((ex - ep) | (hx - hp) | (lx - lp)) == 0)
90 long double p_half = 0.5 * p;
103 *quo = qs ? -cquo : cquo;
109 weak_alias (__remquol, remquol)