SEXP d_asum(SEXP vList, SEXP inc)
{
int n, increment = asInteger(inc);
double * dPtr;
unpackVector(vList, &n, &dPtr);
SEXP out;
PROTECT(out = allocVector(REALSXP, 1));
REAL(out)[0] = cublasDasum(n, dPtr, increment);
checkCublasError("d_asum");
UNPROTECT(1);
return out;
}
void caffe_gpu_asum<double>(const int n, const double* x, double* y) {
CUBLAS_CHECK(cublasDasum(Caffe::cublas_handle(), n, x, 1, y));
}
double asum(const Vector<double> &x)
{
return cublasDasum(x.getSize(), x, x.inc());
}
double magma_dasum(
magma_int_t n,
const double *dx, magma_int_t incx )
{
return cublasDasum( n, dx, incx );
}
