Merge pull request #1 from eddelbuettel/master

gzt · web-flow · commit ae0df9c5c52e · 2021-09-25T20:15:19.000Z
Switch one PI to M_PI to support STRICT_R_HEADERS in Rcpp
diff --git a/src/stats.cpp b/src/stats.cpp
@@ -86,7 +86,7 @@ arma::colvec dmatnorm_calc(arma::cube& x, arma::mat& mean, arma::mat& U,
   arma::colvec logresult(x.n_slices);
   for (int i = 0; i < numslices; i++) {
     arma::mat XM = x.slice(i) - mean;
-    logresult(i) = -0.5 * n * p * log(2 * PI) - 0.5 * p * logdetU -
+    logresult(i) = -0.5 * n * p * log(2 * M_PI) - 0.5 * p * logdetU -
                    0.5 * n * logdetV -
                    0.5 * trace(Vinv * trans(XM) * Uinv * XM);
   }

Original file line number	Diff line number	Diff line change
`@@ -86,7 +86,7 @@ arma::colvec dmatnorm_calc(arma::cube& x, arma::mat& mean, arma::mat& U,`
`86`	`86`	`arma::colvec logresult(x.n_slices);`
`87`	`87`	`for (int i = 0; i < numslices; i++) {`
`88`	`88`	`arma::mat XM = x.slice(i) - mean;`
`89`		`- logresult(i) = -0.5 * n * p * log(2 * PI) - 0.5 * p * logdetU -`
	`89`	`+ logresult(i) = -0.5 * n * p * log(2 * M_PI) - 0.5 * p * logdetU -`
`90`	`90`	`0.5 * n * logdetV -`
`91`	`91`	`0.5 * trace(Vinv * trans(XM) * Uinv * XM);`
`92`	`92`	`}`