lufloat 0.1.1

Fastest FP16 Math and AI Library for AMD APUs
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

use crate::memory::{UnifiedBuffer, hip_check};
use std::ffi::c_int;

unsafe extern "C" {
    fn negative_mask(data: *const u16, size: usize, mask: *mut u16) -> c_int;
}

pub(crate) fn apply<'a>(data: &UnifiedBuffer<'a>, mask: &mut UnifiedBuffer<'a>) {
    assert_eq!(data.len, mask.len);
    let mut remaining = data.len;
    let mut offset = 0;
    while remaining > 0 {
        let current = remaining.min(1 << 34);
        let data_ptr = unsafe { data.ptr.add(offset) };
        let mask_ptr = unsafe { mask.ptr.add(offset) };
        let err = unsafe { negative_mask(data_ptr, current, mask_ptr) };
        hip_check(err, file!(), line!());
        remaining -= current;
        offset += current;
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use crate::memory::{Arena, UnifiedBuffer};

    #[test]
    fn exhaustive_negative_mask() {
        let arena = Arena::new(1 << 17);
        let mut data = UnifiedBuffer::new(&arena, 1 << 16);
        let mut mask = UnifiedBuffer::new(&arena, 1 << 16);
        let input_data = data.slice_mut();
        for i in 0..(1 << 16) {
            input_data[i] = i as u16;
        }
        apply(&data, &mut mask);
        let output_data = mask.slice();
        for i in 0..(1 << 16) {
            let expected = if (i as u16 & 0x8000) == 0 {
                0x0000
            } else {
                0x3C00
            };
            assert_eq!(output_data[i], expected, "Failed at {:016b}", i as u16);
        }
    }
}