tests/testthat/test-simulate_TL.R

In RLumModel: Solving Ordinary Differential Equations to Understand Luminescence

parms <- RLumModel:::.set_pars("Bailey2001")
n <- parms$n$n
test_simulate_TL <- RLumModel:::.simulate_TL(
  temp_begin = 20,
  temp_end = 50,
  heating_rate = 5,
  n = n,
  parms = parms)


test_that("check output",{
  skip_on_cran()
  local_edition(3)

  expect_equal(class(test_simulate_TL)[1], "RLum.Results")

  expect_equal(length(test_simulate_TL$n), length(parms$N) + 2)

  expect_equal(test_simulate_TL$temp, 50)

  ##check concentrations
  expect_equal(length(test_simulate_TL$concentrations), length(parms$N) + 2)

  expect_equal(class(test_simulate_TL$concentrations), "list")

  expect_equal(class(test_simulate_TL$concentrations[[1]])[1], "RLum.Data.Curve")

})

test_that("test controlled crash conditions", {
  skip_on_cran()
  local_edition(3)

  parms <- RLumModel:::.set_pars("Bailey2001")
  n <- parms$n$n

  expect_error(
    RLumModel:::.simulate_TL(
      temp_begin = 10,
      temp_end = 20,
      heating_rate = -1,
      n = n,
      parms = parms),
    regexp = ".simulate_TL()] Heatingrate has the wrong algebraic sign!",
    fixed = TRUE)

  expect_error(
    RLumModel:::.simulate_TL(
      temp_begin = -274,
      temp_end = 20,
      heating_rate = 10,
      n = n,
      parms = parms),
    regexp = "[.simulate_TL()] Argument 'temp' has to be > 0 K!",
    fixed = TRUE)

  expect_error(
    RLumModel:::.simulate_TL(
      temp_begin = 10,
      temp_end = -20,
      heating_rate = -2,
      n = n,
      parms = parms),
    regexp = "[.simulate_TL()] Argument 'heating_rate' has to be a positive number!",
    fixed = TRUE)

})