1 SUBROUTINE strti2( UPLO, DIAG, N, A, LDA, INFO )
65 parameter( one = 1.0e+0 )
77 EXTERNAL sscal, strmv,
xerbla
87 upper = lsame(
uplo,
'U' )
88 nounit = lsame(
diag,
'N' )
89 IF( .NOT.upper .AND. .NOT.lsame(
uplo,
'L' ) )
THEN
91 ELSE IF( .NOT.nounit .AND. .NOT.lsame(
diag,
'U' ) )
THEN
93 ELSE IF( n.LT.0 )
THEN
95 ELSE IF( lda.LT.
max( 1, n ) )
THEN
99 CALL
xerbla(
'STRTI2', -info )
109 a( j, j ) = one / a( j, j )
117 CALL strmv(
'Upper',
'No transpose',
diag, j-1, a, lda,
119 CALL sscal( j-1, ajj, a( 1, j ), 1 )
127 a( j, j ) = one / a( j, j )
136 CALL strmv(
'Lower',
'No transpose',
diag, n-j,
137 $ a( j+1, j+1 ), lda, a( j+1, j ), 1 )
138 CALL sscal( n-j, ajj, a( j+1, j ), 1 )