version 1.7, 2010/12/21 13:53:36
|
version 1.8, 2011/07/22 07:38:10
|
Line 1
|
Line 1
|
SUBROUTINE DPOTRS( UPLO, N, NRHS, A, LDA, B, LDB, INFO ) |
SUBROUTINE DPOTRS( UPLO, N, NRHS, A, LDA, B, LDB, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2) -- |
* -- LAPACK routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 100
|
Line 100
|
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
* |
* |
* Solve A*X = B where A = U'*U. |
* Solve A*X = B where A = U**T *U. |
* |
* |
* Solve U'*X = B, overwriting B with X. |
* Solve U**T *X = B, overwriting B with X. |
* |
* |
CALL DTRSM( 'Left', 'Upper', 'Transpose', 'Non-unit', N, NRHS, |
CALL DTRSM( 'Left', 'Upper', 'Transpose', 'Non-unit', N, NRHS, |
$ ONE, A, LDA, B, LDB ) |
$ ONE, A, LDA, B, LDB ) |
Line 113
|
Line 113
|
$ NRHS, ONE, A, LDA, B, LDB ) |
$ NRHS, ONE, A, LDA, B, LDB ) |
ELSE |
ELSE |
* |
* |
* Solve A*X = B where A = L*L'. |
* Solve A*X = B where A = L*L**T. |
* |
* |
* Solve L*X = B, overwriting B with X. |
* Solve L*X = B, overwriting B with X. |
* |
* |
CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', N, |
CALL DTRSM( 'Left', 'Lower', 'No transpose', 'Non-unit', N, |
$ NRHS, ONE, A, LDA, B, LDB ) |
$ NRHS, ONE, A, LDA, B, LDB ) |
* |
* |
* Solve L'*X = B, overwriting B with X. |
* Solve L**T *X = B, overwriting B with X. |
* |
* |
CALL DTRSM( 'Left', 'Lower', 'Transpose', 'Non-unit', N, NRHS, |
CALL DTRSM( 'Left', 'Lower', 'Transpose', 'Non-unit', N, NRHS, |
$ ONE, A, LDA, B, LDB ) |
$ ONE, A, LDA, B, LDB ) |