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

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

version 1.6, 2010/12/21 13:51:24 version 1.7, 2011/07/22 07:38:01
Line 17 Line 17
 *  *
 *  where  op( X ) is one of  *  where  op( X ) is one of
 *  *
 *     op( X ) = X   or   op( X ) = X',  *     op( X ) = X   or   op( X ) = X**T,
 *  *
 *  alpha and beta are scalars, and A, B and C are matrices, with op( A )  *  alpha and beta are scalars, and A, B and C are matrices, with op( A )
 *  an m by k matrix,  op( B )  a  k by n matrix and  C an m by n matrix.  *  an m by k matrix,  op( B )  a  k by n matrix and  C an m by n matrix.
Line 31 Line 31
 *  *
 *              TRANSA = 'N' or 'n',  op( A ) = A.  *              TRANSA = 'N' or 'n',  op( A ) = A.
 *  *
 *              TRANSA = 'T' or 't',  op( A ) = A'.  *              TRANSA = 'T' or 't',  op( A ) = A**T.
 *  *
 *              TRANSA = 'C' or 'c',  op( A ) = A'.  *              TRANSA = 'C' or 'c',  op( A ) = A**T.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
 *  *
Line 43 Line 43
 *  *
 *              TRANSB = 'N' or 'n',  op( B ) = B.  *              TRANSB = 'N' or 'n',  op( B ) = B.
 *  *
 *              TRANSB = 'T' or 't',  op( B ) = B'.  *              TRANSB = 'T' or 't',  op( B ) = B**T.
 *  *
 *              TRANSB = 'C' or 'c',  op( B ) = B'.  *              TRANSB = 'C' or 'c',  op( B ) = B**T.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
 *  *
Line 249 Line 249
    90         CONTINUE     90         CONTINUE
           ELSE            ELSE
 *  *
 *           Form  C := alpha*A'*B + beta*C  *           Form  C := alpha*A**T*B + beta*C
 *  *
               DO 120 J = 1,N                DO 120 J = 1,N
                   DO 110 I = 1,M                    DO 110 I = 1,M
Line 268 Line 268
       ELSE        ELSE
           IF (NOTA) THEN            IF (NOTA) THEN
 *  *
 *           Form  C := alpha*A*B' + beta*C  *           Form  C := alpha*A*B**T + beta*C
 *  *
               DO 170 J = 1,N                DO 170 J = 1,N
                   IF (BETA.EQ.ZERO) THEN                    IF (BETA.EQ.ZERO) THEN
Line 291 Line 291
   170         CONTINUE    170         CONTINUE
           ELSE            ELSE
 *  *
 *           Form  C := alpha*A'*B' + beta*C  *           Form  C := alpha*A**T*B**T + beta*C
 *  *
               DO 200 J = 1,N                DO 200 J = 1,N
                   DO 190 I = 1,M                    DO 190 I = 1,M
Removed from v.1.6  
changed lines
  Added in v.1.7

CVSweb interface <joel.bertrand@systella.fr>