[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.11 and 1.16

version 1.11, 2014/01/27 09:28:13 version 1.16, 2023/08/07 08:38:44
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE ZDSCAL(N,DA,ZX,INCX)  *       SUBROUTINE ZDSCAL(N,DA,ZX,INCX)
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       DOUBLE PRECISION DA  *       DOUBLE PRECISION DA
 *       INTEGER INCX,N  *       INTEGER INCX,N
Line 17 Line 17
 *       .. Array Arguments ..  *       .. Array Arguments ..
 *       COMPLEX*16 ZX(*)  *       COMPLEX*16 ZX(*)
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 27 Line 27
 *>    ZDSCAL scales a vector by a constant.  *>    ZDSCAL scales a vector by a constant.
 *> \endverbatim  *> \endverbatim
 *  *
   *  Arguments:
   *  ==========
   *
   *> \param[in] N
   *> \verbatim
   *>          N is INTEGER
   *>         number of elements in input vector(s)
   *> \endverbatim
   *>
   *> \param[in] DA
   *> \verbatim
   *>          DA is DOUBLE PRECISION
   *>           On entry, DA specifies the scalar alpha.
   *> \endverbatim
   *>
   *> \param[in,out] ZX
   *> \verbatim
   *>          ZX is COMPLEX*16 array, dimension ( 1 + ( N - 1 )*abs( INCX ) )
   *> \endverbatim
   *>
   *> \param[in] INCX
   *> \verbatim
   *>          INCX is INTEGER
   *>         storage spacing between elements of ZX
   *> \endverbatim
   *
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  
 *> \date November 2011  
 *  *
 *> \ingroup complex16_blas_level1  *> \ingroup complex16_blas_level1
 *  *
Line 52 Line 76
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZDSCAL(N,DA,ZX,INCX)        SUBROUTINE ZDSCAL(N,DA,ZX,INCX)
 *  *
 *  -- Reference BLAS level1 routine (version 3.4.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 2011  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       DOUBLE PRECISION DA        DOUBLE PRECISION DA
Line 69 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 87 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.11  
changed lines
  Added in v.1.16

CVSweb interface <joel.bertrand@systella.fr>