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

Diff for /rpl/lapack/blas/dgemv.f between versions 1.11 and 1.13

version 1.11, 2014/01/27 09:28:12 version 1.13, 2016/08/27 15:37:52
Line 134 Line 134
 *> \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_level2  *> \ingroup double_blas_level2
 *  *
Line 156 Line 156
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)        SUBROUTINE DGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
 *  *
 *  -- Reference BLAS level2 routine (version 3.4.0) --  *  -- Reference BLAS level2 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 278 Line 278
           JX = KX            JX = KX
           IF (INCY.EQ.1) THEN            IF (INCY.EQ.1) THEN
               DO 60 J = 1,N                DO 60 J = 1,N
                   IF (X(JX).NE.ZERO) THEN                    TEMP = ALPHA*X(JX)
                       TEMP = ALPHA*X(JX)                    DO 50 I = 1,M
                       DO 50 I = 1,M                        Y(I) = Y(I) + TEMP*A(I,J)
                           Y(I) = Y(I) + TEMP*A(I,J)     50             CONTINUE
    50                 CONTINUE  
                   END IF  
                   JX = JX + INCX                    JX = JX + INCX
    60         CONTINUE     60         CONTINUE
           ELSE            ELSE
               DO 80 J = 1,N                DO 80 J = 1,N
                   IF (X(JX).NE.ZERO) THEN                    TEMP = ALPHA*X(JX)
                       TEMP = ALPHA*X(JX)                    IY = KY
                       IY = KY                    DO 70 I = 1,M
                       DO 70 I = 1,M                        Y(IY) = Y(IY) + TEMP*A(I,J)
                           Y(IY) = Y(IY) + TEMP*A(I,J)                        IY = IY + INCY
                           IY = IY + INCY     70             CONTINUE
    70                 CONTINUE  
                   END IF  
                   JX = JX + INCX                    JX = JX + INCX
    80         CONTINUE     80         CONTINUE
           END IF            END IF
Removed from v.1.11  
changed lines
  Added in v.1.13

CVSweb interface <joel.bertrand@systella.fr>