200 SUBROUTINE dpftrs( TRANSR, UPLO, N, NRHS, A, B, LDB, INFO )
208 CHARACTER TRANSR, UPLO
209 INTEGER INFO, LDB, N, NRHS
212 DOUBLE PRECISION A( 0: * ), B( ldb, * )
219 parameter ( one = 1.0d+0 )
222 LOGICAL LOWER, NORMALTRANSR
239 normaltransr = lsame( transr,
'N' )
240 lower = lsame( uplo,
'L' )
241 IF( .NOT.normaltransr .AND. .NOT.lsame( transr,
'T' ) )
THEN
243 ELSE IF( .NOT.lower .AND. .NOT.lsame( uplo,
'U' ) )
THEN
245 ELSE IF( n.LT.0 )
THEN
247 ELSE IF( nrhs.LT.0 )
THEN
249 ELSE IF( ldb.LT.max( 1, n ) )
THEN
253 CALL xerbla(
'DPFTRS', -info )
259 IF( n.EQ.0 .OR. nrhs.EQ.0 )
265 CALL dtfsm( transr,
'L', uplo,
'N',
'N', n, nrhs, one, a, b,
267 CALL dtfsm( transr,
'L', uplo,
'T',
'N', n, nrhs, one, a, b,
270 CALL dtfsm( transr,
'L', uplo,
'T',
'N', n, nrhs, one, a, b,
272 CALL dtfsm( transr,
'L', uplo,
'N',
'N', n, nrhs, one, a, b,
subroutine dpftrs(TRANSR, UPLO, N, NRHS, A, B, LDB, INFO)
DPFTRS
subroutine dtfsm(TRANSR, SIDE, UPLO, TRANS, DIAG, M, N, ALPHA, A, B, LDB)
DTFSM solves a matrix equation (one operand is a triangular matrix in RFP format).
subroutine xerbla(SRNAME, INFO)
XERBLA