version 1.2, 2010/04/21 13:45:38
|
version 1.8, 2011/07/22 07:38:19
|
Line 1
|
Line 1
|
SUBROUTINE ZPTCON( N, D, E, ANORM, RCOND, RWORK, INFO ) |
SUBROUTINE ZPTCON( N, D, E, ANORM, RCOND, RWORK, 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, N |
INTEGER INFO, N |
Line 118
|
Line 118
|
* m(i,j) = abs(A(i,j)), i = j, |
* m(i,j) = abs(A(i,j)), i = j, |
* m(i,j) = -abs(A(i,j)), i .ne. j, |
* m(i,j) = -abs(A(i,j)), i .ne. j, |
* |
* |
* and e = [ 1, 1, ..., 1 ]'. Note M(A) = M(L)*D*M(L)'. |
* and e = [ 1, 1, ..., 1 ]**T. Note M(A) = M(L)*D*M(L)**H. |
* |
* |
* Solve M(L) * x = e. |
* Solve M(L) * x = e. |
* |
* |
Line 127
|
Line 127
|
RWORK( I ) = ONE + RWORK( I-1 )*ABS( E( I-1 ) ) |
RWORK( I ) = ONE + RWORK( I-1 )*ABS( E( I-1 ) ) |
20 CONTINUE |
20 CONTINUE |
* |
* |
* Solve D * M(L)' * x = b. |
* Solve D * M(L)**H * x = b. |
* |
* |
RWORK( N ) = RWORK( N ) / D( N ) |
RWORK( N ) = RWORK( N ) / D( N ) |
DO 30 I = N - 1, 1, -1 |
DO 30 I = N - 1, 1, -1 |