Lines Matching +full:x +full:- +full:alpha
1 /* zhpmv.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 zhpmv_(char *uplo, integer *n, doublecomplex *alpha, in zhpmv_() argument
16 doublecomplex *ap, doublecomplex *x, integer *incx, doublecomplex * in zhpmv_() argument
41 /* ZHPMV performs the matrix-vector operation */ in zhpmv_()
43 /* y := alpha*A*x + beta*y, */ in zhpmv_()
45 /* where alpha and beta are scalars, x and y are n element vectors and */ in zhpmv_()
51 /* UPLO - CHARACTER*1. */ in zhpmv_()
64 /* N - INTEGER. */ in zhpmv_()
69 /* ALPHA - COMPLEX*16 . */ in zhpmv_()
70 /* On entry, ALPHA specifies the scalar alpha. */ in zhpmv_()
73 /* AP - COMPLEX*16 array of DIMENSION at least */ in zhpmv_()
89 /* X - COMPLEX*16 array of dimension at least */ in zhpmv_()
90 /* ( 1 + ( n - 1 )*abs( INCX ) ). */ in zhpmv_()
91 /* Before entry, the incremented array X must contain the n */ in zhpmv_()
92 /* element vector x. */ in zhpmv_()
95 /* INCX - INTEGER. */ in zhpmv_()
97 /* X. INCX must not be zero. */ in zhpmv_()
100 /* BETA - COMPLEX*16 . */ in zhpmv_()
105 /* Y - COMPLEX*16 array of dimension at least */ in zhpmv_()
106 /* ( 1 + ( n - 1 )*abs( INCY ) ). */ in zhpmv_()
111 /* INCY - INTEGER. */ in zhpmv_()
121 /* -- Written on 22-October-1986. */ in zhpmv_()
143 --y; in zhpmv_()
144 --x; in zhpmv_()
145 --ap; in zhpmv_()
166 if (*n == 0 || (alpha->r == 0. && alpha->i == 0. && (beta->r == 1. && in zhpmv_()
167 beta->i == 0.))) { in zhpmv_()
171 /* Set up the start points in X and Y. */ in zhpmv_()
176 kx = 1 - (*n - 1) * *incx; in zhpmv_()
181 ky = 1 - (*n - 1) * *incy; in zhpmv_()
189 if (beta->r != 1. || beta->i != 0.) { in zhpmv_()
191 if (beta->r == 0. && beta->i == 0.) { in zhpmv_()
203 z__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in zhpmv_()
204 z__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in zhpmv_()
212 if (beta->r == 0. && beta->i == 0.) { in zhpmv_()
225 z__1.r = beta->r * y[i__3].r - beta->i * y[i__3].i, in zhpmv_()
226 z__1.i = beta->r * y[i__3].i + beta->i * y[i__3] in zhpmv_()
235 if (alpha->r == 0. && alpha->i == 0.) { in zhpmv_()
247 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
248 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
252 i__2 = j - 1; in zhpmv_()
257 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
264 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
265 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
273 i__4 = kk + j - 1; in zhpmv_()
277 z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i = in zhpmv_()
278 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
290 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
291 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
296 i__2 = kk + j - 2; in zhpmv_()
301 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
308 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
309 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
318 i__4 = kk + j - 1; in zhpmv_()
322 z__4.r = alpha->r * temp2.r - alpha->i * temp2.i, z__4.i = in zhpmv_()
323 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
340 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
341 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
357 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
364 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
365 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
373 z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i = in zhpmv_()
374 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
377 kk += *n - j + 1; in zhpmv_()
386 z__1.r = alpha->r * x[i__2].r - alpha->i * x[i__2].i, z__1.i = in zhpmv_()
387 alpha->r * x[i__2].i + alpha->i * x[i__2].r; in zhpmv_()
399 i__2 = kk + *n - j; in zhpmv_()
406 z__2.r = temp1.r * ap[i__5].r - temp1.i * ap[i__5].i, in zhpmv_()
413 z__2.r = z__3.r * x[i__3].r - z__3.i * x[i__3].i, z__2.i = in zhpmv_()
414 z__3.r * x[i__3].i + z__3.i * x[i__3].r; in zhpmv_()
421 z__2.r = alpha->r * temp2.r - alpha->i * temp2.i, z__2.i = in zhpmv_()
422 alpha->r * temp2.i + alpha->i * temp2.r; in zhpmv_()
427 kk += *n - j + 1; in zhpmv_()