1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125

use crate::Float;
use crate::InnerFloat::Finite;
use malachite_base::num::arithmetic::traits::{PowerOf2, RoundToMultipleOfPowerOf2};
use malachite_base::num::basic::integers::PrimitiveInt;
use malachite_base::num::conversion::traits::ExactFrom;
use malachite_base::rounding_modes::RoundingMode;
use malachite_nz::natural::Natural;
use malachite_nz::platform::Limb;

impl Float {
    /// Raises 2 to an integer power, returning a [`Float`] with the specified precision.
    ///
    /// If you need a [`Float`] with precision 1, then the
    /// [`PowerOfTwo`](malachite_base::num::arithmetic::traits::PowerOf2) implementation may be
    /// used instead.
    ///
    /// $f(k) = 2^k$,
    ///
    /// and the result has precision `prec`.
    ///
    /// # Worst-case complexity
    /// $T(n) = O(n)$
    ///
    /// $M(n) = O(n)$
    ///
    /// where $T$ is time, $M$ is additional memory, and $n$ is `prec`.
    ///
    /// # Examples
    /// ```
    /// use malachite_float::Float;
    ///
    /// assert_eq!(Float::power_of_2_prec(0, 1).to_string(), "1.0");
    /// assert_eq!(Float::power_of_2_prec(0, 10).to_string(), "1.0");
    /// assert_eq!(Float::power_of_2_prec(0, 100).to_string(), "1.0");
    ///
    /// assert_eq!(Float::power_of_2_prec(100, 1).to_string(), "1.0e30");
    /// assert_eq!(Float::power_of_2_prec(100, 10).to_string(), "1.268e30");
    /// assert_eq!(
    ///     Float::power_of_2_prec(100, 100).to_string(),
    ///     "1267650600228229401496703205376.0"
    /// );
    ///
    /// assert_eq!(Float::power_of_2_prec(-100, 1).to_string(), "8.0e-31");
    /// assert_eq!(Float::power_of_2_prec(-100, 10).to_string(), "7.89e-31");
    /// assert_eq!(
    ///     Float::power_of_2_prec(-100, 100).to_string(),
    ///     "7.88860905221011805411728565283e-31"
    /// );
    /// ```
    pub fn power_of_2_prec(pow: i64, prec: u64) -> Float {
        assert_ne!(prec, 0);
        Float(Finite {
            sign: true,
            exponent: pow + 1,
            precision: prec,
            significand: Natural::power_of_2(
                prec.round_to_multiple_of_power_of_2(Limb::LOG_WIDTH, RoundingMode::Ceiling)
                    .0
                    - 1,
            ),
        })
    }
}

impl PowerOf2<u64> for Float {
    /// Raises 2 to an integer power, returning a [`Float`] with precision 1.
    ///
    /// To get a [`Float`] with a higher precision, try [`Float::power_of_2_prec`].
    ///
    /// $f(k) = 2^k$.
    ///
    /// # Worst-case complexity
    /// Constant time and additional memory.
    ///
    /// # Examples
    /// ```
    /// use malachite_base::num::arithmetic::traits::PowerOf2;
    /// use malachite_float::Float;
    ///
    /// assert_eq!(Float::power_of_2(0u64).to_string(), "1.0");
    /// assert_eq!(Float::power_of_2(3u64).to_string(), "8.0");
    /// assert_eq!(Float::power_of_2(100u64).to_string(), "1.0e30");
    /// ```
    #[inline]
    fn power_of_2(pow: u64) -> Float {
        Float(Finite {
            sign: true,
            exponent: i64::exact_from(pow + 1),
            precision: 1,
            significand: Natural::HIGH_BIT,
        })
    }
}

impl PowerOf2<i64> for Float {
    /// Raises 2 to an integer power, returning a [`Float`] with precision 1.
    ///
    /// To get a [`Float`] with a higher precision, try [`Float::power_of_2_prec`].
    ///
    /// $f(k) = 2^k$.
    ///
    /// # Worst-case complexity
    /// Constant time and additional memory.
    ///
    /// # Examples
    /// ```
    /// use malachite_base::num::arithmetic::traits::PowerOf2;
    /// use malachite_float::Float;
    ///
    /// assert_eq!(Float::power_of_2(0i64).to_string(), "1.0");
    /// assert_eq!(Float::power_of_2(3i64).to_string(), "8.0");
    /// assert_eq!(Float::power_of_2(100i64).to_string(), "1.0e30");
    /// assert_eq!(Float::power_of_2(-3i64).to_string(), "0.1");
    /// assert_eq!(Float::power_of_2(-100i64).to_string(), "8.0e-31");
    /// ```
    #[inline]
    fn power_of_2(pow: i64) -> Float {
        Float(Finite {
            sign: true,
            exponent: pow + 1,
            precision: 1,
            significand: Natural::HIGH_BIT,
        })
    }
}