|
5 | 5 | ! Online html documentation available at |
6 | 6 | ! http://www.netlib.org/lapack/explore-html/ |
7 | 7 | ! |
8 | | -! Definition: |
9 | | -! =========== |
10 | | -! |
11 | | -! CROTG constructs a plane rotation |
12 | | -! [ c s ] [ a ] = [ r ] |
13 | | -! [ -conjg(s) c ] [ b ] [ 0 ] |
14 | | -! where c is real, s is complex, and c**2 + conjg(s)*s = 1. |
15 | | -! |
16 | 8 | !> \par Purpose: |
17 | 9 | ! ============= |
18 | 10 | !> |
19 | 11 | !> \verbatim |
20 | 12 | !> |
| 13 | +!> CROTG constructs a plane rotation |
| 14 | +!> [ c s ] [ a ] = [ r ] |
| 15 | +!> [ -conjg(s) c ] [ b ] [ 0 ] |
| 16 | +!> where c is real, s is complex, and c**2 + conjg(s)*s = 1. |
| 17 | +!> |
21 | 18 | !> The computation uses the formulas |
22 | 19 | !> |x| = sqrt( Re(x)**2 + Im(x)**2 ) |
23 | 20 | !> sgn(x) = x / |x| if x /= 0 |
|
33 | 30 | !> if the signs of a and b are not the same. |
34 | 31 | !> |
35 | 32 | !> \endverbatim |
| 33 | +!> |
| 34 | +!> @see lartg, @see lartgp |
36 | 35 | ! |
37 | 36 | ! Arguments: |
38 | 37 | ! ========== |
@@ -181,7 +180,7 @@ subroutine CROTG( a, b, c, s ) |
181 | 180 | f2 = ABSSQ( f ) |
182 | 181 | g2 = ABSSQ( g ) |
183 | 182 | h2 = f2 + g2 |
184 | | - ! safmin <= f2 <= h2 <= safmax |
| 183 | + ! safmin <= f2 <= h2 <= safmax |
185 | 184 | if( f2 >= h2 * safmin ) then |
186 | 185 | ! safmin <= f2/h2 <= 1, and h2/f2 is finite |
187 | 186 | c = sqrt( f2 / h2 ) |
@@ -237,7 +236,7 @@ subroutine CROTG( a, b, c, s ) |
237 | 236 | f2 = ABSSQ( fs ) |
238 | 237 | h2 = f2 + g2 |
239 | 238 | end if |
240 | | - ! safmin <= f2 <= h2 <= safmax |
| 239 | + ! safmin <= f2 <= h2 <= safmax |
241 | 240 | if( f2 >= h2 * safmin ) then |
242 | 241 | ! safmin <= f2/h2 <= 1, and h2/f2 is finite |
243 | 242 | c = sqrt( f2 / h2 ) |
|
0 commit comments