56 DOUBLE PRECISION norm,scale
59 INTRINSIC cdabs,dcmplx,dconjg,dsqrt
61 IF (cdabs(ca).EQ.0.0d0)
THEN
66 scale = cdabs(ca) + cdabs(cb)
67 norm = scale*dsqrt((cdabs(ca/dcmplx(scale,0.0d0)))**2+
68 $ (cdabs(cb/dcmplx(scale,0.0d0)))**2)
71 s = alpha*dconjg(cb)/norm