[Rcpp-commits] r4384 - in pkg/Rcpp: . inst/unitTests inst/unitTests/cpp
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Tue Jul 2 17:41:13 CEST 2013
Author: romain
Date: 2013-07-02 17:41:13 +0200 (Tue, 02 Jul 2013)
New Revision: 4384
Added:
pkg/Rcpp/inst/unitTests/cpp/rmath.cpp
Modified:
pkg/Rcpp/ChangeLog
pkg/Rcpp/inst/unitTests/runit.rmath.R
Log:
using sourceCpp in runit.rmath
Modified: pkg/Rcpp/ChangeLog
===================================================================
--- pkg/Rcpp/ChangeLog 2013-07-02 14:45:24 UTC (rev 4383)
+++ pkg/Rcpp/ChangeLog 2013-07-02 15:41:13 UTC (rev 4384)
@@ -6,6 +6,7 @@
* unitTests/runit.misc.R: using sourceCpp
* unitTests/runit.wrap.R: using sourceCpp
* unitTests/runit.support.R: using sourceCpp
+ * unitTests/runit.rmath.R: using sourceCpp
* unitTests/runit.Vector.R: testing List( int, IntegerVector ) which
eventually uses fill__dispatch
* include/Rcpp/traits/r_type_traits.h: support for as<T&> and as<const T&>
Added: pkg/Rcpp/inst/unitTests/cpp/rmath.cpp
===================================================================
--- pkg/Rcpp/inst/unitTests/cpp/rmath.cpp (rev 0)
+++ pkg/Rcpp/inst/unitTests/cpp/rmath.cpp 2013-07-02 15:41:13 UTC (rev 4384)
@@ -0,0 +1,398 @@
+// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 8 -*-
+//
+// rmath.cpp: Rcpp R/C++ interface class library -- rmath unit tests
+//
+// Copyright (C) 2013 Dirk Eddelbuettel and Romain Francois
+//
+// 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/>.
+
+#include <Rcpp.h>
+using namespace Rcpp ;
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnorm( double x, double a, double b ){
+ return NumericVector::create(R::pnorm(x, a, b, 1, 0), R::pnorm(log(x), a, b, 1, 1),
+ R::pnorm(x, a, b, 0, 0), R::pnorm(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnorm( double x, double a, double b ){
+ return NumericVector::create(R::qnorm(x, a, b, 1, 0), R::qnorm(log(x), a, b, 1, 1),
+ R::qnorm(x, a, b, 0, 0), R::qnorm(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_punif( double x, double a, double b ){
+ return NumericVector::create(R::punif(x, a, b, 1, 0), R::punif(log(x), a, b, 1, 1),
+ R::punif(x, a, b, 0, 0), R::punif(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qunif( double x, double a, double b ){
+ return NumericVector::create(R::qunif(x, a, b, 1, 0), R::qunif(log(x), a, b, 1, 1),
+ R::qunif(x, a, b, 0, 0), R::qunif(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pgamma( double x, double a, double b ){
+ return NumericVector::create(R::pgamma(x, a, b, 1, 0), R::pgamma(log(x), a, b, 1, 1),
+ R::pgamma(x, a, b, 0, 0), R::pgamma(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qgamma( double x, double a, double b ){
+ return NumericVector::create(R::qgamma(x, a, b, 1, 0), R::qgamma(log(x), a, b, 1, 1),
+ R::qgamma(x, a, b, 0, 0), R::qgamma(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pbeta( double x, double a, double b ){
+ return NumericVector::create(R::pbeta(x, a, b, 1, 0), R::pbeta(log(x), a, b, 1, 1),
+ R::pbeta(x, a, b, 0, 0), R::pbeta(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qbeta( double x, double a, double b ){
+ return NumericVector::create(R::qbeta(x, a, b, 1, 0), R::qbeta(log(x), a, b, 1, 1),
+ R::qbeta(x, a, b, 0, 0), R::qbeta(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_plnorm( double x, double a, double b ){
+ return NumericVector::create(R::plnorm(x, a, b, 1, 0), R::plnorm(log(x), a, b, 1, 1),
+ R::plnorm(x, a, b, 0, 0), R::plnorm(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qlnorm( double x, double a, double b ){
+ return NumericVector::create(R::qlnorm(x, a, b, 1, 0), R::qlnorm(log(x), a, b, 1, 1),
+ R::qlnorm(x, a, b, 0, 0), R::qlnorm(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_dchisq( double x, double a ){
+ return NumericVector::create(R::dchisq(x, a, 0), R::dchisq(x, a, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pchisq( double x, double a ){
+ return NumericVector::create(R::pchisq(x, a, 1, 0), R::pchisq(log(x), a, 1, 1),
+ R::pchisq(x, a, 0, 0), R::pchisq(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qchisq( double x, double a ){
+ return NumericVector::create(R::qchisq(x, a, 1, 0), R::qchisq(log(x), a, 1, 1),
+ R::qchisq(x, a, 0, 0), R::qchisq(log(x), a, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnchisq( double x, double a, double b ){
+ return NumericVector::create(R::pnchisq(x, a, b, 1, 0), R::pnchisq(log(x), a, b, 1, 1),
+ R::pnchisq(x, a, b, 0, 0), R::pnchisq(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnchisq( double x, double a, double b ){
+ return NumericVector::create(R::qnchisq(x, a, b, 1, 0), R::qnchisq(log(x), a, b, 1, 1),
+ R::qnchisq(x, a, b, 0, 0), R::qnchisq(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pf( double x, double a, double b ){
+ return NumericVector::create(R::pf(x, a, b, 1, 0), R::pf(log(x), a, b, 1, 1),
+ R::pf(x, a, b, 0, 0), R::pf(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qf( double x, double a, double b ){
+ return NumericVector::create(R::qf(x, a, b, 1, 0), R::qf(log(x), a, b, 1, 1),
+ R::qf(x, a, b, 0, 0), R::qf(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_dt( double x, double a ){
+ return NumericVector::create(R::dt(x, a, 0), R::dt(x, a, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pt( double x, double a ){
+ return NumericVector::create(R::pt(x, a, 1, 0), R::pt(log(x), a, 1, 1),
+ R::pt(x, a, 0, 0), R::pt(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qt( double x, double a ){
+ return NumericVector::create(R::qt(x, a, 1, 0), R::qt(log(x), a, 1, 1),
+ R::qt(x, a, 0, 0), R::qt(log(x), a, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pbinom( double x, double a, double b ){
+ return NumericVector::create(R::pbinom(x, a, b, 1, 0), R::pbinom(log(x), a, b, 1, 1),
+ R::pbinom(x, a, b, 0, 0), R::pbinom(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qbinom( double x, double a, double b ){
+ return NumericVector::create(R::qbinom(x, a, b, 1, 0), R::qbinom(log(x), a, b, 1, 1),
+ R::qbinom(x, a, b, 0, 0), R::qbinom(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pcauchy( double x, double a, double b ){
+ return NumericVector::create(R::pcauchy(x, a, b, 1, 0), R::pcauchy(log(x), a, b, 1, 1),
+ R::pcauchy(x, a, b, 0, 0), R::pcauchy(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qcauchy( double x, double a, double b ){
+ return NumericVector::create(R::qcauchy(x, a, b, 1, 0), R::qcauchy(log(x), a, b, 1, 1),
+ R::qcauchy(x, a, b, 0, 0), R::qcauchy(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_dexp( double x, double a ){
+ return NumericVector::create(R::dexp(x, a, 0), R::dexp(x, a, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pexp( double x, double a ){
+ return NumericVector::create(R::pexp(x, a, 1, 0), R::pexp(log(x), a, 1, 1),
+ R::pexp(x, a, 0, 0), R::pexp(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qexp( double x, double a ){
+ return NumericVector::create(R::qexp(x, a, 1, 0), R::qexp(log(x), a, 1, 1),
+ R::qexp(x, a, 0, 0), R::qexp(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_dgeom( double x, double a ){
+ return NumericVector::create(R::dgeom(x, a, 0), R::dgeom(x, a, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pgeom( double x, double a ){
+ return NumericVector::create(R::pgeom(x, a, 1, 0), R::pgeom(log(x), a, 1, 1),
+ R::pgeom(x, a, 0, 0), R::pgeom(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qgeom( double x, double a ){
+ return NumericVector::create(R::qgeom(x, a, 1, 0), R::qgeom(log(x), a, 1, 1),
+ R::qgeom(x, a, 0, 0), R::qgeom(log(x), a, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_phyper( double x, double a, double b, double c ){
+ return NumericVector::create(R::phyper(x, a, b, c, 1, 0), R::phyper(log(x), a, b, c, 1, 1),
+ R::phyper(x, a, b, c, 0, 0), R::phyper(log(x), a, b, c, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qhyper( double x, double a, double b, double c ){
+ return NumericVector::create(R::qhyper(x, a, b, c, 1, 0), R::qhyper(log(x), a, b, c, 1, 1),
+ R::qhyper(x, a, b, c, 0, 0), R::qhyper(log(x), a, b, c, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnbinom( double x, double a, double b ){
+ return NumericVector::create(R::pnbinom(x, a, b, 1, 0), R::pnbinom(log(x), a, b, 1, 1),
+ R::pnbinom(x, a, b, 0, 0), R::pnbinom(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnbinom( double x, double a, double b ){
+ return NumericVector::create(R::qnbinom(x, a, b, 1, 0), R::qnbinom(log(x), a, b, 1, 1),
+ R::qnbinom(x, a, b, 0, 0), R::qnbinom(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_dpois( double x, double a ){
+ return NumericVector::create(R::dpois(x, a, 0), R::dpois(x, a, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_ppois( double x, double a ){
+ return NumericVector::create(R::ppois(x, a, 1, 0), R::ppois(log(x), a, 1, 1),
+ R::ppois(x, a, 0, 0), R::ppois(log(x), a, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qpois( double x, double a ){
+ return NumericVector::create(R::qpois(x, a, 1, 0), R::qpois(log(x), a, 1, 1),
+ R::qpois(x, a, 0, 0), R::qpois(log(x), a, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pweibull( double x, double a, double b ){
+ return NumericVector::create(R::pweibull(x, a, b, 1, 0), R::pweibull(log(x), a, b, 1, 1),
+ R::pweibull(x, a, b, 0, 0), R::pweibull(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qweibull( double x, double a, double b ){
+ return NumericVector::create(R::qweibull(x, a, b, 1, 0), R::qweibull(log(x), a, b, 1, 1),
+ R::qweibull(x, a, b, 0, 0), R::qweibull(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_plogis( double x, double a, double b ){
+ return NumericVector::create(R::plogis(x, a, b, 1, 0), R::plogis(log(x), a, b, 1, 1),
+ R::plogis(x, a, b, 0, 0), R::plogis(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qlogis( double x, double a, double b ){
+ return NumericVector::create(R::qlogis(x, a, b, 1, 0), R::qlogis(log(x), a, b, 1, 1),
+ R::qlogis(x, a, b, 0, 0), R::qlogis(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnbeta( double x, double a, double b, double c ){
+ return NumericVector::create(R::pnbeta(x, a, b, c, 1, 0), R::pnbeta(log(x), a, b, c, 1, 1),
+ R::pnbeta(x, a, b, c, 0, 0), R::pnbeta(log(x), a, b, c, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnbeta( double x, double a, double b, double c ){
+ return NumericVector::create(R::qnbeta(x, a, b, c, 1, 0), R::qnbeta(log(x), a, b, c, 1, 1),
+ R::qnbeta(x, a, b, c, 0, 0), R::qnbeta(log(x), a, b, c, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnf( double x, double a, double b, double c ){
+ return NumericVector::create(R::pnf(x, a, b, c, 1, 0), R::pnf(log(x), a, b, c, 1, 1),
+ R::pnf(x, a, b, c, 0, 0), R::pnf(log(x), a, b, c, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnf( double x, double a, double b, double c ){
+ return NumericVector::create(R::qnf(x, a, b, c, 1, 0), R::qnf(log(x), a, b, c, 1, 1),
+ R::qnf(x, a, b, c, 0, 0), R::qnf(log(x), a, b, c, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pnt( double x, double a, double b ){
+ return NumericVector::create(R::pnt(x, a, b, 1, 0), R::pnt(log(x), a, b, 1, 1),
+ R::pnt(x, a, b, 0, 0), R::pnt(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qnt( double x, double a, double b ){
+ return NumericVector::create(R::qnt(x, a, b, 1, 0), R::qnt(log(x), a, b, 1, 1),
+ R::qnt(x, a, b, 0, 0), R::qnt(log(x), a, b, 0, 1));
+}
+
+// [[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));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_pwilcox( double x, double a, double b ){
+ return NumericVector::create(R::pwilcox(x, a, b, 1, 0), R::pwilcox(log(x), a, b, 1, 1),
+ R::pwilcox(x, a, b, 0, 0), R::pwilcox(log(x), a, b, 0, 1));
+}
+
+// [[Rcpp::export]]
+NumericVector runit_qwilcox( double x, double a, double b ){
+ return NumericVector::create(R::qwilcox(x, a, b, 1, 0), R::qwilcox(log(x), a, b, 1, 1),
+ R::qwilcox(x, a, b, 0, 0), R::qwilcox(log(x), a, b, 0, 1));
+}
+
Modified: pkg/Rcpp/inst/unitTests/runit.rmath.R
===================================================================
--- pkg/Rcpp/inst/unitTests/runit.rmath.R 2013-07-02 14:45:24 UTC (rev 4383)
+++ pkg/Rcpp/inst/unitTests/runit.rmath.R 2013-07-02 15:41:13 UTC (rev 4384)
@@ -1,7 +1,7 @@
#!/usr/bin/r -t
# -*- mode: R; ess-indent-level: 4; tab-width: 4; indent-tabs-mode: nil; -*
#
-# Copyright (C) 2012 Dirk Eddelbuettel and Romain Francois
+# Copyright (C) 2012 - 2013 Dirk Eddelbuettel and Romain Francois
#
# This file is part of Rcpp.
#
@@ -22,318 +22,22 @@
if (.runThisTest) {
-definitions <- function() {
- list("runit_dnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dnorm(x, a, b, 0), R::dnorm(x, a, b, 1));')
- ,"runit_pnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pnorm(x, a, b, 1, 0), R::pnorm(log(x), a, b, 1, 1),
- R::pnorm(x, a, b, 0, 0), R::pnorm(log(x), a, b, 0, 1));')
- ,"runit_qnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qnorm(x, a, b, 1, 0), R::qnorm(log(x), a, b, 1, 1),
- R::qnorm(x, a, b, 0, 0), R::qnorm(log(x), a, b, 0, 1));')
+.setUp <- Rcpp:::unit_test_setup( "rmath.cpp" )
-
- ,"runit_dunif" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dunif(x, a, b, 0), R::dunif(x, a, b, 1));')
- ,"runit_punif" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::punif(x, a, b, 1, 0), R::punif(log(x), a, b, 1, 1),
- R::punif(x, a, b, 0, 0), R::punif(log(x), a, b, 0, 1));')
- ,"runit_qunif" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qunif(x, a, b, 1, 0), R::qunif(log(x), a, b, 1, 1),
- R::qunif(x, a, b, 0, 0), R::qunif(log(x), a, b, 0, 1));')
-
-
- ,"runit_dgamma" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dgamma(x, a, b, 0), R::dgamma(x, a, b, 1));')
- ,"runit_pgamma" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pgamma(x, a, b, 1, 0), R::pgamma(log(x), a, b, 1, 1),
- R::pgamma(x, a, b, 0, 0), R::pgamma(log(x), a, b, 0, 1));')
- ,"runit_qgamma" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qgamma(x, a, b, 1, 0), R::qgamma(log(x), a, b, 1, 1),
- R::qgamma(x, a, b, 0, 0), R::qgamma(log(x), a, b, 0, 1));')
-
-
- ,"runit_dbeta" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dbeta(x, a, b, 0), R::dbeta(x, a, b, 1));')
- ,"runit_pbeta" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pbeta(x, a, b, 1, 0), R::pbeta(log(x), a, b, 1, 1),
- R::pbeta(x, a, b, 0, 0), R::pbeta(log(x), a, b, 0, 1));')
- ,"runit_qbeta" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qbeta(x, a, b, 1, 0), R::qbeta(log(x), a, b, 1, 1),
- R::qbeta(x, a, b, 0, 0), R::qbeta(log(x), a, b, 0, 1));')
-
-
- ,"runit_dlnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dlnorm(x, a, b, 0), R::dlnorm(x, a, b, 1));')
- ,"runit_plnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::plnorm(x, a, b, 1, 0), R::plnorm(log(x), a, b, 1, 1),
- R::plnorm(x, a, b, 0, 0), R::plnorm(log(x), a, b, 0, 1));')
- ,"runit_qlnorm" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qlnorm(x, a, b, 1, 0), R::qlnorm(log(x), a, b, 1, 1),
- R::qlnorm(x, a, b, 0, 0), R::qlnorm(log(x), a, b, 0, 1));')
-
-
- ,"runit_dchisq" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::dchisq(x, a, 0), R::dchisq(x, a, 1));')
- ,"runit_pchisq" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::pchisq(x, a, 1, 0), R::pchisq(log(x), a, 1, 1),
- R::pchisq(x, a, 0, 0), R::pchisq(log(x), a, 0, 1));')
- ,"runit_qchisq" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::qchisq(x, a, 1, 0), R::qchisq(log(x), a, 1, 1),
- R::qchisq(x, a, 0, 0), R::qchisq(log(x), a, 0, 1));')
-
-
- ,"runit_dnchisq" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dnchisq(x, a, b, 0), R::dnchisq(x, a, b, 1));')
- ,"runit_pnchisq" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pnchisq(x, a, b, 1, 0), R::pnchisq(log(x), a, b, 1, 1),
- R::pnchisq(x, a, b, 0, 0), R::pnchisq(log(x), a, b, 0, 1));')
- ,"runit_qnchisq" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qnchisq(x, a, b, 1, 0), R::qnchisq(log(x), a, b, 1, 1),
- R::qnchisq(x, a, b, 0, 0), R::qnchisq(log(x), a, b, 0, 1));')
-
-
- ,"runit_df" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::df(x, a, b, 0), R::df(x, a, b, 1));')
- ,"runit_pf" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pf(x, a, b, 1, 0), R::pf(log(x), a, b, 1, 1),
- R::pf(x, a, b, 0, 0), R::pf(log(x), a, b, 0, 1));')
- ,"runit_qf" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qf(x, a, b, 1, 0), R::qf(log(x), a, b, 1, 1),
- R::qf(x, a, b, 0, 0), R::qf(log(x), a, b, 0, 1));')
-
-
- ,"runit_dt" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::dt(x, a, 0), R::dt(x, a, 1));')
- ,"runit_pt" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::pt(x, a, 1, 0), R::pt(log(x), a, 1, 1),
- R::pt(x, a, 0, 0), R::pt(log(x), a, 0, 1));')
- ,"runit_qt" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::qt(x, a, 1, 0), R::qt(log(x), a, 1, 1),
- R::qt(x, a, 0, 0), R::qt(log(x), a, 0, 1));')
-
-
- ,"runit_dbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dbinom(x, a, b, 0), R::dbinom(x, a, b, 1));')
- ,"runit_pbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pbinom(x, a, b, 1, 0), R::pbinom(log(x), a, b, 1, 1),
- R::pbinom(x, a, b, 0, 0), R::pbinom(log(x), a, b, 0, 1));')
- ,"runit_qbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qbinom(x, a, b, 1, 0), R::qbinom(log(x), a, b, 1, 1),
- R::qbinom(x, a, b, 0, 0), R::qbinom(log(x), a, b, 0, 1));')
-
-
- ,"runit_dcauchy" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dcauchy(x, a, b, 0), R::dcauchy(x, a, b, 1));')
- ,"runit_pcauchy" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pcauchy(x, a, b, 1, 0), R::pcauchy(log(x), a, b, 1, 1),
- R::pcauchy(x, a, b, 0, 0), R::pcauchy(log(x), a, b, 0, 1));')
- ,"runit_qcauchy" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qcauchy(x, a, b, 1, 0), R::qcauchy(log(x), a, b, 1, 1),
- R::qcauchy(x, a, b, 0, 0), R::qcauchy(log(x), a, b, 0, 1));')
-
-
- ,"runit_dexp" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::dexp(x, a, 0), R::dexp(x, a, 1));')
- ,"runit_pexp" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::pexp(x, a, 1, 0), R::pexp(log(x), a, 1, 1),
- R::pexp(x, a, 0, 0), R::pexp(log(x), a, 0, 1));')
- ,"runit_qexp" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::qexp(x, a, 1, 0), R::qexp(log(x), a, 1, 1),
- R::qexp(x, a, 0, 0), R::qexp(log(x), a, 0, 1));')
-
-
- ,"runit_dgeom" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::dgeom(x, a, 0), R::dgeom(x, a, 1));')
- ,"runit_pgeom" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::pgeom(x, a, 1, 0), R::pgeom(log(x), a, 1, 1),
- R::pgeom(x, a, 0, 0), R::pgeom(log(x), a, 0, 1));')
- ,"runit_qgeom" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::qgeom(x, a, 1, 0), R::qgeom(log(x), a, 1, 1),
- R::qgeom(x, a, 0, 0), R::qgeom(log(x), a, 0, 1));')
-
-
- ,"runit_dhyper" = list(signature(x_ = "double", a_ = "double", b_ = "double", c_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_), c = as<double>(c_);
- return NumericVector::create(R::dhyper(x, a, b, c, 0), R::dhyper(x, a, b, c, 1));')
- ,"runit_phyper" = list(signature(x_ = "double", a_ = "double", b_ = "double", c_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_), c = as<double>(c_);
- return NumericVector::create(R::phyper(x, a, b, c, 1, 0), R::phyper(log(x), a, b, c, 1, 1),
- R::phyper(x, a, b, c, 0, 0), R::phyper(log(x), a, b, c, 0, 1));')
- ,"runit_qhyper" = list(signature(x_ = "double", a_ = "double", b_ = "double", c_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_), c = as<double>(c_);
- return NumericVector::create(R::qhyper(x, a, b, c, 1, 0), R::qhyper(log(x), a, b, c, 1, 1),
- R::qhyper(x, a, b, c, 0, 0), R::qhyper(log(x), a, b, c, 0, 1));')
-
-
- ,"runit_dnbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::dnbinom(x, a, b, 0), R::dnbinom(x, a, b, 1));')
- ,"runit_pnbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_) ;
- return NumericVector::create(R::pnbinom(x, a, b, 1, 0), R::pnbinom(log(x), a, b, 1, 1),
- R::pnbinom(x, a, b, 0, 0), R::pnbinom(log(x), a, b, 0, 1));')
- ,"runit_qnbinom" = list(signature(x_ = "double", a_ = "double", b_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_), b = as<double>(b_);
- return NumericVector::create(R::qnbinom(x, a, b, 1, 0), R::qnbinom(log(x), a, b, 1, 1),
- R::qnbinom(x, a, b, 0, 0), R::qnbinom(log(x), a, b, 0, 1));')
-
-
- ,"runit_dpois" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::dpois(x, a, 0), R::dpois(x, a, 1));')
- ,"runit_ppois" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
- return NumericVector::create(R::ppois(x, a, 1, 0), R::ppois(log(x), a, 1, 1),
- R::ppois(x, a, 0, 0), R::ppois(log(x), a, 0, 1));')
- ,"runit_qpois" = list(signature(x_ = "double", a_ = "double"), '
- double x = as<double>(x_), a = as<double>(a_);
[TRUNCATED]
To get the complete diff run:
svnlook diff /svnroot/rcpp -r 4384
More information about the Rcpp-commits
mailing list