blob: f6a4aa127deccbb06b04f95f1672e2b0c1f9259d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
subroutine zrotg(ca,cb,c,s)
double complex ca,cb,s
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)
norm = scale*dsqrt((cdabs(ca/dcmplx(scale,0.0d0)))**2 +
* (cdabs(cb/dcmplx(scale,0.0d0)))**2)
alpha = ca /cdabs(ca)
c = cdabs(ca) / norm
s = alpha * dconjg(cb) / norm
ca = alpha * norm
20 continue
return
end
|
