151 integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2;
153 integer i__, i1, i2, i3, ic, jc, mi, ni, nq;
156 extern logical lsame_(
char *,
char *);
158 int slarf_(
char *,
integer *,
integer *,
real *,
integer *,
real *,
real *,
integer *,
real *), xerbla_(
char *,
integer *);
183 a_offset = 1 + a_dim1;
187 c_offset = 1 + c_dim1;
192 left = lsame_(side,
"L");
193 notran = lsame_(trans,
"N");
203 if (! left && ! lsame_(side,
"R"))
207 else if (! notran && ! lsame_(trans,
"T"))
219 else if (*k < 0 || *k > nq)
223 else if (*lda < max(1,*k))
227 else if (*ldc < max(1,*m))
234 xerbla_(
"SORML2", &i__1);
238 if (*m == 0 || *n == 0 || *k == 0)
242 if (left && notran || ! left && ! notran)
267 i__2 < 0 ? i__ >= i__1 : i__ <= i__1;
283 aii = a[i__ + i__ * a_dim1];
284 a[i__ + i__ * a_dim1] = 1.f;
285 slarf_(side, &mi, &ni, &a[i__ + i__ * a_dim1], lda, &tau[i__], &c__[ ic + jc * c_dim1], ldc, &work[1]);
286 a[i__ + i__ * a_dim1] = aii;
float real
Definition: FLA_f2c.h:30
int logical
Definition: FLA_f2c.h:36
int integer
Definition: FLA_f2c.h:25