1 SUBROUTINE dsyr ( UPLO, N, ALPHA, X, INCX, A, LDA )
7 DOUBLE PRECISION A( LDA, * ), X( * )
88 parameter( zero = 0.0d+0 )
91 INTEGER I, INFO, IX, J, JX, KX
105 IF ( .NOT.lsame( uplo,
'U' ).AND.
106 $ .NOT.lsame( uplo,
'L' ) )
THEN
108 ELSE IF( n.LT.0 )
THEN
110 ELSE IF( incx.EQ.0 )
THEN
112 ELSE IF( lda.LT.max( 1, n ) )
THEN
116 CALL xerbla(
'DSYR ', info )
122 IF( ( n.EQ.0 ).OR.( alpha.EQ.zero ) )
128 kx = 1 - ( n - 1 )*incx
129 ELSE IF( incx.NE.1 )
THEN
137 IF( lsame( uplo,
'U' ) )
THEN
143 IF( x( j ).NE.zero )
THEN
146 a( i, j ) = a( i, j ) + x( i )*temp
153 IF( x( jx ).NE.zero )
THEN
157 a( i, j ) = a( i, j ) + x( ix )*temp
170 IF( x( j ).NE.zero )
THEN
173 a( i, j ) = a( i, j ) + x( i )*temp
180 IF( x( jx ).NE.zero )
THEN
184 a( i, j ) = a( i, j ) + x( ix )*temp
subroutine dsyr(UPLO, N, ALPHA, X, INCX, A, LDA)
subroutine xerbla(SRNAME, INFO)