version 1.6, 2010/12/21 13:51:25
|
version 1.7, 2011/07/22 07:38:01
|
Line 29
|
Line 29
|
INTRINSIC MOD |
INTRINSIC MOD |
* .. |
* .. |
IF (N.LE.0 .OR. INCX.LE.0) RETURN |
IF (N.LE.0 .OR. INCX.LE.0) RETURN |
IF (INCX.EQ.1) GO TO 20 |
IF (INCX.EQ.1) THEN |
* |
|
* code for increment not equal to 1 |
|
* |
|
NINCX = N*INCX |
|
DO 10 I = 1,NINCX,INCX |
|
DX(I) = DA*DX(I) |
|
10 CONTINUE |
|
RETURN |
|
* |
* |
* code for increment equal to 1 |
* code for increment equal to 1 |
* |
* |
* |
* |
* clean-up loop |
* clean-up loop |
* |
* |
20 M = MOD(N,5) |
M = MOD(N,5) |
IF (M.EQ.0) GO TO 40 |
IF (M.NE.0) THEN |
DO 30 I = 1,M |
DO I = 1,M |
DX(I) = DA*DX(I) |
DX(I) = DA*DX(I) |
30 CONTINUE |
END DO |
IF (N.LT.5) RETURN |
IF (N.LT.5) RETURN |
40 MP1 = M + 1 |
END IF |
DO 50 I = MP1,N,5 |
MP1 = M + 1 |
DX(I) = DA*DX(I) |
DO I = MP1,N,5 |
DX(I+1) = DA*DX(I+1) |
DX(I) = DA*DX(I) |
DX(I+2) = DA*DX(I+2) |
DX(I+1) = DA*DX(I+1) |
DX(I+3) = DA*DX(I+3) |
DX(I+2) = DA*DX(I+2) |
DX(I+4) = DA*DX(I+4) |
DX(I+3) = DA*DX(I+3) |
50 CONTINUE |
DX(I+4) = DA*DX(I+4) |
|
END DO |
|
ELSE |
|
* |
|
* code for increment not equal to 1 |
|
* |
|
NINCX = N*INCX |
|
DO I = 1,NINCX,INCX |
|
DX(I) = DA*DX(I) |
|
END DO |
|
END IF |
RETURN |
RETURN |
END |
END |