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

version 1.13, 2017/06/17 10:53:45 version 1.16, 2023/08/07 08:38:44
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:
 *  ========  *  ========
 *  *
Line 35 Line 61
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date December 2016  
 *  
 *> \ingroup complex16_blas_level1  *> \ingroup complex16_blas_level1
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 52 Line 76
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZDSCAL(N,DA,ZX,INCX)        SUBROUTINE ZDSCAL(N,DA,ZX,INCX)
 *  *
 *  -- Reference BLAS level1 routine (version 3.7.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..--
 *     December 2016  
 *  *
 *     .. 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.13  
changed lines
  Added in v.1.16


CVSweb interface <joel.bertrand@systella.fr>