1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
#[inline] pub fn sinc<T: num::Float>(x: T) -> T { if x.is_zero() { T::one() } else { x.sin() / x } } #[cfg(test)] mod tests { use super::*; #[test] fn test_sinc() { assert!(f64::abs(sinc(0.0) - 1.0) < 1e-10_f64); } }