Annotation of OpenXM_contrib/gmp/tests/mpf/t-trunc.c, Revision 1.1.1.1
1.1 ohara 1: /* Test mpf_trunc, mpf_ceil, mpf_floor.
2:
3: Copyright 2001, 2002 Free Software Foundation, Inc.
4:
5: This file is part of the GNU MP Library.
6:
7: The GNU MP Library is free software; you can redistribute it and/or modify
8: it under the terms of the GNU Lesser General Public License as published by
9: the Free Software Foundation; either version 2.1 of the License, or (at your
10: option) any later version.
11:
12: The GNU MP Library is distributed in the hope that it will be useful, but
13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15: License for more details.
16:
17: You should have received a copy of the GNU Lesser General Public License
18: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
20: MA 02111-1307, USA. */
21:
22: #include <stdio.h>
23: #include <stdlib.h>
24: #include "gmp.h"
25: #include "gmp-impl.h"
26: #include "tests.h"
27:
28:
29: void
30: check_print (mpf_srcptr src, mpf_srcptr got, mpf_srcptr want)
31: {
32: mp_trace_base = 16;
33: mpf_trace ("src ", src);
34: mpf_trace ("got ", got);
35: mpf_trace ("want", want);
36:
37: printf ("got size=%d exp=%ld\n", SIZ(got), EXP(got));
38: mpn_trace (" limbs=", PTR(got), ABSIZ(got));
39:
40: printf ("want size=%d exp=%ld\n", SIZ(want), EXP(want));
41: mpn_trace (" limbs=", PTR(want), ABSIZ(want));
42: }
43:
44: void
45: check_one (mpf_srcptr src, mpf_srcptr trunc, mpf_srcptr ceil, mpf_srcptr floor)
46: {
47: mpf_t got;
48:
49: mpf_init2 (got, mpf_get_prec (trunc));
50: ASSERT_ALWAYS (PREC(got) == PREC(trunc));
51: ASSERT_ALWAYS (PREC(got) == PREC(ceil));
52: ASSERT_ALWAYS (PREC(got) == PREC(floor));
53:
54: #define CHECK_SEP(name, fun, want) \
55: mpf_set_ui (got, 54321L); /* initial junk */ \
56: fun (got, src); \
57: MPF_CHECK_FORMAT (got); \
58: if (mpf_cmp (got, want) != 0) \
59: { \
60: printf ("%s wrong\n", name); \
61: check_print (src, got, want); \
62: abort (); \
63: }
64:
65: CHECK_SEP ("mpf_trunc", mpf_trunc, trunc);
66: CHECK_SEP ("mpf_ceil", mpf_ceil, ceil);
67: CHECK_SEP ("mpf_floor", mpf_floor, floor);
68:
69: #define CHECK_INPLACE(name, fun, want) \
70: mpf_set (got, src); \
71: fun (got, got); \
72: MPF_CHECK_FORMAT (got); \
73: if (mpf_cmp (got, want) != 0) \
74: { \
75: printf ("%s wrong\n", name); \
76: check_print (src, got, want); \
77: abort (); \
78: }
79:
80: CHECK_INPLACE ("mpf_trunc", mpf_trunc, trunc);
81:
82: /* Can't do these unconditionally in case truncation by mpf_set strips
83: some low non-zero limbs which would have rounded the result. */
84: if (ABSIZ(src) <= PREC(trunc)+1)
85: {
86: CHECK_INPLACE ("mpf_ceil", mpf_ceil, ceil);
87: CHECK_INPLACE ("mpf_floor", mpf_floor, floor);
88: }
89:
90: mpf_clear (got);
91: }
92:
93: void
94: check_all (mpf_ptr src, mpf_ptr trunc, mpf_ptr ceil, mpf_ptr floor)
95: {
96: /* some of these values are generated with direct field assignments */
97: MPF_CHECK_FORMAT (src);
98: MPF_CHECK_FORMAT (trunc);
99: MPF_CHECK_FORMAT (ceil);
100: MPF_CHECK_FORMAT (floor);
101:
102: check_one (src, trunc, ceil, floor);
103:
104: mpf_neg (src, src);
105: mpf_neg (trunc, trunc);
106: mpf_neg (ceil, ceil);
107: mpf_neg (floor, floor);
108: check_one (src, trunc, floor, ceil);
109: }
110:
111: void
112: check_various (void)
113: {
114: mpf_t src, trunc, ceil, floor;
115: int n, i;
116:
117: mpf_init2 (src, 512L);
118: mpf_init2 (trunc, 256L);
119: mpf_init2 (ceil, 256L);
120: mpf_init2 (floor, 256L);
121:
122: /* 0 */
123: mpf_set_ui (src, 0L);
124: mpf_set_ui (trunc, 0L);
125: mpf_set_ui (ceil, 0L);
126: mpf_set_ui (floor, 0L);
127: check_all (src, trunc, ceil, floor);
128:
129: /* 1 */
130: mpf_set_ui (src, 1L);
131: mpf_set_ui (trunc, 1L);
132: mpf_set_ui (ceil, 1L);
133: mpf_set_ui (floor, 1L);
134: check_all (src, trunc, ceil, floor);
135:
136: /* 2^1024 */
137: mpf_set_ui (src, 1L);
138: mpf_mul_2exp (src, src, 1024L);
139: mpf_set (trunc, src);
140: mpf_set (ceil, src);
141: mpf_set (floor, src);
142: check_all (src, trunc, ceil, floor);
143:
144: /* 1/2^1024, fraction only */
145: mpf_set_ui (src, 1L);
146: mpf_div_2exp (src, src, 1024L);
147: mpf_set_si (trunc, 0L);
148: mpf_set_si (ceil, 1L);
149: mpf_set_si (floor, 0L);
150: check_all (src, trunc, ceil, floor);
151:
152: /* 1/2 */
153: mpf_set_ui (src, 1L);
154: mpf_div_2exp (src, src, 1L);
155: mpf_set_si (trunc, 0L);
156: mpf_set_si (ceil, 1L);
157: mpf_set_si (floor, 0L);
158: check_all (src, trunc, ceil, floor);
159:
160: /* 123+1/2^64 */
161: mpf_set_ui (src, 1L);
162: mpf_div_2exp (src, src, 64L);
163: mpf_add_ui (src, src, 123L);
164: mpf_set_si (trunc, 123L);
165: mpf_set_si (ceil, 124L);
166: mpf_set_si (floor, 123L);
167: check_all (src, trunc, ceil, floor);
168:
169: /* integer of full prec+1 limbs, unchanged */
170: n = PREC(trunc)+1;
171: ASSERT_ALWAYS (n <= PREC(src)+1);
172: EXP(src) = n;
173: SIZ(src) = n;
174: for (i = 0; i < SIZ(src); i++)
175: PTR(src)[i] = i+100;
176: mpf_set (trunc, src);
177: mpf_set (ceil, src);
178: mpf_set (floor, src);
179: check_all (src, trunc, ceil, floor);
180:
181: /* full prec+1 limbs, 1 trimmed for integer */
182: n = PREC(trunc)+1;
183: ASSERT_ALWAYS (n <= PREC(src)+1);
184: EXP(src) = n-1;
185: SIZ(src) = n;
186: for (i = 0; i < SIZ(src); i++)
187: PTR(src)[i] = i+200;
188: EXP(trunc) = n-1;
189: SIZ(trunc) = n-1;
190: for (i = 0; i < SIZ(trunc); i++)
191: PTR(trunc)[i] = i+201;
192: mpf_set (floor, trunc);
193: mpf_add_ui (ceil, trunc, 1L);
194: check_all (src, trunc, ceil, floor);
195:
196: /* prec+3 limbs, 2 trimmed for size */
197: n = PREC(trunc)+3;
198: ASSERT_ALWAYS (n <= PREC(src)+1);
199: EXP(src) = n;
200: SIZ(src) = n;
201: for (i = 0; i < SIZ(src); i++)
202: PTR(src)[i] = i+300;
203: EXP(trunc) = n;
204: SIZ(trunc) = n-2;
205: for (i = 0; i < SIZ(trunc); i++)
206: PTR(trunc)[i] = i+302;
207: mpf_set (floor, trunc);
208: mpf_set (ceil, trunc);
209: PTR(ceil)[0]++;
210: check_all (src, trunc, ceil, floor);
211:
212: /* prec+4 limbs, 2 trimmed for size, 1 trimmed for integer */
213: n = PREC(trunc)+4;
214: ASSERT_ALWAYS (n <= PREC(src)+1);
215: EXP(src) = n-1;
216: SIZ(src) = n;
217: for (i = 0; i < SIZ(src); i++)
218: PTR(src)[i] = i+400;
219: EXP(trunc) = n-1;
220: SIZ(trunc) = n-3;
221: for (i = 0; i < SIZ(trunc); i++)
222: PTR(trunc)[i] = i+403;
223: mpf_set (floor, trunc);
224: mpf_set (ceil, trunc);
225: PTR(ceil)[0]++;
226: check_all (src, trunc, ceil, floor);
227:
228: /* F.F, carry out of ceil */
229: EXP(src) = 1;
230: SIZ(src) = 2;
231: PTR(src)[0] = GMP_NUMB_MAX;
232: PTR(src)[1] = GMP_NUMB_MAX;
233: EXP(trunc) = 1;
234: SIZ(trunc) = 1;
235: PTR(trunc)[0] = GMP_NUMB_MAX;
236: mpf_set (floor, trunc);
237: EXP(ceil) = 2;
238: SIZ(ceil) = 1;
239: PTR(ceil)[0] = 1;
240: check_all (src, trunc, ceil, floor);
241:
242: /* FF.F, carry out of ceil */
243: EXP(src) = 2;
244: SIZ(src) = 3;
245: PTR(src)[0] = GMP_NUMB_MAX;
246: PTR(src)[1] = GMP_NUMB_MAX;
247: PTR(src)[2] = GMP_NUMB_MAX;
248: EXP(trunc) = 2;
249: SIZ(trunc) = 2;
250: PTR(trunc)[0] = GMP_NUMB_MAX;
251: PTR(trunc)[1] = GMP_NUMB_MAX;
252: mpf_set (floor, trunc);
253: EXP(ceil) = 3;
254: SIZ(ceil) = 1;
255: PTR(ceil)[0] = 1;
256: check_all (src, trunc, ceil, floor);
257:
258: mpf_clear (src);
259: mpf_clear (trunc);
260: mpf_clear (ceil);
261: mpf_clear (floor);
262: }
263:
264: int
265: main (void)
266: {
267: tests_start ();
268:
269: check_various ();
270:
271: tests_end ();
272: exit (0);
273: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to t-trunc.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / tests / mpf