func axpy32(alpha float32, X []float32, Y []float32) {
blas.Saxpy(len(X), alpha, X, 1, Y, 1)
}
// Add performs v += a * u (in-place).
func (v Vector) Add(a float32, u Vector) {
blas.Saxpy(len(v), a, u, 1, v, 1)
}
func (y Vector) Add(alpha float32, x Vector) {
blas.Saxpy(len(y), alpha, x, 1, y, 1)
}