37 #ifndef TEMPLATE_BLAS_AXPY_HEADER
38 #define TEMPLATE_BLAS_AXPY_HEADER
64 if (*incx == 1 && *incy == 1) {
72 ix = (-(*n) + 1) * *incx + 1;
75 iy = (-(*n) + 1) * *incy + 1;
78 for (i__ = 1; i__ <= i__1; ++i__) {
79 dy[iy] += *da * dx[ix];
93 for (i__ = 1; i__ <= i__1; ++i__) {
94 dy[i__] += *da * dx[i__];
103 for (i__ = mp1; i__ <= i__1; i__ += 4) {
104 dy[i__] += *da * dx[i__];
105 dy[i__ + 1] += *da * dx[i__ + 1];
106 dy[i__ + 2] += *da * dx[i__ + 2];
107 dy[i__ + 3] += *da * dx[i__ + 3];