126 integer a_dim1, a_offset, i__1, i__2, i__3, i__4;
128 integer i__, j, l, ib, nb, ki, kk, nx, iws, nbmin, iinfo;
133 int zlarfb_(
char *,
char *,
char *,
char *,
integer *,
integer *,
integer *,
doublecomplex *,
integer *,
doublecomplex *,
integer *,
doublecomplex *,
integer *,
doublecomplex *,
integer *);
162 a_offset = 1 + a_dim1;
168 nb = ilaenv_(&c__1,
"ZUNGQR",
" ", m, n, k, &c_n1);
169 lwkopt = max(1,*n) * nb;
172 lquery = *lwork == -1;
177 else if (*n < 0 || *n > *m)
181 else if (*k < 0 || *k > *n)
185 else if (*lda < max(1,*m))
189 else if (*lwork < max(1,*n) && ! lquery)
196 xerbla_(
"ZUNGQR", &i__1);
213 if (nb > 1 && nb < *k)
218 i__2 = ilaenv_(&c__3,
"ZUNGQR",
" ", m, n, k, &c_n1);
229 nb = *lwork / ldwork;
232 i__2 = ilaenv_(&c__2,
"ZUNGQR",
" ", m, n, k, &c_n1);
233 nbmin = max(i__1,i__2);
237 if (nb >= nbmin && nb < *k && nx < *k)
241 ki = (*k - nx - 1) / nb * nb;
257 i__3 = i__ + j * a_dim1;
275 zung2r_fla(&i__1, &i__2, &i__3, &a[kk + 1 + (kk + 1) * a_dim1], lda, & tau[kk + 1], &work[1], &iinfo);
282 i__1 < 0 ? i__ >= 1 : i__ <= 1;
294 zlarft_(
"Forward",
"Columnwise", &i__2, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], &work[1], &ldwork);
297 i__3 = *n - i__ - ib + 1;
298 zlarfb_(
"Left",
"No transpose",
"Forward",
"Columnwise", & i__2, &i__3, &ib, &a[i__ + i__ * a_dim1], lda, &work[ 1], &ldwork, &a[i__ + (i__ + ib) * a_dim1], lda, & work[ib + 1], &ldwork);
302 zung2r_fla(&i__2, &ib, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], & work[1], &iinfo);
314 i__4 = l + j * a_dim1;
int zung2r_fla(integer *m, integer *n, integer *k, doublecomplex *a, integer *lda, doublecomplex *tau, doublecomplex *work, integer *info)
Definition: zung2r.c:105
doublereal r
Definition: FLA_f2c.h:33
double doublereal
Definition: FLA_f2c.h:31
doublereal i
Definition: FLA_f2c.h:33
int logical
Definition: FLA_f2c.h:36
int integer
Definition: FLA_f2c.h:25