version 1.7, 2010/12/21 13:53:36
|
version 1.8, 2011/07/22 07:38:09
|
Line 1
|
Line 1
|
SUBROUTINE DPOCON( UPLO, N, A, LDA, ANORM, RCOND, WORK, IWORK, |
SUBROUTINE DPOCON( UPLO, N, A, LDA, ANORM, RCOND, WORK, IWORK, |
$ INFO ) |
$ 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 -- |
* |
* |
* Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH. |
* Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH. |
* |
* |
Line 129
|
Line 129
|
IF( KASE.NE.0 ) THEN |
IF( KASE.NE.0 ) THEN |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
* |
* |
* Multiply by inv(U'). |
* Multiply by inv(U**T). |
* |
* |
CALL DLATRS( 'Upper', 'Transpose', 'Non-unit', NORMIN, N, A, |
CALL DLATRS( 'Upper', 'Transpose', 'Non-unit', NORMIN, N, A, |
$ LDA, WORK, SCALEL, WORK( 2*N+1 ), INFO ) |
$ LDA, WORK, SCALEL, WORK( 2*N+1 ), INFO ) |
Line 147
|
Line 147
|
$ A, LDA, WORK, SCALEL, WORK( 2*N+1 ), INFO ) |
$ A, LDA, WORK, SCALEL, WORK( 2*N+1 ), INFO ) |
NORMIN = 'Y' |
NORMIN = 'Y' |
* |
* |
* Multiply by inv(L'). |
* Multiply by inv(L**T). |
* |
* |
CALL DLATRS( 'Lower', 'Transpose', 'Non-unit', NORMIN, N, A, |
CALL DLATRS( 'Lower', 'Transpose', 'Non-unit', NORMIN, N, A, |
$ LDA, WORK, SCALEU, WORK( 2*N+1 ), INFO ) |
$ LDA, WORK, SCALEU, WORK( 2*N+1 ), INFO ) |