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];
int integer
Definition: template_blas_common.h:38
int template_blas_axpy(const integer *n, const Treal *da, const Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_axpy.h:41