Example #1
0
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 );
}
Example #2
0
File: ardblas.c Project: rforge/gcb
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");
}