Diff for /rpl/lapack/blas/zherk.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 13 Line 13
 *  *
 *  ZHERK  performs one of the hermitian rank k operations  *  ZHERK  performs one of the hermitian rank k operations
 *  *
 *     C := alpha*A*conjg( A' ) + beta*C,  *     C := alpha*A*A**H + beta*C,
 *  *
 *  or  *  or
 *  *
 *     C := alpha*conjg( A' )*A + beta*C,  *     C := alpha*A**H*A + beta*C,
 *  *
 *  where  alpha and beta  are  real scalars,  C is an  n by n  hermitian  *  where  alpha and beta  are  real scalars,  C is an  n by n  hermitian
 *  matrix and  A  is an  n by k  matrix in the  first case and a  k by n  *  matrix and  A  is an  n by k  matrix in the  first case and a  k by n
Line 43 Line 43
 *           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( A' ) + beta*C.  *              TRANS = 'N' or 'n'   C := alpha*A*A**H + beta*C.
 *  *
 *              TRANS = 'C' or 'c'   C := alpha*conjg( A' )*A + beta*C.  *              TRANS = 'C' or 'c'   C := alpha*A**H*A + beta*C.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
 *  *
Line 219 Line 219
 *  *
       IF (LSAME(TRANS,'N')) THEN        IF (LSAME(TRANS,'N')) THEN
 *  *
 *        Form  C := alpha*A*conjg( A' ) + beta*C.  *        Form  C := alpha*A*A**H + beta*C.
 *  *
           IF (UPPER) THEN            IF (UPPER) THEN
               DO 130 J = 1,N                DO 130 J = 1,N
Line 272 Line 272
           END IF            END IF
       ELSE        ELSE
 *  *
 *        Form  C := alpha*conjg( A' )*A + beta*C.  *        Form  C := alpha*A**H*A + beta*C.
 *  *
           IF (UPPER) THEN            IF (UPPER) THEN
               DO 220 J = 1,N                DO 220 J = 1,N

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


CVSweb interface <joel.bertrand@systella.fr>