version 1.12, 2012/12/14 14:22:52
|
version 1.13, 2014/01/27 09:24:36
|
Line 85
|
Line 85
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date September 2012 |
*> \date November 2013 |
* |
* |
*> \ingroup complex16OTHERauxiliary |
*> \ingroup complex16OTHERauxiliary |
* |
* |
Line 103
|
Line 103
|
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZLARTG( F, G, CS, SN, R ) |
SUBROUTINE ZLARTG( F, G, CS, SN, R ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.4.2) -- |
* -- LAPACK auxiliary routine (version 3.5.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..-- |
* September 2012 |
* November 2013 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
DOUBLE PRECISION CS |
DOUBLE PRECISION CS |
Line 130
|
Line 130
|
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
DOUBLE PRECISION DLAMCH, DLAPY2 |
DOUBLE PRECISION DLAMCH, DLAPY2 |
EXTERNAL DLAMCH, DLAPY2 |
LOGICAL DISNAN |
|
EXTERNAL DLAMCH, DLAPY2, DISNAN |
* .. |
* .. |
* .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
INTRINSIC ABS, DBLE, DCMPLX, DCONJG, DIMAG, INT, LOG, |
INTRINSIC ABS, DBLE, DCMPLX, DCONJG, DIMAG, INT, LOG, |
Line 139
|
Line 140
|
* .. Statement Functions .. |
* .. Statement Functions .. |
DOUBLE PRECISION ABS1, ABSSQ |
DOUBLE PRECISION ABS1, ABSSQ |
* .. |
* .. |
* .. Save statement .. |
|
* SAVE FIRST, SAFMX2, SAFMIN, SAFMN2 |
|
* .. |
|
* .. Data statements .. |
|
* DATA FIRST / .TRUE. / |
|
* .. |
|
* .. Statement Function definitions .. |
* .. Statement Function definitions .. |
ABS1( FF ) = MAX( ABS( DBLE( FF ) ), ABS( DIMAG( FF ) ) ) |
ABS1( FF ) = MAX( ABS( DBLE( FF ) ), ABS( DIMAG( FF ) ) ) |
ABSSQ( FF ) = DBLE( FF )**2 + DIMAG( FF )**2 |
ABSSQ( FF ) = DBLE( FF )**2 + DIMAG( FF )**2 |
* .. |
* .. |
* .. Executable Statements .. |
* .. Executable Statements .. |
* |
* |
* IF( FIRST ) THEN |
SAFMIN = DLAMCH( 'S' ) |
SAFMIN = DLAMCH( 'S' ) |
EPS = DLAMCH( 'E' ) |
EPS = DLAMCH( 'E' ) |
SAFMN2 = DLAMCH( 'B' )**INT( LOG( SAFMIN / EPS ) / |
SAFMN2 = DLAMCH( 'B' )**INT( LOG( SAFMIN / EPS ) / |
$ LOG( DLAMCH( 'B' ) ) / TWO ) |
$ LOG( DLAMCH( 'B' ) ) / TWO ) |
SAFMX2 = ONE / SAFMN2 |
SAFMX2 = ONE / SAFMN2 |
|
* FIRST = .FALSE. |
|
* END IF |
|
SCALE = MAX( ABS1( F ), ABS1( G ) ) |
SCALE = MAX( ABS1( F ), ABS1( G ) ) |
FS = F |
FS = F |
GS = G |
GS = G |
Line 172
|
Line 164
|
IF( SCALE.GE.SAFMX2 ) |
IF( SCALE.GE.SAFMX2 ) |
$ GO TO 10 |
$ GO TO 10 |
ELSE IF( SCALE.LE.SAFMN2 ) THEN |
ELSE IF( SCALE.LE.SAFMN2 ) THEN |
IF( G.EQ.CZERO ) THEN |
IF( G.EQ.CZERO.OR.DISNAN( ABS( G ) ) ) THEN |
CS = ONE |
CS = ONE |
SN = CZERO |
SN = CZERO |
R = F |
R = F |