home *** CD-ROM | disk | FTP | other *** search
- INTEGER FUNCTION IWAMAX (N, XR, XI, INCX)
- IMPLICIT NONE
- C
- C INDEX OF NORMINF(X)
- C
- INTEGER N, INCX
- DOUBLE PRECISION XR(*), XI(*)
- C
- INTEGER I, K, IX
- DOUBLE PRECISION S, P
- C
- C
- K = 0
- IF (N.LE.0) GO TO 20
- K = 1
- S = 0.0D0
- IX = 1
- DO 10 I = 1, N
- P = DABS (XR(IX))+DABS (XI(IX))
- IF (P.GT.S) K = I
- IF (P.GT.S) S = P
- IX = IX+INCX
- 10 CONTINUE
- 20 CONTINUE
- IWAMAX = K
- C
- RETURN
- END
-