version 1.11, 2012/12/14 12:30:33
|
version 1.17, 2017/06/17 11:06:56
|
Line 2
|
Line 2
|
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
* Online html documentation available at |
* Online html documentation available at |
* http://www.netlib.org/lapack/explore-html/ |
* http://www.netlib.org/lapack/explore-html/ |
* |
* |
*> \htmlonly |
*> \htmlonly |
*> Download ZLARTG + dependencies |
*> Download ZLARTG + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zlartg.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zlartg.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zlartg.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zlartg.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zlartg.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zlartg.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
* |
* |
* SUBROUTINE ZLARTG( F, G, CS, SN, R ) |
* SUBROUTINE ZLARTG( F, G, CS, SN, R ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* DOUBLE PRECISION CS |
* DOUBLE PRECISION CS |
* COMPLEX*16 F, G, R, SN |
* COMPLEX*16 F, G, R, SN |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 80
|
Line 80
|
* Authors: |
* Authors: |
* ======== |
* ======== |
* |
* |
*> \author Univ. of Tennessee |
*> \author Univ. of Tennessee |
*> \author Univ. of California Berkeley |
*> \author Univ. of California Berkeley |
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date September 2012 |
*> \date December 2016 |
* |
* |
*> \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.7.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 |
* December 2016 |
* |
* |
* .. 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 |