#!/usr/bin/env r # -*- mode: R; ess-indent-level: 4; tab-width: 4; indent-tabs-mode: nil; -* # # Copyright (C) 2012 - 2016 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 . .runThisTest <- Sys.getenv("RunAllRcppTests") == "yes" if (FALSE && .runThisTest) { #.setUp <- Rcpp:::unitTestSetup("rmath.cpp") sourceCpp("cpp/rmath.cpp") test.rmath.norm <- function() { x <- 0.25 a <- 1.25 b <- 2.50 checkEquals(runit_dnorm(x, a, b), c(dnorm(x, a, b, log=FALSE), dnorm(x, a, b, log=TRUE)), msg = " rmath.dnorm") checkEquals(runit_pnorm(x, a, b), c(pnorm(x, a, b, lower=TRUE, log=FALSE), pnorm(log(x), a, b, lower=TRUE, log=TRUE), pnorm(x, a, b, lower=FALSE, log=FALSE), pnorm(log(x), a, b, lower=FALSE, log=TRUE)), msg = " rmath.pnorm") checkEquals(runit_qnorm(x, a, b), c(qnorm(x, a, b, lower=TRUE, log=FALSE), qnorm(log(x), a, b, lower=TRUE, log=TRUE), qnorm(x, a, b, lower=FALSE, log=FALSE), qnorm(log(x), a, b, lower=FALSE, log=TRUE)), msg = " rmath.qnorm") set.seed(333) rcpp_result <- runit_rnorm(a, b) set.seed(333) rcpp_result_sugar <- runit_rnorm_sugar(a, b) set.seed(333) r_result <- rnorm(5, a, b) checkEquals(rcpp_result, r_result, msg = " rmath.rnorm") checkEquals(rcpp_result_sugar, r_result, msg = " rmath.rnorm.sugar") } }