pub fn kahan_sigma_return_counter<E, I: Iterator<Item = E>, F, Dtype>(
    element_iterator: I,
    op: F
) -> (Dtype, usize)
where F: Fn(E) -> Dtype, Dtype: Float,