CRAN Package Check Results for Package expint

Last updated on 2024-03-28 19:48:30 CET.

Flavor Version Tinstall Tcheck Ttotal Status Flags
r-devel-linux-x86_64-debian-clang 0.1-8 3.12 32.75 35.87 OK
r-devel-linux-x86_64-debian-gcc 0.1-8 2.31 25.79 28.10 OK
r-devel-linux-x86_64-fedora-clang 0.1-8 46.20 OK
r-devel-linux-x86_64-fedora-gcc 0.1-8 56.89 OK
r-devel-windows-x86_64 0.1-8 6.00 230.00 236.00 ERROR
r-patched-linux-x86_64 0.1-8 2.57 32.38 34.95 OK
r-release-linux-x86_64 0.1-8 3.17 31.59 34.76 OK
r-release-macos-arm64 0.1-8 27.00 OK
r-release-macos-x86_64 0.1-8 37.00 OK
r-release-windows-x86_64 0.1-8 9.00 65.00 74.00 OK
r-oldrel-macos-arm64 0.1-8 21.00 OK
r-oldrel-windows-x86_64 0.1-8 13.00 85.00 98.00 OK

Check Details

Version: 0.1-8
Check: tests
Result: ERROR Running 'expint-tests.R' [0s] Running 'gammainc-tests.R' [164s] Running the tests in 'tests/gammainc-tests.R' failed. Complete output: > ### == expint: Exponential Integral and Incomplete Gamma Function == > ### > ### Tests for the incomplete gamma function > ### > ### G(a,x) = int_x^infty t^{a-1} exp(-t) dt > ### > ### for a *real* and x >= 0. > ### > ### AUTHOR: Vincent Goulet <vincent.goulet@act.ulaval.ca> > > ## Load the package > library(expint) > > ## a > 0; direct link to the standard incomplete gamma function > x <- c(0.2, 2.5, 5, 8, 10) > a <- 1.2 > stopifnot(exprs = { + identical(gammainc(a, x), + gamma(a) * pgamma(x, a, 1, lower = FALSE)) + }) > > ## a = 0; direct link to the exponential integral > x <- c(0.2, 2.5, 5, 8, 10) > a <- 0 > stopifnot(exprs = { + identical(gammainc(a, x), expint(x)) + identical(gammainc(a, x), expint_E1(x)) + }) > > ## a < 0; compare to the recursive formula > x <- c(0.2, 2.5, 5, 8, 10) > a <- c(-0.25, -1.2, -2) > stopifnot(exprs = { + all.equal(gammainc(a[1], x), + -(x^a[1] * exp(-x))/a[1] + + gamma(a[1] + 1) * pgamma(x, a[1] + 1, 1, lower = FALSE)/a[1]) + all.equal(gammainc(a[2], x), + -(x^a[2] * exp(-x))/a[2] + + (-(x^(a[2] + 1) * exp(-x))/(a[2] + 1) + + gamma(a[2] + 2) * pgamma(x, a[2] + 2, 1, lower = FALSE)/(a[2] + 1))/a[2]) + all.equal( + gammainc(a[3], x), + -(x^a[3] * exp(-x))/a[3] + + (-(x^(a[3] + 1) * exp(-x))/(a[3] + 1) + expint_E1(x)/(a[3] + 1))/a[3]) + }) > > proc.time() user system elapsed 0.20 0.01 0.20 Flavor: r-devel-windows-x86_64