tests/testthat/test-external-generator.R

In dqrng: Fast Pseudo Random Number Generators

context("external-generators")

Rcpp::sourceCpp("cpp/external-generator.cpp")

n <- 1e2L
rate <- 0.4
use_seed <- 1623

test_that("external RNG (normal)", {
  dqset.seed(use_seed)
  expected <- dqrexp(n, rate)
  dqset.seed(use_seed)
  actual <- dqrexp_extrng(n, rate)
  expect_equal(actual, expected)

  dqset.seed(use_seed)
  expected2 <- expected
  actual2 <- sapply(1:n, function(x) dqrexp_extrng(1, rate))
  expect_equal(actual2, expected2)
})

test_that("external RNG (normal, Xoshiro256++)", {
  dqrng::dqRNGkind("Xoshiro256++")
  dqset.seed(use_seed)
  expected <- dqrexp(n, rate)
  dqset.seed(use_seed)
  actual <- dqrexp_extrng(n, rate)
  expect_equal(actual, expected)

  dqset.seed(use_seed)
  expected2 <- expected
  actual2 <- sapply(1:n, function(x) dqrexp_extrng(1, rate))
  expect_equal(actual2, expected2)
})

test_that("external RNG (parallel, Threefry)", {
  cl <- parallel::makeCluster(2)
  expected3 <- parallel::clusterApply(cl, 1:8, function(stream, seed, N, rate) {
    dqrng::dqRNGkind("Threefry")
    dqrng::dqset.seed(seed, stream)
    dqrng::dqrexp(N, rate)
  }, use_seed, 1e6L, rate)
  parallel::stopCluster(cl)

  cl <- parallel::makeCluster(2)
  actual3 <- parallel::clusterApply(cl, 1:8, function(stream, seed, N, rate) {
    Rcpp::sourceCpp("cpp/external-generator.cpp") ## must be recompiled
    dqrng::dqRNGkind("Threefry")
    dqrng::dqset.seed(seed, stream)
    dqrexp_extrng(N, rate)
  }, use_seed, 1e6L, rate)
  parallel::stopCluster(cl)

  expect_equal(actual3, expected3)
})

test_that("cloned external default RNG gives different result only when a different stream is selected", {
  dqrng::dqRNGkind("default")
  dqset.seed(use_seed)
  expect_true(cloned_calls(stream = 0))
  expect_false(cloned_calls(stream = 1))
})

test_that("cloned external Xoshiro256++ gives different result only when a different stream is selected", {
  dqrng::dqRNGkind("Xoshiro256++")
  dqset.seed(use_seed)
  expect_true(cloned_calls(stream = 0))
  expect_false(cloned_calls(stream = 1))
  dqrng::dqRNGkind("default")
})

test_that("cloned external PCG64 gives different result", {
  skip_if(powerpc_apple)
  dqrng::dqRNGkind("PCG64")
  dqset.seed(use_seed)
  expect_true(cloned_calls(stream = 0))
  expect_false(cloned_calls(stream = 1))
  dqrng::dqRNGkind("default")
})

test_that("cloned external Threefry gives different result only when a different stream is selected", {
  dqrng::dqRNGkind("Threefry")
  dqset.seed(use_seed)
  expect_true(cloned_calls(stream = 0))
  expect_false(cloned_calls(stream = 1))
  dqrng::dqRNGkind("default")
})