Diff for /rpl/lapack/blas/zher2k.f between versions 1.6 and 1.7

version 1.6, 2010/12/21 13:51:26 version 1.7, 2011/07/22 07:38:03
Line 14 Line 14
 *  *
 *  ZHER2K  performs one of the hermitian rank 2k operations  *  ZHER2K  performs one of the hermitian rank 2k operations
 *  *
 *     C := alpha*A*conjg( B' ) + conjg( alpha )*B*conjg( A' ) + beta*C,  *     C := alpha*A*B**H + conjg( alpha )*B*A**H + beta*C,
 *  *
 *  or  *  or
 *  *
 *     C := alpha*conjg( A' )*B + conjg( alpha )*conjg( B' )*A + beta*C,  *     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  *  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  *  hermitian matrix and  A and B  are  n by k matrices in the first case
Line 44 Line 44
 *           On entry,  TRANS  specifies the operation to be performed as  *           On entry,  TRANS  specifies the operation to be performed as
 *           follows:  *           follows:
 *  *
 *              TRANS = 'N' or 'n'    C := alpha*A*conjg( B' )          +  *              TRANS = 'N' or 'n'    C := alpha*A*B**H          +
 *                                         conjg( alpha )*B*conjg( A' ) +  *                                         conjg( alpha )*B*A**H +
 *                                         beta*C.  *                                         beta*C.
 *  *
 *              TRANS = 'C' or 'c'    C := alpha*conjg( A' )*B          +  *              TRANS = 'C' or 'c'    C := alpha*A**H*B          +
 *                                         conjg( alpha )*conjg( B' )*A +  *                                         conjg( alpha )*B**H*A +
 *                                         beta*C.  *                                         beta*C.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
Line 242 Line 242
 *  *
       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 302
           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

Removed from v.1.6  
changed lines
  Added in v.1.7


CVSweb interface <joel.bertrand@systella.fr>