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

use super::{F64Like, LikeF64};

impl<F: F64Like> crate::generic::Tan<LikeF64> for F {
    #[inline]
    fn tan_poly(x2: Self, x3: Self) -> Self {
        // Generated with `./run-generator.sh f64::tan::tan_poly`
        const K3: u64 = 0x3FD5555555555575; // 3.333333333333351e-1
        const K5: u64 = 0x3FC111111110D100; // 1.333333333328781e-1
        const K7: u64 = 0x3FABA1BA1BFA9D6E; // 5.396825400867554e-2
        const K9: u64 = 0x3F9664F469F9BCB5; // 2.1869486778481845e-2
        const K11: u64 = 0x3F8226E9208F9C1A; // 8.863278680975074e-3
        const K13: u64 = 0x3F6D6BE6B0D57F2A; // 3.5914903334485397e-3
        const K15: u64 = 0x3F57F26B97303CD6; // 1.46160609142836e-3
        const K17: u64 = 0x3F4250936AB460E5; // 5.589217957178398e-4
        const K19: u64 = 0x3F35902B995E4FFF; // 3.29027790365588e-4

        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);

        horner!(x3, x2, [k3, k5, k7, k9, k11, k13, k15, k17, k19])
    }
}