rpl/lapack/lapack/zgebal.f - diff

Return to zgebal.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.4, 2010/08/06 15:32:38	version 1.5, 2010/08/07 13:18:08
Line 1	Line 1
SUBROUTINE ZGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )	SUBROUTINE ZGEBAL( 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 118	Line 118
COMPLEX*16 CDUM	COMPLEX*16 CDUM
* ..	* ..
* .. External Functions ..	* .. External Functions ..
LOGICAL LSAME	LOGICAL DISNAN, LSAME
INTEGER IZAMAX	INTEGER IZAMAX
DOUBLE PRECISION DLAMCH	DOUBLE PRECISION DLAMCH
EXTERNAL LSAME, IZAMAX, DLAMCH	EXTERNAL DISNAN, LSAME, IZAMAX, DLAMCH
* ..	* ..
* .. External Subroutines ..	* .. External Subroutines ..
EXTERNAL XERBLA, ZDSCAL, ZSWAP	EXTERNAL XERBLA, ZDSCAL, ZSWAP
Line 274	Line 274
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( 'ZGEBAL', -INFO )
	RETURN
	END IF
F = F*SCLFAC	F = F*SCLFAC
C = C*SCLFAC	C = C*SCLFAC
CA = CA*SCLFAC	CA = CA*SCLFAC