Missed these files

Author: kevinushey

inst/include/Rcpp/traits/NAComparator.h

@@ -0,0 +1,99 @@
+// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 8 -*-
+/* :tabSize=4:indentSize=4:noTabs=false:folding=explicit:collapseFolds=1: */
+//
+// NAComparator.h: Rcpp R/C++ interface class library -- comparator
+//
+// Copyright (C) 2012-2014 Dirk Eddelbuettel, Romain Francois and Kevin Ushey
+//
+// This file is part of Rcpp.
+//
+// Rcpp is free software: you can redistribute it and/or modify it
+// under the terms of the GNU General Public License as published by
+// the Free Software Foundation, either version 2 of the License, or
+// (at your option) any later version.
+//
+// Rcpp is distributed in the hope that it will be useful, but
+// WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with Rcpp.  If not, see <http://www.gnu.org/licenses/>.
+
+#ifndef Rcpp__traits__NAComparator__h
+#define Rcpp__traits__NAComparator__h
+
+namespace Rcpp{
+namespace traits{
+   
+inline bool Rcpp_IsNA(double x) {
+  return memcmp(
+    &x,
+    &NA_REAL,
+    sizeof(double)
+  ) == 0;
+}
+
+inline bool Rcpp_IsNaN(double x) {
+  return memcmp(
+    &x,
+    &R_NaN,
+    sizeof(double)
+  ) == 0;
+}
+
+inline int StrCmp(SEXP x, SEXP y) {
+    if (x == NA_STRING) return (y == NA_STRING ? 0 : 1);
+    if (y == NA_STRING) return -1;
+    if (x == y) return 0;  // same string in cache
+    return strcmp(char_nocheck(x), char_nocheck(y));
+}
+
+template <typename T>
+struct NAComparator {
+  inline bool operator()(T left, T right) const {
+    return left < right;
+  }
+};
+
+template <>
+struct NAComparator<int> {
+  inline bool operator()(int left, int right) const {
+    if (left == NA_INTEGER) return false;
+    if (right == NA_INTEGER) return true;
+    return left < right;
+  }
+};
+
+template <>
+struct NAComparator<double> {  
+  inline bool operator()(double left, double right) const {
+        
+    bool leftNaN = (left != left);
+    bool rightNaN = (right != right);
+    
+    // this branch inspired by data.table: see
+    // https://github.com/arunsrinivasan/datatable/commit/1a3e476d3f746e18261662f484d2afa84ac7a146#commitcomment-4885242
+    if (Rcpp_IsNaN(right) and Rcpp_IsNA(left)) return true;
+    
+    if (leftNaN != rightNaN) {
+      return leftNaN < rightNaN;
+    } else {
+      return left < right;
+    }
+
+  }
+  
+};
+
+template <>
+struct NAComparator<SEXP> {
+  inline bool operator()(SEXP left, SEXP right) const {
+    return StrCmp(left, right) < 0;
+  }
+};
+
+}
+}     
+
+#endif

inst/include/Rcpp/traits/NAEquals.h

@@ -0,0 +1,45 @@
+// NAEquals.h: Rcpp R/C++ interface class library -- for allowing NA == NA
+//
+// Copyright (C) 2014 Kevin UShey
+//
+// This file is part of Rcpp.
+//
+// Rcpp is free software: you can redistribute it and/or modify it
+// under the terms of the GNU General Public License as published by
+// the Free Software Foundation, either version 2 of the License, or
+// (at your option) any later version.
+//
+// Rcpp is distributed in the hope that it will be useful, but
+// WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with Rcpp.  If not, see <http://www.gnu.org/licenses/>.
+
+#ifndef Rcpp__traits__NAEquals__h
+#define Rcpp__traits__NAEquals__h
+
+namespace Rcpp {
+  
+namespace traits {
+  
+template <typename T>
+struct NAEquals {
+  inline bool operator()(T left, T right) const {
+    return left == right;
+  }
+};
+
+template <>
+struct NAEquals<double> {
+  inline bool operator()(double left, double right) const {
+    return memcmp(&left, &right, sizeof(double)) == 0;
+  }
+};
+
+}
+
+}
+
+#endif