tests/testthat/test-simulate_CW_OSL.R

In RLumModel: Solving Ordinary Differential Equations to Understand Luminescence

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

  parms <- RLumModel:::.set_pars("Bailey2001")
  n <- parms$n$n
  test_simulate_CW_OSL <- RLumModel:::.simulate_CW_OSL(
    temp = 125,
    duration = 100,
    n = n,
    parms = parms)



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

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

  expect_equal(test_simulate_CW_OSL$temp, 125)

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

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

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

})


test_that("check output of Bailey 2002/2004",{
  skip_on_cran()
  local_edition(3)


  parms <- RLumModel:::.set_pars("Bailey2002")
  n <- parms$n$n
  test_simulate_CW_OSL <- RLumModel:::.simulate_LM_OSL(
    temp = 125,
    duration = 100,
    n = n,
    parms = parms)

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

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

  expect_equal(test_simulate_CW_OSL$temp, 125)

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

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

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

})

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


  expect_error(
    RLumModel:::.simulate_CW_OSL(
      temp = -274,
      duration = 100,
      n = n,
      parms = parms),
    regexp = "\n [.simulate_CW_OSL()] Argument 'temp' has to be > 0 K!", fixed = TRUE)

  expect_error(
    RLumModel:::.simulate_CW_OSL(
      temp = 125,
      duration = -10,
      n = n,
      parms = parms),
    regexp = "\n [.simulate_CW_OSL()] Argument 'duration' has to be a positive number!", fixed = TRUE)


  expect_error(
    RLumModel:::.simulate_CW_OSL(
      temp = 125,
      duration = 100,
      optical_power = -10,
      n = n,
      parms = parms),
    regexp = "\n [.simulate_CW_OSL()] Argument 'optical_power' has to be a positive number!", fixed = TRUE)

})