187 SUBROUTINE stpttf( TRANSR, UPLO, N, AP, ARF, INFO )
195 CHARACTER TRANSR, UPLO
199 REAL AP( 0: * ), ARF( 0: * )
206 LOGICAL LOWER, NISODD, NORMALTRANSR
207 INTEGER N1, N2, K, NT
209 INTEGER IJP, JP, LDA, JS
226 normaltransr = lsame( transr,
'N' )
227 lower = lsame( uplo,
'L' )
228 IF( .NOT.normaltransr .AND. .NOT.lsame( transr,
'T' ) )
THEN
230 ELSE IF( .NOT.lower .AND. .NOT.lsame( uplo,
'U' ) )
THEN
232 ELSE IF( n.LT.0 )
THEN
236 CALL xerbla(
'STPTTF', -info )
246 IF( normaltransr )
THEN
274 IF( mod( n, 2 ).EQ.0 )
THEN
285 IF( .NOT.normaltransr )
294 IF( normaltransr )
THEN
307 arf( ij ) = ap( ijp )
315 arf( ij ) = ap( ijp )
328 arf( ij ) = ap( ijp )
337 arf( ij ) = ap( ijp )
355 DO ij = i*( lda+1 ), n*lda - 1, lda
356 arf( ij ) = ap( ijp )
362 DO ij = js, js + n2 - j - 1
363 arf( ij ) = ap( ijp )
377 arf( ij ) = ap( ijp )
383 DO ij = i, i + ( n1+i )*lda, lda
384 arf( ij ) = ap( ijp )
397 IF( normaltransr )
THEN
410 arf( ij ) = ap( ijp )
418 arf( ij ) = ap( ijp )
431 arf( ij ) = ap( ijp )
440 arf( ij ) = ap( ijp )
458 DO ij = i + ( i+1 )*lda, ( n+1 )*lda - 1, lda
459 arf( ij ) = ap( ijp )
465 DO ij = js, js + k - j - 1
466 arf( ij ) = ap( ijp )
480 arf( ij ) = ap( ijp )
486 DO ij = i, i + ( k+i )*lda, lda
487 arf( ij ) = ap( ijp )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine stpttf(TRANSR, UPLO, N, AP, ARF, INFO)
STPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...