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];
int integer
Definition: template_blas_common.h:40
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:43