Add BLAS function axpby.
Missing from original implementation.
Note: changed some .reverse() calls in axpy since y.reverse() += alpha * x should be equivalent to y += alpha * x.reverse().
.reverse()
axpy
y.reverse() += alpha * x
y += alpha * x.reverse()