84 INTRINSIC abs,sign,sqrt
86 scale = abs(sa) + abs(sb)
87 IF (scale.EQ.0.0)
THEN
94 IF (abs(sa).GT.abs(sb)) roe = sa
95 r = scale*sqrt((sa/scale)**2+ (sb/scale)**2)
100 IF (abs(sa).GT.abs(sb)) z = s
101 IF (abs(sb).GE.abs(sa) .AND. c.NE.0.0) z = 1.0/c