diff options
| author | jason <jason@8a072113-8704-0410-8d35-dd094bca7971> | 2008-10-28 01:38:50 +0000 |
|---|---|---|
| committer | jason <jason@8a072113-8704-0410-8d35-dd094bca7971> | 2008-10-28 01:38:50 +0000 |
| commit | baba851215b44ac3b60b9248eb02bcce7eb76247 (patch) | |
| tree | 8c0f5c006875532a30d4409f5e94b0f310ff00a7 /SRC/clartv.f | |
Move LAPACK trunk into position.
Diffstat (limited to 'SRC/clartv.f')
| -rw-r--r-- | SRC/clartv.f | 78 |
1 files changed, 78 insertions, 0 deletions
diff --git a/SRC/clartv.f b/SRC/clartv.f new file mode 100644 index 00000000..553a2ded --- /dev/null +++ b/SRC/clartv.f @@ -0,0 +1,78 @@ + SUBROUTINE CLARTV( N, X, INCX, Y, INCY, C, S, INCC ) +* +* -- LAPACK auxiliary routine (version 3.1) -- +* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. +* November 2006 +* +* .. Scalar Arguments .. + INTEGER INCC, INCX, INCY, N +* .. +* .. Array Arguments .. + REAL C( * ) + COMPLEX S( * ), X( * ), Y( * ) +* .. +* +* 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) ) +* +* Arguments +* ========= +* +* N (input) INTEGER +* The number of plane rotations to be applied. +* +* X (input/output) COMPLEX array, dimension (1+(N-1)*INCX) +* The vector x. +* +* INCX (input) INTEGER +* The increment between elements of X. INCX > 0. +* +* Y (input/output) COMPLEX array, dimension (1+(N-1)*INCY) +* The vector y. +* +* INCY (input) INTEGER +* The increment between elements of Y. INCY > 0. +* +* C (input) REAL array, dimension (1+(N-1)*INCC) +* The cosines of the plane rotations. +* +* S (input) COMPLEX array, dimension (1+(N-1)*INCC) +* The sines of the plane rotations. +* +* INCC (input) INTEGER +* The increment between elements of C and S. INCC > 0. +* +* ===================================================================== +* +* .. 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 +* + END |