197 SUBROUTINE stfttr( TRANSR, UPLO, N, ARF, A, LDA, INFO )
205 CHARACTER TRANSR, UPLO
209 REAL A( 0: lda-1, 0: * ), ARF( 0: * )
216 LOGICAL LOWER, NISODD, NORMALTRANSR
217 INTEGER N1, N2, K, NT, NX2, NP1X2
235 normaltransr = lsame( transr,
'N' )
236 lower = lsame( uplo,
'L' )
237 IF( .NOT.normaltransr .AND. .NOT.lsame( transr,
'T' ) )
THEN
239 ELSE IF( .NOT.lower .AND. .NOT.lsame( uplo,
'U' ) )
THEN
241 ELSE IF( n.LT.0 )
THEN
243 ELSE IF( lda.LT.max( 1, n ) )
THEN
247 CALL xerbla(
'STFTTR', -info )
278 IF( mod( n, 2 ).EQ.0 )
THEN
293 IF( normaltransr )
THEN
304 a( n2+j, i ) = arf( ij )
308 a( i, j ) = arf( ij )
320 a( i, j ) = arf( ij )
323 DO l = j - n1, n1 - 1
324 a( j-n1, l ) = arf( ij )
343 a( j, i ) = arf( ij )
347 a( i, n1+j ) = arf( ij )
353 a( j, i ) = arf( ij )
365 a( j, i ) = arf( ij )
371 a( i, j ) = arf( ij )
375 a( n2+j, l ) = arf( ij )
388 IF( normaltransr )
THEN
399 a( k+j, i ) = arf( ij )
403 a( i, j ) = arf( ij )
415 a( i, j ) = arf( ij )
419 a( j-k, l ) = arf( ij )
438 a( i, j ) = arf( ij )
443 a( j, i ) = arf( ij )
446 DO i = k + 1 + j, n - 1
447 a( i, k+1+j ) = arf( ij )
453 a( j, i ) = arf( ij )
465 a( j, i ) = arf( ij )
471 a( i, j ) = arf( ij )
474 DO l = k + 1 + j, n - 1
475 a( k+1+j, l ) = arf( ij )
481 a( i, j ) = arf( ij )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine stfttr(TRANSR, UPLO, N, ARF, A, LDA, INFO)
STFTTR copies a triangular matrix from the rectangular full packed format (TF) to the standard full f...