--- rpl/lapack/blas/zdrot.f	2010/08/13 21:03:41	1.5
+++ rpl/lapack/blas/zdrot.f	2017/06/17 10:53:45	1.13
@@ -1,5 +1,108 @@
+*> \brief \b ZDROT
+*
+*  =========== DOCUMENTATION ===========
+*
+* Online html documentation available at
+*            http://www.netlib.org/lapack/explore-html/
+*
+*  Definition:
+*  ===========
+*
+*       SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
+*
+*       .. Scalar Arguments ..
+*       INTEGER            INCX, INCY, N
+*       DOUBLE PRECISION   C, S
+*       ..
+*       .. Array Arguments ..
+*       COMPLEX*16         CX( * ), CY( * )
+*       ..
+*
+*
+*> \par Purpose:
+*  =============
+*>
+*> \verbatim
+*>
+*> Applies a plane rotation, where the cos and sin (c and s) are real
+*> and the vectors cx and cy are complex.
+*> jack dongarra, linpack, 3/11/78.
+*> \endverbatim
+*
+*  Arguments:
+*  ==========
+*
+*> \param[in] N
+*> \verbatim
+*>          N is INTEGER
+*>           On entry, N specifies the order of the vectors cx and cy.
+*>           N must be at least zero.
+*> \endverbatim
+*>
+*> \param[in,out] CX
+*> \verbatim
+*>          CX 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
+*>           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.
+*> \endverbatim
+*>
+*> \param[in,out] CY
+*> \verbatim
+*>          CY 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
+*>           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.
+*> \endverbatim
+*>
+*> \param[in] C
+*> \verbatim
+*>          C is DOUBLE PRECISION
+*>           On entry, C specifies the cosine, cos.
+*> \endverbatim
+*>
+*> \param[in] S
+*> \verbatim
+*>          S is DOUBLE PRECISION
+*>           On entry, S specifies the sine, sin.
+*> \endverbatim
+*
+*  Authors:
+*  ========
+*
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
+*
+*> \date December 2016
+*
+*> \ingroup complex16_blas_level1
+*
+*  =====================================================================
       SUBROUTINE ZDROT( N, CX, INCX, CY, INCY, C, S )
 *
+*  -- Reference BLAS level1 routine (version 3.7.0) --
+*  -- Reference BLAS is a software package provided by Univ. of Tennessee,    --
+*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
+*     December 2016
+*
 *     .. Scalar Arguments ..
       INTEGER            INCX, INCY, N
       DOUBLE PRECISION   C, S
@@ -8,51 +111,6 @@
       COMPLEX*16         CX( * ), CY( * )
 *     ..
 *
-*  Purpose
-*  =======
-*
-*  Applies a plane rotation, where the cos and sin (c and s) are real
-*  and the vectors cx and cy are complex.
-*  jack dongarra, linpack, 3/11/78.
-*
-*  Arguments
-*  ==========
-*
-*  N        (input) INTEGER
-*           On entry, N specifies the order of the vectors cx and cy.
-*           N must be at least zero.
-*           Unchanged on exit.
-*
-*  CX       (input) 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
-*           vector cx.
-*
-*  INCX     (input) INTEGER
-*           On entry, INCX specifies the increment for the elements of
-*           CX. INCX must not be zero.
-*           Unchanged on exit.
-*
-*  CY       (input) 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
-*           vector cy.
-*
-*  INCY     (input) INTEGER
-*           On entry, INCY specifies the increment for the elements of
-*           CY. INCY must not be zero.
-*           Unchanged on exit.
-*
-*  C        (input) DOUBLE PRECISION
-*           On entry, C specifies the cosine, cos.
-*           Unchanged on exit.
-*
-*  S        (input) DOUBLE PRECISION
-*           On entry, S specifies the sine, sin.
-*           Unchanged on exit.
-*
 * =====================================================================
 *
 *     .. Local Scalars ..
@@ -63,34 +121,33 @@
 *
       IF( N.LE.0 )
      $   RETURN
-      IF( INCX.EQ.1 .AND. INCY.EQ.1 )
-     $   GO TO 20
+      IF( INCX.EQ.1 .AND. INCY.EQ.1 ) THEN
 *
-*        code for unequal increments or equal increments not equal
-*          to 1
+*        code for both increments equal to 1
 *
-      IX = 1
-      IY = 1
-      IF( INCX.LT.0 )
-     $   IX = ( -N+1 )*INCX + 1
-      IF( INCY.LT.0 )
-     $   IY = ( -N+1 )*INCY + 1
-      DO 10 I = 1, N
-         CTEMP = C*CX( IX ) + S*CY( IY )
-         CY( IY ) = C*CY( IY ) - S*CX( IX )
-         CX( IX ) = CTEMP
-         IX = IX + INCX
-         IY = IY + INCY
-   10 CONTINUE
-      RETURN
+         DO I = 1, N
+            CTEMP = C*CX( I ) + S*CY( I )
+            CY( I ) = C*CY( I ) - S*CX( I )
+            CX( I ) = CTEMP
+         END DO
+      ELSE
 *
-*        code for both increments equal to 1
+*        code for unequal increments or equal increments not equal
+*          to 1
 *
-   20 CONTINUE
-      DO 30 I = 1, N
-         CTEMP = C*CX( I ) + S*CY( I )
-         CY( I ) = C*CY( I ) - S*CX( I )
-         CX( I ) = CTEMP
-   30 CONTINUE
+         IX = 1
+         IY = 1
+         IF( INCX.LT.0 )
+     $      IX = ( -N+1 )*INCX + 1
+         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
+            IX = IX + INCX
+            IY = IY + INCY
+         END DO
+      END IF
       RETURN
       END