[BACK]Return to dsyr2k.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / blas

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

version 1.6, 2010/12/21 13:51:25 version 1.7, 2011/07/22 07:38:02
Line 13 Line 13
 *  *
 *  DSYR2K  performs one of the symmetric rank 2k operations  *  DSYR2K  performs one of the symmetric rank 2k operations
 *  *
 *     C := alpha*A*B' + alpha*B*A' + beta*C,  *     C := alpha*A*B**T + alpha*B*A**T + beta*C,
 *  *
 *  or  *  or
 *  *
 *     C := alpha*A'*B + alpha*B'*A + beta*C,  *     C := alpha*A**T*B + alpha*B**T*A + beta*C,
 *  *
 *  where  alpha and beta  are scalars, C is an  n by n  symmetric matrix  *  where  alpha and beta  are scalars, C is an  n by n  symmetric matrix
 *  and  A and B  are  n by k  matrices  in the  first  case  and  k by n  *  and  A and B  are  n by k  matrices  in the  first  case  and  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*B' + alpha*B*A' +  *              TRANS = 'N' or 'n'   C := alpha*A*B**T + alpha*B*A**T +
 *                                        beta*C.  *                                        beta*C.
 *  *
 *              TRANS = 'T' or 't'   C := alpha*A'*B + alpha*B'*A +  *              TRANS = 'T' or 't'   C := alpha*A**T*B + alpha*B**T*A +
 *                                        beta*C.  *                                        beta*C.
 *  *
 *              TRANS = 'C' or 'c'   C := alpha*A'*B + alpha*B'*A +  *              TRANS = 'C' or 'c'   C := alpha*A**T*B + alpha*B**T*A +
 *                                        beta*C.  *                                        beta*C.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
Line 234 Line 234
 *  *
       IF (LSAME(TRANS,'N')) THEN        IF (LSAME(TRANS,'N')) THEN
 *  *
 *        Form  C := alpha*A*B' + alpha*B*A' + C.  *        Form  C := alpha*A*B**T + alpha*B*A**T + C.
 *  *
           IF (UPPER) THEN            IF (UPPER) THEN
               DO 130 J = 1,N                DO 130 J = 1,N
Line 283 Line 283
           END IF            END IF
       ELSE        ELSE
 *  *
 *        Form  C := alpha*A'*B + alpha*B'*A + C.  *        Form  C := alpha*A**T*B + alpha*B**T*A + C.
 *  *
           IF (UPPER) THEN            IF (UPPER) THEN
               DO 210 J = 1,N                DO 210 J = 1,N
Removed from v.1.6  
changed lines
  Added in v.1.7

CVSweb interface <joel.bertrand@systella.fr>