version 1.5, 2010/08/07 13:18:08
|
version 1.9, 2011/07/22 07:38:11
|
Line 1
|
Line 1
|
SUBROUTINE DSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, |
SUBROUTINE DSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, |
$ IWORK, INFO ) |
$ IWORK, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2.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..-- |
* June 2010 |
* -- 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 148
|
Line 148
|
CALL DLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE ) |
CALL DLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE ) |
IF( KASE.NE.0 ) THEN |
IF( KASE.NE.0 ) THEN |
* |
* |
* Multiply by inv(L*D*L') or inv(U*D*U'). |
* Multiply by inv(L*D*L**T) or inv(U*D*U**T). |
* |
* |
CALL DSYTRS( UPLO, N, 1, A, LDA, IPIV, WORK, N, INFO ) |
CALL DSYTRS( UPLO, N, 1, A, LDA, IPIV, WORK, N, INFO ) |
GO TO 30 |
GO TO 30 |