Lines Matching +full:x +full:- +full:alpha
1 /* chpmv.f -- translated by f2c (version 20100827).
4 on Linux or Unix systems, link with .../path/to/libf2c.a -lm
5 or, if you install libf2c.a in a standard place, with -lf2c -lm
6 -- in that order, at the end of the command line, as in
7 cc *.o -lf2c -lm
15 /* Subroutine */ int chpmv_(char *uplo, integer *n, complex *alpha, complex * in chpmv_() argument
16 ap, complex *x, integer *incx, complex *beta, complex *y, integer * in chpmv_() argument
41 /* CHPMV performs the matrix-vector operation */ in chpmv_()
43 /* y := alpha*A*x + beta*y, */ in chpmv_()
45 /* where alpha and beta are scalars, x and y are n element vectors and */ in chpmv_()
51 /* UPLO - CHARACTER*1. */ in chpmv_()
64 /* N - INTEGER. */ in chpmv_()
69 /* ALPHA - COMPLEX . */ in chpmv_()
70 /* On entry, ALPHA specifies the scalar alpha. */ in chpmv_()
73 /* AP - COMPLEX array of DIMENSION at least */ in chpmv_()
89 /* X - COMPLEX array of dimension at least */ in chpmv_()
90 /* ( 1 + ( n - 1 )*abs( INCX ) ). */ in chpmv_()
91 /* Before entry, the incremented array X must contain the n */ in chpmv_()
92 /* element vector x. */ in chpmv_()
95 /* INCX - INTEGER. */ in chpmv_()
97 /* X. INCX must not be zero. */ in chpmv_()
100 /* BETA - COMPLEX . */ in chpmv_()
105 /* Y - COMPLEX array of dimension at least */ in chpmv_()
106 /* ( 1 + ( n - 1 )*abs( INCY ) ). */ in chpmv_()
111 /* INCY - INTEGER. */ in chpmv_()
121 /* -- Written on 22-October-1986. */ in chpmv_()
143 --y; in chpmv_()
144 --x; in chpmv_()
145 --ap; in chpmv_()
166 if (*n == 0 || (alpha->r == 0.f && alpha->i == 0.f && (beta->r == 1.f && in chpmv_()
167 beta->i == 0.f))) { in chpmv_()
171 /* Set up the start points in X and Y. */ in chpmv_()
176 kx = 1 - (*n - 1) * *incx; in chpmv_()
181 ky = 1 - (*n - 1) * *incy; in chpmv_()
189 if (beta->r != 1.f || beta->i != 0.f) { in chpmv_()
191 if (beta->r == 0.f && beta->i == 0.f) { in chpmv_()
203 q__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in chpmv_()
204 q__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in chpmv_()
212 if (beta->r == 0.f && beta->i == 0.f) { in chpmv_()
225 q__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in chpmv_()
226 q__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in chpmv_()
235 if (alpha->r == 0.f && alpha->i == 0.f) { in chpmv_()
247 q__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, q__1.i = in chpmv_()
248 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in chpmv_()
252 i__2 = j - 1; in chpmv_()
257 q__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in chpmv_()
264 q__2.r = q__3.r * x[i__3].r - q__3.i * x[i__3].i, q__2.i = in chpmv_()
265 q__3.r * x[i__3].i + q__3.i * x[i__3].r; in chpmv_()
273 i__4 = kk + j - 1; in chpmv_()
277 q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i = in chpmv_()
278 alpha->r * temp2.i + alpha->i * temp2.r; in chpmv_()
290 q__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, q__1.i = in chpmv_()
291 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in chpmv_()
296 i__2 = kk + j - 2; in chpmv_()
301 q__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in chpmv_()
308 q__2.r = q__3.r * x[i__3].r - q__3.i * x[i__3].i, q__2.i = in chpmv_()
309 q__3.r * x[i__3].i + q__3.i * x[i__3].r; in chpmv_()
318 i__4 = kk + j - 1; in chpmv_()
322 q__4.r = alpha->r * temp2.r - alpha->i * temp2.i, q__4.i = in chpmv_()
323 alpha->r * temp2.i + alpha->i * temp2.r; in chpmv_()
340 q__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, q__1.i = in chpmv_()
341 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in chpmv_()
357 q__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in chpmv_()
364 q__2.r = q__3.r * x[i__3].r - q__3.i * x[i__3].i, q__2.i = in chpmv_()
365 q__3.r * x[i__3].i + q__3.i * x[i__3].r; in chpmv_()
373 q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i = in chpmv_()
374 alpha->r * temp2.i + alpha->i * temp2.r; in chpmv_()
377 kk += *n - j + 1; in chpmv_()
386 q__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, q__1.i = in chpmv_()
387 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in chpmv_()
399 i__2 = kk + *n - j; in chpmv_()
406 q__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in chpmv_()
413 q__2.r = q__3.r * x[i__3].r - q__3.i * x[i__3].i, q__2.i = in chpmv_()
414 q__3.r * x[i__3].i + q__3.i * x[i__3].r; in chpmv_()
421 q__2.r = alpha->r * temp2.r - alpha->i * temp2.i, q__2.i = in chpmv_()
422 alpha->r * temp2.i + alpha->i * temp2.r; in chpmv_()
427 kk += *n - j + 1; in chpmv_()