def daxpy(alpha, x, y): """ This function computes the sum \M{y = S{alpha} x + y} for the vectors x and y. """ yn = array_typed_copy(y, Float) _gslwrap.gsl_blas_daxpy(alpha, x, yn) return yn
def daxpy_cr(alpha, x, y_CR): """ This function computes the sum \M{y = S{alpha} x + y} for the vectors x and y. """ _gslwrap.gsl_blas_daxpy(alpha, x, y_CR)