version 1.4, 2010/08/06 15:32:48
|
version 1.8, 2011/07/22 07:38:19
|
Line 1
|
Line 1
|
SUBROUTINE ZPTSV( N, NRHS, D, E, B, LDB, INFO ) |
SUBROUTINE ZPTSV( N, NRHS, D, E, 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 .. |
INTEGER INFO, LDB, N, NRHS |
INTEGER INFO, LDB, N, NRHS |
Line 45
|
Line 45
|
* A. E can also be regarded as the superdiagonal of the unit |
* A. E can also be regarded as the superdiagonal of the unit |
* bidiagonal factor U from the U**H*D*U factorization of A. |
* bidiagonal factor U from the U**H*D*U factorization of A. |
* |
* |
* B (input/output) COMPLEX*16 array, dimension (LDB,N) |
* B (input/output) COMPLEX*16 array, dimension (LDB,NRHS) |
* On entry, the N-by-NRHS right hand side matrix B. |
* On entry, the N-by-NRHS right hand side matrix B. |
* On exit, if INFO = 0, the N-by-NRHS solution matrix X. |
* On exit, if INFO = 0, the N-by-NRHS solution matrix X. |
* |
* |
Line 85
|
Line 85
|
RETURN |
RETURN |
END IF |
END IF |
* |
* |
* Compute the L*D*L' (or U'*D*U) factorization of A. |
* Compute the L*D*L**H (or U**H*D*U) factorization of A. |
* |
* |
CALL ZPTTRF( N, D, E, INFO ) |
CALL ZPTTRF( N, D, E, INFO ) |
IF( INFO.EQ.0 ) THEN |
IF( INFO.EQ.0 ) THEN |