rpl/lapack/lapack/dhgeqz.f - diff

Return to dhgeqz.f CVS log

Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dhgeqz.f between versions 1.10 and 1.11

version 1.10, 2011/11/21 22:19:29	version 1.11, 2012/07/31 11:06:35
Line 282	Line 282
*> \author Univ. of Colorado Denver	*> \author Univ. of Colorado Denver
*> \author NAG Ltd.	*> \author NAG Ltd.
*	*
*> \date November 2011	*> \date April 2012
*	*
*> \ingroup doubleGEcomputational	*> \ingroup doubleGEcomputational
*	*
Line 304	Line 304
$ ALPHAR, ALPHAI, BETA, Q, LDQ, Z, LDZ, WORK,	$ ALPHAR, ALPHAI, BETA, Q, LDQ, Z, LDZ, WORK,
$ LWORK, INFO )	$ LWORK, INFO )
*	*
* -- LAPACK computational routine (version 3.4.0) --	* -- LAPACK computational routine (version 3.4.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 2011	* April 2012
*	*
* .. Scalar Arguments ..	* .. Scalar Arguments ..
CHARACTER COMPQ, COMPZ, JOB	CHARACTER COMPQ, COMPZ, JOB
Line 741	Line 741
*	*
IF( ( DBLE( MAXIT )SAFMIN )ABS( H( ILAST-1, ILAST ) ).LT.	IF( ( DBLE( MAXIT )SAFMIN )ABS( H( ILAST-1, ILAST ) ).LT.
$ ABS( T( ILAST-1, ILAST-1 ) ) ) THEN	$ ABS( T( ILAST-1, ILAST-1 ) ) ) THEN
ESHIFT = ESHIFT + H( ILAST-1, ILAST ) /	ESHIFT = ESHIFT + H( ILAST, ILAST-1 ) /
$ T( ILAST-1, ILAST-1 )	$ T( ILAST-1, ILAST-1 )
ELSE	ELSE
ESHIFT = ESHIFT + ONE / ( SAFMIN*DBLE( MAXIT ) )	ESHIFT = ESHIFT + ONE / ( SAFMIN*DBLE( MAXIT ) )
Line 920	Line 920
Z( J, ILAST ) = -Z( J, ILAST )	Z( J, ILAST ) = -Z( J, ILAST )
220 CONTINUE	220 CONTINUE
END IF	END IF
	B22 = -B22
END IF	END IF
*	*
* Step 2: Compute ALPHAR, ALPHAI, and BETA (see refs.)	* Step 2: Compute ALPHAR, ALPHAI, and BETA (see refs.)

CVSweb interface <joel.bertrand@systella.fr>

Removed from v.1.10
changed lines
	Added in v.1.11