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 )
real function sdot(N, SX, INCX, SY, INCY)
SDOT
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...
logical function lsame(CA, CB)
LSAME
subroutine sgemv(TRANS, M, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
SGEMV
subroutine sscal(N, SA, SX, INCX)
SSCAL