Subroutine */
| 17504 | } /* zhemv_ */ |
| 17505 | |
| 17506 | /* Subroutine */ int zher2_(char *uplo, integer *n, doublecomplex *alpha, |
| 17507 | doublecomplex *x, integer *incx, doublecomplex *y, integer *incy, |
| 17508 | doublecomplex *a, integer *lda) |
| 17509 | { |
| 17510 | /* System generated locals */ |
| 17511 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5, i__6; |
| 17512 | doublereal d__1; |
| 17513 | doublecomplex z__1, z__2, z__3, z__4; |
| 17514 | |
| 17515 | /* Local variables */ |
| 17516 | static integer i__, j, ix, iy, jx, jy, kx, ky, info; |
| 17517 | static doublecomplex temp1, temp2; |
| 17518 | extern logical lsame_(char *, char *); |
| 17519 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 17520 | |
| 17521 | |
| 17522 | /* |
| 17523 | Purpose |
| 17524 | ======= |
| 17525 | |
| 17526 | ZHER2 performs the hermitian rank 2 operation |
| 17527 | |
| 17528 | A := alpha*x*conjg( y' ) + conjg( alpha )*y*conjg( x' ) + A, |
| 17529 | |
| 17530 | where alpha is a scalar, x and y are n element vectors and A is an n |
| 17531 | by n hermitian matrix. |
| 17532 | |
| 17533 | Arguments |
| 17534 | ========== |
| 17535 | |
| 17536 | UPLO - CHARACTER*1. |
| 17537 | On entry, UPLO specifies whether the upper or lower |
| 17538 | triangular part of the array A is to be referenced as |
| 17539 | follows: |
| 17540 | |
| 17541 | UPLO = 'U' or 'u' Only the upper triangular part of A |
| 17542 | is to be referenced. |
| 17543 | |
| 17544 | UPLO = 'L' or 'l' Only the lower triangular part of A |
| 17545 | is to be referenced. |
| 17546 | |
| 17547 | Unchanged on exit. |
| 17548 | |
| 17549 | N - INTEGER. |
| 17550 | On entry, N specifies the order of the matrix A. |
| 17551 | N must be at least zero. |
| 17552 | Unchanged on exit. |
| 17553 | |
| 17554 | ALPHA - COMPLEX*16 . |
| 17555 | On entry, ALPHA specifies the scalar alpha. |
| 17556 | Unchanged on exit. |
| 17557 | |
| 17558 | X - COMPLEX*16 array of dimension at least |
| 17559 | ( 1 + ( n - 1 )*abs( INCX ) ). |
| 17560 | Before entry, the incremented array X must contain the n |
| 17561 | element vector x. |
| 17562 | Unchanged on exit. |
| 17563 |