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

use super::{F32Like, LikeF32};
use crate::double::SemiDouble;

// Generated with `./run-generator.sh f32::sin_cos::consts`
const FRAC_1_6_HI: u32 = 0x3E2AA000; // 1.6662598e-1
const FRAC_1_6_LO: u32 = 0x382AAAAB; // 4.0690105e-5

impl<F: F32Like> crate::generic::SinCos<LikeF32> for F {
    #[inline]
    fn frac_1_6_ex() -> SemiDouble<Self> {
        SemiDouble::with_parts(Self::from_raw(FRAC_1_6_HI), Self::from_raw(FRAC_1_6_LO))
    }

    #[inline]
    fn sin_poly(x2: Self, x5: Self) -> (Self, Self) {
        // Generated with `./run-generator.sh f32::sin_cos::sin_poly`
        const K3: u32 = 0xBE2AAAA3; // -1.6666655e-1
        const K5: u32 = 0x3C0883AC; // 8.332174e-3
        const K7: u32 = 0xB94CA607; // -1.9516806e-4

        let k3 = Self::from_raw(K3);
        let k5 = Self::from_raw(K5);
        let k7 = Self::from_raw(K7);

        let r = horner!(x5, x2, [k5, k7]);
        (r, k3)
    }

    #[inline]
    fn sin_poly_ex(x2: Self, x5: Self) -> Self {
        // Generated with `./run-generator.sh f32::sin_cos::sin_poly_ex`
        const K5: u32 = 0x3C088602; // 8.332731e-3
        const K7: u32 = 0xB94D49A3; // -1.9577755e-4

        let k5 = Self::from_raw(K5);
        let k7 = Self::from_raw(K7);

        horner!(x5, x2, [k5, k7])
    }

    #[inline]
    fn cos_poly(x2: Self, x4: Self) -> Self {
        // Generated with `./run-generator.sh f32::sin_cos::cos_poly`
        const K4: u32 = 0x3D2AAAA5; // 4.1666646e-2
        const K6: u32 = 0xBAB60642; // -1.3887363e-3
        const K8: u32 = 0x37CCFFFD; // 2.4437899e-5

        let k4 = Self::from_raw(K4);
        let k6 = Self::from_raw(K6);
        let k8 = Self::from_raw(K8);

        horner!(x4, x2, [k4, k6, k8])
    }
}