version 1.6, 2010/08/13 21:03:56
|
version 1.8, 2011/07/22 07:38:10
|
Line 1
|
Line 1
|
SUBROUTINE DPTRFS( N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR, |
SUBROUTINE DPTRFS( N, NRHS, D, E, DF, EF, B, LDB, X, LDX, FERR, |
$ BERR, WORK, INFO ) |
$ BERR, WORK, 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, LDX, N, NRHS |
INTEGER INFO, LDB, LDX, N, NRHS |
Line 263
|
Line 263
|
* 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)**T. |
* |
* |
* Solve M(L) * x = e. |
* Solve M(L) * x = e. |
* |
* |
Line 272
|
Line 272
|
WORK( I ) = ONE + WORK( I-1 )*ABS( EF( I-1 ) ) |
WORK( I ) = ONE + WORK( I-1 )*ABS( EF( I-1 ) ) |
60 CONTINUE |
60 CONTINUE |
* |
* |
* Solve D * M(L)' * x = b. |
* Solve D * M(L)**T * x = b. |
* |
* |
WORK( N ) = WORK( N ) / DF( N ) |
WORK( N ) = WORK( N ) / DF( N ) |
DO 70 I = N - 1, 1, -1 |
DO 70 I = N - 1, 1, -1 |