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 xerbla(SRNAME, INFO)
XERBLA
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...
logical function lsame(CA, CB)
LSAME
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j