inline void gbmv( const Order, const Trans, const int m, const int n,
const int kl, const int ku, const double alpha, const double* a,
const int lda, const double* x, const int incx, const double beta,
double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDgbmv( blas_option< Trans >::value, m, n, kl, ku, alpha, a, lda, x,
incx, beta, y, incy );
}
void d_gbmv(SEXP rtrans, SEXP rkl, SEXP rku, SEXP ralpha, SEXP ra, SEXP rlda,
SEXP rx, SEXP rincx, SEXP rbeta, SEXP ry, SEXP rincy)
{
char
trans = getTranspose(rtrans);
double
alpha = asReal(ralpha), beta = asReal(rbeta),
* a, * x, * y;
int
rowsa, colsa,
lda = asInteger(rlda),
kl = asInteger(rkl), ku = asInteger(rku),
nx, ny,
incx = asInteger(rincx), incy = asInteger(rincy);
unpackVector(rx, &nx, &x);
unpackVector(ry, &ny, &y);
unpackMatrix(ra, &rowsa, &colsa, &a);
cublasDgbmv(trans, rowsa, colsa, kl, ku, alpha, a, lda, x, incx,
beta, y, incy);
checkCublasError("d_gbmv");
}