t-prec.cc revision 1.1.1.2 
1/* Test precision of mpf_class expressions.
2
3Copyright 2001, 2002, 2003, 2008 Free Software Foundation, Inc.
4
5This file is part of the GNU MP Library test suite.
6
7The GNU MP Library test suite is free software; you can redistribute it
8and/or modify it under the terms of the GNU General Public License as
9published by the Free Software Foundation; either version 3 of the License,
10or (at your option) any later version.
11
12The GNU MP Library test suite is distributed in the hope that it will be
13useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
14MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
15Public License for more details.
16
17You should have received a copy of the GNU General Public License along with
18the GNU MP Library test suite.  If not, see http://www.gnu.org/licenses/.  */
19
20#include "config.h"
21
22#include <iostream>
23
24#include "gmp.h"
25#include "gmpxx.h"
26#include "gmp-impl.h"
27#include "tests.h"
28
29using namespace std;
30
31
32const int
33small_prec = 64, medium_prec = 128, large_prec = 192, very_large_prec = 256;
34
35#define ASSERT_ALWAYS_PREC(a, s, prec) \
36{                                      \
37  mpf_srcptr _a = a.get_mpf_t();       \
38  mpf_class _b(s, prec);               \
39  mpf_srcptr _c = _b.get_mpf_t();      \
40  ASSERT_ALWAYS(mpf_eq(_a, _c, prec)); \
41}
42
43
44
45void
46check_mpf (void)
47{
48  mpf_set_default_prec(medium_prec);
49
50  // simple expressions
51  {
52    mpf_class f(3.0, small_prec);
53    mpf_class g(1 / f, very_large_prec);
54    ASSERT_ALWAYS_PREC
55      (g, "0.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
56       "     33333 33333 33333 33333 33333 333", very_large_prec);
57  }
58  {
59    mpf_class f(9.0, medium_prec);
60    mpf_class g(0.0, very_large_prec);
61    g = 1 / f;
62    ASSERT_ALWAYS_PREC
63      (g, "0.11111 11111 11111 11111 11111 11111 11111 11111 11111 11111"
64       "     11111 11111 11111 11111 11111 111", very_large_prec);
65  }
66  {
67    mpf_class f(15.0, large_prec);
68    mpf_class g(0.0, very_large_prec);
69    g = 1 / f;
70    ASSERT_ALWAYS_PREC
71      (g, "0.06666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
72       "     66666 66666 66666 66666 66666 667", very_large_prec);
73  }
74
75  // compound expressions
76  {
77    mpf_class f(3.0, small_prec);
78    mpf_class g(-(-(-1 / f)), very_large_prec);
79    ASSERT_ALWAYS_PREC
80      (g, "-0.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
81       "      33333 33333 33333 33333 33333 333", very_large_prec);
82  }
83  {
84    mpf_class f(3.0, small_prec), g(9.0, medium_prec);
85    mpf_class h(0.0, very_large_prec);
86    h = 1/f + 1/g;
87    ASSERT_ALWAYS_PREC
88      (h, "0.44444 44444 44444 44444 44444 44444 44444 44444 44444 44444"
89       "     44444 44444 44444 44444 44444 444", very_large_prec);
90  }
91  {
92    mpf_class f(3.0, small_prec), g(9.0, medium_prec), h(15.0, large_prec);
93    mpf_class i(0.0, very_large_prec);
94    i = f / g + h;
95    ASSERT_ALWAYS_PREC
96      (i, "15.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
97       "      33333 33333 33333 33333 33333 3", very_large_prec);
98  }
99  {
100    mpf_class f(3.0, small_prec);
101    mpf_class g(-(1 + f) / 3, very_large_prec);
102    ASSERT_ALWAYS_PREC
103      (g, "-1.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
104       "      33333 33333 33333 33333 33333 33", very_large_prec);
105  }
106  {
107    mpf_class f(9.0, medium_prec);
108    mpf_class g(0.0, very_large_prec);
109    g = sqrt(1 / f);
110    ASSERT_ALWAYS_PREC
111      (g, "0.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
112       "     33333 33333 33333 33333 33333 333", very_large_prec);
113  }
114  {
115    mpf_class f(15.0, large_prec);
116    mpf_class g(0.0, very_large_prec);
117    g = hypot(1 + 5 / f, 1.0);
118    ASSERT_ALWAYS_PREC
119      (g, "1.66666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
120       "     66666 66666 66666 66666 66666 67", very_large_prec);
121  }
122
123  // compound assignments
124  {
125    mpf_class f(3.0, small_prec), g(9.0, medium_prec);
126    mpf_class h(1.0, very_large_prec);
127    h -= f / g;
128    ASSERT_ALWAYS_PREC
129      (h, "0.66666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
130       "     66666 66666 66666 66666 66666 667", very_large_prec);
131  }
132
133  // construction from expressions
134  {
135    mpf_class f(3.0, small_prec);
136    mpf_class g(0.0, very_large_prec);
137    g = mpf_class(1 / f);
138    ASSERT_ALWAYS_PREC(g, "0.33333 33333 33333 33333", small_prec);
139  }
140  {
141    mpf_class f(9.0, medium_prec);
142    mpf_class g(0.0, very_large_prec);
143    g = mpf_class(1 / f);
144    ASSERT_ALWAYS_PREC
145      (g, "0.11111 11111 11111 11111 11111 11111 11111 1111", medium_prec);
146  }
147  {
148    mpf_class f(15.0, large_prec);
149    mpf_class g(0.0, very_large_prec);
150    g = mpf_class(1 / f);
151    ASSERT_ALWAYS_PREC
152      (g, "0.06666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
153       "     66666 6667", large_prec);
154  }
155
156  {
157    mpf_class f(3.0, small_prec), g(9.0, medium_prec);
158    mpf_class h(0.0, very_large_prec);
159    h = mpf_class(f / g + 1, large_prec);
160    ASSERT_ALWAYS_PREC
161      (h, "1.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
162       "     33333 333",
163       large_prec);
164  }
165
166  // mixed mpf/mpq expressions
167  {
168    mpf_class f(3.0, small_prec);
169    mpq_class q(1, 3);
170    mpf_class g(0.0, very_large_prec);
171    g = f - q;
172    ASSERT_ALWAYS_PREC
173      (g, "2.66666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
174       "     66666 66666 66666 66666 66666 67", very_large_prec);
175  }
176
177  {
178    mpf_class f(3.0, small_prec);
179    mpq_class q(1, 3);
180    mpf_class g(0.0, very_large_prec);
181    g = mpf_class(f - q, large_prec);
182    ASSERT_ALWAYS_PREC
183      (g, "2.66666 66666 66666 66666 66666 66666 66666 66666 66666 66666"
184       "     66666 667",
185       large_prec);
186  }
187  {
188    mpf_class f(3.0, small_prec);
189    mpq_class q(1, 3);
190    mpf_class g(0.0, very_large_prec);
191    g = mpf_class(f - q);
192    ASSERT_ALWAYS_PREC
193      (g, "2.66666 66666 66666 66666 66666 66666 66666 667", medium_prec);
194  }
195  {
196    mpf_class f(15.0, large_prec);
197    mpq_class q(1, 3);
198    mpf_class g(0.0, very_large_prec);
199    g = mpf_class(f + q);
200    ASSERT_ALWAYS_PREC
201      (g, "15.33333 33333 33333 33333 33333 33333 33333 33333 33333 33333"
202       "      33333 33",
203       large_prec);
204  }
205}
206
207
208int
209main (void)
210{
211  tests_start();
212
213  check_mpf();
214
215  tests_end();
216  return 0;
217}
218