41 lines
1.3 KiB
Text
41 lines
1.3 KiB
Text
*** 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
|