108 integer a_dim1, a_offset, i__1, i__2, i__3;
113 int cscal_(
integer *,
complex *,
complex *,
integer *), clarf_(
char *,
integer *,
integer *,
complex *,
integer *,
complex *,
complex *,
integer *,
complex *), 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_(
"CUNG2R", &i__1);
179 i__3 = l + j * a_dim1;
184 i__2 = j + j * a_dim1;
196 i__1 = i__ + i__ * a_dim1;
201 clarf_(
"Left", &i__1, &i__2, &a[i__ + i__ * a_dim1], &c__1, &tau[ i__], &a[i__ + (i__ + 1) * a_dim1], lda, &work[1]);
207 q__1.
r = -tau[i__2].
r;
208 q__1.
i = -tau[i__2].
i;
209 cscal_(&i__1, &q__1, &a[i__ + 1 + i__ * a_dim1], &c__1);
211 i__1 = i__ + i__ * a_dim1;
213 q__1.
r = 1.f - tau[i__2].
r;
214 q__1.
i = 0.f - tau[i__2].
i;
223 i__2 = l + i__ * a_dim1;
int integer
Definition: FLA_f2c.h:25
real r
Definition: FLA_f2c.h:32
real i
Definition: FLA_f2c.h:32