Merge pull request #452 from timleslie/use-zx-in-zdrot

Use ZX, ZY rather than CX, CY in ZDROT

Author: langou

BLAS/SRC/zdrot.f

@@ -8,14 +8,14 @@
 *  Definition:
 *  ===========
 *
-*       SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
+*       SUBROUTINE ZDROT( N, ZX, INCX, ZY, INCY, C, S )
 *
 *       .. Scalar Arguments ..
 *       INTEGER            INCX, INCY, N
 *       DOUBLE PRECISION   C, S
 *       ..
 *       .. Array Arguments ..
-*       COMPLEX*16         CX( * ), CY( * )
+*       COMPLEX*16         ZX( * ), ZY( * )
 *       ..
 *
 *
@@ -39,36 +39,36 @@
 *>           N must be at least zero.
 *> \endverbatim
 *>
-*> \param[in,out] CX
+*> \param[in,out] ZX
 *> \verbatim
-*>          CX is COMPLEX*16 array, dimension at least
+*>          ZX is COMPLEX*16 array, dimension at least
 *>           ( 1 + ( N - 1 )*abs( INCX ) ).
-*>           Before entry, the incremented array CX must contain the n
-*>           element vector cx. On exit, CX is overwritten by the updated
+*>           Before entry, the incremented array ZX must contain the n
+*>           element vector cx. On exit, ZX is overwritten by the updated
 *>           vector cx.
 *> \endverbatim
 *>
 *> \param[in] INCX
 *> \verbatim
 *>          INCX is INTEGER
 *>           On entry, INCX specifies the increment for the elements of
-*>           CX. INCX must not be zero.
+*>           ZX. INCX must not be zero.
 *> \endverbatim
 *>
-*> \param[in,out] CY
+*> \param[in,out] ZY
 *> \verbatim
-*>          CY is COMPLEX*16 array, dimension at least
+*>          ZY is COMPLEX*16 array, dimension at least
 *>           ( 1 + ( N - 1 )*abs( INCY ) ).
-*>           Before entry, the incremented array CY must contain the n
-*>           element vector cy. On exit, CY is overwritten by the updated
+*>           Before entry, the incremented array ZY must contain the n
+*>           element vector cy. On exit, ZY is overwritten by the updated
 *>           vector cy.
 *> \endverbatim
 *>
 *> \param[in] INCY
 *> \verbatim
 *>          INCY is INTEGER
 *>           On entry, INCY specifies the increment for the elements of
-*>           CY. INCY must not be zero.
+*>           ZY. INCY must not be zero.
 *> \endverbatim
 *>
 *> \param[in] C
@@ -96,7 +96,7 @@
 *> \ingroup complex16_blas_level1
 *
 *  =====================================================================
-      SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
+      SUBROUTINE ZDROT( N, ZX, INCX, ZY, INCY, C, S )
 *
 *  -- Reference BLAS level1 routine (version 3.7.0) --
 *  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
@@ -108,7 +108,7 @@ SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
       DOUBLE PRECISION   C, S
 *     ..
 *     .. Array Arguments ..
-      COMPLEX*16         CX( * ), CY( * )
+      COMPLEX*16         ZX( * ), ZY( * )
 *     ..
 *
 * =====================================================================
@@ -126,9 +126,9 @@ SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
 *        code for both increments equal to 1
 *
          DO I = 1, N
-            CTEMP = C*CX( I ) + S*CY( I )
-            CY( I ) = C*CY( I ) - S*CX( I )
-            CX( I ) = CTEMP
+            CTEMP = C*ZX( I ) + S*ZY( I )
+            ZY( I ) = C*ZY( I ) - S*ZX( I )
+            ZX( I ) = CTEMP
          END DO
       ELSE
 *
@@ -142,9 +142,9 @@ SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
          IF( INCY.LT.0 )
      $      IY = ( -N+1 )*INCY + 1
          DO I = 1, N
-            CTEMP = C*CX( IX ) + S*CY( IY )
-            CY( IY ) = C*CY( IY ) - S*CX( IX )
-            CX( IX ) = CTEMP
+            CTEMP = C*ZX( IX ) + S*ZY( IY )
+            ZY( IY ) = C*ZY( IY ) - S*ZX( IX )
+            ZX( IX ) = CTEMP
             IX = IX + INCX
             IY = IY + INCY
          END DO