{
integer a_dim1, a_offset, c_dim1, c_offset, i__2;
char ch__1[2];
integer i1, i2, nb, mi, nh, ni, nq, nw;
logical left;
extern logical lsame_(char *, char *);
integer iinfo;
extern
int xerbla_(char *, integer *);
extern integer ilaenv_(integer *, char *, char *, integer *, integer *, integer *, integer *);
extern
int dormqr_fla(char *, char *, integer *, integer *, integer *, doublereal *, integer *, doublereal *, doublereal *, integer *, doublereal *, integer *, integer *);
integer lwkopt;
logical lquery;
a_dim1 = *lda;
a_offset = 1 + a_dim1;
a -= a_offset;
--tau;
c_dim1 = *ldc;
c_offset = 1 + c_dim1;
c__ -= c_offset;
--work;
*info = 0;
nh = *ihi - *ilo;
left = lsame_(side, "L");
lquery = *lwork == -1;
if (left)
{
nq = *m;
nw = *n;
}
else
{
nq = *n;
nw = *m;
}
if (! left && ! lsame_(side, "R"))
{
*info = -1;
}
else if (! lsame_(trans, "N") && ! lsame_(trans, "T"))
{
*info = -2;
}
else if (*m < 0)
{
*info = -3;
}
else if (*n < 0)
{
*info = -4;
}
else if (*ilo < 1 || *ilo > max(1,nq))
{
*info = -5;
}
else if (*ihi < min(*ilo,nq) || *ihi > nq)
{
*info = -6;
}
else if (*lda < max(1,nq))
{
*info = -8;
}
else if (*ldc < max(1,*m))
{
*info = -11;
}
else if (*lwork < max(1,nw) && ! lquery)
{
*info = -13;
}
if (*info == 0)
{
if (left)
{
nb = ilaenv_(&c__1, "DORMQR", ch__1, &nh, n, &nh, &c_n1);
}
else
{
nb = ilaenv_(&c__1, "DORMQR", ch__1, m, &nh, &nh, &c_n1);
}
lwkopt = max(1,nw) * nb;
work[1] = (doublereal) lwkopt;
}
if (*info != 0)
{
i__2 = -(*info);
xerbla_("DORMHR", &i__2);
return 0;
}
else if (lquery)
{
return 0;
}
if (*m == 0 || *n == 0 || nh == 0)
{
work[1] = 1.;
return 0;
}
if (left)
{
mi = nh;
ni = *n;
i1 = *ilo + 1;
i2 = 1;
}
else
{
mi = *m;
ni = nh;
i1 = 1;
i2 = *ilo + 1;
}
dormqr_fla(side, trans, &mi, &ni, &nh, &a[*ilo + 1 + *ilo * a_dim1], lda, & tau[*ilo], &c__[i1 + i2 * c_dim1], ldc, &work[1], lwork, &iinfo);
work[1] = (doublereal) lwkopt;
return 0;
}