compare to std::accumulate

Author: kevinushey

examples/boost-simd-accumulate.cpp

@@ -17,17 +17,22 @@ struct plus
 // [[Rcpp::export]]
 double simd_sum(NumericVector x)
 {
-   return boost::simd::accumulate(
-      x.begin(),
-      x.end(),
-      0.0,
-      plus()
-   );
+   return boost::simd::accumulate(x.begin(), x.end(), 0.0, plus());
+}
+
+// [[Rcpp::export]]
+double cpp_sum(NumericVector x)
+{
+   return std::accumulate(x.begin(), x.end(), 0.0, plus());
 }
 
 /***R
 n <- 1024 * 1000
 data <- rnorm(n)
 simd_sum(data)
-microbenchmark(R = sum(data), simd = simd_sum(data))
-*/
+microbenchmark::microbenchmark(
+   R    = sum(data),
+   simd = simd_sum(data),
+   cpp  = cpp_sum(data)
+)
+*/