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
|
// Copyright 2024 the Fearless_SIMD Authors
// SPDX-License-Identifier: Apache-2.0 OR MIT
#![expect(
missing_docs,
reason = "TODO: https://github.com/linebender/fearless_simd/issues/40"
)]
use fearless_simd::{Level, Simd, SimdBase, SimdFloat, dispatch};
#[inline(always)]
fn sigmoid<S: Simd>(simd: S, x: &[f32], out: &mut [f32]) {
let n = S::f32s::N;
for (x, y) in x.chunks_exact(n).zip(out.chunks_exact_mut(n)) {
let a = S::f32s::from_slice(simd, x);
let b = a / (a * a + 1.0).sqrt();
y.copy_from_slice(b.as_slice());
}
}
fn main() {
let level = Level::new();
let inp = [0.1, -0.2, 0.001, 0.4, 1., 2., 3., 4.];
let mut out = [0.; 8];
dispatch!(level, simd => sigmoid(simd, &inp, &mut out));
println!("{out:?}");
}
|
