Subroutine */
| 18578 | } /* zher2k_ */ |
| 18579 | |
| 18580 | /* Subroutine */ int zherk_(char *uplo, char *trans, integer *n, integer *k, |
| 18581 | doublereal *alpha, doublecomplex *a, integer *lda, doublereal *beta, |
| 18582 | doublecomplex *c__, integer *ldc) |
| 18583 | { |
| 18584 | /* System generated locals */ |
| 18585 | integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3, i__4, i__5, |
| 18586 | i__6; |
| 18587 | doublereal d__1; |
| 18588 | doublecomplex z__1, z__2, z__3; |
| 18589 | |
| 18590 | /* Local variables */ |
| 18591 | static integer i__, j, l, info; |
| 18592 | static doublecomplex temp; |
| 18593 | extern logical lsame_(char *, char *); |
| 18594 | static integer nrowa; |
| 18595 | static doublereal rtemp; |
| 18596 | static logical upper; |
| 18597 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 18598 | |
| 18599 | |
| 18600 | /* |
| 18601 | Purpose |
| 18602 | ======= |
| 18603 | |
| 18604 | ZHERK performs one of the hermitian rank k operations |
| 18605 | |
| 18606 | C := alpha*A*conjg( A' ) + beta*C, |
| 18607 | |
| 18608 | or |
| 18609 | |
| 18610 | C := alpha*conjg( A' )*A + beta*C, |
| 18611 | |
| 18612 | where alpha and beta are real scalars, C is an n by n hermitian |
| 18613 | matrix and A is an n by k matrix in the first case and a k by n |
| 18614 | matrix in the second case. |
| 18615 | |
| 18616 | Arguments |
| 18617 | ========== |
| 18618 | |
| 18619 | UPLO - CHARACTER*1. |
| 18620 | On entry, UPLO specifies whether the upper or lower |
| 18621 | triangular part of the array C is to be referenced as |
| 18622 | follows: |
| 18623 | |
| 18624 | UPLO = 'U' or 'u' Only the upper triangular part of C |
| 18625 | is to be referenced. |
| 18626 | |
| 18627 | UPLO = 'L' or 'l' Only the lower triangular part of C |
| 18628 | is to be referenced. |
| 18629 | |
| 18630 | Unchanged on exit. |
| 18631 | |
| 18632 | TRANS - CHARACTER*1. |
| 18633 | On entry, TRANS specifies the operation to be performed as |
| 18634 | follows: |
| 18635 | |
| 18636 | TRANS = 'N' or 'n' C := alpha*A*conjg( A' ) + beta*C. |
| 18637 |