1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
ix!(); #[inline] pub fn sincf(x: f64) -> f64 { if x == 0.0 { return 1.0; } (PI * x).sin() / (PI * x) } #[inline] pub fn sinc(x: f64) -> f64 { if x.abs() < 0.0000000000000000000001 { 1.0 } else { x.sin() / x } }