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