version 1.1.1.1, 2010/01/26 15:22:45
|
version 1.7, 2011/07/22 07:38:02
|
Line 26
|
Line 26
|
ZTEMP = (0.0d0,0.0d0) |
ZTEMP = (0.0d0,0.0d0) |
ZDOTU = (0.0d0,0.0d0) |
ZDOTU = (0.0d0,0.0d0) |
IF (N.LE.0) RETURN |
IF (N.LE.0) RETURN |
IF (INCX.EQ.1 .AND. INCY.EQ.1) GO TO 20 |
IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN |
* |
* |
* code for unequal increments or equal increments |
* code for both increments equal to 1 |
* not equal to 1 |
|
* |
* |
IX = 1 |
DO I = 1,N |
IY = 1 |
ZTEMP = ZTEMP + ZX(I)*ZY(I) |
IF (INCX.LT.0) IX = (-N+1)*INCX + 1 |
END DO |
IF (INCY.LT.0) IY = (-N+1)*INCY + 1 |
ELSE |
DO 10 I = 1,N |
|
ZTEMP = ZTEMP + ZX(IX)*ZY(IY) |
|
IX = IX + INCX |
|
IY = IY + INCY |
|
10 CONTINUE |
|
ZDOTU = ZTEMP |
|
RETURN |
|
* |
* |
* code for both increments equal to 1 |
* code for unequal increments or equal increments |
|
* not equal to 1 |
* |
* |
20 DO 30 I = 1,N |
IX = 1 |
ZTEMP = ZTEMP + ZX(I)*ZY(I) |
IY = 1 |
30 CONTINUE |
IF (INCX.LT.0) IX = (-N+1)*INCX + 1 |
|
IF (INCY.LT.0) IY = (-N+1)*INCY + 1 |
|
DO I = 1,N |
|
ZTEMP = ZTEMP + ZX(IX)*ZY(IY) |
|
IX = IX + INCX |
|
IY = IY + INCY |
|
END DO |
|
END IF |
ZDOTU = ZTEMP |
ZDOTU = ZTEMP |
RETURN |
RETURN |
END |
END |