version 1.11, 2017/06/17 10:53:56
|
version 1.14, 2023/08/07 08:38:58
|
Line 88
|
Line 88
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup OTHERauxiliary |
*> \ingroup OTHERauxiliary |
* |
* |
* ===================================================================== |
* ===================================================================== |
SUBROUTINE DLARTGP( F, G, CS, SN, R ) |
SUBROUTINE DLARTGP( F, G, CS, SN, R ) |
* |
* |
* -- 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 .. |
DOUBLE PRECISION CS, F, G, R, SN |
DOUBLE PRECISION CS, F, G, R, SN |
Line 161
|
Line 158
|
F1 = F1*SAFMN2 |
F1 = F1*SAFMN2 |
G1 = G1*SAFMN2 |
G1 = G1*SAFMN2 |
SCALE = MAX( ABS( F1 ), ABS( G1 ) ) |
SCALE = MAX( ABS( F1 ), ABS( G1 ) ) |
IF( SCALE.GE.SAFMX2 ) |
IF( SCALE.GE.SAFMX2 .AND. COUNT .LT. 20 ) |
$ GO TO 10 |
$ GO TO 10 |
R = SQRT( F1**2+G1**2 ) |
R = SQRT( F1**2+G1**2 ) |
CS = F1 / R |
CS = F1 / R |