version 1.6, 2010/12/21 13:51:26
|
version 1.7, 2011/07/22 07:38:02
|
Line 21
|
Line 21
|
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. Local Scalars .. |
* .. Local Scalars .. |
DOUBLE PRECISION SMAX |
DOUBLE PRECISION DMAX |
INTEGER I,IX |
INTEGER I,IX |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
Line 32
|
Line 32
|
IF (N.LT.1 .OR. INCX.LE.0) RETURN |
IF (N.LT.1 .OR. INCX.LE.0) RETURN |
IZAMAX = 1 |
IZAMAX = 1 |
IF (N.EQ.1) RETURN |
IF (N.EQ.1) RETURN |
IF (INCX.EQ.1) GO TO 20 |
IF (INCX.EQ.1) THEN |
* |
* |
* code for increment not equal to 1 |
* code for increment equal to 1 |
* |
* |
IX = 1 |
DMAX = DCABS1(ZX(1)) |
SMAX = DCABS1(ZX(1)) |
DO I = 2,N |
IX = IX + INCX |
IF (DCABS1(ZX(I)).GT.DMAX) THEN |
DO 10 I = 2,N |
IZAMAX = I |
IF (DCABS1(ZX(IX)).LE.SMAX) GO TO 5 |
DMAX = DCABS1(ZX(I)) |
IZAMAX = I |
END IF |
SMAX = DCABS1(ZX(IX)) |
END DO |
5 IX = IX + INCX |
ELSE |
10 CONTINUE |
|
RETURN |
|
* |
* |
* code for increment equal to 1 |
* code for increment not equal to 1 |
* |
* |
20 SMAX = DCABS1(ZX(1)) |
IX = 1 |
DO 30 I = 2,N |
DMAX = DCABS1(ZX(1)) |
IF (DCABS1(ZX(I)).LE.SMAX) GO TO 30 |
IX = IX + INCX |
IZAMAX = I |
DO I = 2,N |
SMAX = DCABS1(ZX(I)) |
IF (DCABS1(ZX(IX)).GT.DMAX) THEN |
30 CONTINUE |
IZAMAX = I |
|
DMAX = DCABS1(ZX(IX)) |
|
END IF |
|
IX = IX + INCX |
|
END DO |
|
END IF |
RETURN |
RETURN |
END |
END |