Subroutine */
| 16313 | } /* zgemm_ */ |
| 16314 | |
| 16315 | /* Subroutine */ int zgemv_(char *trans, integer *m, integer *n, |
| 16316 | doublecomplex *alpha, doublecomplex *a, integer *lda, doublecomplex * |
| 16317 | x, integer *incx, doublecomplex *beta, doublecomplex *y, integer * |
| 16318 | incy) |
| 16319 | { |
| 16320 | /* System generated locals */ |
| 16321 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5; |
| 16322 | doublecomplex z__1, z__2, z__3; |
| 16323 | |
| 16324 | /* Local variables */ |
| 16325 | static integer i__, j, ix, iy, jx, jy, kx, ky, info; |
| 16326 | static doublecomplex temp; |
| 16327 | static integer lenx, leny; |
| 16328 | extern logical lsame_(char *, char *); |
| 16329 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 16330 | static logical noconj; |
| 16331 | |
| 16332 | |
| 16333 | /* |
| 16334 | Purpose |
| 16335 | ======= |
| 16336 | |
| 16337 | ZGEMV performs one of the matrix-vector operations |
| 16338 | |
| 16339 | y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, or |
| 16340 | |
| 16341 | y := alpha*conjg( A' )*x + beta*y, |
| 16342 | |
| 16343 | where alpha and beta are scalars, x and y are vectors and A is an |
| 16344 | m by n matrix. |
| 16345 | |
| 16346 | Arguments |
| 16347 | ========== |
| 16348 | |
| 16349 | TRANS - CHARACTER*1. |
| 16350 | On entry, TRANS specifies the operation to be performed as |
| 16351 | follows: |
| 16352 | |
| 16353 | TRANS = 'N' or 'n' y := alpha*A*x + beta*y. |
| 16354 | |
| 16355 | TRANS = 'T' or 't' y := alpha*A'*x + beta*y. |
| 16356 | |
| 16357 | TRANS = 'C' or 'c' y := alpha*conjg( A' )*x + beta*y. |
| 16358 | |
| 16359 | Unchanged on exit. |
| 16360 | |
| 16361 | M - INTEGER. |
| 16362 | On entry, M specifies the number of rows of the matrix A. |
| 16363 | M must be at least zero. |
| 16364 | Unchanged on exit. |
| 16365 | |
| 16366 | N - INTEGER. |
| 16367 | On entry, N specifies the number of columns of the matrix A. |
| 16368 | N must be at least zero. |
| 16369 | Unchanged on exit. |
| 16370 | |
| 16371 | ALPHA - COMPLEX*16 . |
| 16372 | On entry, ALPHA specifies the scalar alpha. |