lbfgs::vec_ops
pub fn inplace_vec_add<T>(out: &mut [T], a: &[T], s: T) where T: Float,
Calculates out = out + s * a