|
version 1.7, 2010/12/21 13:53:37
|
version 1.8, 2011/07/22 07:38:10
|
|
Line 1
|
Line 1
|
| SUBROUTINE DPTCON( N, D, E, ANORM, RCOND, WORK, INFO ) |
SUBROUTINE DPTCON( N, D, E, ANORM, RCOND, 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, N |
INTEGER INFO, N |
|
Line 117
|
Line 117
|
| * 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 126
|
Line 126
|
| WORK( I ) = ONE + WORK( I-1 )*ABS( E( I-1 ) ) |
WORK( I ) = ONE + WORK( I-1 )*ABS( E( I-1 ) ) |
| 20 CONTINUE |
20 CONTINUE |
| * |
* |
| * Solve D * M(L)' * x = b. |
* Solve D * M(L)**T * x = b. |
| * |
* |
| WORK( N ) = WORK( N ) / D( N ) |
WORK( N ) = WORK( N ) / D( N ) |
| DO 30 I = N - 1, 1, -1 |
DO 30 I = N - 1, 1, -1 |
![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to dptcon.f CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [local] / rpl / lapack / lapack