35 #ifndef TEMPLATE_BLAS_SPR2_HEADER
36 #define TEMPLATE_BLAS_SPR2_HEADER
42 const Treal *x,
const integer *incx,
const Treal *y,
const integer *incy,
51 integer kk, ix, iy, jx, jy, kx, ky;
122 jx = jy = kx = ky = 0;
129 }
else if (*incx == 0) {
131 }
else if (*incy == 0) {
139 if (*n == 0 || *alpha == 0.) {
144 if (*incx != 1 || *incy != 1) {
148 kx = 1 - (*n - 1) * *incx;
153 ky = 1 - (*n - 1) * *incy;
163 if (*incx == 1 && *incy == 1) {
165 for (j = 1; j <= i__1; ++j) {
166 if (x[j] != 0. || y[j] != 0.) {
167 temp1 = *alpha * y[j];
168 temp2 = *alpha * x[j];
171 for (i__ = 1; i__ <= i__2; ++i__) {
172 ap[k] = ap[k] + x[i__] * temp1 + y[i__] * temp2;
182 for (j = 1; j <= i__1; ++j) {
183 if (x[jx] != 0. || y[jy] != 0.) {
184 temp1 = *alpha * y[jy];
185 temp2 = *alpha * x[jx];
189 for (k = kk; k <= i__2; ++k) {
190 ap[k] = ap[k] + x[ix] * temp1 + y[iy] * temp2;
204 if (*incx == 1 && *incy == 1) {
206 for (j = 1; j <= i__1; ++j) {
207 if (x[j] != 0. || y[j] != 0.) {
208 temp1 = *alpha * y[j];
209 temp2 = *alpha * x[j];
212 for (i__ = j; i__ <= i__2; ++i__) {
213 ap[k] = ap[k] + x[i__] * temp1 + y[i__] * temp2;
218 kk = kk + *n - j + 1;
223 for (j = 1; j <= i__1; ++j) {
224 if (x[jx] != 0. || y[jy] != 0.) {
225 temp1 = *alpha * y[jy];
226 temp2 = *alpha * x[jx];
230 for (k = kk; k <= i__2; ++k) {
231 ap[k] = ap[k] + x[ix] * temp1 + y[iy] * temp2;
239 kk = kk + *n - j + 1;