rpl/lapack/blas/dzasum.f - diff

Return to dzasum.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 22	Line 22
*	*
* .. Local Scalars ..	* .. Local Scalars ..
DOUBLE PRECISION STEMP	DOUBLE PRECISION STEMP
INTEGER I,IX	INTEGER I,NINCX
* ..	* ..
* .. External Functions ..	* .. External Functions ..
DOUBLE PRECISION DCABS1	DOUBLE PRECISION DCABS1
Line 31	Line 31
DZASUM = 0.0d0	DZASUM = 0.0d0
STEMP = 0.0d0	STEMP = 0.0d0
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	* code for increment equal to 1
*	*
IX = 1	DO I = 1,N
DO 10 I = 1,N	STEMP = STEMP + DCABS1(ZX(I))
STEMP = STEMP + DCABS1(ZX(IX))	END DO
IX = IX + INCX	ELSE
10 CONTINUE
DZASUM = STEMP
RETURN
*	*
* code for increment equal to 1	* code for increment not equal to 1
*	*
20 DO 30 I = 1,N	NINCX = N*INCX
STEMP = STEMP + DCABS1(ZX(I))	DO I = 1,NINCX,INCX
30 CONTINUE	STEMP = STEMP + DCABS1(ZX(I))
	END DO
	END IF
DZASUM = STEMP	DZASUM = STEMP
RETURN	RETURN
END	END