lbfgs::vec_ops
pub fn difference_and_save<T>(out: &mut [T], a: &[T], b: &[T]) where T: Float,
Calculates the difference of two slices and saves it in the third: out = a - b