120 integer a_dim1, a_offset, i__1, i__2, i__3;
123 extern logical lsame_(
char *,
char *);
127 int xerbla_(
char *,
integer *);
130 int dorgql_(
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *,
integer *),
dorgqr_fla(
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *,
integer *);
156 a_offset = 1 + a_dim1;
162 lquery = *lwork == -1;
163 upper = lsame_(uplo,
"U");
164 if (! upper && ! lsame_(uplo,
"L"))
172 else if (*lda < max(1,*n))
181 if (*lwork < max(i__1,i__2) && ! lquery)
193 nb = ilaenv_(&c__1,
"DORGQL",
" ", &i__1, &i__2, &i__3, &c_n1);
200 nb = ilaenv_(&c__1,
"DORGQR",
" ", &i__1, &i__2, &i__3, &c_n1);
205 lwkopt = max(i__1,i__2) * nb;
211 xerbla_(
"DORGTR", &i__1);
240 a[i__ + j * a_dim1] = a[i__ + (j + 1) * a_dim1];
243 a[*n + j * a_dim1] = 0.;
251 a[i__ + *n * a_dim1] = 0.;
254 a[*n + *n * a_dim1] = 1.;
259 dorgql_(&i__1, &i__2, &i__3, &a[a_offset], lda, &tau[1], &work[1], lwork, &iinfo);
271 a[j * a_dim1 + 1] = 0.;
277 a[i__ + j * a_dim1] = a[i__ + (j - 1) * a_dim1];
288 a[i__ + a_dim1] = 0.;
297 dorgqr_fla(&i__1, &i__2, &i__3, &a[(a_dim1 << 1) + 2], lda, &tau[1], &work[1], lwork, &iinfo);
int dorgqr_fla(integer *m, integer *n, integer *k, doublereal *a, integer *lda, doublereal *tau, doublereal *work, integer *lwork, integer *info)
Definition: dorgqr.c:123
double doublereal
Definition: FLA_f2c.h:31
int logical
Definition: FLA_f2c.h:36
int integer
Definition: FLA_f2c.h:25