Nothing
test_that("check class and output", {
skip_on_cran()
local_edition(3)
output <- model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(
OSL = c(temp = 20, duration = 1, optical_power = 100)),
plot = FALSE,
verbose = FALSE)
expect_equal(length(output), 12)
expect_equal(is(output), c("RLum.Analysis", "RLum"))
expect_silent(model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(
OSL = c(20, 1, 100)),
plot = TRUE,
verbose = FALSE))
})
test_that("check simulate_sample_history = TRUE", {
skip_on_cran()
local_edition(3)
output <- model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(
RF = c(20, 1, 1)),
plot = FALSE,
simulate_sample_history = TRUE,
show_structure = TRUE,
verbose = FALSE)
expect_equal(length(output), 12)
expect_equal(is(output), c("RLum.Analysis", "RLum"))
})
test_that("check SAR sequence", {
skip_on_cran()
local_edition(3)
sequence <- list(
RegDose = c(0,10,20),
TestDose = 5,
PH = 240,
CH = 200,
OSL_temp = 125,
OSL_duration = 70)
expect_s4_class(model_LuminescenceSignals(
sequence = sequence,
model = "Pagonis2007",
plot = FALSE,
verbose = FALSE
), "RLum.Analysis")
})
test_that("check DRT sequence", {
skip_on_cran()
local_edition(3)
sequence <- list(
Irr_2recover = c(20),
RegDose = c(0,10,20),
TestDose = 5,
PH = 240,
CH = 200,
OSL_temp = 125,
OSL_duration = 70)
expect_s4_class(model_LuminescenceSignals(
sequence = sequence,
model = "Pagonis2007",
plot = FALSE,
verbose = FALSE
), "RLum.Analysis")
})
test_that("check Risoe SEQ", {
skip_on_cran()
local_edition(3)
path <- system.file("extdata", "example_SAR_cycle.SEQ", package="RLumModel")
sequence <- read_SEQ2R(file = path,
txtProgressBar = FALSE)
expect_s4_class(model_LuminescenceSignals(
sequence = sequence,
model = "Bailey2001",
plot = FALSE,
verbose = FALSE
), "RLum.Analysis")
expect_s4_class(model_LuminescenceSignals(
sequence = path,
model = "Bailey2001",
plot = FALSE,
verbose = FALSE
), "RLum.Analysis")
## Peng 2022
expect_s4_class(model_LuminescenceSignals(
model = "Peng2022",
sequence = list(
OSL = c(temp = 20, duration = 1, optical_power = 100)),
plot = FALSE,
verbose = FALSE), "RLum.Analysis")
})
test_that("check custom models", {
skip_on_cran()
local_edition(3)
own_parameters <- list(
N = c(2e15, 2e15, 1e17, 2.4e16),
E = c(0, 0, 0, 0),
s = c(0, 0, 0, 0),
A = c(2e-8, 2e-9, 4e-9, 1e-8),
B = c(0, 0, 5e-11, 4e-8),
Th = c(0, 0),
E_th = c(0, 0),
k_B = 8.617e-5,
W = 0.64,
K = 2.8e7,
model = "customized",
R = 1.7e15
)
own_state_parameters <- c(0, 0, 0, 9.4e15)
own_start_temperature <- 10
sequence <- list(
RF = c(20, 1, 1))
## test routine
temp <- expect_s4_class(model_LuminescenceSignals(
model = "customized",
sequence = sequence,
own_parameters = own_parameters,
own_state_parameters = own_state_parameters,
own_start_temperature = own_start_temperature,
plot = FALSE,
verbose = FALSE), "RLum.Analysis")
## test state case where state parameters are an object
expect_s4_class(model_LuminescenceSignals(
model = "customized",
sequence = sequence,
own_parameters = .set_pars("Bailey2001"),
own_state_parameters = .set_pars("Bailey2001")$n,
plot = FALSE,
verbose = FALSE), "RLum.Analysis")
expect_warning(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(TL = c(20, 30, 10)),
own_parameters = list(),
plot = FALSE,
verbose = FALSE
),
regexp = "[model_LuminescenceSignals()] Argument 'own_parameters' set, but argument 'model' not set to 'customized'. Used 'Bailey2001' as argument for 'model'.",
fixed = TRUE
)
expect_warning(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(TL = c(20, 30, 10)),
own_state_parameters = 10,
plot = FALSE,
verbose = FALSE
),
regexp = "[model_LuminescenceSignals()] Argument 'own_sate_parameters' set, but argument 'model' not set to 'customized'. Ignored argument 'own_state_parameters'.",
fixed = TRUE
)
expect_warning(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(TL = c(20, 30, 10)),
own_start_temperature = 10,
plot = FALSE,
verbose = FALSE
),
regexp = "[model_LuminescenceSignals()] Argument 'own_start_temperature' set, but argument 'model' not set to 'customized'. Ignored argument 'own_start_temperature'.",
fixed = TRUE
)
})
test_that("test controlled crash conditions", {
skip_on_cran()
local_edition(3)
expect_error(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(
OSL = c(20, 1, 100)),
plot = FALSE,
verbose = FALSE,
lab.dose_rate = -1),
regexp = "[model_LuminescenceSignals()] lab.dose_rate has to be a positive number!",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = "test.jpg",
plot = FALSE,
verbose = FALSE),
regexp = "[model_LuminescenceSignals()] Argument 'sequence' is not a *.SEQ file!",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
sequence = list(
OSL = c(20, 1, 100))),
"Argument 'model' not given!")
expect_error(
model_LuminescenceSignals(
model = "Bailey2001"),
"Argument 'sequence' not given!")
expect_error(
model_LuminescenceSignals(
model = "Bailey2000",
sequence = list(
OSL = c(20,1,100))),
regexp = "[model_LuminescenceSignals()] Model not supported. Supported models are: Bailey2001, Bailey2004, Pagonis2008, Pagonis2007, Bailey2002, Friedrich2017, Friedrich2018, Peng2022, customized, customised",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = list(
OSL = c("20",1,100))),
regexp = "[model_LuminescenceSignals()] Sequence comprises non-numeric arguments!",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
model = "Bailey2001",
sequence = matrix(NA)),
regexp = "[model_LuminescenceSignals()] Sequence has to be of class list or a *.seq file",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
model = "customized",
sequence = list(TL = c(20,100,5))),
regexp = "[model_LuminescenceSignals()] Argument 'model' set to 'customized', but no own parameters are given!",
fixed = TRUE)
expect_error(
model_LuminescenceSignals(
model = "customized",
sequence = list(TLL = c(20,100,5))),
regexp = "[model_LuminescenceSignals()] Unknow sequence arguments: Allowed arguments are: IRR, PH, CH, TL, OSL, PAUSE, LM_OSL, RL, RF, ILL, RF_heating",
fixed = TRUE)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.