|
|
1.1 ! root 1: INTEGER FUNCTION IWAMAX(N,XR,XI,INCX) ! 2: DOUBLE PRECISION XR(1),XI(1),S,P ! 3: C INDEX OF NORMINF(X) ! 4: K = 0 ! 5: IF (N .LE. 0) GO TO 20 ! 6: K = 1 ! 7: S = 0.0D0 ! 8: IX = 1 ! 9: DO 10 I = 1, N ! 10: P = DABS(XR(IX)) + DABS(XI(IX)) ! 11: IF (P .GT. S) K = I ! 12: IF (P .GT. S) S = P ! 13: IX = IX + INCX ! 14: 10 CONTINUE ! 15: 20 IWAMAX = K ! 16: RETURN ! 17: END
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.