103 SUBROUTINE slauu2( UPLO, N, A, LDA, INFO )
122 parameter ( one = 1.0e+0 )
145 upper = lsame( uplo,
'U' )
146 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
148 ELSE IF( n.LT.0 )
THEN
150 ELSE IF( lda.LT.max( 1, n ) )
THEN
154 CALL xerbla(
'SLAUU2', -info )
170 a( i, i ) = sdot( n-i+1, a( i, i ), lda, a( i, i ), lda )
171 CALL sgemv(
'No transpose', i-1, n-i, one, a( 1, i+1 ),
172 $ lda, a( i, i+1 ), lda, aii, a( 1, i ), 1 )
174 CALL sscal( i, aii, a( 1, i ), 1 )
185 a( i, i ) = sdot( n-i+1, a( i, i ), 1, a( i, i ), 1 )
186 CALL sgemv(
'Transpose', n-i, i-1, one, a( i+1, 1 ), lda,
187 $ a( i+1, i ), 1, aii, a( i, 1 ), lda )
189 CALL sscal( i, aii, a( i, 1 ), lda )
subroutine sgemv(TRANS, M, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
SGEMV
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine slauu2(UPLO, N, A, LDA, INFO)
SLAUU2 computes the product UUH or LHL, where U and L are upper or lower triangular matrices (unblock...
subroutine sscal(N, SA, SX, INCX)
SSCAL