Mercurial > octave-nkf
diff libcruft/lapack/dtrcon.f @ 7034:68db500cb558
[project @ 2007-10-16 18:54:19 by jwe]
author | jwe |
---|---|
date | Tue, 16 Oct 2007 18:54:23 +0000 |
parents | 6b9cec830d72 |
children |
line wrap: on
line diff
--- a/libcruft/lapack/dtrcon.f Tue Oct 16 17:46:44 2007 +0000 +++ b/libcruft/lapack/dtrcon.f Tue Oct 16 18:54:23 2007 +0000 @@ -1,10 +1,11 @@ SUBROUTINE DTRCON( NORM, UPLO, DIAG, N, A, LDA, RCOND, WORK, $ IWORK, INFO ) * -* -- LAPACK routine (version 3.0) -- -* Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd., -* Courant Institute, Argonne National Lab, and Rice University -* March 31, 1993 +* -- LAPACK routine (version 3.1) -- +* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. +* November 2006 +* +* Modified to call DLACN2 in place of DLACON, 5 Feb 03, SJH. * * .. Scalar Arguments .. CHARACTER DIAG, NORM, UPLO @@ -84,6 +85,9 @@ INTEGER IX, KASE, KASE1 DOUBLE PRECISION AINVNM, ANORM, SCALE, SMLNUM, XNORM * .. +* .. Local Arrays .. + INTEGER ISAVE( 3 ) +* .. * .. External Functions .. LOGICAL LSAME INTEGER IDAMAX @@ -91,7 +95,7 @@ EXTERNAL LSAME, IDAMAX, DLAMCH, DLANTR * .. * .. External Subroutines .. - EXTERNAL DLACON, DLATRS, DRSCL, XERBLA + EXTERNAL DLACN2, DLATRS, DRSCL, XERBLA * .. * .. Intrinsic Functions .. INTRINSIC ABS, DBLE, MAX @@ -150,7 +154,7 @@ END IF KASE = 0 10 CONTINUE - CALL DLACON( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE ) + CALL DLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE ) IF( KASE.NE.0 ) THEN IF( KASE.EQ.KASE1 ) THEN *