|
|
1.1 ! root 1: SUBROUTINE WSWAP(N,XR,XI,INCX,YR,YI,INCY) ! 2: DOUBLE PRECISION XR(1),XI(1),YR(1),YI(1),T ! 3: IF (N .LE. 0) RETURN ! 4: IX = 1 ! 5: IY = 1 ! 6: IF (INCX.LT.0) IX = (-N+1)*INCX + 1 ! 7: IF (INCY.LT.0) IY = (-N+1)*INCY + 1 ! 8: DO 10 I = 1, N ! 9: T = XR(IX) ! 10: XR(IX) = YR(IY) ! 11: YR(IY) = T ! 12: T = XI(IX) ! 13: XI(IX) = YI(IY) ! 14: YI(IY) = T ! 15: IX = IX + INCX ! 16: IY = IY + INCY ! 17: 10 CONTINUE ! 18: RETURN ! 19: END
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.