173 integer a_dim1, a_offset, c_dim1, c_offset, i__2, i__3;
179 integer i1, i2, nb, mi, ni, nq, nw;
181 extern logical lsame_(
char *,
char *);
185 int xerbla_(
char *,
integer *);
188 int sormql_(
char *,
char *,
integer *,
integer *,
integer *,
real *,
integer *,
real *,
real *,
integer *,
real *,
integer *,
integer *);
192 int sormqr_fla(
char *,
char *,
integer *,
integer *,
integer *,
real *,
integer *,
real *,
real *,
integer *,
real *,
integer *,
integer *);
214 a_offset = 1 + a_dim1;
218 c_offset = 1 + c_dim1;
223 left = lsame_(side,
"L");
224 upper = lsame_(uplo,
"U");
225 lquery = *lwork == -1;
237 if (! left && ! lsame_(side,
"R"))
241 else if (! upper && ! lsame_(uplo,
"L"))
245 else if (! lsame_(trans,
"N") && ! lsame_(trans,
"T"))
257 else if (*lda < max(1,nq))
261 else if (*ldc < max(1,*m))
265 else if (*lwork < max(1,nw) && ! lquery)
277 nb = ilaenv_(&c__1,
"SORMQL", ch__1, &i__2, n, &i__3, &c_n1);
283 nb = ilaenv_(&c__1,
"SORMQL", ch__1, m, &i__2, &i__3, &c_n1);
292 nb = ilaenv_(&c__1,
"SORMQR", ch__1, &i__2, n, &i__3, &c_n1);
298 nb = ilaenv_(&c__1,
"SORMQR", ch__1, m, &i__2, &i__3, &c_n1);
301 lwkopt = max(1,nw) * nb;
302 work[1] = (
real) lwkopt;
307 xerbla_(
"SORMTR", &i__2);
315 if (*m == 0 || *n == 0 || nq == 1)
334 sormql_(side, trans, &mi, &ni, &i__2, &a[(a_dim1 << 1) + 1], lda, & tau[1], &c__[c_offset], ldc, &work[1], lwork, &iinfo);
350 sormqr_fla(side, trans, &mi, &ni, &i__2, &a[a_dim1 + 2], lda, &tau[1], & c__[i1 + i2 * c_dim1], ldc, &work[1], lwork, &iinfo);
352 work[1] = (
real) lwkopt;
float real
Definition: FLA_f2c.h:30
int logical
Definition: FLA_f2c.h:36
int integer
Definition: FLA_f2c.h:25
int sormqr_fla(char *side, char *trans, integer *m, integer *n, integer *k, real *a, integer *lda, real *tau, real *c__, integer *ldc, real *work, integer *lwork, integer *info)
Definition: sormqr.c:169