1 /* Compute remainder and a congruent to the quotient.
2 Copyright (C) 1997-2023 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
21 #include <math_private.h>
22 #include <libm-alias-ldouble.h>
25 static const long double zero
= 0.0;
29 __remquol (long double x
, long double p
, int *quo
)
35 GET_LDOUBLE_WORDS (ex
, hx
, lx
, x
);
36 GET_LDOUBLE_WORDS (ep
, hp
, lp
, p
);
38 qs
= (sx
^ (ep
& 0x8000)) >> 15;
42 /* Purge off exception values. */
43 if ((ep
| hp
| lp
) == 0)
44 return (x
* p
) / (x
* p
); /* p = 0 */
45 if ((ex
== 0x7fff) /* x not finite */
46 || ((ep
== 0x7fff) /* p is NaN */
47 && (((hp
& 0x7fffffff) | lp
) != 0)))
48 return (x
* p
) / (x
* p
);
51 x
= __ieee754_fmodl (x
, 8 * p
); /* now x < 8p */
53 if (((ex
- ep
) | (hx
- hp
) | (lx
- lp
)) == 0)
63 if (ep
<= 0x7ffc && x
>= 4 * p
)
68 if (ep
<= 0x7ffd && x
>= 2 * p
)
89 long double p_half
= 0.5 * p
;
102 *quo
= qs
? -cquo
: cquo
;
104 /* Ensure correct sign of zero result in round-downward mode. */
111 libm_alias_ldouble (__remquo
, remquo
)