2 REAL FUNCTION PYTHAG (A, B)
3 C***BEGIN PROLOGUE PYTHAG
5 C***PURPOSE Compute the complex square root of a complex number without
6 C destructive overflow or underflow.
8 C***TYPE SINGLE PRECISION (PYTHAG-S)
12 C Finds sqrt(A**2+B**2) without overflow or destructive underflow
15 C***ROUTINES CALLED (NONE)
16 C***REVISION HISTORY (YYMMDD)
18 C 890531 Changed all specific intrinsics to generic. (WRB)
19 C 891214 Prologue converted to Version 4.0 format. (BAB)
20 C 900402 Added TYPE section. (WRB)
21 C***END PROLOGUE PYTHAG
25 C***FIRST EXECUTABLE STATEMENT PYTHAG
26 P = MAX(ABS(A),ABS(B))
27 Q = MIN(ABS(A),ABS(B))
28 IF (Q .EQ. 0.0E0) GO TO 20
32 IF (T .EQ. 4.0E0) GO TO 20