version 1.11, 2012/08/22 09:48:33
|
version 1.15, 2016/08/27 15:27:13
|
Line 190
|
Line 190
|
*> \param[in,out] Q |
*> \param[in,out] Q |
*> \verbatim |
*> \verbatim |
*> Q is COMPLEX*16 array, dimension (LDQ, N) |
*> Q is COMPLEX*16 array, dimension (LDQ, N) |
*> On entry, if COMPZ = 'V', the unitary matrix Q1 used in the |
*> On entry, if COMPQ = 'V', the unitary matrix Q1 used in the |
*> reduction of (A,B) to generalized Hessenberg form. |
*> reduction of (A,B) to generalized Hessenberg form. |
*> On exit, if COMPZ = 'I', the unitary matrix of left Schur |
*> On exit, if COMPQ = 'I', the unitary matrix of left Schur |
*> vectors of (H,T), and if COMPZ = 'V', the unitary matrix of |
*> vectors of (H,T), and if COMPQ = 'V', the unitary matrix of |
*> left Schur vectors of (A,B). |
*> left Schur vectors of (A,B). |
*> Not referenced if COMPZ = 'N'. |
*> Not referenced if COMPQ = 'N'. |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in] LDQ |
*> \param[in] LDQ |
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.1) -- |
* -- LAPACK computational routine (version 3.6.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..-- |
* April 2012 |
* April 2012 |
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 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 ) |