20 integer a_dim1, a_offset, i__1, i__2, i__3, i__4;
170 a_offset = 1 + a_dim1;
184 }
else if (*
lda < *k + 1) {
186 }
else if (*
incx == 0) {
188 }
else if (*
incy == 0) {
198 if (*
n == 0 || (*
alpha == 0. && *
beta == 1.)) {
207 kx = 1 - (*
n - 1) * *
incx;
212 ky = 1 - (*
n - 1) * *
incy;
224 for (i__ = 1; i__ <= i__1; ++i__) {
230 for (i__ = 1; i__ <= i__1; ++i__) {
239 for (i__ = 1; i__ <= i__1; ++i__) {
246 for (i__ = 1; i__ <= i__1; ++i__) {
264 for (
j = 1;
j <= i__1; ++
j) {
269 i__2 = 1, i__3 =
j - *k;
271 for (i__ =
max(i__2,i__3); i__ <= i__4; ++i__) {
272 y[i__] += temp1 *
a[
l + i__ +
j * a_dim1];
273 temp2 +=
a[
l + i__ +
j * a_dim1] *
x[i__];
276 y[
j] =
y[
j] + temp1 *
a[kplus1 +
j * a_dim1] + *
alpha * temp2;
283 for (
j = 1;
j <= i__1; ++
j) {
290 i__4 = 1, i__2 =
j - *k;
292 for (i__ =
max(i__4,i__2); i__ <= i__3; ++i__) {
293 y[iy] += temp1 *
a[
l + i__ +
j * a_dim1];
294 temp2 +=
a[
l + i__ +
j * a_dim1] *
x[ix];
299 y[jy] =
y[jy] + temp1 *
a[kplus1 +
j * a_dim1] + *
alpha *
316 for (
j = 1;
j <= i__1; ++
j) {
319 y[
j] += temp1 *
a[
j * a_dim1 + 1];
322 i__4 = *
n, i__2 =
j + *k;
323 i__3 =
min(i__4,i__2);
324 for (i__ =
j + 1; i__ <= i__3; ++i__) {
325 y[i__] += temp1 *
a[
l + i__ +
j * a_dim1];
326 temp2 +=
a[
l + i__ +
j * a_dim1] *
x[i__];
336 for (
j = 1;
j <= i__1; ++
j) {
339 y[jy] += temp1 *
a[
j * a_dim1 + 1];
344 i__4 = *
n, i__2 =
j + *k;
345 i__3 =
min(i__4,i__2);
346 for (i__ =
j + 1; i__ <= i__3; ++i__) {
349 y[iy] += temp1 *
a[
l + i__ +
j * a_dim1];
350 temp2 +=
a[
l + i__ +
j * a_dim1] *
x[ix];