pub fn reciprocal<T, D>(input: &Array<T, D>) -> Result<Array<T, D>, FerrayError>where T: Element + Float, D: Dimension,
Elementwise reciprocal: 1/x.
Uses hardware SIMD for contiguous f64 arrays.