217 SUBROUTINE ctrttf( TRANSR, UPLO, N, A, LDA, ARF, INFO )
225 CHARACTER TRANSR, UPLO
229 COMPLEX A( 0: lda-1, 0: * ), ARF( 0: * )
237 LOGICAL LOWER, NISODD, NORMALTRANSR
238 INTEGER I, IJ, J, K, L, N1, N2, NT, NX2, NP1X2
248 INTRINSIC conjg, max, mod
255 normaltransr = lsame( transr,
'N' )
256 lower = lsame( uplo,
'L' )
257 IF( .NOT.normaltransr .AND. .NOT.lsame( transr,
'C' ) )
THEN
259 ELSE IF( .NOT.lower .AND. .NOT.lsame( uplo,
'U' ) )
THEN
261 ELSE IF( n.LT.0 )
THEN
263 ELSE IF( lda.LT.max( 1, n ) )
THEN
267 CALL xerbla(
'CTRTTF', -info )
275 IF( normaltransr )
THEN
278 arf( 0 ) = conjg( a( 0, 0 ) )
302 IF( mod( n, 2 ).EQ.0 )
THEN
317 IF( normaltransr )
THEN
330 arf( ij ) = conjg( a( n2+j, i ) )
334 arf( ij ) = a( i, j )
348 arf( ij ) = a( i, j )
351 DO l = j - n1, n1 - 1
352 arf( ij ) = conjg( a( j-n1, l ) )
373 arf( ij ) = conjg( a( j, i ) )
377 arf( ij ) = a( i, n1+j )
383 arf( ij ) = conjg( a( j, i ) )
397 arf( ij ) = conjg( a( j, i ) )
403 arf( ij ) = a( i, j )
407 arf( ij ) = conjg( a( n2+j, l ) )
420 IF( normaltransr )
THEN
433 arf( ij ) = conjg( a( k+j, i ) )
437 arf( ij ) = a( i, j )
451 arf( ij ) = a( i, j )
455 arf( ij ) = conjg( a( j-k, l ) )
476 arf( ij ) = a( i, j )
481 arf( ij ) = conjg( a( j, i ) )
484 DO i = k + 1 + j, n - 1
485 arf( ij ) = a( i, k+1+j )
491 arf( ij ) = conjg( a( j, i ) )
505 arf( ij ) = conjg( a( j, i ) )
511 arf( ij ) = a( i, j )
514 DO l = k + 1 + j, n - 1
515 arf( ij ) = conjg( a( k+1+j, l ) )
523 arf( ij ) = a( i, j )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ctrttf(TRANSR, UPLO, N, A, LDA, ARF, INFO)
CTRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...