fpmath 0.1.1

A floating point math library
Documentation 
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

use super::{F64Like, LikeF64};

// Generated with `./run-generator.sh f64::atan::consts`
const FRAC_PI_2_HI: u64 = 0x3FF921FB54442D18; // 1.5707963267948966e0
const FRAC_PI_2_LO: u64 = 0x3C91A62633145C07; // 6.123233995736766e-17
const FRAC_3PI_4: u64 = 0x4002D97C7F3321D2; // 2.356194490192345e0

impl<F: F64Like> crate::generic::Atan<LikeF64> for F {
    #[inline]
    fn frac_pi_2_hi() -> Self {
        Self::from_raw(FRAC_PI_2_HI)
    }

    #[inline]
    fn frac_pi_2_lo() -> Self {
        Self::from_raw(FRAC_PI_2_LO)
    }

    #[inline]
    fn frac_3pi_4() -> Self {
        Self::from_raw(FRAC_3PI_4)
    }

    #[inline]
    fn atan_poly(x2: Self) -> (Self, Self) {
        // Generated with `./run-generator.sh f64::atan::atan_poly`
        const K3: u64 = 0xBFD55555555554A8; // -3.333333333333237e-1
        const K5: u64 = 0x3FC999999998DDDB; // 1.99999999998666e-1
        const K7: u64 = 0xBFC2492492222261; // -1.428571427861867e-1
        const K9: u64 = 0x3FBC71C713BABB10; // 1.1111110908198563e-1
        const K11: u64 = 0xBFB745D0D8DBF8A8; // -9.090905470290755e-2
        const K13: u64 = 0x3FB3B133B63E4546; // 7.692263793275797e-2
        const K15: u64 = 0xBFB110D0F16E6DDA; // -6.666284460373309e-2
        const K17: u64 = 0x3FAE1ADCFBD54AF4; // 5.87987000859779e-2
        const K19: u64 = 0xBFAAE24E0347F882; // -5.2507818127183206e-2
        const K21: u64 = 0x3FA822298173BE98; // 4.713563637949142e-2
        const K23: u64 = 0xBFA57D989E95BD38; // -4.1973847731123504e-2
        const K25: u64 = 0x3FA28B28F449C44D; // 3.621795637162304e-2
        const K27: u64 = 0xBF9DF6081CE612E8; // -2.9258848918760544e-2
        const K29: u64 = 0x3F95BAFB6FDFA68F; // 2.122109289476853e-2
        const K31: u64 = 0xBF8B0B07B02A8E04; // -1.3204631866137191e-2
        const K33: u64 = 0x3F7B84A0C8C39E9E; // 6.718280852716111e-3
        const K35: u64 = 0xBF65ABAC56B5B133; // -2.6453367224147594e-3
        const K37: u64 = 0x3F48893B30409082; // 7.487811094546578e-4
        const K39: u64 = 0xBF21A73A6FE96DF7; // -1.346834978502712e-4
        const K41: u64 = 0x3EE823E8FE5305A0; // 1.151097962759606e-5

        let k3 = Self::from_raw(K3);
        let k5 = Self::from_raw(K5);
        let k7 = Self::from_raw(K7);
        let k9 = Self::from_raw(K9);
        let k11 = Self::from_raw(K11);
        let k13 = Self::from_raw(K13);
        let k15 = Self::from_raw(K15);
        let k17 = Self::from_raw(K17);
        let k19 = Self::from_raw(K19);
        let k21 = Self::from_raw(K21);
        let k23 = Self::from_raw(K23);
        let k25 = Self::from_raw(K25);
        let k27 = Self::from_raw(K27);
        let k29 = Self::from_raw(K29);
        let k31 = Self::from_raw(K31);
        let k33 = Self::from_raw(K33);
        let k35 = Self::from_raw(K35);
        let k37 = Self::from_raw(K37);
        let k39 = Self::from_raw(K39);
        let k41 = Self::from_raw(K41);

        let t = horner!(
            x2,
            x2,
            [
                k5, k7, k9, k11, k13, k15, k17, k19, k21, k23, k25, k27, k29, k31, k33, k35, k37,
                k39, k41
            ]
        );
        (k3, t)
    }
}