version 1.15, 2017/06/17 11:06:56
|
version 1.18, 2023/08/07 08:39:31
|
Line 97
|
Line 97
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16OTHERauxiliary |
*> \ingroup complex16OTHERauxiliary |
* |
* |
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZLARFGP( N, ALPHA, X, INCX, TAU ) |
SUBROUTINE ZLARFGP( N, ALPHA, X, INCX, TAU ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.7.0) -- |
* -- LAPACK auxiliary 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 .. |
INTEGER INCX, N |
INTEGER INCX, N |
Line 197
|
Line 194
|
BETA = BETA*BIGNUM |
BETA = BETA*BIGNUM |
ALPHI = ALPHI*BIGNUM |
ALPHI = ALPHI*BIGNUM |
ALPHR = ALPHR*BIGNUM |
ALPHR = ALPHR*BIGNUM |
IF( ABS( BETA ).LT.SMLNUM ) |
IF( (ABS( BETA ).LT.SMLNUM) .AND. (KNT .LT. 20) ) |
$ GO TO 10 |
$ GO TO 10 |
* |
* |
* New BETA is at most 1, at least SMLNUM |
* New BETA is at most 1, at least SMLNUM |
Line 238
|
Line 235
|
DO J = 1, N-1 |
DO J = 1, N-1 |
X( 1 + (J-1)*INCX ) = ZERO |
X( 1 + (J-1)*INCX ) = ZERO |
END DO |
END DO |
BETA = -SAVEALPHA |
BETA = DBLE( -SAVEALPHA ) |
END IF |
END IF |
ELSE |
ELSE |
XNORM = DLAPY2( ALPHR, ALPHI ) |
XNORM = DLAPY2( ALPHR, ALPHI ) |