version 1.9, 2012/08/22 09:36:40
|
version 1.14, 2017/06/17 10:53:43
|
Line 2
|
Line 2
|
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
* Online html documentation available at |
* Online html documentation available at |
* http://www.netlib.org/lapack/explore-html/ |
* http://www.netlib.org/lapack/explore-html/ |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* 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) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* DOUBLE PRECISION ALPHA,BETA |
* DOUBLE PRECISION ALPHA,BETA |
* INTEGER INCX,INCY,LDA,M,N |
* INTEGER INCX,INCY,LDA,M,N |
Line 18
|
Line 18
|
* .. Array Arguments .. |
* .. Array Arguments .. |
* DOUBLE PRECISION A(LDA,*),X(*),Y(*) |
* DOUBLE PRECISION A(LDA,*),X(*),Y(*) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 129
|
Line 129
|
* Authors: |
* Authors: |
* ======== |
* ======== |
* |
* |
*> \author Univ. of Tennessee |
*> \author Univ. of Tennessee |
*> \author Univ. of California Berkeley |
*> \author Univ. of California Berkeley |
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date December 2016 |
* |
* |
*> \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.7.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 |
* December 2016 |
* |
* |
* .. 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 |