125 integer a_dim1, a_offset, i__1, i__2, i__3;
127 integer i__, j, l, ib, nb, ki, kk, nx, iws, nbmin, iinfo;
129 int dorgl2_fla(
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *), dlarfb_(
char *,
char *,
char *,
char *,
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
integer *,
doublereal *,
integer *,
doublereal *,
integer *), dlarft_(
char *,
char *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *), xerbla_(
char *,
integer *);
156 a_offset = 1 + a_dim1;
162 nb = ilaenv_(&c__1,
"DORGLQ",
" ", m, n, k, &c_n1);
163 lwkopt = max(1,*m) * nb;
165 lquery = *lwork == -1;
174 else if (*k < 0 || *k > *m)
178 else if (*lda < max(1,*m))
182 else if (*lwork < max(1,*m) && ! lquery)
189 xerbla_(
"DORGLQ", &i__1);
205 if (nb > 1 && nb < *k)
210 i__2 = ilaenv_(&c__3,
"DORGLQ",
" ", m, n, k, &c_n1);
221 nb = *lwork / ldwork;
224 i__2 = ilaenv_(&c__2,
"DORGLQ",
" ", m, n, k, &c_n1);
225 nbmin = max(i__1,i__2);
229 if (nb >= nbmin && nb < *k && nx < *k)
233 ki = (*k - nx - 1) / nb * nb;
249 a[i__ + j * a_dim1] = 0.;
265 dorgl2_fla(&i__1, &i__2, &i__3, &a[kk + 1 + (kk + 1) * a_dim1], lda, & tau[kk + 1], &work[1], &iinfo);
272 i__1 < 0 ? i__ >= 1 : i__ <= 1;
284 dlarft_(
"Forward",
"Rowwise", &i__2, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], &work[1], &ldwork);
286 i__2 = *m - i__ - ib + 1;
288 dlarfb_(
"Right",
"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);
292 dorgl2_fla(&ib, &i__2, &ib, &a[i__ + i__ * a_dim1], lda, &tau[i__], & work[1], &iinfo);
304 a[l + j * a_dim1] = 0.;
int dorgl2_fla(integer *m, integer *n, integer *k, doublereal *a, integer *lda, doublereal *tau, doublereal *work, integer *info)
Definition: dorgl2.c:102
double doublereal
Definition: FLA_f2c.h:31
int logical
Definition: FLA_f2c.h:36
int integer
Definition: FLA_f2c.h:25