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;
int template_blas_spr(const char *uplo, integer *n, Treal *alpha, Treal *x, integer *incx, Treal *ap)
Definition: template_blas_spr.h:43
int integer
Definition: template_blas_common.h:40
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46