195 SUBROUTINE dtrttf( TRANSR, UPLO, N, A, LDA, ARF, INFO )
203 CHARACTER TRANSR, UPLO
207 DOUBLE PRECISION A( 0: lda-1, 0: * ), ARF( 0: * )
214 LOGICAL LOWER, NISODD, NORMALTRANSR
215 INTEGER I, IJ, J, K, L, N1, N2, NT, NX2, NP1X2
232 normaltransr = lsame( transr,
'N' )
233 lower = lsame( uplo,
'L' )
234 IF( .NOT.normaltransr .AND. .NOT.lsame( transr,
'T' ) )
THEN
236 ELSE IF( .NOT.lower .AND. .NOT.lsame( uplo,
'U' ) )
THEN
238 ELSE IF( n.LT.0 )
THEN
240 ELSE IF( lda.LT.max( 1, n ) )
THEN
244 CALL xerbla(
'DTRTTF', -info )
275 IF( mod( n, 2 ).EQ.0 )
THEN
290 IF( normaltransr )
THEN
301 arf( ij ) = a( n2+j, i )
305 arf( ij ) = a( i, j )
317 arf( ij ) = a( i, j )
320 DO l = j - n1, n1 - 1
321 arf( ij ) = a( j-n1, l )
340 arf( ij ) = a( j, i )
344 arf( ij ) = a( i, n1+j )
350 arf( ij ) = a( j, i )
362 arf( ij ) = a( j, i )
368 arf( ij ) = a( i, j )
372 arf( ij ) = a( n2+j, l )
385 IF( normaltransr )
THEN
396 arf( ij ) = a( k+j, i )
400 arf( ij ) = a( i, j )
412 arf( ij ) = a( i, j )
416 arf( ij ) = a( j-k, l )
435 arf( ij ) = a( i, j )
440 arf( ij ) = a( j, i )
443 DO i = k + 1 + j, n - 1
444 arf( ij ) = a( i, k+1+j )
450 arf( ij ) = a( j, i )
462 arf( ij ) = a( j, i )
468 arf( ij ) = a( i, j )
471 DO l = k + 1 + j, n - 1
472 arf( ij ) = a( k+1+j, l )
478 arf( ij ) = a( i, j )
subroutine dtrttf(TRANSR, UPLO, N, A, LDA, ARF, INFO)
DTRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...
subroutine xerbla(SRNAME, INFO)
XERBLA