helpers for constructing bitwise 1s and 0s; more BH includes

kevinushey · kevinushey · commit 352f344e2712 · 2016-01-29T12:28:04.000-08:00
diff --git a/R/plugin.R b/R/plugin.R
@@ -7,6 +7,7 @@ inlineCxxPlugin <- function() {
         "-include", system.file("include/RcppNT2.h", package = "RcppNT2")
       )
     ),
+    LinkingTo = "BH",
     includes = "#include <RcppNT2.h>"
   )
 }
diff --git a/R/precompile.R b/R/precompile.R
@@ -25,6 +25,7 @@ precompileRcppNT2 <- function() {
     CXX,
     "-x c++-header", header,
     paste("-I", include, sep = ""),
+    paste("-I", system.file("include", package = "BH")),
     paste("-I", R.home("include"), sep = ""),
     "-std=c++11",
     cxxFlags
diff --git a/inst/examples/example-na-handling-variance.cpp b/inst/examples/example-na-handling-variance.cpp
@@ -9,23 +9,12 @@ using namespace RcppNT2;
 #include <Rcpp.h>
 using namespace Rcpp;
 
-namespace detail {
-
-union DoublePunner {
-  uint64_t i;
-  double d;
-};
-
-static const DoublePunner punner = {-1ULL};
-
-} // namespace detail
-
 struct IsNaN
 {
   template <typename T>
   T operator()(const T& data)
   {
-    return nt2::if_else(data == data, detail::punner.d, (double) 0x0);
+    return nt2::if_else(data == data, bitwise::ones(), bitwise::zeroes());
   }
 };
 
diff --git a/inst/include/RcppNT2.h b/inst/include/RcppNT2.h
@@ -41,6 +41,7 @@
 
 #include <RcppNT2/core/core.h>
 #include <RcppNT2/algorithm/algorithm.h>
+#include <RcppNT2/bitwise/bitwise.h>
 #include <RcppNT2/convert/convert.h>
 #include <RcppNT2/functor/functor.h>
 #include <RcppNT2/traits/traits.h>
diff --git a/inst/include/RcppNT2/bitwise/bitwise.h b/inst/include/RcppNT2/bitwise/bitwise.h
@@ -0,0 +1,52 @@
+#ifndef RCPP_NT2_BITWISE_BITWISE_H
+#define RCPP_NT2_BITWISE_BITWISE_H
+
+#include <utility>
+#include <cstdint>
+
+namespace RcppNT2 {
+namespace bitwise {
+
+namespace detail {
+
+// Used for constructing 'double's with all bits set to 1.
+union DoublePunner {
+  uint64_t i;
+  double d;
+};
+
+} // namespace detail
+
+// A class which can be (lazily) converted to an integral value
+// with all bits equal to 1. Primarily useful for 'double's.
+// Intended use case:
+//
+//    double ones = bitwise::ones();
+//    foo(static_cast<double>(bitwise::ones()));
+//
+class ones
+{
+public:
+
+  operator double() const {
+    detail::DoublePunner punner { -1ULL };
+    return punner.d;
+  }
+
+  operator int() const {
+    return ~0;
+  }
+
+};
+
+class zeroes
+{
+public:
+  template <typename T>
+  operator T() const { return T(); }
+};
+
+} // namespace bitwise
+} // namespace RcppNT2
+
+#endif /* RCPP_NT2_BITWISE_BITWISE_H */

Original file line number	Diff line number	Diff line change
`@@ -7,6 +7,7 @@ inlineCxxPlugin <- function() {`
`7`	`7`	`"-include", system.file("include/RcppNT2.h", package = "RcppNT2")`
`8`	`8`	`)`
`9`	`9`	`),`
	`10`	`+ LinkingTo = "BH",`
`10`	`11`	`includes = "#include <RcppNT2.h>"`
`11`	`12`	`)`
`12`	`13`	`}`