text: /* ../netlib/clarfx.f -- translated by f2c (version 20100827). 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 "FLA_f2c.h" /* Table of constant values */
  static integer c__1 = 1;
  /* > \brief \b CLARFX applies an elementary reflector to a general rectangular matrix, with loop unrolling whe n the reflector has order ≤ 10. */
  /* =========== DOCUMENTATION =========== */
  /* Online html documentation available at */
  /* http://www.netlib.org/lapack/explore-html/ */
  /* > \htmlonly */
  /* > Download CLARFX + dependencies */
  /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/clarfx. f"> */
  /* > [TGZ]</a> */
  /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/clarfx. f"> */
  /* > [ZIP]</a> */
  /* > <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/clarfx. f"> */
  /* > [TXT]</a> */
  /* > \endhtmlonly */
  /* Definition: */
  /* =========== */
  /* SUBROUTINE CLARFX( SIDE, M, N, V, TAU, C, LDC, WORK ) */
  /* .. Scalar Arguments .. */
  /* CHARACTER SIDE */
  /* INTEGER LDC, M, N */
  /* COMPLEX TAU */
  /* .. */
  /* .. Array Arguments .. */
  /* COMPLEX C( LDC, * ), V( * ), WORK( * ) */
  /* .. */
  /* > \par Purpose: */
  /* ============= */
  /* > */
  /* > \verbatim */
  /* > */
  /* > CLARFX applies a complex elementary reflector H to a complex m by n */
  /* > matrix C, from either the left or the right. H is represented in the */
  /* > form */
  /* > */
  /* > H = I - tau * v * v**H */
  /* > */
  /* > where tau is a complex scalar and v is a complex vector. */
  /* > */
  /* > If tau = 0, then H is taken to be the unit matrix */
  /* > */
  /* > This version uses inline code if H has order < 11. */
  /* > \endverbatim */
  /* Arguments: */
  /* ========== */
  /* > \param[in] SIDE */
  /* > \verbatim */
  /* > SIDE is CHARACTER*1 */
  /* > = 'L': form H * C */
  /* > = 'R': form C * H */
  /* > \endverbatim */
  /* > */
  /* > \param[in] M */
  /* > \verbatim */
  /* > M is INTEGER */
  /* > The number of rows of the matrix C. */
  /* > \endverbatim */
  /* > */
  /* > \param[in] N */
  /* > \verbatim */
  /* > N is INTEGER */
  /* > The number of columns of the matrix C. */
  /* > \endverbatim */
  /* > */
  /* > \param[in] V */
  /* > \verbatim */
  /* > V is COMPLEX array, dimension (M) if SIDE = 'L' */
  /* > or (N) if SIDE = 'R' */
  /* > The vector v in the representation of H. */
  /* > \endverbatim */
  /* > */
  /* > \param[in] TAU */
  /* > \verbatim */
  /* > TAU is COMPLEX */
  /* > The value tau in the representation of H. */
  /* > \endverbatim */
  /* > */
  /* > \param[in,out] C */
  /* > \verbatim */
  /* > C is COMPLEX array, dimension (LDC,N) */
  /* > On entry, the m by n matrix C. */
  /* > On exit, C is overwritten by the matrix H * C if SIDE = 'L', */
  /* > or C * H if SIDE = 'R'. */
  /* > \endverbatim */
  /* > */
  /* > \param[in] LDC */
  /* > \verbatim */
  /* > LDC is INTEGER */
  /* > The leading dimension of the array C. LDA >= max(1,M). */
  /* > \endverbatim */
  /* > */
  /* > \param[out] WORK */
  /* > \verbatim */
  /* > WORK is COMPLEX array, dimension (N) if SIDE = 'L' */
  /* > or (M) if SIDE = 'R' */
  /* > WORK is not referenced if H has order < 11. */
  /* > \endverbatim */
  /* Authors: */
  /* ======== */
  /* > \author Univ. of Tennessee */
  /* > \author Univ. of California Berkeley */
  /* > \author Univ. of Colorado Denver */
  /* > \author NAG Ltd. */
  /* > \date September 2012 */
  /* > \ingroup complexOTHERauxiliary */
  /* ===================================================================== */
  /* Subroutine */
  int clarfx_(char *side, integer *m, integer *n, complex *v, complex *tau, complex *c__, integer *ldc, complex *work)
  {
      AOCL_DTL_TRACE_ENTRY(AOCL_DTL_LEVEL_TRACE_5);
  #if AOCL_DTL_LOG_ENABLE 
      char buffer[256]; 
  #if FLA_ENABLE_ILP64 
      snprintf(buffer, 256,"clarfx inputs: side %c, m %lld, n %lld, ldc %lld",*side, *m, *n, *ldc);
  #else 
      snprintf(buffer, 256,"clarfx inputs: side %c, m %d, n %d, ldc %d",*side, *m, *n, *ldc);
  #endif
      AOCL_DTL_LOG(AOCL_DTL_LEVEL_TRACE_5, buffer);
  #endif
      /* System generated locals */
      integer c_dim1, c_offset, i__1, i__2, i__3, i__4, i__5, i__6, i__7, i__8, i__9, i__10, i__11;
      complex q__1, q__2, q__3, q__4, q__5, q__6, q__7, q__8, q__9, q__10, q__11, q__12, q__13, q__14, q__15, q__16, q__17, q__18, q__19;
      /* Builtin functions */
      void r_cnjg(complex *, complex *);
      /* Local variables */
      integer j;
      complex t1, t2, t3, t4, t5, t6, t7, t8, t9, v1, v2, v3, v4, v5, v6, v7, v8, v9, t10, v10, sum;
      extern /* Subroutine */
      int clarf_(char *, integer *, integer *, complex * , integer *, complex *, complex *, integer *, complex *);
      extern logical lsame_(char *, char *);
      /* -- LAPACK auxiliary routine (version 3.4.2) -- */
      /* -- LAPACK is a software package provided by Univ. of Tennessee, -- */
      /* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- */
      /* September 2012 */
      /* .. Scalar Arguments .. */
      /* .. */
      /* .. Array Arguments .. */
      /* .. */
      /* ===================================================================== */
      /* .. Parameters .. */
      /* .. */
      /* .. Local Scalars .. */
      /* .. */
      /* .. External Functions .. */
      /* .. */
      /* .. External Subroutines .. */
      /* .. */
      /* .. Intrinsic Functions .. */
      /* .. */
      /* .. Executable Statements .. */
      /* Parameter adjustments */
      --v;
      c_dim1 = *ldc;
      c_offset = 1 + c_dim1;
      c__ -= c_offset;
      --work;
      /* Function Body */
      if (tau->r == 0.f && tau->i == 0.f)
      {
          AOCL_DTL_TRACE_EXIT(AOCL_DTL_LEVEL_TRACE_5);
          return 0;
      }
      if (lsame_(side, "L"))
      {
          /* Form H * C, where H has order m. */
          switch (*m)
          {
          case 1:
              goto L10;
          case 2:
              goto L30;
          case 3:
              goto L50;
          case 4:
              goto L70;
          case 5:
              goto L90;
          case 6:
              goto L110;
          case 7:
              goto L130;
          case 8:
              goto L150;
          case 9:
              goto L170;
          case 10:
              goto L190;
          }
          /* Code for general M */
          clarf_(side, m, n, &v[1], &c__1, tau, &c__[c_offset], ldc, &work[1]);
          goto L410;
  L10: /* Special code for 1 x 1 Householder */
          q__3.r = tau->r * v[1].r - tau->i * v[1].i;
          q__3.i = tau->r * v[1].i + tau->i * v[1].r; // , expr subst
          r_cnjg(&q__4, &v[1]);
          q__2.r = q__3.r * q__4.r - q__3.i * q__4.i;
          q__2.i = q__3.r * q__4.i + q__3.i * q__4.r; // , expr subst
          q__1.r = 1.f - q__2.r;
          q__1.i = 0.f - q__2.i; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__1.r = t1.r * c__[i__3].r - t1.i * c__[i__3].i;
              q__1.i = t1.r * c__[i__3].i + t1.i * c__[i__3].r; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L20: */
          }
          goto L410;
  L30: /* Special code for 2 x 2 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__2.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__2.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__3.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__3.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__1.r = q__2.r + q__3.r;
              q__1.i = q__2.i + q__3.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L40: */
          }
          goto L410;
  L50: /* Special code for 3 x 3 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__3.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__3.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__4.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__4.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__2.r = q__3.r + q__4.r;
              q__2.i = q__3.i + q__4.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__5.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__5.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__1.r = q__2.r + q__5.r;
              q__1.i = q__2.i + q__5.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L60: */
          }
          goto L410;
  L70: /* Special code for 4 x 4 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__4.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__4.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__5.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__5.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__3.r = q__4.r + q__5.r;
              q__3.i = q__4.i + q__5.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__6.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__6.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__2.r = q__3.r + q__6.r;
              q__2.i = q__3.i + q__6.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__7.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__7.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__1.r = q__2.r + q__7.r;
              q__1.i = q__2.i + q__7.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L80: */
          }
          goto L410;
  L90: /* Special code for 5 x 5 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__5.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__5.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__6.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__6.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__4.r = q__5.r + q__6.r;
              q__4.i = q__5.i + q__6.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__7.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__7.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__3.r = q__4.r + q__7.r;
              q__3.i = q__4.i + q__7.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__8.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__8.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__2.r = q__3.r + q__8.r;
              q__2.i = q__3.i + q__8.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__9.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__9.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__1.r = q__2.r + q__9.r;
              q__1.i = q__2.i + q__9.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L100: */
          }
          goto L410;
  L110: /* Special code for 6 x 6 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[6]);
          v6.r = q__1.r;
          v6.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__6.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__6.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__7.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__7.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__5.r = q__6.r + q__7.r;
              q__5.i = q__6.i + q__7.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__8.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__8.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__4.r = q__5.r + q__8.r;
              q__4.i = q__5.i + q__8.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__9.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__9.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__3.r = q__4.r + q__9.r;
              q__3.i = q__4.i + q__9.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__10.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__10.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__2.r = q__3.r + q__10.r;
              q__2.i = q__3.i + q__10.i; // , expr subst
              i__7 = j * c_dim1 + 6;
              q__11.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__11.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__1.r = q__2.r + q__11.r;
              q__1.i = q__2.i + q__11.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 6;
              i__3 = j * c_dim1 + 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L120: */
          }
          goto L410;
  L130: /* Special code for 7 x 7 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[6]);
          v6.r = q__1.r;
          v6.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[7]);
          v7.r = q__1.r;
          v7.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__7.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__7.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__8.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__8.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__6.r = q__7.r + q__8.r;
              q__6.i = q__7.i + q__8.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__9.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__9.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__5.r = q__6.r + q__9.r;
              q__5.i = q__6.i + q__9.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__10.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__10.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__4.r = q__5.r + q__10.r;
              q__4.i = q__5.i + q__10.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__11.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__11.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__3.r = q__4.r + q__11.r;
              q__3.i = q__4.i + q__11.i; // , expr subst
              i__7 = j * c_dim1 + 6;
              q__12.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__12.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__2.r = q__3.r + q__12.r;
              q__2.i = q__3.i + q__12.i; // , expr subst
              i__8 = j * c_dim1 + 7;
              q__13.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__13.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__1.r = q__2.r + q__13.r;
              q__1.i = q__2.i + q__13.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 6;
              i__3 = j * c_dim1 + 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 7;
              i__3 = j * c_dim1 + 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L140: */
          }
          goto L410;
  L150: /* Special code for 8 x 8 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[6]);
          v6.r = q__1.r;
          v6.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[7]);
          v7.r = q__1.r;
          v7.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[8]);
          v8.r = q__1.r;
          v8.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__8.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__8.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__9.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__9.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__7.r = q__8.r + q__9.r;
              q__7.i = q__8.i + q__9.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__10.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__10.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__6.r = q__7.r + q__10.r;
              q__6.i = q__7.i + q__10.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__11.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__11.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__5.r = q__6.r + q__11.r;
              q__5.i = q__6.i + q__11.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__12.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__12.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__4.r = q__5.r + q__12.r;
              q__4.i = q__5.i + q__12.i; // , expr subst
              i__7 = j * c_dim1 + 6;
              q__13.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__13.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__3.r = q__4.r + q__13.r;
              q__3.i = q__4.i + q__13.i; // , expr subst
              i__8 = j * c_dim1 + 7;
              q__14.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__14.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__2.r = q__3.r + q__14.r;
              q__2.i = q__3.i + q__14.i; // , expr subst
              i__9 = j * c_dim1 + 8;
              q__15.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__15.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__1.r = q__2.r + q__15.r;
              q__1.i = q__2.i + q__15.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 6;
              i__3 = j * c_dim1 + 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 7;
              i__3 = j * c_dim1 + 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 8;
              i__3 = j * c_dim1 + 8;
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L160: */
          }
          goto L410;
  L170: /* Special code for 9 x 9 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[6]);
          v6.r = q__1.r;
          v6.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[7]);
          v7.r = q__1.r;
          v7.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[8]);
          v8.r = q__1.r;
          v8.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[9]);
          v9.r = q__1.r;
          v9.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v9);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t9.r = q__1.r;
          t9.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__9.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__9.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__10.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__10.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__8.r = q__9.r + q__10.r;
              q__8.i = q__9.i + q__10.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__11.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__11.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__7.r = q__8.r + q__11.r;
              q__7.i = q__8.i + q__11.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__12.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__12.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__6.r = q__7.r + q__12.r;
              q__6.i = q__7.i + q__12.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__13.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__13.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__5.r = q__6.r + q__13.r;
              q__5.i = q__6.i + q__13.i; // , expr subst
              i__7 = j * c_dim1 + 6;
              q__14.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__14.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__4.r = q__5.r + q__14.r;
              q__4.i = q__5.i + q__14.i; // , expr subst
              i__8 = j * c_dim1 + 7;
              q__15.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__15.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__3.r = q__4.r + q__15.r;
              q__3.i = q__4.i + q__15.i; // , expr subst
              i__9 = j * c_dim1 + 8;
              q__16.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__16.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__2.r = q__3.r + q__16.r;
              q__2.i = q__3.i + q__16.i; // , expr subst
              i__10 = j * c_dim1 + 9;
              q__17.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i;
              q__17.i = v9.r * c__[i__10].i + v9.i * c__[i__10].r; // , expr subst
              q__1.r = q__2.r + q__17.r;
              q__1.i = q__2.i + q__17.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 6;
              i__3 = j * c_dim1 + 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 7;
              i__3 = j * c_dim1 + 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 8;
              i__3 = j * c_dim1 + 8;
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 9;
              i__3 = j * c_dim1 + 9;
              q__2.r = sum.r * t9.r - sum.i * t9.i;
              q__2.i = sum.r * t9.i + sum.i * t9.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L180: */
          }
          goto L410;
  L190: /* Special code for 10 x 10 Householder */
          r_cnjg(&q__1, &v[1]);
          v1.r = q__1.r;
          v1.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[2]);
          v2.r = q__1.r;
          v2.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[3]);
          v3.r = q__1.r;
          v3.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[4]);
          v4.r = q__1.r;
          v4.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[5]);
          v5.r = q__1.r;
          v5.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[6]);
          v6.r = q__1.r;
          v6.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[7]);
          v7.r = q__1.r;
          v7.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[8]);
          v8.r = q__1.r;
          v8.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[9]);
          v9.r = q__1.r;
          v9.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v9);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t9.r = q__1.r;
          t9.i = q__1.i; // , expr subst
          r_cnjg(&q__1, &v[10]);
          v10.r = q__1.r;
          v10.i = q__1.i; // , expr subst
          r_cnjg(&q__2, &v10);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t10.r = q__1.r;
          t10.i = q__1.i; // , expr subst
          i__1 = *n;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j * c_dim1 + 1;
              q__10.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__10.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j * c_dim1 + 2;
              q__11.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__11.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__9.r = q__10.r + q__11.r;
              q__9.i = q__10.i + q__11.i; // , expr subst
              i__4 = j * c_dim1 + 3;
              q__12.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__12.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__8.r = q__9.r + q__12.r;
              q__8.i = q__9.i + q__12.i; // , expr subst
              i__5 = j * c_dim1 + 4;
              q__13.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__13.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__7.r = q__8.r + q__13.r;
              q__7.i = q__8.i + q__13.i; // , expr subst
              i__6 = j * c_dim1 + 5;
              q__14.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__14.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__6.r = q__7.r + q__14.r;
              q__6.i = q__7.i + q__14.i; // , expr subst
              i__7 = j * c_dim1 + 6;
              q__15.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__15.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__5.r = q__6.r + q__15.r;
              q__5.i = q__6.i + q__15.i; // , expr subst
              i__8 = j * c_dim1 + 7;
              q__16.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__16.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__4.r = q__5.r + q__16.r;
              q__4.i = q__5.i + q__16.i; // , expr subst
              i__9 = j * c_dim1 + 8;
              q__17.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__17.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__3.r = q__4.r + q__17.r;
              q__3.i = q__4.i + q__17.i; // , expr subst
              i__10 = j * c_dim1 + 9;
              q__18.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i;
              q__18.i = v9.r * c__[i__10].i + v9.i * c__[i__10].r; // , expr subst
              q__2.r = q__3.r + q__18.r;
              q__2.i = q__3.i + q__18.i; // , expr subst
              i__11 = j * c_dim1 + 10;
              q__19.r = v10.r * c__[i__11].r - v10.i * c__[i__11].i;
              q__19.i = v10.r * c__[i__11].i + v10.i * c__[i__11].r; // , expr subst
              q__1.r = q__2.r + q__19.r;
              q__1.i = q__2.i + q__19.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 1;
              i__3 = j * c_dim1 + 1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 2;
              i__3 = j * c_dim1 + 2;
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 3;
              i__3 = j * c_dim1 + 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 4;
              i__3 = j * c_dim1 + 4;
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 5;
              i__3 = j * c_dim1 + 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 6;
              i__3 = j * c_dim1 + 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 7;
              i__3 = j * c_dim1 + 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 8;
              i__3 = j * c_dim1 + 8;
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 9;
              i__3 = j * c_dim1 + 9;
              q__2.r = sum.r * t9.r - sum.i * t9.i;
              q__2.i = sum.r * t9.i + sum.i * t9.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j * c_dim1 + 10;
              i__3 = j * c_dim1 + 10;
              q__2.r = sum.r * t10.r - sum.i * t10.i;
              q__2.i = sum.r * t10.i + sum.i * t10.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L200: */
          }
          goto L410;
      }
      else
      {
          /* Form C * H, where H has order n. */
          switch (*n)
          {
          case 1:
              goto L210;
          case 2:
              goto L230;
          case 3:
              goto L250;
          case 4:
              goto L270;
          case 5:
              goto L290;
          case 6:
              goto L310;
          case 7:
              goto L330;
          case 8:
              goto L350;
          case 9:
              goto L370;
          case 10:
              goto L390;
          }
          /* Code for general N */
          clarf_(side, m, n, &v[1], &c__1, tau, &c__[c_offset], ldc, &work[1]);
          goto L410;
  L210: /* Special code for 1 x 1 Householder */
          q__3.r = tau->r * v[1].r - tau->i * v[1].i;
          q__3.i = tau->r * v[1].i + tau->i * v[1].r; // , expr subst
          r_cnjg(&q__4, &v[1]);
          q__2.r = q__3.r * q__4.r - q__3.i * q__4.i;
          q__2.i = q__3.r * q__4.i + q__3.i * q__4.r; // , expr subst
          q__1.r = 1.f - q__2.r;
          q__1.i = 0.f - q__2.i; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__1.r = t1.r * c__[i__3].r - t1.i * c__[i__3].i;
              q__1.i = t1.r * c__[i__3].i + t1.i * c__[i__3].r; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L220: */
          }
          goto L410;
  L230: /* Special code for 2 x 2 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__2.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__2.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__3.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__3.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__1.r = q__2.r + q__3.r;
              q__1.i = q__2.i + q__3.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L240: */
          }
          goto L410;
  L250: /* Special code for 3 x 3 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__3.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__3.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__4.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__4.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__2.r = q__3.r + q__4.r;
              q__2.i = q__3.i + q__4.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__5.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__5.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__1.r = q__2.r + q__5.r;
              q__1.i = q__2.i + q__5.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L260: */
          }
          goto L410;
  L270: /* Special code for 4 x 4 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__4.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__4.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__5.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__5.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__3.r = q__4.r + q__5.r;
              q__3.i = q__4.i + q__5.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__6.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__6.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__2.r = q__3.r + q__6.r;
              q__2.i = q__3.i + q__6.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__7.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__7.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__1.r = q__2.r + q__7.r;
              q__1.i = q__2.i + q__7.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L280: */
          }
          goto L410;
  L290: /* Special code for 5 x 5 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__5.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__5.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__6.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__6.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__4.r = q__5.r + q__6.r;
              q__4.i = q__5.i + q__6.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__7.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__7.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__3.r = q__4.r + q__7.r;
              q__3.i = q__4.i + q__7.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__8.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__8.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__2.r = q__3.r + q__8.r;
              q__2.i = q__3.i + q__8.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__9.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__9.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__1.r = q__2.r + q__9.r;
              q__1.i = q__2.i + q__9.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L300: */
          }
          goto L410;
  L310: /* Special code for 6 x 6 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          v6.r = v[6].r;
          v6.i = v[6].i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__6.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__6.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__7.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__7.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__5.r = q__6.r + q__7.r;
              q__5.i = q__6.i + q__7.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__8.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__8.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__4.r = q__5.r + q__8.r;
              q__4.i = q__5.i + q__8.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__9.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__9.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__3.r = q__4.r + q__9.r;
              q__3.i = q__4.i + q__9.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__10.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__10.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__2.r = q__3.r + q__10.r;
              q__2.i = q__3.i + q__10.i; // , expr subst
              i__7 = j + c_dim1 * 6;
              q__11.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__11.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__1.r = q__2.r + q__11.r;
              q__1.i = q__2.i + q__11.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 6;
              i__3 = j + c_dim1 * 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L320: */
          }
          goto L410;
  L330: /* Special code for 7 x 7 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          v6.r = v[6].r;
          v6.i = v[6].i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          v7.r = v[7].r;
          v7.i = v[7].i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__7.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__7.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__8.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__8.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__6.r = q__7.r + q__8.r;
              q__6.i = q__7.i + q__8.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__9.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__9.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__5.r = q__6.r + q__9.r;
              q__5.i = q__6.i + q__9.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__10.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__10.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__4.r = q__5.r + q__10.r;
              q__4.i = q__5.i + q__10.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__11.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__11.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__3.r = q__4.r + q__11.r;
              q__3.i = q__4.i + q__11.i; // , expr subst
              i__7 = j + c_dim1 * 6;
              q__12.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__12.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__2.r = q__3.r + q__12.r;
              q__2.i = q__3.i + q__12.i; // , expr subst
              i__8 = j + c_dim1 * 7;
              q__13.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__13.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__1.r = q__2.r + q__13.r;
              q__1.i = q__2.i + q__13.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 6;
              i__3 = j + c_dim1 * 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 7;
              i__3 = j + c_dim1 * 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L340: */
          }
          goto L410;
  L350: /* Special code for 8 x 8 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          v6.r = v[6].r;
          v6.i = v[6].i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          v7.r = v[7].r;
          v7.i = v[7].i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          v8.r = v[8].r;
          v8.i = v[8].i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__8.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__8.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__9.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__9.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__7.r = q__8.r + q__9.r;
              q__7.i = q__8.i + q__9.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__10.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__10.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__6.r = q__7.r + q__10.r;
              q__6.i = q__7.i + q__10.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__11.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__11.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__5.r = q__6.r + q__11.r;
              q__5.i = q__6.i + q__11.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__12.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__12.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__4.r = q__5.r + q__12.r;
              q__4.i = q__5.i + q__12.i; // , expr subst
              i__7 = j + c_dim1 * 6;
              q__13.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__13.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__3.r = q__4.r + q__13.r;
              q__3.i = q__4.i + q__13.i; // , expr subst
              i__8 = j + c_dim1 * 7;
              q__14.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__14.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__2.r = q__3.r + q__14.r;
              q__2.i = q__3.i + q__14.i; // , expr subst
              i__9 = j + (c_dim1 << 3);
              q__15.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__15.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__1.r = q__2.r + q__15.r;
              q__1.i = q__2.i + q__15.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 6;
              i__3 = j + c_dim1 * 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 7;
              i__3 = j + c_dim1 * 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 3);
              i__3 = j + (c_dim1 << 3);
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L360: */
          }
          goto L410;
  L370: /* Special code for 9 x 9 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          v6.r = v[6].r;
          v6.i = v[6].i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          v7.r = v[7].r;
          v7.i = v[7].i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          v8.r = v[8].r;
          v8.i = v[8].i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          v9.r = v[9].r;
          v9.i = v[9].i; // , expr subst
          r_cnjg(&q__2, &v9);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t9.r = q__1.r;
          t9.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__9.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__9.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__10.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__10.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__8.r = q__9.r + q__10.r;
              q__8.i = q__9.i + q__10.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__11.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__11.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__7.r = q__8.r + q__11.r;
              q__7.i = q__8.i + q__11.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__12.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__12.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__6.r = q__7.r + q__12.r;
              q__6.i = q__7.i + q__12.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__13.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__13.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__5.r = q__6.r + q__13.r;
              q__5.i = q__6.i + q__13.i; // , expr subst
              i__7 = j + c_dim1 * 6;
              q__14.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__14.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__4.r = q__5.r + q__14.r;
              q__4.i = q__5.i + q__14.i; // , expr subst
              i__8 = j + c_dim1 * 7;
              q__15.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__15.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__3.r = q__4.r + q__15.r;
              q__3.i = q__4.i + q__15.i; // , expr subst
              i__9 = j + (c_dim1 << 3);
              q__16.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__16.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__2.r = q__3.r + q__16.r;
              q__2.i = q__3.i + q__16.i; // , expr subst
              i__10 = j + c_dim1 * 9;
              q__17.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i;
              q__17.i = v9.r * c__[i__10].i + v9.i * c__[i__10].r; // , expr subst
              q__1.r = q__2.r + q__17.r;
              q__1.i = q__2.i + q__17.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 6;
              i__3 = j + c_dim1 * 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 7;
              i__3 = j + c_dim1 * 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 3);
              i__3 = j + (c_dim1 << 3);
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 9;
              i__3 = j + c_dim1 * 9;
              q__2.r = sum.r * t9.r - sum.i * t9.i;
              q__2.i = sum.r * t9.i + sum.i * t9.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L380: */
          }
          goto L410;
  L390: /* Special code for 10 x 10 Householder */
          v1.r = v[1].r;
          v1.i = v[1].i; // , expr subst
          r_cnjg(&q__2, &v1);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t1.r = q__1.r;
          t1.i = q__1.i; // , expr subst
          v2.r = v[2].r;
          v2.i = v[2].i; // , expr subst
          r_cnjg(&q__2, &v2);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t2.r = q__1.r;
          t2.i = q__1.i; // , expr subst
          v3.r = v[3].r;
          v3.i = v[3].i; // , expr subst
          r_cnjg(&q__2, &v3);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t3.r = q__1.r;
          t3.i = q__1.i; // , expr subst
          v4.r = v[4].r;
          v4.i = v[4].i; // , expr subst
          r_cnjg(&q__2, &v4);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t4.r = q__1.r;
          t4.i = q__1.i; // , expr subst
          v5.r = v[5].r;
          v5.i = v[5].i; // , expr subst
          r_cnjg(&q__2, &v5);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t5.r = q__1.r;
          t5.i = q__1.i; // , expr subst
          v6.r = v[6].r;
          v6.i = v[6].i; // , expr subst
          r_cnjg(&q__2, &v6);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t6.r = q__1.r;
          t6.i = q__1.i; // , expr subst
          v7.r = v[7].r;
          v7.i = v[7].i; // , expr subst
          r_cnjg(&q__2, &v7);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t7.r = q__1.r;
          t7.i = q__1.i; // , expr subst
          v8.r = v[8].r;
          v8.i = v[8].i; // , expr subst
          r_cnjg(&q__2, &v8);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t8.r = q__1.r;
          t8.i = q__1.i; // , expr subst
          v9.r = v[9].r;
          v9.i = v[9].i; // , expr subst
          r_cnjg(&q__2, &v9);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t9.r = q__1.r;
          t9.i = q__1.i; // , expr subst
          v10.r = v[10].r;
          v10.i = v[10].i; // , expr subst
          r_cnjg(&q__2, &v10);
          q__1.r = tau->r * q__2.r - tau->i * q__2.i;
          q__1.i = tau->r * q__2.i + tau->i * q__2.r; // , expr subst
          t10.r = q__1.r;
          t10.i = q__1.i; // , expr subst
          i__1 = *m;
          for (j = 1;
                  j <= i__1;
                  ++j)
          {
              i__2 = j + c_dim1;
              q__10.r = v1.r * c__[i__2].r - v1.i * c__[i__2].i;
              q__10.i = v1.r * c__[i__2].i + v1.i * c__[i__2].r; // , expr subst
              i__3 = j + (c_dim1 << 1);
              q__11.r = v2.r * c__[i__3].r - v2.i * c__[i__3].i;
              q__11.i = v2.r * c__[i__3].i + v2.i * c__[i__3].r; // , expr subst
              q__9.r = q__10.r + q__11.r;
              q__9.i = q__10.i + q__11.i; // , expr subst
              i__4 = j + c_dim1 * 3;
              q__12.r = v3.r * c__[i__4].r - v3.i * c__[i__4].i;
              q__12.i = v3.r * c__[i__4].i + v3.i * c__[i__4].r; // , expr subst
              q__8.r = q__9.r + q__12.r;
              q__8.i = q__9.i + q__12.i; // , expr subst
              i__5 = j + (c_dim1 << 2);
              q__13.r = v4.r * c__[i__5].r - v4.i * c__[i__5].i;
              q__13.i = v4.r * c__[i__5].i + v4.i * c__[i__5].r; // , expr subst
              q__7.r = q__8.r + q__13.r;
              q__7.i = q__8.i + q__13.i; // , expr subst
              i__6 = j + c_dim1 * 5;
              q__14.r = v5.r * c__[i__6].r - v5.i * c__[i__6].i;
              q__14.i = v5.r * c__[i__6].i + v5.i * c__[i__6].r; // , expr subst
              q__6.r = q__7.r + q__14.r;
              q__6.i = q__7.i + q__14.i; // , expr subst
              i__7 = j + c_dim1 * 6;
              q__15.r = v6.r * c__[i__7].r - v6.i * c__[i__7].i;
              q__15.i = v6.r * c__[i__7].i + v6.i * c__[i__7].r; // , expr subst
              q__5.r = q__6.r + q__15.r;
              q__5.i = q__6.i + q__15.i; // , expr subst
              i__8 = j + c_dim1 * 7;
              q__16.r = v7.r * c__[i__8].r - v7.i * c__[i__8].i;
              q__16.i = v7.r * c__[i__8].i + v7.i * c__[i__8].r; // , expr subst
              q__4.r = q__5.r + q__16.r;
              q__4.i = q__5.i + q__16.i; // , expr subst
              i__9 = j + (c_dim1 << 3);
              q__17.r = v8.r * c__[i__9].r - v8.i * c__[i__9].i;
              q__17.i = v8.r * c__[i__9].i + v8.i * c__[i__9].r; // , expr subst
              q__3.r = q__4.r + q__17.r;
              q__3.i = q__4.i + q__17.i; // , expr subst
              i__10 = j + c_dim1 * 9;
              q__18.r = v9.r * c__[i__10].r - v9.i * c__[i__10].i;
              q__18.i = v9.r * c__[i__10].i + v9.i * c__[i__10].r; // , expr subst
              q__2.r = q__3.r + q__18.r;
              q__2.i = q__3.i + q__18.i; // , expr subst
              i__11 = j + c_dim1 * 10;
              q__19.r = v10.r * c__[i__11].r - v10.i * c__[i__11].i;
              q__19.i = v10.r * c__[i__11].i + v10.i * c__[i__11].r; // , expr subst
              q__1.r = q__2.r + q__19.r;
              q__1.i = q__2.i + q__19.i; // , expr subst
              sum.r = q__1.r;
              sum.i = q__1.i; // , expr subst
              i__2 = j + c_dim1;
              i__3 = j + c_dim1;
              q__2.r = sum.r * t1.r - sum.i * t1.i;
              q__2.i = sum.r * t1.i + sum.i * t1.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 1);
              i__3 = j + (c_dim1 << 1);
              q__2.r = sum.r * t2.r - sum.i * t2.i;
              q__2.i = sum.r * t2.i + sum.i * t2.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 3;
              i__3 = j + c_dim1 * 3;
              q__2.r = sum.r * t3.r - sum.i * t3.i;
              q__2.i = sum.r * t3.i + sum.i * t3.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 2);
              i__3 = j + (c_dim1 << 2);
              q__2.r = sum.r * t4.r - sum.i * t4.i;
              q__2.i = sum.r * t4.i + sum.i * t4.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 5;
              i__3 = j + c_dim1 * 5;
              q__2.r = sum.r * t5.r - sum.i * t5.i;
              q__2.i = sum.r * t5.i + sum.i * t5.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 6;
              i__3 = j + c_dim1 * 6;
              q__2.r = sum.r * t6.r - sum.i * t6.i;
              q__2.i = sum.r * t6.i + sum.i * t6.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 7;
              i__3 = j + c_dim1 * 7;
              q__2.r = sum.r * t7.r - sum.i * t7.i;
              q__2.i = sum.r * t7.i + sum.i * t7.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + (c_dim1 << 3);
              i__3 = j + (c_dim1 << 3);
              q__2.r = sum.r * t8.r - sum.i * t8.i;
              q__2.i = sum.r * t8.i + sum.i * t8.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 9;
              i__3 = j + c_dim1 * 9;
              q__2.r = sum.r * t9.r - sum.i * t9.i;
              q__2.i = sum.r * t9.i + sum.i * t9.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              i__2 = j + c_dim1 * 10;
              i__3 = j + c_dim1 * 10;
              q__2.r = sum.r * t10.r - sum.i * t10.i;
              q__2.i = sum.r * t10.i + sum.i * t10.r; // , expr subst
              q__1.r = c__[i__3].r - q__2.r;
              q__1.i = c__[i__3].i - q__2.i; // , expr subst
              c__[i__2].r = q__1.r;
              c__[i__2].i = q__1.i; // , expr subst
              /* L400: */
          }
          goto L410;
      }
  L410:
      AOCL_DTL_TRACE_EXIT(AOCL_DTL_LEVEL_TRACE_5);
      return 0;
      /* End of CLARFX */
  }
  /* clarfx_ */
  