1 SUBROUTINE dpotrs( UPLO, N, NRHS, A, LDA, B, LDB, INFO )
10 INTEGER INFO, LDA, LDB, N, NRHS
13 DOUBLE PRECISION A( LDA, * ), B( LDB, * )
59 parameter( one = 1.0d+0 )
79 upper = lsame( uplo,
'U' )
80 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
82 ELSE IF( n.LT.0 )
THEN
84 ELSE IF( nrhs.LT.0 )
THEN
86 ELSE IF( lda.LT.max( 1, n ) )
THEN
88 ELSE IF( ldb.LT.max( 1, n ) )
THEN
92 CALL xerbla(
'DPOTRS', -info )
98 IF( n.EQ.0 .OR. nrhs.EQ.0 )
107 CALL dtrsm(
'Left',
'Upper',
'Transpose',
'Non-unit', n, nrhs,
108 $ one, a, lda, b, ldb )
112 CALL dtrsm(
'Left',
'Upper',
'No transpose',
'Non-unit', n,
113 $ nrhs, one, a, lda, b, ldb )
120 CALL dtrsm(
'Left',
'Lower',
'No transpose',
'Non-unit', n,
121 $ nrhs, one, a, lda, b, ldb )
125 CALL dtrsm(
'Left',
'Lower',
'Transpose',
'Non-unit', n, nrhs,
126 $ one, a, lda, b, ldb )
subroutine dpotrs(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
subroutine dtrsm(SIDE, UPLO, TRANSA, DIAG, M, N, ALPHA, A, LDA, B, LDB)
subroutine xerbla(SRNAME, INFO)