{
integer a_dim1, a_offset, i__1, i__2, i__3, i__4;
integer i__, j, l, ib, nb, ki, kk, nx, iws, nbmin, iinfo;
extern
int cungl2_fla(integer *, integer *, integer *, complex *, integer *, complex *, complex *, integer *), clarfb_( char *, char *, char *, char *, integer *, integer *, integer *, complex *, integer *, complex *, integer *, complex *, integer *, complex *, integer *), clarft_( char *, char *, integer *, integer *, complex *, integer *, complex *, complex *, integer *), xerbla_(char *, integer *);
extern integer ilaenv_(integer *, char *, char *, integer *, integer *, integer *, integer *);
integer ldwork, lwkopt;
logical lquery;
a_dim1 = *lda;
a_offset = 1 + a_dim1;
a -= a_offset;
--tau;
--work;
*info = 0;
nb = ilaenv_(&c__1, "CUNGLQ", " ", m, n, k, &c_n1);
lwkopt = max(1,*m) * nb;
work[1].r = (real) lwkopt;
work[1].i = 0.f;
lquery = *lwork == -1;
if (*m < 0)
{
*info = -1;
}
else if (*n < *m)
{
*info = -2;
}
else if (*k < 0 || *k > *m)
{
*info = -3;
}
else if (*lda < max(1,*m))
{
*info = -5;
}
else if (*lwork < max(1,*m) && ! lquery)
{
*info = -8;
}
if (*info != 0)
{
i__1 = -(*info);
xerbla_("CUNGLQ", &i__1);
return 0;
}
else if (lquery)
{
return 0;
}
if (*m <= 0)
{
work[1].r = 1.f;
work[1].i = 0.f;
return 0;
}
nbmin = 2;
nx = 0;
iws = *m;
if (nb > 1 && nb < *k)
{
i__1 = 0;
i__2 = ilaenv_(&c__3, "CUNGLQ", " ", m, n, k, &c_n1);
nx = max(i__1,i__2);
if (nx < *k)
{
ldwork = *m;
iws = ldwork * nb;
if (*lwork < iws)
{
nb = *lwork / ldwork;
i__1 = 2;
i__2 = ilaenv_(&c__2, "CUNGLQ", " ", m, n, k, &c_n1);
nbmin = max(i__1,i__2);
}
}
}
if (nb >= nbmin && nb < *k && nx < *k)
{
ki = (*k - nx - 1) / nb * nb;
i__1 = *k;
i__2 = ki + nb;
kk = min(i__1,i__2);
i__1 = kk;
for (j = 1;
j <= i__1;
++j)
{
i__2 = *m;
for (i__ = kk + 1;
i__ <= i__2;
++i__)
{
i__3 = i__ + j * a_dim1;
a[i__3].r = 0.f;
a[i__3].i = 0.f;
}
}
}
else
{
kk = 0;
}
if (kk < *m)
{
i__1 = *m - kk;
i__2 = *n - kk;
i__3 = *k - kk;
cungl2_fla(&i__1, &i__2, &i__3, &a[kk + 1 + (kk + 1) * a_dim1], lda, & tau[kk + 1], &work[1], &iinfo);
}
if (kk > 0)
{
i__1 = -nb;
for (i__ = ki + 1;
i__1 < 0 ? i__ >= 1 : i__ <= 1;
i__ += i__1)
{
i__2 = nb;
i__3 = *k - i__ + 1;
ib = min(i__2,i__3);
if (i__ + ib <= *m)
{
i__2 = *n - i__ + 1;
clarft_("Forward", "Rowwise", &i__2, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], &work[1], &ldwork);
i__2 = *m - i__ - ib + 1;
i__3 = *n - i__ + 1;
clarfb_("Right", "Conjugate transpose", "Forward", "Rowwise", &i__2, &i__3, &ib, &a[i__ + i__ * a_dim1], lda, &work[ 1], &ldwork, &a[i__ + ib + i__ * a_dim1], lda, &work[ ib + 1], &ldwork);
}
i__2 = *n - i__ + 1;
cungl2_fla(&ib, &i__2, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], & work[1], &iinfo);
i__2 = i__ - 1;
for (j = 1;
j <= i__2;
++j)
{
i__3 = i__ + ib - 1;
for (l = i__;
l <= i__3;
++l)
{
i__4 = l + j * a_dim1;
a[i__4].r = 0.f;
a[i__4].i = 0.f;
}
}
}
}
work[1].r = (real) iws;
work[1].i = 0.f;
return 0;
}