14 typedef long long BLASLONG;
15 typedef unsigned long long BLASULONG;
17 typedef long BLASLONG;
18 typedef unsigned long BLASULONG;
22 typedef BLASLONG blasint;
24 #define blasabs(x) llabs(x)
26 #define blasabs(x) labs(x)
30 #define blasabs(x) abs(x)
33 typedef blasint integer;
35 typedef unsigned int uinteger;
36 typedef char *address;
37 typedef short int shortint;
39 typedef double doublereal;
40 typedef struct { real r, i; } complex;
41 typedef struct { doublereal r, i; } doublecomplex;
43 static inline _Fcomplex Cf(complex *z) {_Fcomplex zz={z->r , z->i}; return zz;}
44 static inline _Dcomplex Cd(doublecomplex *z) {_Dcomplex zz={z->r , z->i};return zz;}
45 static inline _Fcomplex * _pCf(complex *z) {return (_Fcomplex*)z;}
46 static inline _Dcomplex * _pCd(doublecomplex *z) {return (_Dcomplex*)z;}
48 static inline _Complex float Cf(complex *z) {return z->r + z->i*_Complex_I;}
49 static inline _Complex double Cd(doublecomplex *z) {return z->r + z->i*_Complex_I;}
50 static inline _Complex float * _pCf(complex *z) {return (_Complex float*)z;}
51 static inline _Complex double * _pCd(doublecomplex *z) {return (_Complex double*)z;}
53 #define pCf(z) (*_pCf(z))
54 #define pCd(z) (*_pCd(z))
56 typedef short int shortlogical;
57 typedef char logical1;
58 typedef char integer1;
63 /* Extern is for use with -E */
74 /*external read, write*/
83 /*internal read, write*/
113 /*rewind, backspace, endfile*/
125 ftnint *inex; /*parameters in standard's order*/
151 union Multitype { /* for multiple entry points */
162 typedef union Multitype Multitype;
164 struct Vardesc { /* for Namelist */
170 typedef struct Vardesc Vardesc;
177 typedef struct Namelist Namelist;
179 #define abs(x) ((x) >= 0 ? (x) : -(x))
180 #define dabs(x) (fabs(x))
181 #define f2cmin(a,b) ((a) <= (b) ? (a) : (b))
182 #define f2cmax(a,b) ((a) >= (b) ? (a) : (b))
183 #define dmin(a,b) (f2cmin(a,b))
184 #define dmax(a,b) (f2cmax(a,b))
185 #define bit_test(a,b) ((a) >> (b) & 1)
186 #define bit_clear(a,b) ((a) & ~((uinteger)1 << (b)))
187 #define bit_set(a,b) ((a) | ((uinteger)1 << (b)))
189 #define abort_() { sig_die("Fortran abort routine called", 1); }
190 #define c_abs(z) (cabsf(Cf(z)))
191 #define c_cos(R,Z) { pCf(R)=ccos(Cf(Z)); }
193 #define c_div(c, a, b) {Cf(c)._Val[0] = (Cf(a)._Val[0]/Cf(b)._Val[0]); Cf(c)._Val[1]=(Cf(a)._Val[1]/Cf(b)._Val[1]);}
194 #define z_div(c, a, b) {Cd(c)._Val[0] = (Cd(a)._Val[0]/Cd(b)._Val[0]); Cd(c)._Val[1]=(Cd(a)._Val[1]/Cd(b)._Val[1]);}
196 #define c_div(c, a, b) {pCf(c) = Cf(a)/Cf(b);}
197 #define z_div(c, a, b) {pCd(c) = Cd(a)/Cd(b);}
199 #define c_exp(R, Z) {pCf(R) = cexpf(Cf(Z));}
200 #define c_log(R, Z) {pCf(R) = clogf(Cf(Z));}
201 #define c_sin(R, Z) {pCf(R) = csinf(Cf(Z));}
202 //#define c_sqrt(R, Z) {*(R) = csqrtf(Cf(Z));}
203 #define c_sqrt(R, Z) {pCf(R) = csqrtf(Cf(Z));}
204 #define d_abs(x) (fabs(*(x)))
205 #define d_acos(x) (acos(*(x)))
206 #define d_asin(x) (asin(*(x)))
207 #define d_atan(x) (atan(*(x)))
208 #define d_atn2(x, y) (atan2(*(x),*(y)))
209 #define d_cnjg(R, Z) { pCd(R) = conj(Cd(Z)); }
210 #define r_cnjg(R, Z) { pCf(R) = conjf(Cf(Z)); }
211 #define d_cos(x) (cos(*(x)))
212 #define d_cosh(x) (cosh(*(x)))
213 #define d_dim(__a, __b) ( *(__a) > *(__b) ? *(__a) - *(__b) : 0.0 )
214 #define d_exp(x) (exp(*(x)))
215 #define d_imag(z) (cimag(Cd(z)))
216 #define r_imag(z) (cimagf(Cf(z)))
217 #define d_int(__x) (*(__x)>0 ? floor(*(__x)) : -floor(- *(__x)))
218 #define r_int(__x) (*(__x)>0 ? floor(*(__x)) : -floor(- *(__x)))
219 #define d_lg10(x) ( 0.43429448190325182765 * log(*(x)) )
220 #define r_lg10(x) ( 0.43429448190325182765 * log(*(x)) )
221 #define d_log(x) (log(*(x)))
222 #define d_mod(x, y) (fmod(*(x), *(y)))
223 #define u_nint(__x) ((__x)>=0 ? floor((__x) + .5) : -floor(.5 - (__x)))
224 #define d_nint(x) u_nint(*(x))
225 #define u_sign(__a,__b) ((__b) >= 0 ? ((__a) >= 0 ? (__a) : -(__a)) : -((__a) >= 0 ? (__a) : -(__a)))
226 #define d_sign(a,b) u_sign(*(a),*(b))
227 #define r_sign(a,b) u_sign(*(a),*(b))
228 #define d_sin(x) (sin(*(x)))
229 #define d_sinh(x) (sinh(*(x)))
230 #define d_sqrt(x) (sqrt(*(x)))
231 #define d_tan(x) (tan(*(x)))
232 #define d_tanh(x) (tanh(*(x)))
233 #define i_abs(x) abs(*(x))
234 #define i_dnnt(x) ((integer)u_nint(*(x)))
235 #define i_len(s, n) (n)
236 #define i_nint(x) ((integer)u_nint(*(x)))
237 #define i_sign(a,b) ((integer)u_sign((integer)*(a),(integer)*(b)))
238 #define pow_dd(ap, bp) ( pow(*(ap), *(bp)))
239 #define pow_si(B,E) spow_ui(*(B),*(E))
240 #define pow_ri(B,E) spow_ui(*(B),*(E))
241 #define pow_di(B,E) dpow_ui(*(B),*(E))
242 #define pow_zi(p, a, b) {pCd(p) = zpow_ui(Cd(a), *(b));}
243 #define pow_ci(p, a, b) {pCf(p) = cpow_ui(Cf(a), *(b));}
244 #define pow_zz(R,A,B) {pCd(R) = cpow(Cd(A),*(B));}
245 #define s_cat(lpp, rpp, rnp, np, llp) { ftnlen i, nc, ll; char *f__rp, *lp; ll = (llp); lp = (lpp); for(i=0; i < (int)*(np); ++i) { nc = ll; if((rnp)[i] < nc) nc = (rnp)[i]; ll -= nc; f__rp = (rpp)[i]; while(--nc >= 0) *lp++ = *(f__rp)++; } while(--ll >= 0) *lp++ = ' '; }
246 #define s_cmp(a,b,c,d) ((integer)strncmp((a),(b),f2cmin((c),(d))))
247 #define s_copy(A,B,C,D) { int __i,__m; for (__i=0, __m=f2cmin((C),(D)); __i<__m && (B)[__i] != 0; ++__i) (A)[__i] = (B)[__i]; }
248 #define sig_die(s, kill) { exit(1); }
249 #define s_stop(s, n) {exit(0);}
250 static char junk[] = "\n@(#)LIBF77 VERSION 19990503\n";
251 #define z_abs(z) (cabs(Cd(z)))
252 #define z_exp(R, Z) {pCd(R) = cexp(Cd(Z));}
253 #define z_sqrt(R, Z) {pCd(R) = csqrt(Cd(Z));}
254 #define myexit_() break;
255 #define mycycle() continue;
256 #define myceiling(w) {ceil(w)}
257 #define myhuge(w) {HUGE_VAL}
258 //#define mymaxloc_(w,s,e,n) {if (sizeof(*(w)) == sizeof(double)) dmaxloc_((w),*(s),*(e),n); else dmaxloc_((w),*(s),*(e),n);}
259 #define mymaxloc(w,s,e,n) {dmaxloc_(w,*(s),*(e),n)}
261 /* procedure parameter types for -A and -C++ */
263 #define F2C_proc_par_types 1
265 typedef logical (*L_fp)(...);
267 typedef logical (*L_fp)();
270 static float spow_ui(float x, integer n) {
271 float pow=1.0; unsigned long int u;
273 if(n < 0) n = -n, x = 1/x;
282 static double dpow_ui(double x, integer n) {
283 double pow=1.0; unsigned long int u;
285 if(n < 0) n = -n, x = 1/x;
295 static _Fcomplex cpow_ui(complex x, integer n) {
296 complex pow={1.0,0.0}; unsigned long int u;
298 if(n < 0) n = -n, x.r = 1/x.r, x.i=1/x.i;
300 if(u & 01) pow.r *= x.r, pow.i *= x.i;
301 if(u >>= 1) x.r *= x.r, x.i *= x.i;
305 _Fcomplex p={pow.r, pow.i};
309 static _Complex float cpow_ui(_Complex float x, integer n) {
310 _Complex float pow=1.0; unsigned long int u;
312 if(n < 0) n = -n, x = 1/x;
323 static _Dcomplex zpow_ui(_Dcomplex x, integer n) {
324 _Dcomplex pow={1.0,0.0}; unsigned long int u;
326 if(n < 0) n = -n, x._Val[0] = 1/x._Val[0], x._Val[1] =1/x._Val[1];
328 if(u & 01) pow._Val[0] *= x._Val[0], pow._Val[1] *= x._Val[1];
329 if(u >>= 1) x._Val[0] *= x._Val[0], x._Val[1] *= x._Val[1];
333 _Dcomplex p = {pow._Val[0], pow._Val[1]};
337 static _Complex double zpow_ui(_Complex double x, integer n) {
338 _Complex double pow=1.0; unsigned long int u;
340 if(n < 0) n = -n, x = 1/x;
350 static integer pow_ii(integer x, integer n) {
351 integer pow; unsigned long int u;
353 if (n == 0 || x == 1) pow = 1;
354 else if (x != -1) pow = x == 0 ? 1/x : 0;
357 if ((n > 0) || !(n == 0 || x == 1 || x != -1)) {
367 static integer dmaxloc_(double *w, integer s, integer e, integer *n)
369 double m; integer i, mi;
370 for(m=w[s-1], mi=s, i=s+1; i<=e; i++)
371 if (w[i-1]>m) mi=i ,m=w[i-1];
374 static integer smaxloc_(float *w, integer s, integer e, integer *n)
376 float m; integer i, mi;
377 for(m=w[s-1], mi=s, i=s+1; i<=e; i++)
378 if (w[i-1]>m) mi=i ,m=w[i-1];
381 static inline void cdotc_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
382 integer n = *n_, incx = *incx_, incy = *incy_, i;
384 _Fcomplex zdotc = {0.0, 0.0};
385 if (incx == 1 && incy == 1) {
386 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
387 zdotc._Val[0] += conjf(Cf(&x[i]))._Val[0] * Cf(&y[i])._Val[0];
388 zdotc._Val[1] += conjf(Cf(&x[i]))._Val[1] * Cf(&y[i])._Val[1];
391 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
392 zdotc._Val[0] += conjf(Cf(&x[i*incx]))._Val[0] * Cf(&y[i*incy])._Val[0];
393 zdotc._Val[1] += conjf(Cf(&x[i*incx]))._Val[1] * Cf(&y[i*incy])._Val[1];
399 _Complex float zdotc = 0.0;
400 if (incx == 1 && incy == 1) {
401 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
402 zdotc += conjf(Cf(&x[i])) * Cf(&y[i]);
405 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
406 zdotc += conjf(Cf(&x[i*incx])) * Cf(&y[i*incy]);
412 static inline void zdotc_(doublecomplex *z, integer *n_, doublecomplex *x, integer *incx_, doublecomplex *y, integer *incy_) {
413 integer n = *n_, incx = *incx_, incy = *incy_, i;
415 _Dcomplex zdotc = {0.0, 0.0};
416 if (incx == 1 && incy == 1) {
417 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
418 zdotc._Val[0] += conj(Cd(&x[i]))._Val[0] * Cd(&y[i])._Val[0];
419 zdotc._Val[1] += conj(Cd(&x[i]))._Val[1] * Cd(&y[i])._Val[1];
422 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
423 zdotc._Val[0] += conj(Cd(&x[i*incx]))._Val[0] * Cd(&y[i*incy])._Val[0];
424 zdotc._Val[1] += conj(Cd(&x[i*incx]))._Val[1] * Cd(&y[i*incy])._Val[1];
430 _Complex double zdotc = 0.0;
431 if (incx == 1 && incy == 1) {
432 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
433 zdotc += conj(Cd(&x[i])) * Cd(&y[i]);
436 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
437 zdotc += conj(Cd(&x[i*incx])) * Cd(&y[i*incy]);
443 static inline void cdotu_(complex *z, integer *n_, complex *x, integer *incx_, complex *y, integer *incy_) {
444 integer n = *n_, incx = *incx_, incy = *incy_, i;
446 _Fcomplex zdotc = {0.0, 0.0};
447 if (incx == 1 && incy == 1) {
448 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
449 zdotc._Val[0] += Cf(&x[i])._Val[0] * Cf(&y[i])._Val[0];
450 zdotc._Val[1] += Cf(&x[i])._Val[1] * Cf(&y[i])._Val[1];
453 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
454 zdotc._Val[0] += Cf(&x[i*incx])._Val[0] * Cf(&y[i*incy])._Val[0];
455 zdotc._Val[1] += Cf(&x[i*incx])._Val[1] * Cf(&y[i*incy])._Val[1];
461 _Complex float zdotc = 0.0;
462 if (incx == 1 && incy == 1) {
463 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
464 zdotc += Cf(&x[i]) * Cf(&y[i]);
467 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
468 zdotc += Cf(&x[i*incx]) * Cf(&y[i*incy]);
474 static inline void zdotu_(doublecomplex *z, integer *n_, doublecomplex *x, integer *incx_, doublecomplex *y, integer *incy_) {
475 integer n = *n_, incx = *incx_, incy = *incy_, i;
477 _Dcomplex zdotc = {0.0, 0.0};
478 if (incx == 1 && incy == 1) {
479 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
480 zdotc._Val[0] += Cd(&x[i])._Val[0] * Cd(&y[i])._Val[0];
481 zdotc._Val[1] += Cd(&x[i])._Val[1] * Cd(&y[i])._Val[1];
484 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
485 zdotc._Val[0] += Cd(&x[i*incx])._Val[0] * Cd(&y[i*incy])._Val[0];
486 zdotc._Val[1] += Cd(&x[i*incx])._Val[1] * Cd(&y[i*incy])._Val[1];
492 _Complex double zdotc = 0.0;
493 if (incx == 1 && incy == 1) {
494 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
495 zdotc += Cd(&x[i]) * Cd(&y[i]);
498 for (i=0;i<n;i++) { /* zdotc = zdotc + dconjg(x(i))* y(i) */
499 zdotc += Cd(&x[i*incx]) * Cd(&y[i*incy]);
505 /* -- translated by f2c (version 20000121).
506 You must link the resulting object file with the libraries:
507 -lf2c -lm (in that order)
513 /* Table of constant values */
515 static integer c__1 = 1;
517 /* > \brief \b ZLARFX applies an elementary reflector to a general rectangular matrix, with loop unrolling whe
518 n the reflector has order ≤ 10. */
520 /* =========== DOCUMENTATION =========== */
522 /* Online html documentation available at */
523 /* http://www.netlib.org/lapack/explore-html/ */
526 /* > Download ZLARFX + dependencies */
527 /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zlarfx.
530 /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zlarfx.
533 /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zlarfx.
541 /* SUBROUTINE ZLARFX( SIDE, M, N, V, TAU, C, LDC, WORK ) */
544 /* INTEGER LDC, M, N */
546 /* COMPLEX*16 C( LDC, * ), V( * ), WORK( * ) */
549 /* > \par Purpose: */
554 /* > ZLARFX applies a complex elementary reflector H to a complex m by n */
555 /* > matrix C, from either the left or the right. H is represented in the */
558 /* > H = I - tau * v * v**H */
560 /* > where tau is a complex scalar and v is a complex vector. */
562 /* > If tau = 0, then H is taken to be the unit matrix */
564 /* > This version uses inline code if H has order < 11. */
570 /* > \param[in] SIDE */
572 /* > SIDE is CHARACTER*1 */
573 /* > = 'L': form H * C */
574 /* > = 'R': form C * H */
580 /* > The number of rows of the matrix C. */
586 /* > The number of columns of the matrix C. */
591 /* > V is COMPLEX*16 array, dimension (M) if SIDE = 'L' */
592 /* > or (N) if SIDE = 'R' */
593 /* > The vector v in the representation of H. */
596 /* > \param[in] TAU */
598 /* > TAU is COMPLEX*16 */
599 /* > The value tau in the representation of H. */
602 /* > \param[in,out] C */
604 /* > C is COMPLEX*16 array, dimension (LDC,N) */
605 /* > On entry, the m by n matrix C. */
606 /* > On exit, C is overwritten by the matrix H * C if SIDE = 'L', */
607 /* > or C * H if SIDE = 'R'. */
610 /* > \param[in] LDC */
612 /* > LDC is INTEGER */
613 /* > The leading dimension of the array C. LDC >= f2cmax(1,M). */
616 /* > \param[out] WORK */
618 /* > WORK is COMPLEX*16 array, dimension (N) if SIDE = 'L' */
619 /* > or (M) if SIDE = 'R' */
620 /* > WORK is not referenced if H has order < 11. */
626 /* > \author Univ. of Tennessee */
627 /* > \author Univ. of California Berkeley */
628 /* > \author Univ. of Colorado Denver */
629 /* > \author NAG Ltd. */
631 /* > \date December 2016 */
633 /* > \ingroup complex16OTHERauxiliary */
635 /* ===================================================================== */
636 /* Subroutine */ int zlarfx_(char *side, integer *m, integer *n,
637 doublecomplex *v, doublecomplex *tau, doublecomplex *c__, integer *
638 ldc, doublecomplex *work)
640 /* System generated locals */
641 integer c_dim1, c_offset, i__1, i__2, i__3, i__4, i__5, i__6, i__7, i__8,
643 doublecomplex z__1, z__2, z__3, z__4, z__5, z__6, z__7, z__8, z__9, z__10,
644 z__11, z__12, z__13, z__14, z__15, z__16, z__17, z__18, z__19;
646 /* Local variables */
648 extern logical lsame_(char *, char *);
649 extern /* Subroutine */ int zlarf_(char *, integer *, integer *,
650 doublecomplex *, integer *, doublecomplex *, doublecomplex *,
651 integer *, doublecomplex *);
652 doublecomplex t1, t2, t3, t4, t5, t6, t7, t8, t9, v1, v2, v3, v4, v5, v6,
653 v7, v8, v9, t10, v10, sum;
656 /* -- LAPACK auxiliary routine (version 3.7.0) -- */
657 /* -- LAPACK is a software package provided by Univ. of Tennessee, -- */
658 /* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- */
662 /* ===================================================================== */
665 /* Parameter adjustments */
668 c_offset = 1 + c_dim1 * 1;
673 if (tau->r == 0. && tau->i == 0.) {
676 if (lsame_(side, "L")) {
678 /* Form H * C, where H has order m. */
693 /* Code for general M */
695 zlarf_(side, m, n, &v[1], &c__1, tau, &c__[c_offset], ldc, &work[1]);
699 /* Special code for 1 x 1 Householder */
701 z__3.r = tau->r * v[1].r - tau->i * v[1].i, z__3.i = tau->r * v[1].i
703 d_cnjg(&z__4, &v[1]);
704 z__2.r = z__3.r * z__4.r - z__3.i * z__4.i, z__2.i = z__3.r * z__4.i
706 z__1.r = 1. - z__2.r, z__1.i = 0. - z__2.i;
707 t1.r = z__1.r, t1.i = z__1.i;
709 for (j = 1; j <= i__1; ++j) {
710 i__2 = j * c_dim1 + 1;
711 i__3 = j * c_dim1 + 1;
712 z__1.r = t1.r * c__[i__3].r - t1.i * c__[i__3].i, z__1.i = t1.r *
713 c__[i__3].i + t1.i * c__[i__3].r;
714 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
720 /* Special code for 2 x 2 Householder */
722 d_cnjg(&z__1, &v[1]);
723 v1.r = z__1.r, v1.i = z__1.i;
725 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
727 t1.r = z__1.r, t1.i = z__1.i;
728 d_cnjg(&z__1, &v[2]);
729 v2.r = z__1.r, v2.i = z__1.i;
731 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
733 t2.r = z__1.r, t2.i = z__1.i;
735 for (j = 1; j <= i__1; ++j) {
736 i__2 = j * c_dim1 + 1;
737 z__2.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__2.i = v1.r *
738 c__[i__2].i + v1.i * c__[i__2].r;
739 i__3 = j * c_dim1 + 2;
740 z__3.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__3.i = v2.r *
741 c__[i__3].i + v2.i * c__[i__3].r;
742 z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
743 sum.r = z__1.r, sum.i = z__1.i;
744 i__2 = j * c_dim1 + 1;
745 i__3 = j * c_dim1 + 1;
746 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
748 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
749 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
750 i__2 = j * c_dim1 + 2;
751 i__3 = j * c_dim1 + 2;
752 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
754 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
755 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
761 /* Special code for 3 x 3 Householder */
763 d_cnjg(&z__1, &v[1]);
764 v1.r = z__1.r, v1.i = z__1.i;
766 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
768 t1.r = z__1.r, t1.i = z__1.i;
769 d_cnjg(&z__1, &v[2]);
770 v2.r = z__1.r, v2.i = z__1.i;
772 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
774 t2.r = z__1.r, t2.i = z__1.i;
775 d_cnjg(&z__1, &v[3]);
776 v3.r = z__1.r, v3.i = z__1.i;
778 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
780 t3.r = z__1.r, t3.i = z__1.i;
782 for (j = 1; j <= i__1; ++j) {
783 i__2 = j * c_dim1 + 1;
784 z__3.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__3.i = v1.r *
785 c__[i__2].i + v1.i * c__[i__2].r;
786 i__3 = j * c_dim1 + 2;
787 z__4.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__4.i = v2.r *
788 c__[i__3].i + v2.i * c__[i__3].r;
789 z__2.r = z__3.r + z__4.r, z__2.i = z__3.i + z__4.i;
790 i__4 = j * c_dim1 + 3;
791 z__5.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__5.i = v3.r *
792 c__[i__4].i + v3.i * c__[i__4].r;
793 z__1.r = z__2.r + z__5.r, z__1.i = z__2.i + z__5.i;
794 sum.r = z__1.r, sum.i = z__1.i;
795 i__2 = j * c_dim1 + 1;
796 i__3 = j * c_dim1 + 1;
797 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
799 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
800 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
801 i__2 = j * c_dim1 + 2;
802 i__3 = j * c_dim1 + 2;
803 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
805 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
806 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
807 i__2 = j * c_dim1 + 3;
808 i__3 = j * c_dim1 + 3;
809 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
811 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
812 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
818 /* Special code for 4 x 4 Householder */
820 d_cnjg(&z__1, &v[1]);
821 v1.r = z__1.r, v1.i = z__1.i;
823 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
825 t1.r = z__1.r, t1.i = z__1.i;
826 d_cnjg(&z__1, &v[2]);
827 v2.r = z__1.r, v2.i = z__1.i;
829 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
831 t2.r = z__1.r, t2.i = z__1.i;
832 d_cnjg(&z__1, &v[3]);
833 v3.r = z__1.r, v3.i = z__1.i;
835 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
837 t3.r = z__1.r, t3.i = z__1.i;
838 d_cnjg(&z__1, &v[4]);
839 v4.r = z__1.r, v4.i = z__1.i;
841 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
843 t4.r = z__1.r, t4.i = z__1.i;
845 for (j = 1; j <= i__1; ++j) {
846 i__2 = j * c_dim1 + 1;
847 z__4.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__4.i = v1.r *
848 c__[i__2].i + v1.i * c__[i__2].r;
849 i__3 = j * c_dim1 + 2;
850 z__5.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__5.i = v2.r *
851 c__[i__3].i + v2.i * c__[i__3].r;
852 z__3.r = z__4.r + z__5.r, z__3.i = z__4.i + z__5.i;
853 i__4 = j * c_dim1 + 3;
854 z__6.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__6.i = v3.r *
855 c__[i__4].i + v3.i * c__[i__4].r;
856 z__2.r = z__3.r + z__6.r, z__2.i = z__3.i + z__6.i;
857 i__5 = j * c_dim1 + 4;
858 z__7.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__7.i = v4.r *
859 c__[i__5].i + v4.i * c__[i__5].r;
860 z__1.r = z__2.r + z__7.r, z__1.i = z__2.i + z__7.i;
861 sum.r = z__1.r, sum.i = z__1.i;
862 i__2 = j * c_dim1 + 1;
863 i__3 = j * c_dim1 + 1;
864 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
866 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
867 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
868 i__2 = j * c_dim1 + 2;
869 i__3 = j * c_dim1 + 2;
870 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
872 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
873 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
874 i__2 = j * c_dim1 + 3;
875 i__3 = j * c_dim1 + 3;
876 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
878 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
879 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
880 i__2 = j * c_dim1 + 4;
881 i__3 = j * c_dim1 + 4;
882 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
884 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
885 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
891 /* Special code for 5 x 5 Householder */
893 d_cnjg(&z__1, &v[1]);
894 v1.r = z__1.r, v1.i = z__1.i;
896 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
898 t1.r = z__1.r, t1.i = z__1.i;
899 d_cnjg(&z__1, &v[2]);
900 v2.r = z__1.r, v2.i = z__1.i;
902 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
904 t2.r = z__1.r, t2.i = z__1.i;
905 d_cnjg(&z__1, &v[3]);
906 v3.r = z__1.r, v3.i = z__1.i;
908 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
910 t3.r = z__1.r, t3.i = z__1.i;
911 d_cnjg(&z__1, &v[4]);
912 v4.r = z__1.r, v4.i = z__1.i;
914 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
916 t4.r = z__1.r, t4.i = z__1.i;
917 d_cnjg(&z__1, &v[5]);
918 v5.r = z__1.r, v5.i = z__1.i;
920 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
922 t5.r = z__1.r, t5.i = z__1.i;
924 for (j = 1; j <= i__1; ++j) {
925 i__2 = j * c_dim1 + 1;
926 z__5.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__5.i = v1.r *
927 c__[i__2].i + v1.i * c__[i__2].r;
928 i__3 = j * c_dim1 + 2;
929 z__6.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__6.i = v2.r *
930 c__[i__3].i + v2.i * c__[i__3].r;
931 z__4.r = z__5.r + z__6.r, z__4.i = z__5.i + z__6.i;
932 i__4 = j * c_dim1 + 3;
933 z__7.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__7.i = v3.r *
934 c__[i__4].i + v3.i * c__[i__4].r;
935 z__3.r = z__4.r + z__7.r, z__3.i = z__4.i + z__7.i;
936 i__5 = j * c_dim1 + 4;
937 z__8.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__8.i = v4.r *
938 c__[i__5].i + v4.i * c__[i__5].r;
939 z__2.r = z__3.r + z__8.r, z__2.i = z__3.i + z__8.i;
940 i__6 = j * c_dim1 + 5;
941 z__9.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__9.i = v5.r *
942 c__[i__6].i + v5.i * c__[i__6].r;
943 z__1.r = z__2.r + z__9.r, z__1.i = z__2.i + z__9.i;
944 sum.r = z__1.r, sum.i = z__1.i;
945 i__2 = j * c_dim1 + 1;
946 i__3 = j * c_dim1 + 1;
947 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
949 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
950 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
951 i__2 = j * c_dim1 + 2;
952 i__3 = j * c_dim1 + 2;
953 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
955 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
956 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
957 i__2 = j * c_dim1 + 3;
958 i__3 = j * c_dim1 + 3;
959 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
961 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
962 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
963 i__2 = j * c_dim1 + 4;
964 i__3 = j * c_dim1 + 4;
965 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
967 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
968 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
969 i__2 = j * c_dim1 + 5;
970 i__3 = j * c_dim1 + 5;
971 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
973 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
974 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
980 /* Special code for 6 x 6 Householder */
982 d_cnjg(&z__1, &v[1]);
983 v1.r = z__1.r, v1.i = z__1.i;
985 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
987 t1.r = z__1.r, t1.i = z__1.i;
988 d_cnjg(&z__1, &v[2]);
989 v2.r = z__1.r, v2.i = z__1.i;
991 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
993 t2.r = z__1.r, t2.i = z__1.i;
994 d_cnjg(&z__1, &v[3]);
995 v3.r = z__1.r, v3.i = z__1.i;
997 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
999 t3.r = z__1.r, t3.i = z__1.i;
1000 d_cnjg(&z__1, &v[4]);
1001 v4.r = z__1.r, v4.i = z__1.i;
1003 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1005 t4.r = z__1.r, t4.i = z__1.i;
1006 d_cnjg(&z__1, &v[5]);
1007 v5.r = z__1.r, v5.i = z__1.i;
1009 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1011 t5.r = z__1.r, t5.i = z__1.i;
1012 d_cnjg(&z__1, &v[6]);
1013 v6.r = z__1.r, v6.i = z__1.i;
1015 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1017 t6.r = z__1.r, t6.i = z__1.i;
1019 for (j = 1; j <= i__1; ++j) {
1020 i__2 = j * c_dim1 + 1;
1021 z__6.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__6.i = v1.r *
1022 c__[i__2].i + v1.i * c__[i__2].r;
1023 i__3 = j * c_dim1 + 2;
1024 z__7.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__7.i = v2.r *
1025 c__[i__3].i + v2.i * c__[i__3].r;
1026 z__5.r = z__6.r + z__7.r, z__5.i = z__6.i + z__7.i;
1027 i__4 = j * c_dim1 + 3;
1028 z__8.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__8.i = v3.r *
1029 c__[i__4].i + v3.i * c__[i__4].r;
1030 z__4.r = z__5.r + z__8.r, z__4.i = z__5.i + z__8.i;
1031 i__5 = j * c_dim1 + 4;
1032 z__9.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__9.i = v4.r *
1033 c__[i__5].i + v4.i * c__[i__5].r;
1034 z__3.r = z__4.r + z__9.r, z__3.i = z__4.i + z__9.i;
1035 i__6 = j * c_dim1 + 5;
1036 z__10.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__10.i = v5.r
1037 * c__[i__6].i + v5.i * c__[i__6].r;
1038 z__2.r = z__3.r + z__10.r, z__2.i = z__3.i + z__10.i;
1039 i__7 = j * c_dim1 + 6;
1040 z__11.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__11.i = v6.r
1041 * c__[i__7].i + v6.i * c__[i__7].r;
1042 z__1.r = z__2.r + z__11.r, z__1.i = z__2.i + z__11.i;
1043 sum.r = z__1.r, sum.i = z__1.i;
1044 i__2 = j * c_dim1 + 1;
1045 i__3 = j * c_dim1 + 1;
1046 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1048 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1049 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1050 i__2 = j * c_dim1 + 2;
1051 i__3 = j * c_dim1 + 2;
1052 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1054 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1055 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1056 i__2 = j * c_dim1 + 3;
1057 i__3 = j * c_dim1 + 3;
1058 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1060 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1061 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1062 i__2 = j * c_dim1 + 4;
1063 i__3 = j * c_dim1 + 4;
1064 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1066 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1067 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1068 i__2 = j * c_dim1 + 5;
1069 i__3 = j * c_dim1 + 5;
1070 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1072 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1073 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1074 i__2 = j * c_dim1 + 6;
1075 i__3 = j * c_dim1 + 6;
1076 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
1078 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1079 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1085 /* Special code for 7 x 7 Householder */
1087 d_cnjg(&z__1, &v[1]);
1088 v1.r = z__1.r, v1.i = z__1.i;
1090 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1092 t1.r = z__1.r, t1.i = z__1.i;
1093 d_cnjg(&z__1, &v[2]);
1094 v2.r = z__1.r, v2.i = z__1.i;
1096 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1098 t2.r = z__1.r, t2.i = z__1.i;
1099 d_cnjg(&z__1, &v[3]);
1100 v3.r = z__1.r, v3.i = z__1.i;
1102 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1104 t3.r = z__1.r, t3.i = z__1.i;
1105 d_cnjg(&z__1, &v[4]);
1106 v4.r = z__1.r, v4.i = z__1.i;
1108 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1110 t4.r = z__1.r, t4.i = z__1.i;
1111 d_cnjg(&z__1, &v[5]);
1112 v5.r = z__1.r, v5.i = z__1.i;
1114 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1116 t5.r = z__1.r, t5.i = z__1.i;
1117 d_cnjg(&z__1, &v[6]);
1118 v6.r = z__1.r, v6.i = z__1.i;
1120 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1122 t6.r = z__1.r, t6.i = z__1.i;
1123 d_cnjg(&z__1, &v[7]);
1124 v7.r = z__1.r, v7.i = z__1.i;
1126 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1128 t7.r = z__1.r, t7.i = z__1.i;
1130 for (j = 1; j <= i__1; ++j) {
1131 i__2 = j * c_dim1 + 1;
1132 z__7.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__7.i = v1.r *
1133 c__[i__2].i + v1.i * c__[i__2].r;
1134 i__3 = j * c_dim1 + 2;
1135 z__8.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__8.i = v2.r *
1136 c__[i__3].i + v2.i * c__[i__3].r;
1137 z__6.r = z__7.r + z__8.r, z__6.i = z__7.i + z__8.i;
1138 i__4 = j * c_dim1 + 3;
1139 z__9.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__9.i = v3.r *
1140 c__[i__4].i + v3.i * c__[i__4].r;
1141 z__5.r = z__6.r + z__9.r, z__5.i = z__6.i + z__9.i;
1142 i__5 = j * c_dim1 + 4;
1143 z__10.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__10.i = v4.r
1144 * c__[i__5].i + v4.i * c__[i__5].r;
1145 z__4.r = z__5.r + z__10.r, z__4.i = z__5.i + z__10.i;
1146 i__6 = j * c_dim1 + 5;
1147 z__11.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__11.i = v5.r
1148 * c__[i__6].i + v5.i * c__[i__6].r;
1149 z__3.r = z__4.r + z__11.r, z__3.i = z__4.i + z__11.i;
1150 i__7 = j * c_dim1 + 6;
1151 z__12.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__12.i = v6.r
1152 * c__[i__7].i + v6.i * c__[i__7].r;
1153 z__2.r = z__3.r + z__12.r, z__2.i = z__3.i + z__12.i;
1154 i__8 = j * c_dim1 + 7;
1155 z__13.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__13.i = v7.r
1156 * c__[i__8].i + v7.i * c__[i__8].r;
1157 z__1.r = z__2.r + z__13.r, z__1.i = z__2.i + z__13.i;
1158 sum.r = z__1.r, sum.i = z__1.i;
1159 i__2 = j * c_dim1 + 1;
1160 i__3 = j * c_dim1 + 1;
1161 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1163 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1164 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1165 i__2 = j * c_dim1 + 2;
1166 i__3 = j * c_dim1 + 2;
1167 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1169 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1170 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1171 i__2 = j * c_dim1 + 3;
1172 i__3 = j * c_dim1 + 3;
1173 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1175 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1176 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1177 i__2 = j * c_dim1 + 4;
1178 i__3 = j * c_dim1 + 4;
1179 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1181 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1182 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1183 i__2 = j * c_dim1 + 5;
1184 i__3 = j * c_dim1 + 5;
1185 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1187 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1188 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1189 i__2 = j * c_dim1 + 6;
1190 i__3 = j * c_dim1 + 6;
1191 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
1193 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1194 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1195 i__2 = j * c_dim1 + 7;
1196 i__3 = j * c_dim1 + 7;
1197 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
1199 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1200 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1206 /* Special code for 8 x 8 Householder */
1208 d_cnjg(&z__1, &v[1]);
1209 v1.r = z__1.r, v1.i = z__1.i;
1211 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1213 t1.r = z__1.r, t1.i = z__1.i;
1214 d_cnjg(&z__1, &v[2]);
1215 v2.r = z__1.r, v2.i = z__1.i;
1217 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1219 t2.r = z__1.r, t2.i = z__1.i;
1220 d_cnjg(&z__1, &v[3]);
1221 v3.r = z__1.r, v3.i = z__1.i;
1223 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1225 t3.r = z__1.r, t3.i = z__1.i;
1226 d_cnjg(&z__1, &v[4]);
1227 v4.r = z__1.r, v4.i = z__1.i;
1229 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1231 t4.r = z__1.r, t4.i = z__1.i;
1232 d_cnjg(&z__1, &v[5]);
1233 v5.r = z__1.r, v5.i = z__1.i;
1235 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1237 t5.r = z__1.r, t5.i = z__1.i;
1238 d_cnjg(&z__1, &v[6]);
1239 v6.r = z__1.r, v6.i = z__1.i;
1241 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1243 t6.r = z__1.r, t6.i = z__1.i;
1244 d_cnjg(&z__1, &v[7]);
1245 v7.r = z__1.r, v7.i = z__1.i;
1247 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1249 t7.r = z__1.r, t7.i = z__1.i;
1250 d_cnjg(&z__1, &v[8]);
1251 v8.r = z__1.r, v8.i = z__1.i;
1253 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1255 t8.r = z__1.r, t8.i = z__1.i;
1257 for (j = 1; j <= i__1; ++j) {
1258 i__2 = j * c_dim1 + 1;
1259 z__8.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__8.i = v1.r *
1260 c__[i__2].i + v1.i * c__[i__2].r;
1261 i__3 = j * c_dim1 + 2;
1262 z__9.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__9.i = v2.r *
1263 c__[i__3].i + v2.i * c__[i__3].r;
1264 z__7.r = z__8.r + z__9.r, z__7.i = z__8.i + z__9.i;
1265 i__4 = j * c_dim1 + 3;
1266 z__10.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__10.i = v3.r
1267 * c__[i__4].i + v3.i * c__[i__4].r;
1268 z__6.r = z__7.r + z__10.r, z__6.i = z__7.i + z__10.i;
1269 i__5 = j * c_dim1 + 4;
1270 z__11.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__11.i = v4.r
1271 * c__[i__5].i + v4.i * c__[i__5].r;
1272 z__5.r = z__6.r + z__11.r, z__5.i = z__6.i + z__11.i;
1273 i__6 = j * c_dim1 + 5;
1274 z__12.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__12.i = v5.r
1275 * c__[i__6].i + v5.i * c__[i__6].r;
1276 z__4.r = z__5.r + z__12.r, z__4.i = z__5.i + z__12.i;
1277 i__7 = j * c_dim1 + 6;
1278 z__13.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__13.i = v6.r
1279 * c__[i__7].i + v6.i * c__[i__7].r;
1280 z__3.r = z__4.r + z__13.r, z__3.i = z__4.i + z__13.i;
1281 i__8 = j * c_dim1 + 7;
1282 z__14.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__14.i = v7.r
1283 * c__[i__8].i + v7.i * c__[i__8].r;
1284 z__2.r = z__3.r + z__14.r, z__2.i = z__3.i + z__14.i;
1285 i__9 = j * c_dim1 + 8;
1286 z__15.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__15.i = v8.r
1287 * c__[i__9].i + v8.i * c__[i__9].r;
1288 z__1.r = z__2.r + z__15.r, z__1.i = z__2.i + z__15.i;
1289 sum.r = z__1.r, sum.i = z__1.i;
1290 i__2 = j * c_dim1 + 1;
1291 i__3 = j * c_dim1 + 1;
1292 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1294 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1295 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1296 i__2 = j * c_dim1 + 2;
1297 i__3 = j * c_dim1 + 2;
1298 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1300 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1301 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1302 i__2 = j * c_dim1 + 3;
1303 i__3 = j * c_dim1 + 3;
1304 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1306 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1307 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1308 i__2 = j * c_dim1 + 4;
1309 i__3 = j * c_dim1 + 4;
1310 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1312 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1313 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1314 i__2 = j * c_dim1 + 5;
1315 i__3 = j * c_dim1 + 5;
1316 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1318 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1319 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1320 i__2 = j * c_dim1 + 6;
1321 i__3 = j * c_dim1 + 6;
1322 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
1324 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1325 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1326 i__2 = j * c_dim1 + 7;
1327 i__3 = j * c_dim1 + 7;
1328 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
1330 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1331 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1332 i__2 = j * c_dim1 + 8;
1333 i__3 = j * c_dim1 + 8;
1334 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
1336 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1337 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1343 /* Special code for 9 x 9 Householder */
1345 d_cnjg(&z__1, &v[1]);
1346 v1.r = z__1.r, v1.i = z__1.i;
1348 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1350 t1.r = z__1.r, t1.i = z__1.i;
1351 d_cnjg(&z__1, &v[2]);
1352 v2.r = z__1.r, v2.i = z__1.i;
1354 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1356 t2.r = z__1.r, t2.i = z__1.i;
1357 d_cnjg(&z__1, &v[3]);
1358 v3.r = z__1.r, v3.i = z__1.i;
1360 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1362 t3.r = z__1.r, t3.i = z__1.i;
1363 d_cnjg(&z__1, &v[4]);
1364 v4.r = z__1.r, v4.i = z__1.i;
1366 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1368 t4.r = z__1.r, t4.i = z__1.i;
1369 d_cnjg(&z__1, &v[5]);
1370 v5.r = z__1.r, v5.i = z__1.i;
1372 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1374 t5.r = z__1.r, t5.i = z__1.i;
1375 d_cnjg(&z__1, &v[6]);
1376 v6.r = z__1.r, v6.i = z__1.i;
1378 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1380 t6.r = z__1.r, t6.i = z__1.i;
1381 d_cnjg(&z__1, &v[7]);
1382 v7.r = z__1.r, v7.i = z__1.i;
1384 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1386 t7.r = z__1.r, t7.i = z__1.i;
1387 d_cnjg(&z__1, &v[8]);
1388 v8.r = z__1.r, v8.i = z__1.i;
1390 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1392 t8.r = z__1.r, t8.i = z__1.i;
1393 d_cnjg(&z__1, &v[9]);
1394 v9.r = z__1.r, v9.i = z__1.i;
1396 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1398 t9.r = z__1.r, t9.i = z__1.i;
1400 for (j = 1; j <= i__1; ++j) {
1401 i__2 = j * c_dim1 + 1;
1402 z__9.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__9.i = v1.r *
1403 c__[i__2].i + v1.i * c__[i__2].r;
1404 i__3 = j * c_dim1 + 2;
1405 z__10.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__10.i = v2.r
1406 * c__[i__3].i + v2.i * c__[i__3].r;
1407 z__8.r = z__9.r + z__10.r, z__8.i = z__9.i + z__10.i;
1408 i__4 = j * c_dim1 + 3;
1409 z__11.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__11.i = v3.r
1410 * c__[i__4].i + v3.i * c__[i__4].r;
1411 z__7.r = z__8.r + z__11.r, z__7.i = z__8.i + z__11.i;
1412 i__5 = j * c_dim1 + 4;
1413 z__12.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__12.i = v4.r
1414 * c__[i__5].i + v4.i * c__[i__5].r;
1415 z__6.r = z__7.r + z__12.r, z__6.i = z__7.i + z__12.i;
1416 i__6 = j * c_dim1 + 5;
1417 z__13.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__13.i = v5.r
1418 * c__[i__6].i + v5.i * c__[i__6].r;
1419 z__5.r = z__6.r + z__13.r, z__5.i = z__6.i + z__13.i;
1420 i__7 = j * c_dim1 + 6;
1421 z__14.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__14.i = v6.r
1422 * c__[i__7].i + v6.i * c__[i__7].r;
1423 z__4.r = z__5.r + z__14.r, z__4.i = z__5.i + z__14.i;
1424 i__8 = j * c_dim1 + 7;
1425 z__15.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__15.i = v7.r
1426 * c__[i__8].i + v7.i * c__[i__8].r;
1427 z__3.r = z__4.r + z__15.r, z__3.i = z__4.i + z__15.i;
1428 i__9 = j * c_dim1 + 8;
1429 z__16.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__16.i = v8.r
1430 * c__[i__9].i + v8.i * c__[i__9].r;
1431 z__2.r = z__3.r + z__16.r, z__2.i = z__3.i + z__16.i;
1432 i__10 = j * c_dim1 + 9;
1433 z__17.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i, z__17.i =
1434 v9.r * c__[i__10].i + v9.i * c__[i__10].r;
1435 z__1.r = z__2.r + z__17.r, z__1.i = z__2.i + z__17.i;
1436 sum.r = z__1.r, sum.i = z__1.i;
1437 i__2 = j * c_dim1 + 1;
1438 i__3 = j * c_dim1 + 1;
1439 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1441 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1442 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1443 i__2 = j * c_dim1 + 2;
1444 i__3 = j * c_dim1 + 2;
1445 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1447 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1448 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1449 i__2 = j * c_dim1 + 3;
1450 i__3 = j * c_dim1 + 3;
1451 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1453 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1454 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1455 i__2 = j * c_dim1 + 4;
1456 i__3 = j * c_dim1 + 4;
1457 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1459 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1460 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1461 i__2 = j * c_dim1 + 5;
1462 i__3 = j * c_dim1 + 5;
1463 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1465 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1466 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1467 i__2 = j * c_dim1 + 6;
1468 i__3 = j * c_dim1 + 6;
1469 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
1471 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1472 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1473 i__2 = j * c_dim1 + 7;
1474 i__3 = j * c_dim1 + 7;
1475 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
1477 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1478 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1479 i__2 = j * c_dim1 + 8;
1480 i__3 = j * c_dim1 + 8;
1481 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
1483 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1484 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1485 i__2 = j * c_dim1 + 9;
1486 i__3 = j * c_dim1 + 9;
1487 z__2.r = sum.r * t9.r - sum.i * t9.i, z__2.i = sum.r * t9.i +
1489 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1490 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1496 /* Special code for 10 x 10 Householder */
1498 d_cnjg(&z__1, &v[1]);
1499 v1.r = z__1.r, v1.i = z__1.i;
1501 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1503 t1.r = z__1.r, t1.i = z__1.i;
1504 d_cnjg(&z__1, &v[2]);
1505 v2.r = z__1.r, v2.i = z__1.i;
1507 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1509 t2.r = z__1.r, t2.i = z__1.i;
1510 d_cnjg(&z__1, &v[3]);
1511 v3.r = z__1.r, v3.i = z__1.i;
1513 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1515 t3.r = z__1.r, t3.i = z__1.i;
1516 d_cnjg(&z__1, &v[4]);
1517 v4.r = z__1.r, v4.i = z__1.i;
1519 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1521 t4.r = z__1.r, t4.i = z__1.i;
1522 d_cnjg(&z__1, &v[5]);
1523 v5.r = z__1.r, v5.i = z__1.i;
1525 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1527 t5.r = z__1.r, t5.i = z__1.i;
1528 d_cnjg(&z__1, &v[6]);
1529 v6.r = z__1.r, v6.i = z__1.i;
1531 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1533 t6.r = z__1.r, t6.i = z__1.i;
1534 d_cnjg(&z__1, &v[7]);
1535 v7.r = z__1.r, v7.i = z__1.i;
1537 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1539 t7.r = z__1.r, t7.i = z__1.i;
1540 d_cnjg(&z__1, &v[8]);
1541 v8.r = z__1.r, v8.i = z__1.i;
1543 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1545 t8.r = z__1.r, t8.i = z__1.i;
1546 d_cnjg(&z__1, &v[9]);
1547 v9.r = z__1.r, v9.i = z__1.i;
1549 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1551 t9.r = z__1.r, t9.i = z__1.i;
1552 d_cnjg(&z__1, &v[10]);
1553 v10.r = z__1.r, v10.i = z__1.i;
1554 d_cnjg(&z__2, &v10);
1555 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1557 t10.r = z__1.r, t10.i = z__1.i;
1559 for (j = 1; j <= i__1; ++j) {
1560 i__2 = j * c_dim1 + 1;
1561 z__10.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__10.i = v1.r
1562 * c__[i__2].i + v1.i * c__[i__2].r;
1563 i__3 = j * c_dim1 + 2;
1564 z__11.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__11.i = v2.r
1565 * c__[i__3].i + v2.i * c__[i__3].r;
1566 z__9.r = z__10.r + z__11.r, z__9.i = z__10.i + z__11.i;
1567 i__4 = j * c_dim1 + 3;
1568 z__12.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__12.i = v3.r
1569 * c__[i__4].i + v3.i * c__[i__4].r;
1570 z__8.r = z__9.r + z__12.r, z__8.i = z__9.i + z__12.i;
1571 i__5 = j * c_dim1 + 4;
1572 z__13.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__13.i = v4.r
1573 * c__[i__5].i + v4.i * c__[i__5].r;
1574 z__7.r = z__8.r + z__13.r, z__7.i = z__8.i + z__13.i;
1575 i__6 = j * c_dim1 + 5;
1576 z__14.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__14.i = v5.r
1577 * c__[i__6].i + v5.i * c__[i__6].r;
1578 z__6.r = z__7.r + z__14.r, z__6.i = z__7.i + z__14.i;
1579 i__7 = j * c_dim1 + 6;
1580 z__15.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__15.i = v6.r
1581 * c__[i__7].i + v6.i * c__[i__7].r;
1582 z__5.r = z__6.r + z__15.r, z__5.i = z__6.i + z__15.i;
1583 i__8 = j * c_dim1 + 7;
1584 z__16.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__16.i = v7.r
1585 * c__[i__8].i + v7.i * c__[i__8].r;
1586 z__4.r = z__5.r + z__16.r, z__4.i = z__5.i + z__16.i;
1587 i__9 = j * c_dim1 + 8;
1588 z__17.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__17.i = v8.r
1589 * c__[i__9].i + v8.i * c__[i__9].r;
1590 z__3.r = z__4.r + z__17.r, z__3.i = z__4.i + z__17.i;
1591 i__10 = j * c_dim1 + 9;
1592 z__18.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i, z__18.i =
1593 v9.r * c__[i__10].i + v9.i * c__[i__10].r;
1594 z__2.r = z__3.r + z__18.r, z__2.i = z__3.i + z__18.i;
1595 i__11 = j * c_dim1 + 10;
1596 z__19.r = v10.r * c__[i__11].r - v10.i * c__[i__11].i, z__19.i =
1597 v10.r * c__[i__11].i + v10.i * c__[i__11].r;
1598 z__1.r = z__2.r + z__19.r, z__1.i = z__2.i + z__19.i;
1599 sum.r = z__1.r, sum.i = z__1.i;
1600 i__2 = j * c_dim1 + 1;
1601 i__3 = j * c_dim1 + 1;
1602 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1604 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1605 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1606 i__2 = j * c_dim1 + 2;
1607 i__3 = j * c_dim1 + 2;
1608 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1610 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1611 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1612 i__2 = j * c_dim1 + 3;
1613 i__3 = j * c_dim1 + 3;
1614 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1616 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1617 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1618 i__2 = j * c_dim1 + 4;
1619 i__3 = j * c_dim1 + 4;
1620 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1622 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1623 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1624 i__2 = j * c_dim1 + 5;
1625 i__3 = j * c_dim1 + 5;
1626 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1628 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1629 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1630 i__2 = j * c_dim1 + 6;
1631 i__3 = j * c_dim1 + 6;
1632 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
1634 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1635 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1636 i__2 = j * c_dim1 + 7;
1637 i__3 = j * c_dim1 + 7;
1638 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
1640 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1641 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1642 i__2 = j * c_dim1 + 8;
1643 i__3 = j * c_dim1 + 8;
1644 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
1646 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1647 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1648 i__2 = j * c_dim1 + 9;
1649 i__3 = j * c_dim1 + 9;
1650 z__2.r = sum.r * t9.r - sum.i * t9.i, z__2.i = sum.r * t9.i +
1652 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1653 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1654 i__2 = j * c_dim1 + 10;
1655 i__3 = j * c_dim1 + 10;
1656 z__2.r = sum.r * t10.r - sum.i * t10.i, z__2.i = sum.r * t10.i +
1658 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1659 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1665 /* Form C * H, where H has order n. */
1680 /* Code for general N */
1682 zlarf_(side, m, n, &v[1], &c__1, tau, &c__[c_offset], ldc, &work[1]);
1686 /* Special code for 1 x 1 Householder */
1688 z__3.r = tau->r * v[1].r - tau->i * v[1].i, z__3.i = tau->r * v[1].i
1690 d_cnjg(&z__4, &v[1]);
1691 z__2.r = z__3.r * z__4.r - z__3.i * z__4.i, z__2.i = z__3.r * z__4.i
1693 z__1.r = 1. - z__2.r, z__1.i = 0. - z__2.i;
1694 t1.r = z__1.r, t1.i = z__1.i;
1696 for (j = 1; j <= i__1; ++j) {
1699 z__1.r = t1.r * c__[i__3].r - t1.i * c__[i__3].i, z__1.i = t1.r *
1700 c__[i__3].i + t1.i * c__[i__3].r;
1701 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1707 /* Special code for 2 x 2 Householder */
1709 v1.r = v[1].r, v1.i = v[1].i;
1711 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1713 t1.r = z__1.r, t1.i = z__1.i;
1714 v2.r = v[2].r, v2.i = v[2].i;
1716 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1718 t2.r = z__1.r, t2.i = z__1.i;
1720 for (j = 1; j <= i__1; ++j) {
1722 z__2.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__2.i = v1.r *
1723 c__[i__2].i + v1.i * c__[i__2].r;
1724 i__3 = j + (c_dim1 << 1);
1725 z__3.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__3.i = v2.r *
1726 c__[i__3].i + v2.i * c__[i__3].r;
1727 z__1.r = z__2.r + z__3.r, z__1.i = z__2.i + z__3.i;
1728 sum.r = z__1.r, sum.i = z__1.i;
1731 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1733 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1734 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1735 i__2 = j + (c_dim1 << 1);
1736 i__3 = j + (c_dim1 << 1);
1737 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1739 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1740 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1746 /* Special code for 3 x 3 Householder */
1748 v1.r = v[1].r, v1.i = v[1].i;
1750 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1752 t1.r = z__1.r, t1.i = z__1.i;
1753 v2.r = v[2].r, v2.i = v[2].i;
1755 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1757 t2.r = z__1.r, t2.i = z__1.i;
1758 v3.r = v[3].r, v3.i = v[3].i;
1760 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1762 t3.r = z__1.r, t3.i = z__1.i;
1764 for (j = 1; j <= i__1; ++j) {
1766 z__3.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__3.i = v1.r *
1767 c__[i__2].i + v1.i * c__[i__2].r;
1768 i__3 = j + (c_dim1 << 1);
1769 z__4.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__4.i = v2.r *
1770 c__[i__3].i + v2.i * c__[i__3].r;
1771 z__2.r = z__3.r + z__4.r, z__2.i = z__3.i + z__4.i;
1772 i__4 = j + c_dim1 * 3;
1773 z__5.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__5.i = v3.r *
1774 c__[i__4].i + v3.i * c__[i__4].r;
1775 z__1.r = z__2.r + z__5.r, z__1.i = z__2.i + z__5.i;
1776 sum.r = z__1.r, sum.i = z__1.i;
1779 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1781 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1782 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1783 i__2 = j + (c_dim1 << 1);
1784 i__3 = j + (c_dim1 << 1);
1785 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1787 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1788 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1789 i__2 = j + c_dim1 * 3;
1790 i__3 = j + c_dim1 * 3;
1791 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1793 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1794 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1800 /* Special code for 4 x 4 Householder */
1802 v1.r = v[1].r, v1.i = v[1].i;
1804 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1806 t1.r = z__1.r, t1.i = z__1.i;
1807 v2.r = v[2].r, v2.i = v[2].i;
1809 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1811 t2.r = z__1.r, t2.i = z__1.i;
1812 v3.r = v[3].r, v3.i = v[3].i;
1814 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1816 t3.r = z__1.r, t3.i = z__1.i;
1817 v4.r = v[4].r, v4.i = v[4].i;
1819 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1821 t4.r = z__1.r, t4.i = z__1.i;
1823 for (j = 1; j <= i__1; ++j) {
1825 z__4.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__4.i = v1.r *
1826 c__[i__2].i + v1.i * c__[i__2].r;
1827 i__3 = j + (c_dim1 << 1);
1828 z__5.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__5.i = v2.r *
1829 c__[i__3].i + v2.i * c__[i__3].r;
1830 z__3.r = z__4.r + z__5.r, z__3.i = z__4.i + z__5.i;
1831 i__4 = j + c_dim1 * 3;
1832 z__6.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__6.i = v3.r *
1833 c__[i__4].i + v3.i * c__[i__4].r;
1834 z__2.r = z__3.r + z__6.r, z__2.i = z__3.i + z__6.i;
1835 i__5 = j + (c_dim1 << 2);
1836 z__7.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__7.i = v4.r *
1837 c__[i__5].i + v4.i * c__[i__5].r;
1838 z__1.r = z__2.r + z__7.r, z__1.i = z__2.i + z__7.i;
1839 sum.r = z__1.r, sum.i = z__1.i;
1842 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1844 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1845 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1846 i__2 = j + (c_dim1 << 1);
1847 i__3 = j + (c_dim1 << 1);
1848 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1850 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1851 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1852 i__2 = j + c_dim1 * 3;
1853 i__3 = j + c_dim1 * 3;
1854 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1856 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1857 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1858 i__2 = j + (c_dim1 << 2);
1859 i__3 = j + (c_dim1 << 2);
1860 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1862 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1863 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1869 /* Special code for 5 x 5 Householder */
1871 v1.r = v[1].r, v1.i = v[1].i;
1873 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1875 t1.r = z__1.r, t1.i = z__1.i;
1876 v2.r = v[2].r, v2.i = v[2].i;
1878 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1880 t2.r = z__1.r, t2.i = z__1.i;
1881 v3.r = v[3].r, v3.i = v[3].i;
1883 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1885 t3.r = z__1.r, t3.i = z__1.i;
1886 v4.r = v[4].r, v4.i = v[4].i;
1888 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1890 t4.r = z__1.r, t4.i = z__1.i;
1891 v5.r = v[5].r, v5.i = v[5].i;
1893 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1895 t5.r = z__1.r, t5.i = z__1.i;
1897 for (j = 1; j <= i__1; ++j) {
1899 z__5.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__5.i = v1.r *
1900 c__[i__2].i + v1.i * c__[i__2].r;
1901 i__3 = j + (c_dim1 << 1);
1902 z__6.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__6.i = v2.r *
1903 c__[i__3].i + v2.i * c__[i__3].r;
1904 z__4.r = z__5.r + z__6.r, z__4.i = z__5.i + z__6.i;
1905 i__4 = j + c_dim1 * 3;
1906 z__7.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__7.i = v3.r *
1907 c__[i__4].i + v3.i * c__[i__4].r;
1908 z__3.r = z__4.r + z__7.r, z__3.i = z__4.i + z__7.i;
1909 i__5 = j + (c_dim1 << 2);
1910 z__8.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__8.i = v4.r *
1911 c__[i__5].i + v4.i * c__[i__5].r;
1912 z__2.r = z__3.r + z__8.r, z__2.i = z__3.i + z__8.i;
1913 i__6 = j + c_dim1 * 5;
1914 z__9.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__9.i = v5.r *
1915 c__[i__6].i + v5.i * c__[i__6].r;
1916 z__1.r = z__2.r + z__9.r, z__1.i = z__2.i + z__9.i;
1917 sum.r = z__1.r, sum.i = z__1.i;
1920 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
1922 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1923 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1924 i__2 = j + (c_dim1 << 1);
1925 i__3 = j + (c_dim1 << 1);
1926 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
1928 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1929 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1930 i__2 = j + c_dim1 * 3;
1931 i__3 = j + c_dim1 * 3;
1932 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
1934 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1935 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1936 i__2 = j + (c_dim1 << 2);
1937 i__3 = j + (c_dim1 << 2);
1938 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
1940 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1941 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1942 i__2 = j + c_dim1 * 5;
1943 i__3 = j + c_dim1 * 5;
1944 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
1946 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
1947 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
1953 /* Special code for 6 x 6 Householder */
1955 v1.r = v[1].r, v1.i = v[1].i;
1957 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1959 t1.r = z__1.r, t1.i = z__1.i;
1960 v2.r = v[2].r, v2.i = v[2].i;
1962 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1964 t2.r = z__1.r, t2.i = z__1.i;
1965 v3.r = v[3].r, v3.i = v[3].i;
1967 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1969 t3.r = z__1.r, t3.i = z__1.i;
1970 v4.r = v[4].r, v4.i = v[4].i;
1972 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1974 t4.r = z__1.r, t4.i = z__1.i;
1975 v5.r = v[5].r, v5.i = v[5].i;
1977 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1979 t5.r = z__1.r, t5.i = z__1.i;
1980 v6.r = v[6].r, v6.i = v[6].i;
1982 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
1984 t6.r = z__1.r, t6.i = z__1.i;
1986 for (j = 1; j <= i__1; ++j) {
1988 z__6.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__6.i = v1.r *
1989 c__[i__2].i + v1.i * c__[i__2].r;
1990 i__3 = j + (c_dim1 << 1);
1991 z__7.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__7.i = v2.r *
1992 c__[i__3].i + v2.i * c__[i__3].r;
1993 z__5.r = z__6.r + z__7.r, z__5.i = z__6.i + z__7.i;
1994 i__4 = j + c_dim1 * 3;
1995 z__8.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__8.i = v3.r *
1996 c__[i__4].i + v3.i * c__[i__4].r;
1997 z__4.r = z__5.r + z__8.r, z__4.i = z__5.i + z__8.i;
1998 i__5 = j + (c_dim1 << 2);
1999 z__9.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__9.i = v4.r *
2000 c__[i__5].i + v4.i * c__[i__5].r;
2001 z__3.r = z__4.r + z__9.r, z__3.i = z__4.i + z__9.i;
2002 i__6 = j + c_dim1 * 5;
2003 z__10.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__10.i = v5.r
2004 * c__[i__6].i + v5.i * c__[i__6].r;
2005 z__2.r = z__3.r + z__10.r, z__2.i = z__3.i + z__10.i;
2006 i__7 = j + c_dim1 * 6;
2007 z__11.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__11.i = v6.r
2008 * c__[i__7].i + v6.i * c__[i__7].r;
2009 z__1.r = z__2.r + z__11.r, z__1.i = z__2.i + z__11.i;
2010 sum.r = z__1.r, sum.i = z__1.i;
2013 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
2015 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2016 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2017 i__2 = j + (c_dim1 << 1);
2018 i__3 = j + (c_dim1 << 1);
2019 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
2021 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2022 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2023 i__2 = j + c_dim1 * 3;
2024 i__3 = j + c_dim1 * 3;
2025 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
2027 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2028 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2029 i__2 = j + (c_dim1 << 2);
2030 i__3 = j + (c_dim1 << 2);
2031 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
2033 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2034 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2035 i__2 = j + c_dim1 * 5;
2036 i__3 = j + c_dim1 * 5;
2037 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
2039 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2040 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2041 i__2 = j + c_dim1 * 6;
2042 i__3 = j + c_dim1 * 6;
2043 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
2045 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2046 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2052 /* Special code for 7 x 7 Householder */
2054 v1.r = v[1].r, v1.i = v[1].i;
2056 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2058 t1.r = z__1.r, t1.i = z__1.i;
2059 v2.r = v[2].r, v2.i = v[2].i;
2061 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2063 t2.r = z__1.r, t2.i = z__1.i;
2064 v3.r = v[3].r, v3.i = v[3].i;
2066 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2068 t3.r = z__1.r, t3.i = z__1.i;
2069 v4.r = v[4].r, v4.i = v[4].i;
2071 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2073 t4.r = z__1.r, t4.i = z__1.i;
2074 v5.r = v[5].r, v5.i = v[5].i;
2076 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2078 t5.r = z__1.r, t5.i = z__1.i;
2079 v6.r = v[6].r, v6.i = v[6].i;
2081 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2083 t6.r = z__1.r, t6.i = z__1.i;
2084 v7.r = v[7].r, v7.i = v[7].i;
2086 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2088 t7.r = z__1.r, t7.i = z__1.i;
2090 for (j = 1; j <= i__1; ++j) {
2092 z__7.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__7.i = v1.r *
2093 c__[i__2].i + v1.i * c__[i__2].r;
2094 i__3 = j + (c_dim1 << 1);
2095 z__8.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__8.i = v2.r *
2096 c__[i__3].i + v2.i * c__[i__3].r;
2097 z__6.r = z__7.r + z__8.r, z__6.i = z__7.i + z__8.i;
2098 i__4 = j + c_dim1 * 3;
2099 z__9.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__9.i = v3.r *
2100 c__[i__4].i + v3.i * c__[i__4].r;
2101 z__5.r = z__6.r + z__9.r, z__5.i = z__6.i + z__9.i;
2102 i__5 = j + (c_dim1 << 2);
2103 z__10.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__10.i = v4.r
2104 * c__[i__5].i + v4.i * c__[i__5].r;
2105 z__4.r = z__5.r + z__10.r, z__4.i = z__5.i + z__10.i;
2106 i__6 = j + c_dim1 * 5;
2107 z__11.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__11.i = v5.r
2108 * c__[i__6].i + v5.i * c__[i__6].r;
2109 z__3.r = z__4.r + z__11.r, z__3.i = z__4.i + z__11.i;
2110 i__7 = j + c_dim1 * 6;
2111 z__12.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__12.i = v6.r
2112 * c__[i__7].i + v6.i * c__[i__7].r;
2113 z__2.r = z__3.r + z__12.r, z__2.i = z__3.i + z__12.i;
2114 i__8 = j + c_dim1 * 7;
2115 z__13.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__13.i = v7.r
2116 * c__[i__8].i + v7.i * c__[i__8].r;
2117 z__1.r = z__2.r + z__13.r, z__1.i = z__2.i + z__13.i;
2118 sum.r = z__1.r, sum.i = z__1.i;
2121 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
2123 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2124 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2125 i__2 = j + (c_dim1 << 1);
2126 i__3 = j + (c_dim1 << 1);
2127 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
2129 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2130 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2131 i__2 = j + c_dim1 * 3;
2132 i__3 = j + c_dim1 * 3;
2133 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
2135 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2136 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2137 i__2 = j + (c_dim1 << 2);
2138 i__3 = j + (c_dim1 << 2);
2139 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
2141 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2142 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2143 i__2 = j + c_dim1 * 5;
2144 i__3 = j + c_dim1 * 5;
2145 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
2147 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2148 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2149 i__2 = j + c_dim1 * 6;
2150 i__3 = j + c_dim1 * 6;
2151 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
2153 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2154 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2155 i__2 = j + c_dim1 * 7;
2156 i__3 = j + c_dim1 * 7;
2157 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
2159 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2160 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2166 /* Special code for 8 x 8 Householder */
2168 v1.r = v[1].r, v1.i = v[1].i;
2170 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2172 t1.r = z__1.r, t1.i = z__1.i;
2173 v2.r = v[2].r, v2.i = v[2].i;
2175 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2177 t2.r = z__1.r, t2.i = z__1.i;
2178 v3.r = v[3].r, v3.i = v[3].i;
2180 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2182 t3.r = z__1.r, t3.i = z__1.i;
2183 v4.r = v[4].r, v4.i = v[4].i;
2185 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2187 t4.r = z__1.r, t4.i = z__1.i;
2188 v5.r = v[5].r, v5.i = v[5].i;
2190 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2192 t5.r = z__1.r, t5.i = z__1.i;
2193 v6.r = v[6].r, v6.i = v[6].i;
2195 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2197 t6.r = z__1.r, t6.i = z__1.i;
2198 v7.r = v[7].r, v7.i = v[7].i;
2200 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2202 t7.r = z__1.r, t7.i = z__1.i;
2203 v8.r = v[8].r, v8.i = v[8].i;
2205 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2207 t8.r = z__1.r, t8.i = z__1.i;
2209 for (j = 1; j <= i__1; ++j) {
2211 z__8.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__8.i = v1.r *
2212 c__[i__2].i + v1.i * c__[i__2].r;
2213 i__3 = j + (c_dim1 << 1);
2214 z__9.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__9.i = v2.r *
2215 c__[i__3].i + v2.i * c__[i__3].r;
2216 z__7.r = z__8.r + z__9.r, z__7.i = z__8.i + z__9.i;
2217 i__4 = j + c_dim1 * 3;
2218 z__10.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__10.i = v3.r
2219 * c__[i__4].i + v3.i * c__[i__4].r;
2220 z__6.r = z__7.r + z__10.r, z__6.i = z__7.i + z__10.i;
2221 i__5 = j + (c_dim1 << 2);
2222 z__11.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__11.i = v4.r
2223 * c__[i__5].i + v4.i * c__[i__5].r;
2224 z__5.r = z__6.r + z__11.r, z__5.i = z__6.i + z__11.i;
2225 i__6 = j + c_dim1 * 5;
2226 z__12.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__12.i = v5.r
2227 * c__[i__6].i + v5.i * c__[i__6].r;
2228 z__4.r = z__5.r + z__12.r, z__4.i = z__5.i + z__12.i;
2229 i__7 = j + c_dim1 * 6;
2230 z__13.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__13.i = v6.r
2231 * c__[i__7].i + v6.i * c__[i__7].r;
2232 z__3.r = z__4.r + z__13.r, z__3.i = z__4.i + z__13.i;
2233 i__8 = j + c_dim1 * 7;
2234 z__14.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__14.i = v7.r
2235 * c__[i__8].i + v7.i * c__[i__8].r;
2236 z__2.r = z__3.r + z__14.r, z__2.i = z__3.i + z__14.i;
2237 i__9 = j + (c_dim1 << 3);
2238 z__15.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__15.i = v8.r
2239 * c__[i__9].i + v8.i * c__[i__9].r;
2240 z__1.r = z__2.r + z__15.r, z__1.i = z__2.i + z__15.i;
2241 sum.r = z__1.r, sum.i = z__1.i;
2244 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
2246 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2247 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2248 i__2 = j + (c_dim1 << 1);
2249 i__3 = j + (c_dim1 << 1);
2250 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
2252 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2253 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2254 i__2 = j + c_dim1 * 3;
2255 i__3 = j + c_dim1 * 3;
2256 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
2258 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2259 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2260 i__2 = j + (c_dim1 << 2);
2261 i__3 = j + (c_dim1 << 2);
2262 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
2264 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2265 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2266 i__2 = j + c_dim1 * 5;
2267 i__3 = j + c_dim1 * 5;
2268 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
2270 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2271 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2272 i__2 = j + c_dim1 * 6;
2273 i__3 = j + c_dim1 * 6;
2274 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
2276 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2277 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2278 i__2 = j + c_dim1 * 7;
2279 i__3 = j + c_dim1 * 7;
2280 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
2282 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2283 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2284 i__2 = j + (c_dim1 << 3);
2285 i__3 = j + (c_dim1 << 3);
2286 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
2288 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2289 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2295 /* Special code for 9 x 9 Householder */
2297 v1.r = v[1].r, v1.i = v[1].i;
2299 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2301 t1.r = z__1.r, t1.i = z__1.i;
2302 v2.r = v[2].r, v2.i = v[2].i;
2304 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2306 t2.r = z__1.r, t2.i = z__1.i;
2307 v3.r = v[3].r, v3.i = v[3].i;
2309 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2311 t3.r = z__1.r, t3.i = z__1.i;
2312 v4.r = v[4].r, v4.i = v[4].i;
2314 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2316 t4.r = z__1.r, t4.i = z__1.i;
2317 v5.r = v[5].r, v5.i = v[5].i;
2319 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2321 t5.r = z__1.r, t5.i = z__1.i;
2322 v6.r = v[6].r, v6.i = v[6].i;
2324 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2326 t6.r = z__1.r, t6.i = z__1.i;
2327 v7.r = v[7].r, v7.i = v[7].i;
2329 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2331 t7.r = z__1.r, t7.i = z__1.i;
2332 v8.r = v[8].r, v8.i = v[8].i;
2334 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2336 t8.r = z__1.r, t8.i = z__1.i;
2337 v9.r = v[9].r, v9.i = v[9].i;
2339 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2341 t9.r = z__1.r, t9.i = z__1.i;
2343 for (j = 1; j <= i__1; ++j) {
2345 z__9.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__9.i = v1.r *
2346 c__[i__2].i + v1.i * c__[i__2].r;
2347 i__3 = j + (c_dim1 << 1);
2348 z__10.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__10.i = v2.r
2349 * c__[i__3].i + v2.i * c__[i__3].r;
2350 z__8.r = z__9.r + z__10.r, z__8.i = z__9.i + z__10.i;
2351 i__4 = j + c_dim1 * 3;
2352 z__11.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__11.i = v3.r
2353 * c__[i__4].i + v3.i * c__[i__4].r;
2354 z__7.r = z__8.r + z__11.r, z__7.i = z__8.i + z__11.i;
2355 i__5 = j + (c_dim1 << 2);
2356 z__12.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__12.i = v4.r
2357 * c__[i__5].i + v4.i * c__[i__5].r;
2358 z__6.r = z__7.r + z__12.r, z__6.i = z__7.i + z__12.i;
2359 i__6 = j + c_dim1 * 5;
2360 z__13.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__13.i = v5.r
2361 * c__[i__6].i + v5.i * c__[i__6].r;
2362 z__5.r = z__6.r + z__13.r, z__5.i = z__6.i + z__13.i;
2363 i__7 = j + c_dim1 * 6;
2364 z__14.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__14.i = v6.r
2365 * c__[i__7].i + v6.i * c__[i__7].r;
2366 z__4.r = z__5.r + z__14.r, z__4.i = z__5.i + z__14.i;
2367 i__8 = j + c_dim1 * 7;
2368 z__15.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__15.i = v7.r
2369 * c__[i__8].i + v7.i * c__[i__8].r;
2370 z__3.r = z__4.r + z__15.r, z__3.i = z__4.i + z__15.i;
2371 i__9 = j + (c_dim1 << 3);
2372 z__16.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__16.i = v8.r
2373 * c__[i__9].i + v8.i * c__[i__9].r;
2374 z__2.r = z__3.r + z__16.r, z__2.i = z__3.i + z__16.i;
2375 i__10 = j + c_dim1 * 9;
2376 z__17.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i, z__17.i =
2377 v9.r * c__[i__10].i + v9.i * c__[i__10].r;
2378 z__1.r = z__2.r + z__17.r, z__1.i = z__2.i + z__17.i;
2379 sum.r = z__1.r, sum.i = z__1.i;
2382 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
2384 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2385 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2386 i__2 = j + (c_dim1 << 1);
2387 i__3 = j + (c_dim1 << 1);
2388 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
2390 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2391 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2392 i__2 = j + c_dim1 * 3;
2393 i__3 = j + c_dim1 * 3;
2394 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
2396 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2397 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2398 i__2 = j + (c_dim1 << 2);
2399 i__3 = j + (c_dim1 << 2);
2400 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
2402 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2403 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2404 i__2 = j + c_dim1 * 5;
2405 i__3 = j + c_dim1 * 5;
2406 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
2408 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2409 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2410 i__2 = j + c_dim1 * 6;
2411 i__3 = j + c_dim1 * 6;
2412 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
2414 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2415 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2416 i__2 = j + c_dim1 * 7;
2417 i__3 = j + c_dim1 * 7;
2418 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
2420 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2421 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2422 i__2 = j + (c_dim1 << 3);
2423 i__3 = j + (c_dim1 << 3);
2424 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
2426 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2427 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2428 i__2 = j + c_dim1 * 9;
2429 i__3 = j + c_dim1 * 9;
2430 z__2.r = sum.r * t9.r - sum.i * t9.i, z__2.i = sum.r * t9.i +
2432 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2433 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2439 /* Special code for 10 x 10 Householder */
2441 v1.r = v[1].r, v1.i = v[1].i;
2443 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2445 t1.r = z__1.r, t1.i = z__1.i;
2446 v2.r = v[2].r, v2.i = v[2].i;
2448 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2450 t2.r = z__1.r, t2.i = z__1.i;
2451 v3.r = v[3].r, v3.i = v[3].i;
2453 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2455 t3.r = z__1.r, t3.i = z__1.i;
2456 v4.r = v[4].r, v4.i = v[4].i;
2458 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2460 t4.r = z__1.r, t4.i = z__1.i;
2461 v5.r = v[5].r, v5.i = v[5].i;
2463 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2465 t5.r = z__1.r, t5.i = z__1.i;
2466 v6.r = v[6].r, v6.i = v[6].i;
2468 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2470 t6.r = z__1.r, t6.i = z__1.i;
2471 v7.r = v[7].r, v7.i = v[7].i;
2473 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2475 t7.r = z__1.r, t7.i = z__1.i;
2476 v8.r = v[8].r, v8.i = v[8].i;
2478 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2480 t8.r = z__1.r, t8.i = z__1.i;
2481 v9.r = v[9].r, v9.i = v[9].i;
2483 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2485 t9.r = z__1.r, t9.i = z__1.i;
2486 v10.r = v[10].r, v10.i = v[10].i;
2487 d_cnjg(&z__2, &v10);
2488 z__1.r = tau->r * z__2.r - tau->i * z__2.i, z__1.i = tau->r * z__2.i
2490 t10.r = z__1.r, t10.i = z__1.i;
2492 for (j = 1; j <= i__1; ++j) {
2494 z__10.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i, z__10.i = v1.r
2495 * c__[i__2].i + v1.i * c__[i__2].r;
2496 i__3 = j + (c_dim1 << 1);
2497 z__11.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i, z__11.i = v2.r
2498 * c__[i__3].i + v2.i * c__[i__3].r;
2499 z__9.r = z__10.r + z__11.r, z__9.i = z__10.i + z__11.i;
2500 i__4 = j + c_dim1 * 3;
2501 z__12.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i, z__12.i = v3.r
2502 * c__[i__4].i + v3.i * c__[i__4].r;
2503 z__8.r = z__9.r + z__12.r, z__8.i = z__9.i + z__12.i;
2504 i__5 = j + (c_dim1 << 2);
2505 z__13.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i, z__13.i = v4.r
2506 * c__[i__5].i + v4.i * c__[i__5].r;
2507 z__7.r = z__8.r + z__13.r, z__7.i = z__8.i + z__13.i;
2508 i__6 = j + c_dim1 * 5;
2509 z__14.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i, z__14.i = v5.r
2510 * c__[i__6].i + v5.i * c__[i__6].r;
2511 z__6.r = z__7.r + z__14.r, z__6.i = z__7.i + z__14.i;
2512 i__7 = j + c_dim1 * 6;
2513 z__15.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i, z__15.i = v6.r
2514 * c__[i__7].i + v6.i * c__[i__7].r;
2515 z__5.r = z__6.r + z__15.r, z__5.i = z__6.i + z__15.i;
2516 i__8 = j + c_dim1 * 7;
2517 z__16.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i, z__16.i = v7.r
2518 * c__[i__8].i + v7.i * c__[i__8].r;
2519 z__4.r = z__5.r + z__16.r, z__4.i = z__5.i + z__16.i;
2520 i__9 = j + (c_dim1 << 3);
2521 z__17.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i, z__17.i = v8.r
2522 * c__[i__9].i + v8.i * c__[i__9].r;
2523 z__3.r = z__4.r + z__17.r, z__3.i = z__4.i + z__17.i;
2524 i__10 = j + c_dim1 * 9;
2525 z__18.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i, z__18.i =
2526 v9.r * c__[i__10].i + v9.i * c__[i__10].r;
2527 z__2.r = z__3.r + z__18.r, z__2.i = z__3.i + z__18.i;
2528 i__11 = j + c_dim1 * 10;
2529 z__19.r = v10.r * c__[i__11].r - v10.i * c__[i__11].i, z__19.i =
2530 v10.r * c__[i__11].i + v10.i * c__[i__11].r;
2531 z__1.r = z__2.r + z__19.r, z__1.i = z__2.i + z__19.i;
2532 sum.r = z__1.r, sum.i = z__1.i;
2535 z__2.r = sum.r * t1.r - sum.i * t1.i, z__2.i = sum.r * t1.i +
2537 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2538 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2539 i__2 = j + (c_dim1 << 1);
2540 i__3 = j + (c_dim1 << 1);
2541 z__2.r = sum.r * t2.r - sum.i * t2.i, z__2.i = sum.r * t2.i +
2543 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2544 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2545 i__2 = j + c_dim1 * 3;
2546 i__3 = j + c_dim1 * 3;
2547 z__2.r = sum.r * t3.r - sum.i * t3.i, z__2.i = sum.r * t3.i +
2549 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2550 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2551 i__2 = j + (c_dim1 << 2);
2552 i__3 = j + (c_dim1 << 2);
2553 z__2.r = sum.r * t4.r - sum.i * t4.i, z__2.i = sum.r * t4.i +
2555 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2556 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2557 i__2 = j + c_dim1 * 5;
2558 i__3 = j + c_dim1 * 5;
2559 z__2.r = sum.r * t5.r - sum.i * t5.i, z__2.i = sum.r * t5.i +
2561 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2562 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2563 i__2 = j + c_dim1 * 6;
2564 i__3 = j + c_dim1 * 6;
2565 z__2.r = sum.r * t6.r - sum.i * t6.i, z__2.i = sum.r * t6.i +
2567 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2568 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2569 i__2 = j + c_dim1 * 7;
2570 i__3 = j + c_dim1 * 7;
2571 z__2.r = sum.r * t7.r - sum.i * t7.i, z__2.i = sum.r * t7.i +
2573 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2574 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2575 i__2 = j + (c_dim1 << 3);
2576 i__3 = j + (c_dim1 << 3);
2577 z__2.r = sum.r * t8.r - sum.i * t8.i, z__2.i = sum.r * t8.i +
2579 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2580 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2581 i__2 = j + c_dim1 * 9;
2582 i__3 = j + c_dim1 * 9;
2583 z__2.r = sum.r * t9.r - sum.i * t9.i, z__2.i = sum.r * t9.i +
2585 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2586 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;
2587 i__2 = j + c_dim1 * 10;
2588 i__3 = j + c_dim1 * 10;
2589 z__2.r = sum.r * t10.r - sum.i * t10.i, z__2.i = sum.r * t10.i +
2591 z__1.r = c__[i__3].r - z__2.r, z__1.i = c__[i__3].i - z__2.i;
2592 c__[i__2].r = z__1.r, c__[i__2].i = z__1.i;