[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.11 and 1.12

version 1.11, 2014/01/27 09:28:12 version 1.12, 2015/11/26 11:44:14
Line 166 Line 166
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date November 2015
 *  *
 *> \ingroup double_blas_level3  *> \ingroup double_blas_level3
 *  *
Line 187 Line 187
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)        SUBROUTINE DGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *  *
 *  -- Reference BLAS level3 routine (version 3.4.0) --  *  -- Reference BLAS level3 routine (version 3.6.0) --
 *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --  *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     November 2011  *     November 2015
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       DOUBLE PRECISION ALPHA,BETA        DOUBLE PRECISION ALPHA,BETA
Line 311 Line 311
    60                 CONTINUE     60                 CONTINUE
                   END IF                    END IF
                   DO 80 L = 1,K                    DO 80 L = 1,K
                       IF (B(L,J).NE.ZERO) THEN                        TEMP = ALPHA*B(L,J)
                           TEMP = ALPHA*B(L,J)                        DO 70 I = 1,M
                           DO 70 I = 1,M                            C(I,J) = C(I,J) + TEMP*A(I,L)
                               C(I,J) = C(I,J) + TEMP*A(I,L)     70                 CONTINUE
    70                     CONTINUE  
                       END IF  
    80             CONTINUE     80             CONTINUE
    90         CONTINUE     90         CONTINUE
           ELSE            ELSE
Line 353 Line 351
   140                 CONTINUE    140                 CONTINUE
                   END IF                    END IF
                   DO 160 L = 1,K                    DO 160 L = 1,K
                       IF (B(J,L).NE.ZERO) THEN                        TEMP = ALPHA*B(J,L)
                           TEMP = ALPHA*B(J,L)                        DO 150 I = 1,M
                           DO 150 I = 1,M                            C(I,J) = C(I,J) + TEMP*A(I,L)
                               C(I,J) = C(I,J) + TEMP*A(I,L)    150                 CONTINUE
   150                     CONTINUE  
                       END IF  
   160             CONTINUE    160             CONTINUE
   170         CONTINUE    170         CONTINUE
           ELSE            ELSE
Removed from v.1.11  
changed lines
  Added in v.1.12

CVSweb interface <joel.bertrand@systella.fr>