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