1 /*++
2 Copyright (c) 2012 Microsoft Corporation
3
4 Module Name:
5
6 mpfx.cpp
7
8 Abstract:
9
10 Multi precision fixed point numbers.
11
12 Author:
13
14 Leonardo de Moura (leonardo) 2012-09-19
15
16 Revision History:
17
18 --*/
19 #include "util/mpfx.h"
20
tst1()21 static void tst1() {
22 mpfx_manager m;
23 scoped_mpfx a(m), b(m), c(m);
24 m.set(a, 1);
25 m.set(b, 2);
26 std::cout << a << " + " << b << " == " << (a+b) << "\n";
27 m.set(a, 5);
28 m.set(c, 3);
29 m.display_raw(std::cout, (a*a*b)/c); std::cout << "\n";
30 m.display_decimal(std::cout, (a*a*b)/c); std::cout << "\n";
31 m.display_decimal(std::cout, (a*a*b)/c, 10); std::cout << "\n";
32 m.round_to_plus_inf();
33 m.display_decimal(std::cout, (a*a*b)/c); std::cout << "\n";
34 m.set(a, -1, 4);
35 m.display_decimal(std::cout, a); std::cout << "\n";
36 }
37
tst_prev_power_2(int64_t n,uint64_t d,unsigned expected)38 static void tst_prev_power_2(int64_t n, uint64_t d, unsigned expected) {
39 mpfx_manager m;
40 scoped_mpfx a(m);
41 m.set(a, n, d);
42 ENSURE(m.prev_power_of_two(a) == expected);
43 }
44
tst_prev_power_2()45 static void tst_prev_power_2() {
46 tst_prev_power_2(-10, 1, 0);
47 tst_prev_power_2(0, 1, 0);
48 tst_prev_power_2(1, 1, 0);
49 tst_prev_power_2(2, 1, 1);
50 tst_prev_power_2(3, 1, 1);
51 tst_prev_power_2(4, 1, 2);
52 tst_prev_power_2(5, 1, 2);
53 tst_prev_power_2(8, 1, 3);
54 tst_prev_power_2(9, 1, 3);
55 tst_prev_power_2(9, 2, 2);
56 tst_prev_power_2(9, 4, 1);
57 tst_prev_power_2(9, 5, 0);
58 tst_prev_power_2((1ll << 60) + 1, 1, 60);
59 tst_prev_power_2((1ll << 60), 1, 60);
60 tst_prev_power_2((1ll << 60) - 1, 1, 59);
61 tst_prev_power_2((1ll << 60), 3, 58);
62 }
63
tst_mpfx()64 void tst_mpfx() {
65 tst_prev_power_2();
66 tst1();
67 }
68