void CNAME(BLASLONG m, BLASLONG n, BLASLONG ku, BLASLONG kl, FLOAT alpha_r, FLOAT alpha_i,
FLOAT *a, BLASLONG lda,
FLOAT *x, BLASLONG incx, FLOAT *y, BLASLONG incy, void *buffer){
BLASLONG i, offset_u, offset_l, start, end, length;
FLOAT *X = x;
FLOAT *Y = y;
FLOAT *gemvbuffer = (FLOAT *)buffer;
FLOAT *bufferY = gemvbuffer;
FLOAT *bufferX = gemvbuffer;
#ifdef TRANS
FLOAT _Complex temp;
#endif
if (incy != 1) {
Y = bufferY;
bufferX = (FLOAT *)(((BLASLONG)bufferY + M * sizeof(FLOAT) * 2 + 4095) & ~4095);
gemvbuffer = bufferX;
COPY_K(M, y, incy, Y, 1);
}
if (incx != 1) {
X = bufferX;
gemvbuffer = (FLOAT *)(((BLASLONG)bufferX + N * sizeof(FLOAT) * 2 + 4095) & ~4095);
COPY_K(N, x, incx, X, 1);
}
offset_u = ku;
offset_l = ku + m;
for (i = 0; i < MIN(n, m + ku); i++) {
start = MAX(offset_u, 0);
end = MIN(offset_l, ku + kl + 1);
length = end - start;
#ifndef TRANS
ZAXPY(length, 0, 0,
#ifndef XCONJ
alpha_r * X[i * 2 + 0] - alpha_i * X[i * 2 + 1],
alpha_i * X[i * 2 + 0] + alpha_r * X[i * 2 + 1],
#else
alpha_r * X[i * 2 + 0] + alpha_i * X[i * 2 + 1],
alpha_i * X[i * 2 + 0] - alpha_r * X[i * 2 + 1],
#endif
a + start * 2, 1, Y + (start - offset_u) * 2, 1, NULL, 0);
#else
#ifndef XCONJ
temp = ZDOT(length, a + start * 2, 1, X + (start - offset_u) * 2, 1);
#else
temp = ZDOT(length, X + (start - offset_u) * 2, 1, a + start * 2, 1);
#endif
#if !defined(XCONJ) || !defined(CONJ)
Y[i * 2 + 0] += alpha_r * CREAL(temp) - alpha_i * CIMAG(temp);
Y[i * 2 + 1] += alpha_i * CREAL(temp) + alpha_r * CIMAG(temp);
#else
Y[i * 2 + 0] += alpha_r * CREAL(temp) + alpha_i * CIMAG(temp);
Y[i * 2 + 1] += alpha_i * CREAL(temp) - alpha_r * CIMAG(temp);
#endif
#endif
offset_u --;
offset_l --;
a += lda * 2;
}
if (incy != 1) {
COPY_K(M, Y, 1, y, incy);
}
return;
}
void Num_axpy_zprimme(int n, Complex_Z alpha, Complex_Z *x, int incx,
Complex_Z *y, int incy) {
ZAXPY(&n, &alpha, x, &incx, y, &incy);
}
