105 integer a_dim1, a_offset, i__1, i__2;
110 int sscal_(
integer *,
real *,
real *,
integer *), slarf_(
char *,
integer *,
integer *,
real *,
integer *,
real *,
real *,
integer *,
real *), 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_(
"SORGL2", &i__1);
178 a[l + j * a_dim1] = 0.f;
181 if (j > *k && j <= *m)
183 a[j + j * a_dim1] = 1.f;
197 a[i__ + i__ * a_dim1] = 1.f;
200 slarf_(
"Right", &i__1, &i__2, &a[i__ + i__ * a_dim1], lda, & tau[i__], &a[i__ + 1 + i__ * a_dim1], lda, &work[1]);
204 sscal_(&i__1, &r__1, &a[i__ + (i__ + 1) * a_dim1], lda);
206 a[i__ + i__ * a_dim1] = 1.f - tau[i__];
213 a[i__ + l * a_dim1] = 0.f;
int integer
Definition: FLA_f2c.h:25
float real
Definition: FLA_f2c.h:30