/* clar2v.f -- translated by f2c (version 20061008).
You must link the resulting object file with libf2c:
on Microsoft Windows system, link with libf2c.lib;
on Linux or Unix systems, link with .../path/to/libf2c.a -lm
or, if you install libf2c.a in a standard place, with -lf2c -lm
-- in that order, at the end of the command line, as in
cc *.o -lf2c -lm
Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
http://www.netlib.org/f2c/libf2c.zip
*/
#include "f2c.h"
#include "blaswrap.h"
/* Subroutine */ int clar2v_(integer *n, complex *x, complex *y, complex *z__,
integer *incx, real *c__, complex *s, integer *incc)
{
/* System generated locals */
integer i__1, i__2;
real r__1;
complex q__1, q__2, q__3, q__4, q__5;
/* Builtin functions */
double r_imag(complex *);
void r_cnjg(complex *, complex *);
/* Local variables */
integer i__;
complex t2, t3, t4;
real t5, t6;
integer ic;
real ci;
complex si;
integer ix;
real xi, yi;
complex zi;
real t1i, t1r, sii, zii, sir, zir;
/* -- LAPACK auxiliary routine (version 3.2) -- */
/* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
/* November 2006 */
/* .. Scalar Arguments .. */
/* .. */
/* .. Array Arguments .. */
/* .. */
/* Purpose */
/* ======= */
/* CLAR2V applies a vector of complex plane rotations with real cosines */
/* from both sides to a sequence of 2-by-2 complex Hermitian matrices, */
/* defined by the elements of the vectors x, y and z. For i = 1,2,...,n */
/* ( x(i) z(i) ) := */
/* ( conjg(z(i)) y(i) ) */
/* ( c(i) conjg(s(i)) ) ( x(i) z(i) ) ( c(i) -conjg(s(i)) ) */
/* ( -s(i) c(i) ) ( conjg(z(i)) y(i) ) ( s(i) c(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; the elements of x are assumed to be real. */
/* Y (input/output) COMPLEX array, dimension (1+(N-1)*INCX) */
/* The vector y; the elements of y are assumed to be real. */
/* Z (input/output) COMPLEX array, dimension (1+(N-1)*INCX) */
/* The vector z. */
/* INCX (input) INTEGER */
/* The increment between elements of X, Y and Z. INCX > 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 .. */
/* .. */
/* .. Intrinsic Functions .. */
/* .. */
/* .. Executable Statements .. */
/* Parameter adjustments */
--s;
--c__;
--z__;
--y;
--x;
/* Function Body */
ix = 1;
ic = 1;
i__1 = *n;
for (i__ = 1; i__ <= i__1; ++i__) {
i__2 = ix;
xi = x[i__2].r;
i__2 = ix;
yi = y[i__2].r;
i__2 = ix;
zi.r = z__[i__2].r, zi.i = z__[i__2].i;
zir = zi.r;
zii = r_imag(&zi);
ci = c__[ic];
i__2 = ic;
si.r = s[i__2].r, si.i = s[i__2].i;
sir = si.r;
sii = r_imag(&si);
t1r = sir * zir - sii * zii;
t1i = sir * zii + sii * zir;
q__1.r = ci * zi.r, q__1.i = ci * zi.i;
t2.r = q__1.r, t2.i = q__1.i;
r_cnjg(&q__3, &si);
q__2.r = xi * q__3.r, q__2.i = xi * q__3.i;
q__1.r = t2.r - q__2.r, q__1.i = t2.i - q__2.i;
t3.r = q__1.r, t3.i = q__1.i;
r_cnjg(&q__2, &t2);
q__3.r = yi * si.r, q__3.i = yi * si.i;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
t4.r = q__1.r, t4.i = q__1.i;
t5 = ci * xi + t1r;
t6 = ci * yi - t1r;
i__2 = ix;
r__1 = ci * t5 + (sir * t4.r + sii * r_imag(&t4));
x[i__2].r = r__1, x[i__2].i = 0.f;
i__2 = ix;
r__1 = ci * t6 - (sir * t3.r - sii * r_imag(&t3));
y[i__2].r = r__1, y[i__2].i = 0.f;
i__2 = ix;
q__2.r = ci * t3.r, q__2.i = ci * t3.i;
r_cnjg(&q__4, &si);
q__5.r = t6, q__5.i = t1i;
q__3.r = q__4.r * q__5.r - q__4.i * q__5.i, q__3.i = q__4.r * q__5.i
+ q__4.i * q__5.r;
q__1.r = q__2.r + q__3.r, q__1.i = q__2.i + q__3.i;
z__[i__2].r = q__1.r, z__[i__2].i = q__1.i;
ix += *incx;
ic += *incc;
/* L10: */
}
return 0;
/* End of CLAR2V */
} /* clar2v_ */
