define STRICT_R_HEADERS, change PI to M_PI

eddelbuettel · eddelbuettel · commit 907536cb4431 · 2021-05-02T16:02:24.000+02:00
diff --git a/src/compute_groupcond_logl.cpp b/src/compute_groupcond_logl.cpp
@@ -1,3 +1,4 @@
+#define STRICT_R_HEADERS
 #include <Rcpp.h>
 using namespace Rcpp;
 
@@ -41,7 +42,7 @@ double compute_groupcond_logl(NumericMatrix x,
   int wn = dims[3];
   
   // constants
-  double sqrtpi = std::sqrt(2*PI);
+  double sqrtpi = std::sqrt(2*M_PI);
   double cz = std::pow(sqrtpi, zn);
   double cw = std::pow(sqrtpi, wn);
   
@@ -128,4 +129,4 @@ double compute_groupcond_logl(NumericMatrix x,
   }
   
   return out - N*std::log(cz*std::sqrt(detveps)*cw*std::sqrt(detvarz)*gamisize);
-}
+}
diff --git a/src/creg_group_logl.cpp b/src/creg_group_logl.cpp
@@ -1,3 +1,4 @@
+#define STRICT_R_HEADERS
 #include <Rcpp.h>
 using namespace Rcpp;
 
@@ -54,7 +55,7 @@ NumericVector creg_group_logl_cpp(List datalist, List modellist){
         int wn = dims[3];
         
         // constants
-        double sqrtpi = std::sqrt(2*PI);
+        double sqrtpi = std::sqrt(2*M_PI);
         double cz = std::pow(sqrtpi, zn);
         double cw = std::pow(sqrtpi, wn);
         
@@ -144,4 +145,4 @@ NumericVector creg_group_logl_cpp(List datalist, List modellist){
         
   }
   return res;
-}
+}
