version 1.3, 2010/08/06 15:28:37
|
version 1.8, 2011/07/22 07:38:05
|
Line 1
|
Line 1
|
SUBROUTINE DGETRS( TRANS, N, NRHS, A, LDA, IPIV, B, LDB, INFO ) |
SUBROUTINE DGETRS( TRANS, N, NRHS, A, LDA, IPIV, 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 TRANS |
CHARACTER TRANS |
Line 18
|
Line 18
|
* ======= |
* ======= |
* |
* |
* DGETRS solves a system of linear equations |
* DGETRS solves a system of linear equations |
* A * X = B or A' * X = B |
* A * X = B or A**T * X = B |
* with a general N-by-N matrix A using the LU factorization computed |
* with a general N-by-N matrix A using the LU factorization computed |
* by DGETRF. |
* by DGETRF. |
* |
* |
Line 28
|
Line 28
|
* TRANS (input) CHARACTER*1 |
* TRANS (input) CHARACTER*1 |
* Specifies the form of the system of equations: |
* Specifies the form of the system of equations: |
* = 'N': A * X = B (No transpose) |
* = 'N': A * X = B (No transpose) |
* = 'T': A'* X = B (Transpose) |
* = 'T': A**T* X = B (Transpose) |
* = 'C': A'* X = B (Conjugate transpose = Transpose) |
* = 'C': A**T* X = B (Conjugate transpose = Transpose) |
* |
* |
* N (input) INTEGER |
* N (input) INTEGER |
* The order of the matrix A. N >= 0. |
* The order of the matrix A. N >= 0. |
Line 126
|
Line 126
|
$ NRHS, ONE, A, LDA, B, LDB ) |
$ NRHS, ONE, A, LDA, B, LDB ) |
ELSE |
ELSE |
* |
* |
* Solve A' * X = B. |
* Solve A**T * X = B. |
* |
* |
* 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 ) |
* |
* |
* Solve L'*X = B, overwriting B with X. |
* Solve L**T *X = B, overwriting B with X. |
* |
* |
CALL DTRSM( 'Left', 'Lower', 'Transpose', 'Unit', N, NRHS, ONE, |
CALL DTRSM( 'Left', 'Lower', 'Transpose', 'Unit', N, NRHS, ONE, |
$ A, LDA, B, LDB ) |
$ A, LDA, B, LDB ) |