clartv()

subroutine clartv	(	integer	n,
		complex, dimension( * )	x,
		integer	incx,
		complex, dimension( * )	y,
		integer	incy,
		real, dimension( * )	c,
		complex, dimension( * )	s,
		integer	incc )

CLARTV applies a vector of plane rotations with real cosines and complex sines to the elements of a pair of vectors.

Download CLARTV + dependencies [TGZ] [ZIP] [TXT]

Purpose:

!>
!> CLARTV applies a vector of complex plane rotations with real cosines
!> to elements of the complex vectors x and y. For i = 1,2,...,n
!>
!>    ( x(i) ) := (        c(i)   s(i) ) ( x(i) )
!>    ( y(i) )    ( -conjg(s(i))  c(i) ) ( y(i) )
!> 

Parameters

[in]	N	!> N is INTEGER !> The number of plane rotations to be applied. !>
[in,out]	X	!> X is COMPLEX array, dimension (1+(N-1)*INCX) !> The vector x. !>
[in]	INCX	!> INCX is INTEGER !> The increment between elements of X. INCX > 0. !>
[in,out]	Y	!> Y is COMPLEX array, dimension (1+(N-1)*INCY) !> The vector y. !>
[in]	INCY	!> INCY is INTEGER !> The increment between elements of Y. INCY > 0. !>
[in]	C	!> C is REAL array, dimension (1+(N-1)*INCC) !> The cosines of the plane rotations. !>
[in]	S	!> S is COMPLEX array, dimension (1+(N-1)*INCC) !> The sines of the plane rotations. !>
[in]	INCC	!> INCC is INTEGER !> The increment between elements of C and S. INCC > 0. !>

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Definition at line 104 of file clartv.f.

*
*  -- LAPACK auxiliary routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      INTEGER            INCC, INCX, INCY, N
*     ..
*     .. Array Arguments ..
      REAL               C( * )
      COMPLEX            S( * ), X( * ), Y( * )
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      INTEGER            I, IC, IX, IY
      COMPLEX            XI, YI
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          conjg
*     ..
*     .. Executable Statements ..
*
      ix = 1
      iy = 1
      ic = 1
      DO 10 i = 1, n
         xi = x( ix )
         yi = y( iy )
         x( ix ) = c( ic )*xi + s( ic )*yi
         y( iy ) = c( ic )*yi - conjg( s( ic ) )*xi
         ix = ix + incx
         iy = iy + incy
         ic = ic + incc
   10 CONTINUE
      RETURN
*
*     End of CLARTV
*

Here is the caller graph for this function: