version 1.6, 2010/12/21 13:51:26
|
version 1.7, 2011/07/22 07:38:02
|
Line 31
|
Line 31
|
IF (N.LT.1 .OR. INCX.LE.0) RETURN |
IF (N.LT.1 .OR. INCX.LE.0) RETURN |
IDAMAX = 1 |
IDAMAX = 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 = DABS(DX(1)) |
DMAX = DABS(DX(1)) |
DO I = 2,N |
IX = IX + INCX |
IF (DABS(DX(I)).GT.DMAX) THEN |
DO 10 I = 2,N |
IDAMAX = I |
IF (DABS(DX(IX)).LE.DMAX) GO TO 5 |
DMAX = DABS(DX(I)) |
IDAMAX = I |
END IF |
DMAX = DABS(DX(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 DMAX = DABS(DX(1)) |
IX = 1 |
DO 30 I = 2,N |
DMAX = DABS(DX(1)) |
IF (DABS(DX(I)).LE.DMAX) GO TO 30 |
IX = IX + INCX |
IDAMAX = I |
DO I = 2,N |
DMAX = DABS(DX(I)) |
IF (DABS(DX(IX)).GT.DMAX) THEN |
30 CONTINUE |
IDAMAX = I |
|
DMAX = DABS(DX(IX)) |
|
END IF |
|
IX = IX + INCX |
|
END DO |
|
END IF |
RETURN |
RETURN |
END |
END |