Added random number generator unit tests (closes #28)

Adds unit test for random generation for the R:: namespace where Rmath
is kept.

Added section breaks in C++ code for each distribution

Author: coatless

ChangeLog

@@ -1,3 +1,9 @@
+2016-07-22  James J Balamuta  <balamut2@illinois.edu>
+
+        * inst/unitTests/cpp/rmath.cpp: Added unit test functions for Rmath
+        random number generators.
+        * inst/unitTests/runit.rmath.R: idem
+
 2016-07-21  Dirk Eddelbuettel  <edd@debian.org>
 
         * DESCRIPTION (Version): Roll minor release

inst/NEWS.Rd

@@ -10,6 +10,12 @@
       \item The \code{NORET} macro is now defined if it was not already defined
       by R (Kevin fixing issue \ghit{512}). 
     }
+    \item Changes in Rcpp unit tests
+    \itemize{
+      \item New unit tests for random number generators the R namespace which 
+      call the standalone Rmath library. (James Balamuta in \ghpr{514} 
+      addressing issue \ghit{28}).
+    }
   }
 }
 

inst/unitTests/cpp/rmath.cpp

@@ -22,6 +22,9 @@
 #include <Rcpp.h>
 using namespace Rcpp ;
 
+
+// ------------------- Normal Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnorm( double x, double a, double b ){
 	     return NumericVector::create(R::dnorm(x, a, b, 0), R::dnorm(x, a, b, 1));
@@ -39,6 +42,13 @@ NumericVector runit_qnorm( double x, double a, double b ){
                                           R::qnorm(x, a, b, 0, 0), R::qnorm(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rnorm( double a, double b ){
+    return NumericVector::create(R::rnorm(a, b), R::rnorm(a, b));
+}
+
+// ------------------- Uniform Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dunif( double x, double a, double b ){
 	     return NumericVector::create(R::dunif(x, a, b, 0), R::dunif(x, a, b, 1));
@@ -56,6 +66,14 @@ NumericVector runit_qunif( double x, double a, double b ){
                                           R::qunif(x, a, b, 0, 0), R::qunif(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_runif( double a, double b ){
+    return NumericVector::create(R::runif(a, b), R::runif(a, b));
+}
+
+
+// ------------------- Gamma Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dgamma( double x, double a, double b ){
 	     return NumericVector::create(R::dgamma(x, a, b, 0), R::dgamma(x, a, b, 1));
@@ -73,6 +91,14 @@ NumericVector runit_qgamma( double x, double a, double b ){
                                           R::qgamma(x, a, b, 0, 0), R::qgamma(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rgamma( double a, double b ){
+    return NumericVector::create(R::rgamma(a, b), R::rgamma(a, b));
+}
+
+
+// ------------------- Beta Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dbeta( double x, double a, double b ){
 	     return NumericVector::create(R::dbeta(x, a, b, 0), R::dbeta(x, a, b, 1));
@@ -90,6 +116,14 @@ NumericVector runit_qbeta( double x, double a, double b ){
                                           R::qbeta(x, a, b, 0, 0), R::qbeta(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rbeta( double a, double b ){
+    return NumericVector::create(R::rbeta(a, b), R::rbeta(a, b));
+}
+
+
+// ------------------- Log Normal Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dlnorm( double x, double a, double b ){
 	     return NumericVector::create(R::dlnorm(x, a, b, 0), R::dlnorm(x, a, b, 1));
@@ -107,6 +141,13 @@ NumericVector runit_qlnorm( double x, double a, double b ){
                                           R::qlnorm(x, a, b, 0, 0), R::qlnorm(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rlnorm( double a, double b ){
+    return NumericVector::create(R::rlnorm(a, b), R::rlnorm(a, b));
+}
+
+// ------------------- Chi-Squared Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dchisq( double x, double a ){
 	     return NumericVector::create(R::dchisq(x, a, 0), R::dchisq(x, a, 1));
@@ -124,6 +165,13 @@ NumericVector runit_qchisq( double x, double a ){
                                           R::qchisq(x, a, 0, 0), R::qchisq(log(x), a, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rchisq( double a ){
+    return NumericVector::create(R::rchisq(a), R::rchisq(a));
+}
+
+// ------------------- Non-central Chi-Squared Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnchisq( double x, double a, double b ){
 	     return NumericVector::create(R::dnchisq(x, a, b, 0), R::dnchisq(x, a, b, 1));
@@ -141,6 +189,8 @@ NumericVector runit_qnchisq( double x, double a, double b ){
                                           R::qnchisq(x, a, b, 0, 0), R::qnchisq(log(x), a, b, 0, 1));
 }
 
+// -------------------  F Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_df( double x, double a, double b ){
 	     return NumericVector::create(R::df(x, a, b, 0), R::df(x, a, b, 1));
@@ -158,6 +208,13 @@ NumericVector runit_qf( double x, double a, double b ){
                                           R::qf(x, a, b, 0, 0), R::qf(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rf( double a, double b ){
+    return NumericVector::create(R::rf(a, b), R::rf(a, b));
+}
+
+// -------------------  Student t Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dt( double x, double a ){
 	     return NumericVector::create(R::dt(x, a, 0), R::dt(x, a, 1));
@@ -175,6 +232,13 @@ NumericVector runit_qt( double x, double a ){
                                           R::qt(x, a, 0, 0), R::qt(log(x), a, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rt( double a ){
+    return NumericVector::create(R::rt(a), R::rt(a));
+}
+
+// -------------------  Binomial Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dbinom( double x, double a, double b ){
 	     return NumericVector::create(R::dbinom(x, a, b, 0), R::dbinom(x, a, b, 1));
@@ -192,6 +256,13 @@ NumericVector runit_qbinom( double x, double a, double b ){
                                           R::qbinom(x, a, b, 0, 0), R::qbinom(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rbinom( double a, double b ){
+    return NumericVector::create(R::rbinom(a, b), R::rbinom(a, b));
+}
+
+// -------------------  Cauchy Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dcauchy( double x, double a, double b ){
 	     return NumericVector::create(R::dcauchy(x, a, b, 0), R::dcauchy(x, a, b, 1));
@@ -209,6 +280,13 @@ NumericVector runit_qcauchy( double x, double a, double b ){
                                           R::qcauchy(x, a, b, 0, 0), R::qcauchy(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rcauchy( double a, double b ){
+    return NumericVector::create(R::rcauchy(a, b), R::rcauchy(a, b));
+}
+
+// -------------------  Exponential Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dexp( double x, double a ){
 	     return NumericVector::create(R::dexp(x, a, 0), R::dexp(x, a, 1));
@@ -226,6 +304,13 @@ NumericVector runit_qexp( double x, double a ){
                                           R::qexp(x, a, 0, 0), R::qexp(log(x), a, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rexp( double a ){
+    return NumericVector::create(R::rexp(a), R::rexp(a));
+}
+
+// -------------------  Geometric Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dgeom( double x, double a ){
 	     return NumericVector::create(R::dgeom(x, a, 0), R::dgeom(x, a, 1));
@@ -243,6 +328,13 @@ NumericVector runit_qgeom( double x, double a ){
                                           R::qgeom(x, a, 0, 0), R::qgeom(log(x), a, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rgeom( double a ){
+    return NumericVector::create(R::rgeom(a), R::rgeom(a));
+}
+
+// -------------------  Hypergeometric Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dhyper( double x, double a, double b, double c ){
 	     return NumericVector::create(R::dhyper(x, a, b, c, 0), R::dhyper(x, a, b, c, 1));
@@ -260,6 +352,13 @@ NumericVector runit_qhyper( double x, double a, double b, double c ){
                                           R::qhyper(x, a, b, c, 0, 0), R::qhyper(log(x), a, b, c, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rhyper( double a, double b, double c ){
+    return NumericVector::create(R::rhyper(a, b, c), R::rhyper(a, b, c));
+}
+
+// -------------------  Negative Binomial Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnbinom( double x, double a, double b ){
 	     return NumericVector::create(R::dnbinom(x, a, b, 0), R::dnbinom(x, a, b, 1));
@@ -277,6 +376,13 @@ NumericVector runit_qnbinom( double x, double a, double b ){
                                           R::qnbinom(x, a, b, 0, 0), R::qnbinom(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rnbinom( double a, double b){
+    return NumericVector::create(R::rnbinom(a, b), R::rnbinom(a, b));
+}
+
+// -------------------  Poisson Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dpois( double x, double a ){
 	     return NumericVector::create(R::dpois(x, a, 0), R::dpois(x, a, 1));
@@ -294,6 +400,13 @@ NumericVector runit_qpois( double x, double a ){
                                           R::qpois(x, a, 0, 0), R::qpois(log(x), a, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rpois( double a ){
+    return NumericVector::create(R::rpois(a), R::rpois(a));
+}
+
+// -------------------  Weibull Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dweibull( double x, double a, double b ){
 	     return NumericVector::create(R::dweibull(x, a, b, 0), R::dweibull(x, a, b, 1));
@@ -311,6 +424,13 @@ NumericVector runit_qweibull( double x, double a, double b ){
                                           R::qweibull(x, a, b, 0, 0), R::qweibull(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rweibull( double a, double b ){
+    return NumericVector::create(R::rweibull(a, b), R::rweibull(a, b));
+}
+
+// -------------------  Logistic Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dlogis( double x, double a, double b ){
 	     return NumericVector::create(R::dlogis(x, a, b, 0), R::dlogis(x, a, b, 1));
@@ -328,6 +448,13 @@ NumericVector runit_qlogis( double x, double a, double b ){
                                           R::qlogis(x, a, b, 0, 0), R::qlogis(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rlogis( double a, double b ){
+    return NumericVector::create(R::rlogis(a, b), R::rlogis(a, b));
+}
+
+// -------------------  Non-central Beta Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnbeta( double x, double a, double b, double c ){
 	     return NumericVector::create(R::dnbeta(x, a, b, c, 0), R::dnbeta(x, a, b, c, 1));
@@ -345,6 +472,8 @@ NumericVector runit_qnbeta( double x, double a, double b, double c ){
                                           R::qnbeta(x, a, b, c, 0, 0), R::qnbeta(log(x), a, b, c, 0, 1));
 }
 
+// -------------------  Non-central F Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnf( double x, double a, double b, double c ){
 	     return NumericVector::create(R::dnf(x, a, b, c, 0), R::dnf(x, a, b, c, 1));
@@ -362,6 +491,8 @@ NumericVector runit_qnf( double x, double a, double b, double c ){
                                           R::qnf(x, a, b, c, 0, 0), R::qnf(log(x), a, b, c, 0, 1));
 }
 
+// -------------------  Non-central Student t Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dnt( double x, double a, double b ){
 	     return NumericVector::create(R::dnt(x, a, b, 0), R::dnt(x, a, b, 1));
@@ -379,6 +510,8 @@ NumericVector runit_qnt( double x, double a, double b ){
                                           R::qnt(x, a, b, 0, 0), R::qnt(log(x), a, b, 0, 1));
 }
 
+// -------------------  Wilcoxon Rank Sum Statistic Distribution
+
 // [[Rcpp::export]]
 NumericVector runit_dwilcox( double x, double a, double b ){
 	     return NumericVector::create(R::dwilcox(x, a, b, 0), R::dwilcox(x, a, b, 1));
@@ -396,3 +529,7 @@ NumericVector runit_qwilcox( double x, double a, double b ){
                                           R::qwilcox(x, a, b, 0, 0), R::qwilcox(log(x), a, b, 0, 1));
 }
 
+// [[Rcpp::export]]
+NumericVector runit_rwilcox( double a, double b ){
+    return NumericVector::create(R::rwilcox(a, b), R::rwilcox(a, b));
+}