[BACK]Return to zdscal.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / blas

Diff for /rpl/lapack/blas/zdscal.f between versions 1.15 and 1.16

version 1.15, 2018/05/29 07:19:42 version 1.16, 2023/08/07 08:38:44
Line 61 Line 61
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date November 2017  
 *  
 *> \ingroup complex16_blas_level1  *> \ingroup complex16_blas_level1
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 78 Line 76
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZDSCAL(N,DA,ZX,INCX)        SUBROUTINE ZDSCAL(N,DA,ZX,INCX)
 *  *
 *  -- Reference BLAS level1 routine (version 3.8.0) --  *  -- Reference BLAS level1 routine --
 *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --  *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     November 2017  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       DOUBLE PRECISION DA        DOUBLE PRECISION DA
Line 95 Line 92
 *  *
 *     .. Local Scalars ..  *     .. Local Scalars ..
       INTEGER I,NINCX        INTEGER I,NINCX
   *     .. Parameters ..
         DOUBLE PRECISION ONE
         PARAMETER (ONE=1.0D+0)
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC DCMPLX        INTRINSIC DBLE, DCMPLX, DIMAG
 *     ..  *     ..
       IF (N.LE.0 .OR. INCX.LE.0) RETURN        IF (N.LE.0 .OR. INCX.LE.0 .OR. DA.EQ.ONE) RETURN
       IF (INCX.EQ.1) THEN        IF (INCX.EQ.1) THEN
 *  *
 *        code for increment equal to 1  *        code for increment equal to 1
 *  *
          DO I = 1,N           DO I = 1,N
             ZX(I) = DCMPLX(DA,0.0d0)*ZX(I)              ZX(I) = DCMPLX(DA*DBLE(ZX(I)),DA*DIMAG(ZX(I)))
          END DO           END DO
       ELSE        ELSE
 *  *
Line 113 Line 113
 *  *
          NINCX = N*INCX           NINCX = N*INCX
          DO I = 1,NINCX,INCX           DO I = 1,NINCX,INCX
             ZX(I) = DCMPLX(DA,0.0d0)*ZX(I)              ZX(I) = DCMPLX(DA*DBLE(ZX(I)),DA*DIMAG(ZX(I)))
          END DO           END DO
       END IF        END IF
       RETURN        RETURN
   *
   *     End of ZDSCAL
   *
       END        END
Removed from v.1.15  
changed lines
  Added in v.1.16

CVSweb interface <joel.bertrand@systella.fr>