pub fn norm2_squared_diff<T>(a: &[T], b: &[T]) -> Twhere T: Float + Sum<T> + Mul<T, Output = T> + AddAssign,
Calculate the squared 2-norm of the difference of two vectors