1 /* tremquo -- test file for mpfr_remquo, mpfr_remainder and mpfr_fmodquo.
2
3 Copyright 2007-2023 Free Software Foundation, Inc.
4 Contributed by the AriC and Caramba projects, INRIA.
5
6 This file is part of the GNU MPFR Library.
7
8 The GNU MPFR Library is free software; you can redistribute it and/or modify
9 it under the terms of the GNU Lesser General Public License as published by
10 the Free Software Foundation; either version 3 of the License, or (at your
11 option) any later version.
12
13 The GNU MPFR Library is distributed in the hope that it will be useful, but
14 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
16 License for more details.
17
18 You should have received a copy of the GNU Lesser General Public License
19 along with the GNU MPFR Library; see the file COPYING.LESSER. If not, see
20 https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
21 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */
22
23 #include "mpfr-test.h"
24
25 static void
26 bug20090227 (void)
27 {
28 mpfr_t x, y, r1, r2;
29
30 mpfr_init2 (x, 118);
31 mpfr_init2 (y, 181);
32 mpfr_init2 (r1, 140);
33 mpfr_init2 (r2, 140);
34 mpfr_set_si (x, -1, MPFR_RNDN);
35 mpfr_set_str_binary (y, "1.100100100001111110110101010001000100001011010001100001000110100110001001100011001100010100010111000000011011100000111001101000100101001000000100100111000001000100010100110011111010");
36 mpfr_remainder (r1, x, y, MPFR_RNDU);
37 /* since the quotient is -1, r1 is the rounding of x+y */
38 mpfr_add (r2, x, y, MPFR_RNDU);
39 if (mpfr_cmp (r1, r2))
40 {
41 printf ("Error in mpfr_remainder (bug20090227)\n");
42 printf ("Expected ");
43 mpfr_dump (r2);
44 printf ("Got ");
45 mpfr_dump (r1);
46 exit (1);
47 }
48 mpfr_clear (x);
49 mpfr_clear (y);
50 mpfr_clear (r1);
51 mpfr_clear (r2);
52 }
53
54 int
55 main (int argc, char *argv[])
56 {
57 mpfr_t x, y, r;
58 long q[1];
59 int inex;
60
61 tests_start_mpfr ();
62
63 if (argc == 3) /* usage: tremquo x y (rnd=MPFR_RNDN implicit) */
64 {
65 mpfr_init2 (x, GMP_NUMB_BITS);
66 mpfr_init2 (y, GMP_NUMB_BITS);
67 mpfr_init2 (r, GMP_NUMB_BITS);
68 mpfr_set_str (x, argv[1], 10, MPFR_RNDN);
69 mpfr_set_str (y, argv[2], 10, MPFR_RNDN);
70 mpfr_remquo (r, q, x, y, MPFR_RNDN);
71 printf ("r=");
72 mpfr_out_str (stdout, 10, 0, r, MPFR_RNDN);
73 printf (" q=%ld\n", q[0]);
74 mpfr_clear (x);
75 mpfr_clear (y);
76 mpfr_clear (r);
77 return 0;
78 }
79
80 bug20090227 ();
81
82 mpfr_init (x);
83 mpfr_init (y);
84 mpfr_init (r);
85
86 /* special values */
87 mpfr_set_nan (x);
88 mpfr_set_ui (y, 1, MPFR_RNDN);
89 mpfr_remquo (r, q, x, y, MPFR_RNDN);
90 MPFR_ASSERTN(mpfr_nan_p (r));
91 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
92 MPFR_ASSERTN(mpfr_nan_p (r));
93
94 mpfr_set_ui (x, 1, MPFR_RNDN);
95 mpfr_set_nan (y);
96 mpfr_remquo (r, q, x, y, MPFR_RNDN);
97 MPFR_ASSERTN(mpfr_nan_p (r));
98 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
99 MPFR_ASSERTN(mpfr_nan_p (r));
100
101 mpfr_set_inf (x, 1); /* +Inf */
102 mpfr_set_ui (y, 1, MPFR_RNDN);
103 mpfr_remquo (r, q, x, y, MPFR_RNDN);
104 MPFR_ASSERTN (mpfr_nan_p (r));
105 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
106 MPFR_ASSERTN (mpfr_nan_p (r));
107
108 mpfr_set_inf (x, 1); /* +Inf */
109 mpfr_set_ui (y, 0, MPFR_RNDN);
110 mpfr_remquo (r, q, x, y, MPFR_RNDN);
111 MPFR_ASSERTN (mpfr_nan_p (r));
112 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
113 MPFR_ASSERTN (mpfr_nan_p (r));
114
115 mpfr_set_inf (x, 1); /* +Inf */
116 mpfr_set_inf (y, 1);
117 mpfr_remquo (r, q, x, y, MPFR_RNDN);
118 MPFR_ASSERTN (mpfr_nan_p (r));
119 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
120 MPFR_ASSERTN (mpfr_nan_p (r));
121
122 mpfr_set_ui (x, 0, MPFR_RNDN);
123 mpfr_set_inf (y, 1);
124 mpfr_remquo (r, q, x, y, MPFR_RNDN);
125 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_POS (r));
126 MPFR_ASSERTN (q[0] == (long) 0);
127 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
128 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_POS (r));
129 MPFR_ASSERTN (q[0] == (long) 0);
130
131 mpfr_set_ui (x, 0, MPFR_RNDN);
132 mpfr_neg (x, x, MPFR_RNDN); /* -0 */
133 mpfr_set_inf (y, 1);
134 mpfr_remquo (r, q, x, y, MPFR_RNDN);
135 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_NEG (r));
136 MPFR_ASSERTN (q[0] == (long) 0);
137 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
138 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_NEG (r));
139 MPFR_ASSERTN (q[0] == (long) 0);
140
141 mpfr_set_ui (x, 17, MPFR_RNDN);
142 mpfr_set_inf (y, 1);
143 mpfr_remquo (r, q, x, y, MPFR_RNDN);
144 MPFR_ASSERTN (mpfr_cmp (r, x) == 0);
145 MPFR_ASSERTN (q[0] == (long) 0);
146 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
147 MPFR_ASSERTN (mpfr_cmp (r, x) == 0);
148 MPFR_ASSERTN (q[0] == (long) 0);
149
150 mpfr_set_ui (x, 17, MPFR_RNDN);
151 mpfr_set_ui (y, 0, MPFR_RNDN);
152 mpfr_remquo (r, q, x, y, MPFR_RNDN);
153 MPFR_ASSERTN (mpfr_nan_p (r));
154 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
155 MPFR_ASSERTN (mpfr_nan_p (r));
156
157 mpfr_set_ui (x, 0, MPFR_RNDN);
158 mpfr_set_ui (y, 17, MPFR_RNDN);
159 mpfr_remquo (r, q, x, y, MPFR_RNDN);
160 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_POS (r));
161 MPFR_ASSERTN (q[0] == (long) 0);
162 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
163 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_POS (r));
164 MPFR_ASSERTN (q[0] == (long) 0);
165
166 mpfr_set_ui (x, 0, MPFR_RNDN);
167 mpfr_neg (x, x, MPFR_RNDN);
168 mpfr_set_ui (y, 17, MPFR_RNDN);
169 mpfr_remquo (r, q, x, y, MPFR_RNDN);
170 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_NEG (r));
171 MPFR_ASSERTN (q[0] == (long) 0);
172 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
173 MPFR_ASSERTN (mpfr_cmp_ui (r, 0) == 0 && MPFR_IS_NEG (r));
174 MPFR_ASSERTN (q[0] == (long) 0);
175
176 mpfr_set_prec (x, 53);
177 mpfr_set_prec (y, 53);
178
179 /* check four possible sign combinations for 42/17 */
180 mpfr_set_ui (x, 42, MPFR_RNDN);
181 mpfr_set_ui (y, 17, MPFR_RNDN);
182 mpfr_remquo (r, q, x, y, MPFR_RNDN);
183 MPFR_ASSERTN (mpfr_cmp_ui (r, 8) == 0);
184 MPFR_ASSERTN (q[0] == (long) 2);
185 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
186 MPFR_ASSERTN (mpfr_cmp_ui (r, 8) == 0);
187 MPFR_ASSERTN (q[0] == (long) 2);
188
189 mpfr_set_si (x, -42, MPFR_RNDN);
190 mpfr_set_ui (y, 17, MPFR_RNDN);
191 mpfr_remquo (r, q, x, y, MPFR_RNDN);
192 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
193 MPFR_ASSERTN (q[0] == (long) -2);
194 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
195 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
196 MPFR_ASSERTN (q[0] == (long) -2);
197
198 mpfr_set_si (x, -42, MPFR_RNDN);
199 mpfr_set_si (y, -17, MPFR_RNDN);
200 mpfr_remquo (r, q, x, y, MPFR_RNDN);
201 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
202 MPFR_ASSERTN (q[0] == (long) 2);
203 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
204 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
205 MPFR_ASSERTN (q[0] == (long) 2);
206
207 mpfr_set_ui (x, 42, MPFR_RNDN);
208 mpfr_set_si (y, -17, MPFR_RNDN);
209 mpfr_remquo (r, q, x, y, MPFR_RNDN);
210 MPFR_ASSERTN (mpfr_cmp_ui (r, 8) == 0);
211 MPFR_ASSERTN (q[0] == (long) -2);
212 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
213 MPFR_ASSERTN (mpfr_cmp_ui (r, 8) == 0);
214 MPFR_ASSERTN (q[0] == (long) -2);
215
216 /* same tests for 43/17, rounded to 3 to nearest, and to 2 to zero */
217 mpfr_set_ui (x, 43, MPFR_RNDN);
218 mpfr_set_ui (y, 17, MPFR_RNDN);
219 mpfr_remquo (r, q, x, y, MPFR_RNDN);
220 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
221 MPFR_ASSERTN (q[0] == (long) 3);
222 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
223 MPFR_ASSERTN (mpfr_cmp_ui (r, 9) == 0);
224 MPFR_ASSERTN (q[0] == (long) 2);
225
226 mpfr_set_si (x, -43, MPFR_RNDN);
227 mpfr_set_ui (y, 17, MPFR_RNDN);
228 mpfr_remquo (r, q, x, y, MPFR_RNDN);
229 MPFR_ASSERTN (mpfr_cmp_si (r, 8) == 0);
230 MPFR_ASSERTN (q[0] == (long) -3);
231 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
232 MPFR_ASSERTN (mpfr_cmp_si (r, -9) == 0);
233 MPFR_ASSERTN (q[0] == (long) -2);
234
235 mpfr_set_si (x, -43, MPFR_RNDN);
236 mpfr_set_si (y, -17, MPFR_RNDN);
237 mpfr_remquo (r, q, x, y, MPFR_RNDN);
238 MPFR_ASSERTN (mpfr_cmp_si (r, 8) == 0);
239 MPFR_ASSERTN (q[0] == (long) 3);
240 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
241 MPFR_ASSERTN (mpfr_cmp_si (r, -9) == 0);
242 MPFR_ASSERTN (q[0] == (long) 2);
243
244 mpfr_set_ui (x, 43, MPFR_RNDN);
245 mpfr_set_si (y, -17, MPFR_RNDN);
246 mpfr_remquo (r, q, x, y, MPFR_RNDN);
247 MPFR_ASSERTN (mpfr_cmp_si (r, -8) == 0);
248 MPFR_ASSERTN (q[0] == (long) -3);
249 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
250 MPFR_ASSERTN (mpfr_cmp_ui (r, 9) == 0);
251 MPFR_ASSERTN (q[0] == (long) -2);
252
253 /* other tests */
254 mpfr_set_prec (x, 100);
255 mpfr_set_prec (y, 50);
256 mpfr_set_ui (x, 42, MPFR_RNDN);
257 mpfr_nextabove (x); /* 42 + 2^(-94) */
258 mpfr_set_ui (y, 21, MPFR_RNDN);
259 mpfr_remquo (r, q, x, y, MPFR_RNDN);
260 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -94) == 0);
261 MPFR_ASSERTN (q[0] == (long) 2);
262 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
263 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -94) == 0);
264 MPFR_ASSERTN (q[0] == (long) 2);
265
266 mpfr_set_prec (x, 50);
267 mpfr_set_prec (y, 100);
268 mpfr_set_ui (x, 42, MPFR_RNDN);
269 mpfr_nextabove (x); /* 42 + 2^(-44) */
270 mpfr_set_ui (y, 21, MPFR_RNDN);
271 mpfr_remquo (r, q, x, y, MPFR_RNDN);
272 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -44) == 0);
273 MPFR_ASSERTN (q[0] == (long) 2);
274 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
275 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -44) == 0);
276 MPFR_ASSERTN (q[0] == (long) 2);
277
278 mpfr_set_prec (x, 100);
279 mpfr_set_prec (y, 50);
280 mpfr_set_ui (x, 42, MPFR_RNDN);
281 mpfr_set_ui (y, 21, MPFR_RNDN);
282 mpfr_nextabove (y); /* 21 + 2^(-45) */
283 mpfr_remquo (r, q, x, y, MPFR_RNDN);
284 /* r should be 42 - 2*(21 + 2^(-45)) = -2^(-44) */
285 MPFR_ASSERTN (mpfr_cmp_si_2exp (r, -1, -44) == 0);
286 MPFR_ASSERTN (q[0] == (long) 2);
287 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
288 /* r should be 42 - (21 + 2^(-45) = 21 - 2^(-45) */
289 mpfr_sub_ui (r, r, 21, MPFR_RNDN);
290 MPFR_ASSERTN (mpfr_cmp_si_2exp (r, -1, -45) == 0);
291 MPFR_ASSERTN (q[0] == (long) 1);
292
293 mpfr_set_prec (x, 50);
294 mpfr_set_prec (y, 100);
295 mpfr_set_ui (x, 42, MPFR_RNDN);
296 mpfr_set_ui (y, 21, MPFR_RNDN);
297 mpfr_nextabove (y); /* 21 + 2^(-95) */
298 mpfr_remquo (r, q, x, y, MPFR_RNDN);
299 /* r should be 42 - 2*(21 + 2^(-95)) = -2^(-94) */
300 MPFR_ASSERTN (mpfr_cmp_si_2exp (r, -1, -94) == 0);
301 MPFR_ASSERTN (q[0] == (long) 2);
302 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
303 /* r should be 42 - (21 + 2^(-95) = 21 - 2^(-94), which rounded to 53 bits
304 should give 21 */
305 MPFR_ASSERTN (mpfr_cmp_ui (r, 21) == 0);
306 MPFR_ASSERTN (q[0] == (long) 1);
307
308 /* exercise large quotient */
309 mpfr_set_ui_2exp (x, 1, 65, MPFR_RNDN);
310 mpfr_set_ui (y, 1, MPFR_RNDN);
311 /* quotient is 2^65 */
312 mpfr_remquo (r, q, x, y, MPFR_RNDN);
313 MPFR_ASSERTN (mpfr_cmp_si (r, 0) == 0);
314 MPFR_ASSERTN (q[0] % 1073741824L == 0L);
315 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
316 MPFR_ASSERTN (mpfr_cmp_si (r, 0) == 0);
317 MPFR_ASSERTN (q[0] % 1073741824L == 0L);
318
319 /* another large quotient */
320 mpfr_set_prec (x, 65);
321 mpfr_set_prec (y, 65);
322 mpfr_const_pi (x, MPFR_RNDN);
323 mpfr_mul_2ui (x, x, 63, MPFR_RNDN);
324 mpfr_const_log2 (y, MPFR_RNDN);
325 mpfr_set_prec (r, 10);
326 mpfr_remquo (r, q, x, y, MPFR_RNDN);
327 /* q should be 41803643793084085130, r should be 605/2048 */
328 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 605, -11) == 0);
329 MPFR_ASSERTN ((q[0] > 0) && ((q[0] % 1073741824L) == 733836170L));
330 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
331 /* q should be 41803643793084085130, r should be 605/2048 */
332 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 605, -11) == 0);
333 MPFR_ASSERTN ((q[0] > 0) && ((q[0] % 1073741824L) == 733836170L));
334
335 /* check cases where quotient is 1.5 +/- eps */
336 mpfr_set_prec (x, 65);
337 mpfr_set_prec (y, 65);
338 mpfr_set_prec (r, 63);
339 mpfr_set_ui (x, 3, MPFR_RNDN);
340 mpfr_set_ui (y, 2, MPFR_RNDN);
341 mpfr_remquo (r, q, x, y, MPFR_RNDN);
342 /* x/y = 1.5, quotient should be 2 (even rule), remainder should be -1 */
343 MPFR_ASSERTN (mpfr_cmp_si (r, -1) == 0);
344 MPFR_ASSERTN (q[0] == 2L);
345 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
346 /* x/y = 1.5, quotient should be 1, remainder should be 1 */
347 MPFR_ASSERTN (mpfr_cmp_si (r, 1) == 0);
348 MPFR_ASSERTN (q[0] == 1L);
349
350 mpfr_set_ui (x, 3, MPFR_RNDN);
351 mpfr_nextabove (x); /* 3 + 2^(-63) */
352 mpfr_set_ui (y, 2, MPFR_RNDN);
353 mpfr_remquo (r, q, x, y, MPFR_RNDN);
354 /* x/y = 1.5 + 2^(-64), quo should be 2, r should be -1 + 2^(-63) */
355 MPFR_ASSERTN (mpfr_add_ui (r, r, 1, MPFR_RNDN) == 0);
356 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -63) == 0);
357 MPFR_ASSERTN (q[0] == 2L);
358 mpfr_set_prec (r, 64);
359 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
360 /* x/y = 1.5 + 2^(-64), quo should be 1, r should be 1 + 2^(-63) */
361 MPFR_ASSERTN (mpfr_sub_ui (r, r, 1, MPFR_RNDN) == 0);
362 MPFR_ASSERTN (mpfr_cmp_ui_2exp (r, 1, -63) == 0);
363 MPFR_ASSERTN (q[0] == 1L);
364
365 mpfr_set_prec (r, 63);
366 mpfr_set_ui (x, 3, MPFR_RNDN);
367 mpfr_set_ui (y, 2, MPFR_RNDN);
368 mpfr_nextabove (y); /* 2 + 2^(-63) */
369 mpfr_remquo (r, q, x, y, MPFR_RNDN);
370 /* x/y = 1.5 - eps, quo should be 1, r should be 1 - 2^(-63) */
371 MPFR_ASSERTN (mpfr_sub_ui (r, r, 1, MPFR_RNDN) == 0);
372 MPFR_ASSERTN (mpfr_cmp_si_2exp (r, -1, -63) == 0);
373 MPFR_ASSERTN (q[0] == 1L);
374 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
375 /* x/y = 1.5 - eps, quo should be 1, r should be 1 - 2^(-63) */
376 MPFR_ASSERTN (mpfr_sub_ui (r, r, 1, MPFR_RNDN) == 0);
377 MPFR_ASSERTN (mpfr_cmp_si_2exp (r, -1, -63) == 0);
378 MPFR_ASSERTN (q[0] == 1L);
379
380 /* bug founds by Kaveh Ghazi, 3 May 2007 */
381 mpfr_set_ui (x, 2, MPFR_RNDN);
382 mpfr_set_ui (y, 3, MPFR_RNDN);
383 /* quotient rounded to nearest is 1, thus remainder is -1 */
384 mpfr_remainder (r, x, y, MPFR_RNDN);
385 MPFR_ASSERTN (mpfr_cmp_si (r, -1) == 0);
386
387 mpfr_set_si (x, -1, MPFR_RNDN);
388 mpfr_set_ui (y, 1, MPFR_RNDN);
389 mpfr_remainder (r, x, y, MPFR_RNDN);
390 MPFR_ASSERTN (mpfr_cmp_si (r, 0) == 0 && MPFR_IS_NEG (r));
391
392 /* check argument reuse */
393 mpfr_set_si (x, -1, MPFR_RNDN);
394 mpfr_set_ui (y, 1, MPFR_RNDN);
395 mpfr_remainder (x, x, y, MPFR_RNDN);
396 MPFR_ASSERTN (mpfr_cmp_si (x, 0) == 0 && MPFR_IS_NEG (x));
397
398 mpfr_set_ui_2exp (x, 1, mpfr_get_emax () - 1, MPFR_RNDN);
399 mpfr_set_ui_2exp (y, 1, mpfr_get_emin (), MPFR_RNDN);
400 mpfr_remquo (r, q, x, y, MPFR_RNDN);
401 MPFR_ASSERTN (mpfr_zero_p (r) && MPFR_IS_POS (r));
402 MPFR_ASSERTN (q[0] == 0);
403 mpfr_fmodquo (r, q, x, y, MPFR_RNDN);
404 MPFR_ASSERTN (mpfr_zero_p (r) && MPFR_IS_POS (r));
405 MPFR_ASSERTN (q[0] == 0);
406
407 mpfr_set_prec (x, 380);
408 mpfr_set_prec (y, 385);
409 mpfr_set_str_binary (x, "0.11011010010110011101011000100100101100101011010001011100110001100101111001010100001011111110111100101110101010110011010101000100000100011101101100001011101110100111101111111010001001000010000110010110011100111000001110111010000100101001010111100100010001101001110100011110010000000001110001111001101100111011001000110110011100100011111110010100011001000001001011010111010000000000E-2");
410 mpfr_set_str_binary (y, "0.1100011000011101011010001100010111001110110111001101010010111100111100011010010011011101111101111001010111111110001001100001111101001000000010100101111001001110010110000111001000101010111001001000100101011111000010100110001111000110011011010101111101100110010101011010011101100001011101001000101111110110110110000001001101110111110110111110111111001001011110001110011111100000000000000E-1");
411 mpfr_set_prec (r, 2);
412 inex = mpfr_remainder (r, x, y, MPFR_RNDA);
413 MPFR_ASSERTN(mpfr_cmp_si_2exp (r, -3, -4) == 0);
414 MPFR_ASSERTN(inex < 0);
415
416 mpfr_clear (x);
417 mpfr_clear (y);
418 mpfr_clear (r);
419
420 tests_end_mpfr ();
421
422 return 0;
423 }