blob: d88b6f107fd73063d35cc24d1d3ad1c690689909 (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
*** BLAS/SRC/zrotg.f~ Fri Oct 30 21:09:39 1992
--- BLAS/SRC/zrotg.f Mon Jul 3 01:28:54 1995
***************
*** 10,17 ****
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
--- 10,18 ----
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
|
