blob: 7c7ca2d70d96544e62fa51afc55d080e77e4ac7b (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
*** BLAS/SRC/zrotg.f~ Mon Jul 3 19:23:21 1995
--- BLAS/SRC/zrotg.f Tue Jul 4 13:11:24 1995
***************
*** 3,20 ****
double precision c
double precision norm,scale
double complex alpha
! if (cdabs(ca) .ne. 0.0d0) go to 10
c = 0.0d0
s = (1.0d0,0.0d0)
ca = cb
go to 20
10 continue
! scale = cdabs(ca) + cdabs(cb)
! c = (cdabs(ca/dcmplx(scale,0.0d0)))**2 +
! * (cdabs(cb/dcmplx(scale,0.0d0)))**2
norm = scale*dsqrt(c)
! alpha = ca /cdabs(ca)
! c = cdabs(ca) / norm
s = alpha * dconjg(cb) / norm
ca = alpha * norm
20 continue
--- 3,20 ----
double precision c
double precision norm,scale
double complex alpha
! if (abs(ca) .ne. 0.0d0) go to 10
c = 0.0d0
s = (1.0d0,0.0d0)
ca = cb
go to 20
10 continue
! scale = abs(ca) + abs(cb)
! c = (abs(ca/dcmplx(scale,0.0d0)))**2 +
! * (abs(cb/dcmplx(scale,0.0d0)))**2
norm = scale*dsqrt(c)
! alpha = ca /abs(ca)
! c = abs(ca) / norm
s = alpha * dconjg(cb) / norm
ca = alpha * norm
20 continue
|
