]>
git.proxmox.com Git - ceph.git/blob - ceph/src/boost/libs/multiprecision/test/test_mpfr_mpc_precisions.cpp
1 ///////////////////////////////////////////////////////////////
2 // Copyright Christopher Kormanyos 2002 - 2011.
3 // Copyright 2011 John Maddock. Distributed under the Boost
4 // Software License, Version 1.0. (See accompanying file
5 // LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt
7 // This work is based on an earlier work:
8 // "Algorithm 910: A Portable C++ Multiple-Precision System for Special-Function Calculations",
9 // in ACM TOMS, {VOL 37, ISSUE 4, (February 2011)} (C) ACM, 2011. http://doi.acm.org/10.1145/1916461.1916469
12 #define _SCL_SECURE_NO_WARNINGS
15 #include <boost/detail/lightweight_test.hpp>
18 #include <boost/multiprecision/mpfr.hpp>
19 #include <boost/multiprecision/mpc.hpp>
22 T
make_rvalue_copy(const T a
)
29 using namespace boost::multiprecision
;
31 // Test change of default precision:
33 mpfr_float::default_precision(100);
35 BOOST_CHECK_EQUAL(a
.precision(), 100);
36 mpfr_float::default_precision(20);
38 // test assignment from lvalue:
40 BOOST_CHECK_EQUAL(b
.precision(), 20);
42 BOOST_CHECK_EQUAL(b
.precision(), a
.precision());
45 // test assignment from rvalue:
47 BOOST_CHECK_EQUAL(b
.precision(), 20);
48 b
= make_rvalue_copy(a
);
49 BOOST_CHECK_EQUAL(b
.precision(), a
.precision());
51 mpfr_float::default_precision(20);
53 // test construct from lvalue:
55 BOOST_CHECK_EQUAL(b
.precision(), 100);
58 // test construct from rvalue:
59 mpfr_float
b(make_rvalue_copy(a
));
60 BOOST_CHECK_EQUAL(b
.precision(), 100);
62 mpc_complex::default_precision(100);
63 mpc_complex
ca("0.1");
64 BOOST_CHECK_EQUAL(ca
.precision(), 100);
65 mpc_complex::default_precision(20);
67 // test assignment from lvalue:
69 BOOST_CHECK_EQUAL(b
.precision(), 20);
71 BOOST_CHECK_EQUAL(b
.precision(), ca
.precision());
74 // test assignment from rvalue:
76 BOOST_CHECK_EQUAL(b
.precision(), 20);
77 b
= make_rvalue_copy(ca
);
78 BOOST_CHECK_EQUAL(b
.precision(), ca
.precision());
81 // test construct from lvalue:
83 BOOST_CHECK_EQUAL(b
.precision(), ca
.precision());
86 // test construct from rvalue:
87 mpc_complex
b(make_rvalue_copy(ca
));
88 BOOST_CHECK_EQUAL(b
.precision(), ca
.precision());
90 // real and imaginary:
91 BOOST_CHECK_EQUAL(ca
.real().precision(), 100);
92 BOOST_CHECK_EQUAL(ca
.imag().precision(), 100);
93 BOOST_CHECK_EQUAL(real(ca
).precision(), 100);
94 BOOST_CHECK_EQUAL(imag(ca
).precision(), 100);
97 // Construction at specific precision:
100 mpfr_float
f150(mpfr_float(), 150u);
101 BOOST_CHECK_EQUAL(f150
.precision(), 150);
102 mpc_complex
f150c(mpc_complex(), 150u);
103 BOOST_CHECK_EQUAL(f150c
.precision(), 150);
104 mpc_complex
f150cc(mpfr_float(), mpfr_float(), 150u);
105 BOOST_CHECK_EQUAL(f150cc
.precision(), 150);
108 mpfr_float
f150(2, 150);
109 BOOST_CHECK_EQUAL(f150
.precision(), 150);
112 mpfr_float
f150("1.2", 150);
113 BOOST_CHECK_EQUAL(f150
.precision(), 150);
116 // Copying precision:
119 mpc_complex
c(ca
.backend().data());
120 BOOST_CHECK_EQUAL(c
.precision(), 100);
121 mpc_complex_100
c100(2);
123 BOOST_CHECK_EQUAL(d
.precision(), 100);
124 mpfr_float_100
f100(2);
126 BOOST_CHECK_EQUAL(d
.precision(), 100);
129 // Check that the overloads for precision don't mess up 2-arg
133 mpc_complex
c(2, 3u);
134 BOOST_CHECK_EQUAL(c
.real(), 2);
135 BOOST_CHECK_EQUAL(c
.imag(), 3);
138 // 3-arg complex number construction with 3rd arg a precision:
141 mpc_complex
c(2, 3, 100);
142 BOOST_CHECK_EQUAL(c
.precision(), 100);
143 mpfr_float_50
x(2), y(3);
144 mpc_complex
z(x
, y
, 100);
145 BOOST_CHECK_EQUAL(c
.precision(), 100);
148 // From https://github.com/boostorg/multiprecision/issues/65
153 BOOST_CHECK_EQUAL(a
, 2);
154 BOOST_CHECK_EQUAL(a
.precision(), 100);
159 BOOST_CHECK_EQUAL(a
.real(), 2);
160 BOOST_CHECK_EQUAL(a
.imag(), 3);
161 BOOST_CHECK_EQUAL(a
.precision(), 100);
164 mpc_complex::default_precision(1000);
165 mpfr_float::default_precision(1000);
166 mpc_complex
a("1.324719827394086120398419082734980126734089612309871092830981236748901273498071240986123094861246981263481263489016238947147129807419028748901273409127349087124612576129076541203975704195690418570914657910465091256016501650916509165097164509164509761409561097561097650791650971465097165097162059761209561029756019265019726509126509172650971625097162450971309756104975610274650917825018740981274098127409182375701465172340923847120836540491320467127043127893281461230951097260126309812374091265091824981231236409851274",
167 "-0.80743891267394610982659071452346156102764312401571972642394120395608291471029347812645125986123123904123471209381289471230512983491286102875870192091283712396550981723409812740981263471230498715096104897123094710923879065981740928740981271801391209238470129560941870129387409812883437894183883841283700483832883218128438938184289148239164079329657861209381892037483468937489237419236509823723705612893489712412306531274812364980127304981648712483248732");
168 mpc_complex::default_precision(40);
169 mpfr_float::default_precision(40);
170 BOOST_CHECK_EQUAL(a
, a
);
174 // string_view with explicit precision:
176 #ifndef BOOST_NO_CXX17_HDR_STRING_VIEW
178 std::string
s("222");
179 std::string_view
v(s
.c_str(), 1);
180 mpfr_float
f(v
, 100);
181 BOOST_CHECK_EQUAL(f
, 2);
182 BOOST_CHECK_EQUAL(f
.precision(), 100);
185 std::string
x("222"), y("333");
186 std::string_view
vx(x
.c_str(), 1), vy(y
.c_str(), 1);
187 mpc_complex
c(vx
, vy
, 100);
188 BOOST_CHECK_EQUAL(c
.real(), 2);
189 BOOST_CHECK_EQUAL(c
.imag(), 3);
190 BOOST_CHECK_EQUAL(c
.precision(), 100);
194 mpc_complex::default_precision(100);
195 mpfr_float::default_precision(100);
199 mpc_complex::default_precision(50);
200 mpfr_float::default_precision(50);
204 BOOST_CHECK_EQUAL(z
.precision(), 100);
208 mpfr_float
x(2, 100); // 100 digits precision.
209 mpfr_float
y(3, 50); // 50 digits precision.
211 BOOST_CHECK_EQUAL(x
, 3);
212 BOOST_CHECK_EQUAL(y
, 2);
213 BOOST_CHECK_EQUAL(x
.precision(), 50);
214 BOOST_CHECK_EQUAL(y
.precision(), 100);
216 BOOST_CHECK_EQUAL(x
, 2);
217 BOOST_CHECK_EQUAL(y
, 3);
218 BOOST_CHECK_EQUAL(x
.precision(), 100);
219 BOOST_CHECK_EQUAL(y
.precision(), 50);
220 x
= std::move(mpfr_float(y
));
221 BOOST_CHECK_EQUAL(x
, y
);
222 BOOST_CHECK_EQUAL(x
.precision(), y
.precision());
225 mpc_complex
x(2, 3, 100); // 100 digits precision.
226 mpc_complex
y(3, 4, 50); // 50 digits precision.
228 BOOST_CHECK_EQUAL(x
.real(), 3);
229 BOOST_CHECK_EQUAL(x
.imag(), 4);
230 BOOST_CHECK_EQUAL(y
.real(), 2);
231 BOOST_CHECK_EQUAL(y
.imag(), 3);
232 BOOST_CHECK_EQUAL(x
.precision(), 50);
233 BOOST_CHECK_EQUAL(y
.precision(), 100);
235 BOOST_CHECK_EQUAL(x
.real(), 2);
236 BOOST_CHECK_EQUAL(x
.imag(), 3);
237 BOOST_CHECK_EQUAL(y
.real(), 3);
238 BOOST_CHECK_EQUAL(y
.imag(), 4);
239 BOOST_CHECK_EQUAL(x
.precision(), 100);
240 BOOST_CHECK_EQUAL(y
.precision(), 50);
241 x
= std::move(mpc_complex(y
));
242 BOOST_CHECK_EQUAL(x
, y
);
243 BOOST_CHECK_EQUAL(x
.precision(), y
.precision());
246 mpfr_float
c(4), d(8), e(9), f
;
248 mpfr_float
g((c
+ d
) * d
/ e
);
251 mpfr_float::default_precision(100);
254 BOOST_CHECK_EQUAL(f1
.precision(), 100);
256 BOOST_CHECK_EQUAL(f1
.precision(), 100);
258 BOOST_CHECK_EQUAL(f2
.precision(), 100);
259 mpfr_float
f3("5.1");
260 BOOST_CHECK_EQUAL(f3
.precision(), 100);
262 mpfr_float::default_precision(50);
263 mpfr_float
f4(f3
, 50);
264 BOOST_CHECK_EQUAL(f4
.precision(), 50);
265 f4
.assign(f1
, f4
.precision());
266 BOOST_CHECK_EQUAL(f4
.precision(), 50);
270 // Overloads of Math lib functions, discovered while fixing
271 // https://github.com/boostorg/multiprecision/issues/91
273 mpfr_float::default_precision(100);
276 BOOST_CHECK_EQUAL(f1
.precision(), 100);
277 mpfr_float::default_precision(20);
278 BOOST_CHECK_EQUAL(asinh(f1
).precision(), 100);
279 BOOST_CHECK_EQUAL(acosh(f1
).precision(), 100);
280 BOOST_CHECK_EQUAL(atanh(f1
).precision(), 100);
281 BOOST_CHECK_EQUAL(cbrt(f1
).precision(), 100);
282 BOOST_CHECK_EQUAL(erf(f1
).precision(), 100);
283 BOOST_CHECK_EQUAL(erfc(f1
).precision(), 100);
284 BOOST_CHECK_EQUAL(expm1(f1
).precision(), 100);
285 BOOST_CHECK_EQUAL(lgamma(f1
).precision(), 100);
286 BOOST_CHECK_EQUAL(tgamma(f1
).precision(), 100);
287 BOOST_CHECK_EQUAL(log1p(f1
).precision(), 100);
290 return boost::report_errors();