version 1.15, 2017/06/17 10:54:30
|
version 1.18, 2023/08/07 08:39:40
|
Line 351
|
Line 351
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16OTHERcomputational |
*> \ingroup complex16OTHERcomputational |
* |
* |
*> \par Further Details: |
*> \par Further Details: |
Line 379
|
Line 377
|
$ LDB, TOLA, TOLB, ALPHA, BETA, U, LDU, V, LDV, |
$ LDB, TOLA, TOLB, ALPHA, BETA, U, LDU, V, LDV, |
$ Q, LDQ, WORK, NCYCLE, INFO ) |
$ Q, LDQ, WORK, NCYCLE, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.7.0) -- |
* -- LAPACK computational routine -- |
* -- 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..-- |
* December 2016 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOBQ, JOBU, JOBV |
CHARACTER JOBQ, JOBU, JOBV |
Line 401
|
Line 398
|
* .. Parameters .. |
* .. Parameters .. |
INTEGER MAXIT |
INTEGER MAXIT |
PARAMETER ( MAXIT = 40 ) |
PARAMETER ( MAXIT = 40 ) |
DOUBLE PRECISION ZERO, ONE |
DOUBLE PRECISION ZERO, ONE, HUGENUM |
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 ) |
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 ) |
COMPLEX*16 CZERO, CONE |
COMPLEX*16 CZERO, CONE |
PARAMETER ( CZERO = ( 0.0D+0, 0.0D+0 ), |
PARAMETER ( CZERO = ( 0.0D+0, 0.0D+0 ), |
Line 424
|
Line 421
|
$ ZLASET, ZROT |
$ ZLASET, ZROT |
* .. |
* .. |
* .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
INTRINSIC ABS, DBLE, DCONJG, MAX, MIN |
INTRINSIC ABS, DBLE, DCONJG, MAX, MIN, HUGE |
|
PARAMETER ( HUGENUM = HUGE(ZERO) ) |
* .. |
* .. |
* .. Executable Statements .. |
* .. Executable Statements .. |
* |
* |
Line 610
|
Line 608
|
* |
* |
A1 = DBLE( A( K+I, N-L+I ) ) |
A1 = DBLE( A( K+I, N-L+I ) ) |
B1 = DBLE( B( I, N-L+I ) ) |
B1 = DBLE( B( I, N-L+I ) ) |
|
GAMMA = B1 / A1 |
* |
* |
IF( A1.NE.ZERO ) THEN |
IF( (GAMMA.LE.HUGENUM).AND.(GAMMA.GE.-HUGENUM) ) THEN |
GAMMA = B1 / A1 |
|
* |
* |
IF( GAMMA.LT.ZERO ) THEN |
IF( GAMMA.LT.ZERO ) THEN |
CALL ZDSCAL( L-I+1, -ONE, B( I, N-L+I ), LDB ) |
CALL ZDSCAL( L-I+1, -ONE, B( I, N-L+I ), LDB ) |