version 1.9, 2011/11/21 22:19:48
|
version 1.14, 2015/11/26 11:44:24
|
Line 266
|
Line 266
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date April 2012 |
* |
* |
*> \ingroup complex16GEcomputational |
*> \ingroup complex16GEcomputational |
* |
* |
Line 284
|
Line 284
|
$ ALPHA, BETA, Q, LDQ, Z, LDZ, WORK, LWORK, |
$ ALPHA, BETA, Q, LDQ, Z, LDZ, WORK, LWORK, |
$ RWORK, INFO ) |
$ RWORK, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.6.0) -- |
* -- 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 454
|
Line 454
|
CALL ZSCAL( J-1, SIGNBC, T( 1, J ), 1 ) |
CALL ZSCAL( J-1, SIGNBC, T( 1, J ), 1 ) |
CALL ZSCAL( J, SIGNBC, H( 1, J ), 1 ) |
CALL ZSCAL( J, SIGNBC, H( 1, J ), 1 ) |
ELSE |
ELSE |
H( J, J ) = H( J, J )*SIGNBC |
CALL ZSCAL( 1, SIGNBC, H( J, J ), 1 ) |
END IF |
END IF |
IF( ILZ ) |
IF( ILZ ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, J ), 1 ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, J ), 1 ) |
Line 666
|
Line 666
|
CALL ZSCAL( ILAST+1-IFRSTM, SIGNBC, H( IFRSTM, ILAST ), |
CALL ZSCAL( ILAST+1-IFRSTM, SIGNBC, H( IFRSTM, ILAST ), |
$ 1 ) |
$ 1 ) |
ELSE |
ELSE |
H( ILAST, ILAST ) = H( ILAST, ILAST )*SIGNBC |
CALL ZSCAL( 1, SIGNBC, H( ILAST, ILAST ), 1 ) |
END IF |
END IF |
IF( ILZ ) |
IF( ILZ ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, ILAST ), 1 ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, ILAST ), 1 ) |
Line 744
|
Line 744
|
* |
* |
* Exceptional shift. Chosen for no particularly good reason. |
* Exceptional shift. Chosen for no particularly good reason. |
* |
* |
ESHIFT = ESHIFT + DCONJG( ( ASCALE*H( ILAST-1, ILAST ) ) / |
ESHIFT = ESHIFT + (ASCALE*H(ILAST,ILAST-1))/ |
$ ( BSCALE*T( ILAST-1, ILAST-1 ) ) ) |
$ (BSCALE*T(ILAST-1,ILAST-1)) |
SHIFT = ESHIFT |
SHIFT = ESHIFT |
END IF |
END IF |
* |
* |
Line 850
|
Line 850
|
CALL ZSCAL( J-1, SIGNBC, T( 1, J ), 1 ) |
CALL ZSCAL( J-1, SIGNBC, T( 1, J ), 1 ) |
CALL ZSCAL( J, SIGNBC, H( 1, J ), 1 ) |
CALL ZSCAL( J, SIGNBC, H( 1, J ), 1 ) |
ELSE |
ELSE |
H( J, J ) = H( J, J )*SIGNBC |
CALL ZSCAL( 1, SIGNBC, H( J, J ), 1 ) |
END IF |
END IF |
IF( ILZ ) |
IF( ILZ ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, J ), 1 ) |
$ CALL ZSCAL( N, SIGNBC, Z( 1, J ), 1 ) |