rpl/lapack/blas/zdscal.f - diff

Return to zdscal.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 ..
INTEGER I,IX	INTEGER I,NINCX
* ..	* ..
* .. Intrinsic Functions ..	* .. Intrinsic Functions ..
INTRINSIC DCMPLX	INTRINSIC DCMPLX
* ..	* ..
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	ZX(I) = DCMPLX(DA,0.0d0)*ZX(I)
ZX(IX) = DCMPLX(DA,0.0d0)*ZX(IX)	END DO
IX = IX + INCX	ELSE
10 CONTINUE
RETURN
*	*
* code for increment equal to 1	* code for increment not equal to 1
*	*
20 DO 30 I = 1,N	NINCX = N*INCX
ZX(I) = DCMPLX(DA,0.0d0)*ZX(I)	DO I = 1,NINCX,INCX
30 CONTINUE	ZX(I) = DCMPLX(DA,0.0d0)*ZX(I)
	END DO
	END IF
RETURN	RETURN
END	END