37 #ifndef TEMPLATE_BLAS_SPR_HEADER
38 #define TEMPLATE_BLAS_SPR_HEADER
44 Treal *x,
integer *incx, Treal *ap)
120 }
else if (*incx == 0) {
128 if (*n == 0 || *alpha == 0.) {
133 kx = 1 - (*n - 1) * *incx;
134 }
else if (*incx != 1) {
144 for (j = 1; j <= i__1; ++j) {
146 temp = *alpha * x[j];
149 for (i__ = 1; i__ <= i__2; ++i__) {
150 ap[k] += x[i__] * temp;
161 for (j = 1; j <= i__1; ++j) {
163 temp = *alpha * x[jx];
166 for (k = kk; k <= i__2; ++k) {
167 ap[k] += x[ix] * temp;
181 for (j = 1; j <= i__1; ++j) {
183 temp = *alpha * x[j];
186 for (i__ = j; i__ <= i__2; ++i__) {
187 ap[k] += x[i__] * temp;
192 kk = kk + *n - j + 1;
198 for (j = 1; j <= i__1; ++j) {
200 temp = *alpha * x[jx];
203 for (k = kk; k <= i__2; ++k) {
204 ap[k] += x[ix] * temp;
210 kk = kk + *n - j + 1;