108 integer a_dim1, a_offset, i__1, i__2;
113 int dscal_(
integer *,
doublereal *,
doublereal *,
integer *), dlarf_(
char *,
integer *,
integer *,
doublereal *,
integer *,
doublereal *,
doublereal *,
integer *,
doublereal *), xerbla_(
char *,
integer *);
135 a_offset = 1 + a_dim1;
145 else if (*n < 0 || *n > *m)
149 else if (*k < 0 || *k > *n)
153 else if (*lda < max(1,*m))
160 xerbla_(
"DORG2R", &i__1);
179 a[l + j * a_dim1] = 0.;
182 a[j + j * a_dim1] = 1.;
192 a[i__ + i__ * a_dim1] = 1.;
195 dlarf_(
"Left", &i__1, &i__2, &a[i__ + i__ * a_dim1], &c__1, &tau[ i__], &a[i__ + (i__ + 1) * a_dim1], lda, &work[1]);
201 dscal_(&i__1, &d__1, &a[i__ + 1 + i__ * a_dim1], &c__1);
203 a[i__ + i__ * a_dim1] = 1. - tau[i__];
210 a[l + i__ * a_dim1] = 0.;
double doublereal
Definition: FLA_f2c.h:31
int integer
Definition: FLA_f2c.h:25