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

//! Fallback implementations of the SIMD traits.

use std::ptr;

use traits::{SimdF32, SimdMask32};

impl SimdF32 for f32 {
    type Raw = f32;

    type Mask = u32;

    #[inline]
    fn width(self) -> usize { 1 }

    #[inline]
    fn floor(self) -> f32 { f32::floor(self) }

    #[inline]
    fn ceil(self) -> f32 { f32::ceil(self) }

    #[inline]
    // See https://github.com/rust-lang/rust/issues/55107 for some discussion
    // of this choice. Basically, this implements round up on even semantics
    // for all values other than -0.5.
    fn round(self) -> f32 { f32::floor(self + (0.5 - 0.25 * ::std::f32::EPSILON)) }

    #[inline]
    fn abs(self) -> f32 { f32::abs(self) }

    #[inline]
    fn min(self, b: f32) -> f32 { f32::min(self, b) }

    #[inline]
    fn max(self, b: f32) -> f32 { f32::max(self, b) }

    #[inline]
    fn recip(self) -> f32 { f32::recip(self) }

    #[inline]
    fn rsqrt(self) -> f32 { f32::recip(f32::sqrt(self)) }

    #[inline]
    fn splat(self, x: f32) -> f32 { x }

    #[inline]
    fn steps(self) -> f32 { 0.0 }

    #[inline]
    unsafe fn from_raw(raw: f32) -> f32 { raw }

    #[inline]
    unsafe fn load(p: *const f32) -> f32 { ptr::read(p) }

    #[inline]
    unsafe fn store(self, p: *mut f32) { ptr::write(p, self); }

    #[inline]
    unsafe fn create() -> f32 { 0.0 }

    #[inline]
    fn eq(self, other: f32) -> u32 {
        if self == other { !0 } else { 0 }
    }
}

impl SimdMask32 for u32 {
    type Raw = u32;
    type F32 = f32;

    #[inline]
    fn select(self, a: f32, b: f32) -> f32 {
        if self & 0x80000000 != 0 { a } else { b }
    }
}