示例#1
0
文件: ardblas.c 项目: rforge/gcb
void d_spr(SEXP ruplo, SEXP ralpha, SEXP rx, SEXP rincx, SEXP ra)
{
	char
		uplo = getSymLoc(ruplo);
	double
		alpha = asReal(ralpha),
		* a, * x;
	int
		rowsa, colsa,
		nx, incx = asInteger(rincx);

	unpackVector(rx, &nx, &x);
	unpackMatrix(ra, &rowsa, &colsa, &a);

	cublasDspr(uplo, rowsa, alpha, x, incx, a);
	checkCublasError("d_spr");
}
inline void hpr( const Order, const UpLo, const int n, const double alpha,
        const double* x, const int incx, double* ap ) {
    BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
    cublasDspr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}