version 1.4, 2010/08/06 15:32:48
|
version 1.8, 2011/07/22 07:38:19
|
Line 1
|
Line 1
|
SUBROUTINE ZPTRFS( UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, |
SUBROUTINE ZPTRFS( UPLO, N, NRHS, D, E, DF, EF, B, LDB, X, LDX, |
$ FERR, BERR, WORK, RWORK, INFO ) |
$ FERR, BERR, WORK, 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 .. |
CHARACTER UPLO |
CHARACTER UPLO |
Line 327
|
Line 327
|
* 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 336
|
Line 336
|
RWORK( I ) = ONE + RWORK( I-1 )*ABS( EF( I-1 ) ) |
RWORK( I ) = ONE + RWORK( I-1 )*ABS( EF( I-1 ) ) |
70 CONTINUE |
70 CONTINUE |
* |
* |
* Solve D * M(L)' * x = b. |
* Solve D * M(L)**H * x = b. |
* |
* |
RWORK( N ) = RWORK( N ) / DF( N ) |
RWORK( N ) = RWORK( N ) / DF( N ) |
DO 80 I = N - 1, 1, -1 |
DO 80 I = N - 1, 1, -1 |