105 integer a_dim1, a_offset, i__1, i__2;
110 int dscal_(
integer *,
doublereal *,
doublereal *,
integer *), dlarf_(
char *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *,
doublereal *), xerbla_(
char *,
integer *);
132 a_offset = 1 + a_dim1;
146 else if (*k < 0 || *k > *m)
150 else if (*lda < max(1,*m))
157 xerbla_(
"DORGL2", &i__1);
178 a[l + j * a_dim1] = 0.;
181 if (j > *k && j <= *m)
183 a[j + j * a_dim1] = 1.;
197 a[i__ + i__ * a_dim1] = 1.;
200 dlarf_(
"Right", &i__1, &i__2, &a[i__ + i__ * a_dim1], lda, & tau[i__], &a[i__ + 1 + i__ * a_dim1], lda, &work[1]);
204 dscal_(&i__1, &d__1, &a[i__ + (i__ + 1) * a_dim1], lda);
206 a[i__ + i__ * a_dim1] = 1. - tau[i__];
213 a[i__ + l * a_dim1] = 0.;
double doublereal
Definition: FLA_f2c.h:31
int integer
Definition: FLA_f2c.h:25