23 lines
790 B
Plaintext
23 lines
790 B
Plaintext
*** 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
|