version 1.1, 2010/01/26 15:22:45
|
version 1.9, 2012/08/22 09:36:41
|
Line 1
|
Line 1
|
|
*> \brief \b ZHER2K |
|
* |
|
* =========== DOCUMENTATION =========== |
|
* |
|
* Online html documentation available at |
|
* http://www.netlib.org/lapack/explore-html/ |
|
* |
|
* Definition: |
|
* =========== |
|
* |
|
* SUBROUTINE ZHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC) |
|
* |
|
* .. Scalar Arguments .. |
|
* COMPLEX*16 ALPHA |
|
* DOUBLE PRECISION BETA |
|
* INTEGER K,LDA,LDB,LDC,N |
|
* CHARACTER TRANS,UPLO |
|
* .. |
|
* .. Array Arguments .. |
|
* COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*) |
|
* .. |
|
* |
|
* |
|
*> \par Purpose: |
|
* ============= |
|
*> |
|
*> \verbatim |
|
*> |
|
*> ZHER2K performs one of the hermitian rank 2k operations |
|
*> |
|
*> C := alpha*A*B**H + conjg( alpha )*B*A**H + beta*C, |
|
*> |
|
*> or |
|
*> |
|
*> C := alpha*A**H*B + conjg( alpha )*B**H*A + beta*C, |
|
*> |
|
*> where alpha and beta are scalars with beta real, C is an n by n |
|
*> hermitian matrix and A and B are n by k matrices in the first case |
|
*> and k by n matrices in the second case. |
|
*> \endverbatim |
|
* |
|
* Arguments: |
|
* ========== |
|
* |
|
*> \param[in] UPLO |
|
*> \verbatim |
|
*> UPLO is CHARACTER*1 |
|
*> On entry, UPLO specifies whether the upper or lower |
|
*> triangular part of the array C is to be referenced as |
|
*> follows: |
|
*> |
|
*> UPLO = 'U' or 'u' Only the upper triangular part of C |
|
*> is to be referenced. |
|
*> |
|
*> UPLO = 'L' or 'l' Only the lower triangular part of C |
|
*> is to be referenced. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] TRANS |
|
*> \verbatim |
|
*> TRANS is CHARACTER*1 |
|
*> On entry, TRANS specifies the operation to be performed as |
|
*> follows: |
|
*> |
|
*> TRANS = 'N' or 'n' C := alpha*A*B**H + |
|
*> conjg( alpha )*B*A**H + |
|
*> beta*C. |
|
*> |
|
*> TRANS = 'C' or 'c' C := alpha*A**H*B + |
|
*> conjg( alpha )*B**H*A + |
|
*> beta*C. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] N |
|
*> \verbatim |
|
*> N is INTEGER |
|
*> On entry, N specifies the order of the matrix C. N must be |
|
*> at least zero. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] K |
|
*> \verbatim |
|
*> K is INTEGER |
|
*> On entry with TRANS = 'N' or 'n', K specifies the number |
|
*> of columns of the matrices A and B, and on entry with |
|
*> TRANS = 'C' or 'c', K specifies the number of rows of the |
|
*> matrices A and B. K must be at least zero. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] ALPHA |
|
*> \verbatim |
|
*> ALPHA is COMPLEX*16 . |
|
*> On entry, ALPHA specifies the scalar alpha. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] A |
|
*> \verbatim |
|
*> A is COMPLEX*16 array of DIMENSION ( LDA, ka ), where ka is |
|
*> k when TRANS = 'N' or 'n', and is n otherwise. |
|
*> Before entry with TRANS = 'N' or 'n', the leading n by k |
|
*> part of the array A must contain the matrix A, otherwise |
|
*> the leading k by n part of the array A must contain the |
|
*> matrix A. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] LDA |
|
*> \verbatim |
|
*> LDA is INTEGER |
|
*> On entry, LDA specifies the first dimension of A as declared |
|
*> in the calling (sub) program. When TRANS = 'N' or 'n' |
|
*> then LDA must be at least max( 1, n ), otherwise LDA must |
|
*> be at least max( 1, k ). |
|
*> \endverbatim |
|
*> |
|
*> \param[in] B |
|
*> \verbatim |
|
*> B is COMPLEX*16 array of DIMENSION ( LDB, kb ), where kb is |
|
*> k when TRANS = 'N' or 'n', and is n otherwise. |
|
*> Before entry with TRANS = 'N' or 'n', the leading n by k |
|
*> part of the array B must contain the matrix B, otherwise |
|
*> the leading k by n part of the array B must contain the |
|
*> matrix B. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] LDB |
|
*> \verbatim |
|
*> LDB is INTEGER |
|
*> On entry, LDB specifies the first dimension of B as declared |
|
*> in the calling (sub) program. When TRANS = 'N' or 'n' |
|
*> then LDB must be at least max( 1, n ), otherwise LDB must |
|
*> be at least max( 1, k ). |
|
*> Unchanged on exit. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] BETA |
|
*> \verbatim |
|
*> BETA is DOUBLE PRECISION . |
|
*> On entry, BETA specifies the scalar beta. |
|
*> \endverbatim |
|
*> |
|
*> \param[in,out] C |
|
*> \verbatim |
|
*> C is COMPLEX*16 array of DIMENSION ( LDC, n ). |
|
*> Before entry with UPLO = 'U' or 'u', the leading n by n |
|
*> upper triangular part of the array C must contain the upper |
|
*> triangular part of the hermitian matrix and the strictly |
|
*> lower triangular part of C is not referenced. On exit, the |
|
*> upper triangular part of the array C is overwritten by the |
|
*> upper triangular part of the updated matrix. |
|
*> Before entry with UPLO = 'L' or 'l', the leading n by n |
|
*> lower triangular part of the array C must contain the lower |
|
*> triangular part of the hermitian matrix and the strictly |
|
*> upper triangular part of C is not referenced. On exit, the |
|
*> lower triangular part of the array C is overwritten by the |
|
*> lower triangular part of the updated matrix. |
|
*> Note that the imaginary parts of the diagonal elements need |
|
*> not be set, they are assumed to be zero, and on exit they |
|
*> are set to zero. |
|
*> \endverbatim |
|
*> |
|
*> \param[in] LDC |
|
*> \verbatim |
|
*> LDC is INTEGER |
|
*> On entry, LDC specifies the first dimension of C as declared |
|
*> in the calling (sub) program. LDC must be at least |
|
*> max( 1, n ). |
|
*> \endverbatim |
|
* |
|
* Authors: |
|
* ======== |
|
* |
|
*> \author Univ. of Tennessee |
|
*> \author Univ. of California Berkeley |
|
*> \author Univ. of Colorado Denver |
|
*> \author NAG Ltd. |
|
* |
|
*> \date November 2011 |
|
* |
|
*> \ingroup complex16_blas_level3 |
|
* |
|
*> \par Further Details: |
|
* ===================== |
|
*> |
|
*> \verbatim |
|
*> |
|
*> Level 3 Blas routine. |
|
*> |
|
*> -- Written on 8-February-1989. |
|
*> Jack Dongarra, Argonne National Laboratory. |
|
*> Iain Duff, AERE Harwell. |
|
*> Jeremy Du Croz, Numerical Algorithms Group Ltd. |
|
*> Sven Hammarling, Numerical Algorithms Group Ltd. |
|
*> |
|
*> -- Modified 8-Nov-93 to set C(J,J) to DBLE( C(J,J) ) when BETA = 1. |
|
*> Ed Anderson, Cray Research Inc. |
|
*> \endverbatim |
|
*> |
|
* ===================================================================== |
SUBROUTINE ZHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC) |
SUBROUTINE ZHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC) |
|
* |
|
* -- Reference BLAS level3 routine (version 3.4.0) -- |
|
* -- Reference BLAS is a software package provided by Univ. of Tennessee, -- |
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
|
* November 2011 |
|
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
DOUBLE COMPLEX ALPHA |
COMPLEX*16 ALPHA |
DOUBLE PRECISION BETA |
DOUBLE PRECISION BETA |
INTEGER K,LDA,LDB,LDC,N |
INTEGER K,LDA,LDB,LDC,N |
CHARACTER TRANS,UPLO |
CHARACTER TRANS,UPLO |
* .. |
* .. |
* .. Array Arguments .. |
* .. Array Arguments .. |
DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*) |
COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*) |
* .. |
* .. |
* |
* |
* Purpose |
|
* ======= |
|
* |
|
* ZHER2K performs one of the hermitian rank 2k operations |
|
* |
|
* C := alpha*A*conjg( B' ) + conjg( alpha )*B*conjg( A' ) + beta*C, |
|
* |
|
* or |
|
* |
|
* C := alpha*conjg( A' )*B + conjg( alpha )*conjg( B' )*A + beta*C, |
|
* |
|
* where alpha and beta are scalars with beta real, C is an n by n |
|
* hermitian matrix and A and B are n by k matrices in the first case |
|
* and k by n matrices in the second case. |
|
* |
|
* Arguments |
|
* ========== |
|
* |
|
* UPLO - CHARACTER*1. |
|
* On entry, UPLO specifies whether the upper or lower |
|
* triangular part of the array C is to be referenced as |
|
* follows: |
|
* |
|
* UPLO = 'U' or 'u' Only the upper triangular part of C |
|
* is to be referenced. |
|
* |
|
* UPLO = 'L' or 'l' Only the lower triangular part of C |
|
* is to be referenced. |
|
* |
|
* Unchanged on exit. |
|
* |
|
* TRANS - CHARACTER*1. |
|
* On entry, TRANS specifies the operation to be performed as |
|
* follows: |
|
* |
|
* TRANS = 'N' or 'n' C := alpha*A*conjg( B' ) + |
|
* conjg( alpha )*B*conjg( A' ) + |
|
* beta*C. |
|
* |
|
* TRANS = 'C' or 'c' C := alpha*conjg( A' )*B + |
|
* conjg( alpha )*conjg( B' )*A + |
|
* beta*C. |
|
* |
|
* Unchanged on exit. |
|
* |
|
* N - INTEGER. |
|
* On entry, N specifies the order of the matrix C. N must be |
|
* at least zero. |
|
* Unchanged on exit. |
|
* |
|
* K - INTEGER. |
|
* On entry with TRANS = 'N' or 'n', K specifies the number |
|
* of columns of the matrices A and B, and on entry with |
|
* TRANS = 'C' or 'c', K specifies the number of rows of the |
|
* matrices A and B. K must be at least zero. |
|
* Unchanged on exit. |
|
* |
|
* ALPHA - COMPLEX*16 . |
|
* On entry, ALPHA specifies the scalar alpha. |
|
* Unchanged on exit. |
|
* |
|
* A - COMPLEX*16 array of DIMENSION ( LDA, ka ), where ka is |
|
* k when TRANS = 'N' or 'n', and is n otherwise. |
|
* Before entry with TRANS = 'N' or 'n', the leading n by k |
|
* part of the array A must contain the matrix A, otherwise |
|
* the leading k by n part of the array A must contain the |
|
* matrix A. |
|
* Unchanged on exit. |
|
* |
|
* LDA - INTEGER. |
|
* On entry, LDA specifies the first dimension of A as declared |
|
* in the calling (sub) program. When TRANS = 'N' or 'n' |
|
* then LDA must be at least max( 1, n ), otherwise LDA must |
|
* be at least max( 1, k ). |
|
* Unchanged on exit. |
|
* |
|
* B - COMPLEX*16 array of DIMENSION ( LDB, kb ), where kb is |
|
* k when TRANS = 'N' or 'n', and is n otherwise. |
|
* Before entry with TRANS = 'N' or 'n', the leading n by k |
|
* part of the array B must contain the matrix B, otherwise |
|
* the leading k by n part of the array B must contain the |
|
* matrix B. |
|
* Unchanged on exit. |
|
* |
|
* LDB - INTEGER. |
|
* On entry, LDB specifies the first dimension of B as declared |
|
* in the calling (sub) program. When TRANS = 'N' or 'n' |
|
* then LDB must be at least max( 1, n ), otherwise LDB must |
|
* be at least max( 1, k ). |
|
* Unchanged on exit. |
|
* |
|
* BETA - DOUBLE PRECISION . |
|
* On entry, BETA specifies the scalar beta. |
|
* Unchanged on exit. |
|
* |
|
* C - COMPLEX*16 array of DIMENSION ( LDC, n ). |
|
* Before entry with UPLO = 'U' or 'u', the leading n by n |
|
* upper triangular part of the array C must contain the upper |
|
* triangular part of the hermitian matrix and the strictly |
|
* lower triangular part of C is not referenced. On exit, the |
|
* upper triangular part of the array C is overwritten by the |
|
* upper triangular part of the updated matrix. |
|
* Before entry with UPLO = 'L' or 'l', the leading n by n |
|
* lower triangular part of the array C must contain the lower |
|
* triangular part of the hermitian matrix and the strictly |
|
* upper triangular part of C is not referenced. On exit, the |
|
* lower triangular part of the array C is overwritten by the |
|
* lower triangular part of the updated matrix. |
|
* Note that the imaginary parts of the diagonal elements need |
|
* not be set, they are assumed to be zero, and on exit they |
|
* are set to zero. |
|
* |
|
* LDC - INTEGER. |
|
* On entry, LDC specifies the first dimension of C as declared |
|
* in the calling (sub) program. LDC must be at least |
|
* max( 1, n ). |
|
* Unchanged on exit. |
|
* |
|
* Further Details |
|
* =============== |
|
* |
|
* Level 3 Blas routine. |
|
* |
|
* -- Written on 8-February-1989. |
|
* Jack Dongarra, Argonne National Laboratory. |
|
* Iain Duff, AERE Harwell. |
|
* Jeremy Du Croz, Numerical Algorithms Group Ltd. |
|
* Sven Hammarling, Numerical Algorithms Group Ltd. |
|
* |
|
* -- Modified 8-Nov-93 to set C(J,J) to DBLE( C(J,J) ) when BETA = 1. |
|
* Ed Anderson, Cray Research Inc. |
|
* |
|
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. External Functions .. |
* .. External Functions .. |
Line 154
|
Line 226
|
INTRINSIC DBLE,DCONJG,MAX |
INTRINSIC DBLE,DCONJG,MAX |
* .. |
* .. |
* .. Local Scalars .. |
* .. Local Scalars .. |
DOUBLE COMPLEX TEMP1,TEMP2 |
COMPLEX*16 TEMP1,TEMP2 |
INTEGER I,INFO,J,L,NROWA |
INTEGER I,INFO,J,L,NROWA |
LOGICAL UPPER |
LOGICAL UPPER |
* .. |
* .. |
* .. Parameters .. |
* .. Parameters .. |
DOUBLE PRECISION ONE |
DOUBLE PRECISION ONE |
PARAMETER (ONE=1.0D+0) |
PARAMETER (ONE=1.0D+0) |
DOUBLE COMPLEX ZERO |
COMPLEX*16 ZERO |
PARAMETER (ZERO= (0.0D+0,0.0D+0)) |
PARAMETER (ZERO= (0.0D+0,0.0D+0)) |
* .. |
* .. |
* |
* |
Line 242
|
Line 314
|
* |
* |
IF (LSAME(TRANS,'N')) THEN |
IF (LSAME(TRANS,'N')) THEN |
* |
* |
* Form C := alpha*A*conjg( B' ) + conjg( alpha )*B*conjg( A' ) + |
* Form C := alpha*A*B**H + conjg( alpha )*B*A**H + |
* C. |
* C. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
Line 302
|
Line 374
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form C := alpha*conjg( A' )*B + conjg( alpha )*conjg( B' )*A + |
* Form C := alpha*A**H*B + conjg( alpha )*B**H*A + |
* C. |
* C. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |