rpl/lapack/lapack/dgebal.f - diff

Return to dgebal.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.4, 2010/08/06 15:32:22	version 1.5, 2010/08/07 13:18:06
Line 1	Line 1
SUBROUTINE DGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )	SUBROUTINE DGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
*	*
* -- LAPACK routine (version 3.2) --	* -- LAPACK routine (version 3.2.2) --
* -- 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	* June 2010
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
CHARACTER JOB	CHARACTER JOB
Line 116	Line 116
$ SFMIN2	$ SFMIN2
* ..	* ..
* .. External Functions ..	* .. External Functions ..
LOGICAL LSAME	LOGICAL DISNAN, LSAME
INTEGER IDAMAX	INTEGER IDAMAX
DOUBLE PRECISION DLAMCH	DOUBLE PRECISION DLAMCH
EXTERNAL LSAME, IDAMAX, DLAMCH	EXTERNAL DISNAN, LSAME, IDAMAX, DLAMCH
* ..	* ..
* .. External Subroutines ..	* .. External Subroutines ..
EXTERNAL DSCAL, DSWAP, XERBLA	EXTERNAL DSCAL, DSWAP, XERBLA
Line 266	Line 266
160 CONTINUE	160 CONTINUE
IF( C.GE.G .OR. MAX( F, C, CA ).GE.SFMAX2 .OR.	IF( C.GE.G .OR. MAX( F, C, CA ).GE.SFMAX2 .OR.
$ MIN( R, G, RA ).LE.SFMIN2 )GO TO 170	$ MIN( R, G, RA ).LE.SFMIN2 )GO TO 170
	IF( DISNAN( C+F+CA+R+G+RA ) ) THEN
	*
	* Exit if NaN to avoid infinite loop
	*
	INFO = -3
	CALL XERBLA( 'DGEBAL', -INFO )
	RETURN
	END IF
F = F*SCLFAC	F = F*SCLFAC
C = C*SCLFAC	C = C*SCLFAC
CA = CA*SCLFAC	CA = CA*SCLFAC