NAME
ZROTG - construct givens plane rotation
SYNOPSIS
subroutine zrotg(ca,cb,c,s)
double
precision norm,scale
if
(cdabs(ca) .ne. 0.0d0) go to 10
scale
= cdabs(ca) + cdabs(cb)
norm
= scale*dsqrt((cdabs(ca/dcmplx(scale,0.0d0)))**2 +
*
(cdabs(cb/dcmplx(scale,0.0d0)))**2)
s
= alpha * dconjg(cb) / norm
PURPOSE
ČLÁNKY DO MAILU