35 #ifndef TEMPLATE_BLAS_SPR_HEADER
36 #define TEMPLATE_BLAS_SPR_HEADER
42 Treal *x,
integer *incx, Treal *ap)
118 }
else if (*incx == 0) {
126 if (*n == 0 || *alpha == 0.) {
131 kx = 1 - (*n - 1) * *incx;
132 }
else if (*incx != 1) {
142 for (j = 1; j <= i__1; ++j) {
144 temp = *alpha * x[j];
147 for (i__ = 1; i__ <= i__2; ++i__) {
148 ap[k] += x[i__] * temp;
159 for (j = 1; j <= i__1; ++j) {
161 temp = *alpha * x[jx];
164 for (k = kk; k <= i__2; ++k) {
165 ap[k] += x[ix] * temp;
179 for (j = 1; j <= i__1; ++j) {
181 temp = *alpha * x[j];
184 for (i__ = j; i__ <= i__2; ++i__) {
185 ap[k] += x[i__] * temp;
190 kk = kk + *n - j + 1;
196 for (j = 1; j <= i__1; ++j) {
198 temp = *alpha * x[jx];
201 for (k = kk; k <= i__2; ++k) {
202 ap[k] += x[ix] * temp;
208 kk = kk + *n - j + 1;