Function axpy 

pub fn axpy(alpha: f64, x: &[f64], y: &mut [f64]) -> Result<(), SolversError>

y[i] += α·x[i] — the BLAS axpy.