35 #ifndef TEMPLATE_BLAS_AXPY_HEADER
36 #define TEMPLATE_BLAS_AXPY_HEADER
62 if (*incx == 1 && *incy == 1) {
70 ix = (-(*n) + 1) * *incx + 1;
73 iy = (-(*n) + 1) * *incy + 1;
76 for (i__ = 1; i__ <= i__1; ++i__) {
77 dy[iy] += *da * dx[ix];
91 for (i__ = 1; i__ <= i__1; ++i__) {
92 dy[i__] += *da * dx[i__];
101 for (i__ = mp1; i__ <= i__1; i__ += 4) {
102 dy[i__] += *da * dx[i__];
103 dy[i__ + 1] += *da * dx[i__ + 1];
104 dy[i__ + 2] += *da * dx[i__ + 2];
105 dy[i__ + 3] += *da * dx[i__ + 3];