107 integer q_dim1, q_offset, i__1, i__2, i__3;
110 extern logical lsame_(
char *,
char *);
114 int dorg2l_(
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *),
dorg2r_fla(
integer *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *), xerbla_(
char *,
integer *);
140 q_offset = 1 + q_dim1;
145 upper = lsame_(uplo,
"U");
146 if (! upper && ! lsame_(uplo,
"L"))
154 else if (*ldq < max(1,*n))
161 xerbla_(
"DOPGTR", &i__1);
186 q[i__ + j * q_dim1] = ap[ij];
191 q[*n + j * q_dim1] = 0.;
199 q[i__ + *n * q_dim1] = 0.;
202 q[*n + *n * q_dim1] = 1.;
207 dorg2l_(&i__1, &i__2, &i__3, &q[q_offset], ldq, &tau[1], &work[1], & iinfo);
221 q[i__ + q_dim1] = 0.;
230 q[j * q_dim1 + 1] = 0.;
236 q[i__ + j * q_dim1] = ap[ij];
249 dorg2r_fla(&i__1, &i__2, &i__3, &q[(q_dim1 << 1) + 2], ldq, &tau[1], &work[1], &iinfo);
double doublereal
Definition: FLA_f2c.h:31
int integer
Definition: FLA_f2c.h:25
int logical
Definition: FLA_f2c.h:36
int dorg2r_fla(integer *m, integer *n, integer *k, doublereal *a, integer *lda, doublereal *tau, doublereal *work, integer *info)
Definition: dorg2r.c:105