rpl/lapack/blas/izamax.f - diff

Return to izamax.f CVS log

Up to [local] / rpl / lapack / blas

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