tests/testthat/test-realms.R

In rice: Radiocarbon Equations

test_that("calBPtoBCAD returns expected BC/AD value", {
  result <- calBPtoBCAD(2000)
  expect_length(result, 1)
  expect_equal(result, -50, tolerance = 1e-5)
})

test_that("calBPtoBCAD returns expected BC/AD values (multiple)", {
  result <- calBPtoBCAD((-10):10)
  expect_length(result, 21)
  expect_equal(result, 1960:1940, tolerance = 1e-5)
})

test_that("calBPtob2k returns expected b2k value", {
  result <- calBPtob2k(-50)
  expect_length(result, 1)
  expect_equal(result, 0, tolerance = 1e-5)
})

test_that("calBPtoC14 returns expected C14 values", {
  result <- as.vector(calBPtoC14(100))
  expect_length(result, 2)
  expect_equal(result, c(124,10), tolerance = 1e-5)
})

test_that("calBPtoF14C returns expected F14C values", {
  result <- as.vector(calBPtoF14C(100))
  expect_length(result, 2)
  expect_equal(result, c(0.98468220, 0.00122656), tolerance = 1e-5)
})

test_that("calBPtopMC returns expected pMC values", {
  result <- as.vector(calBPtopMC(100))
  expect_length(result, 2)
  expect_equal(result, c(98.468220, 0.122656), tolerance = 1e-5)
})

test_that("calBPtoD14C returns expected D14C values", {
  result <- as.vector(calBPtoD14C(1000))
  expect_length(result, 2)
  expect_equal(result, c(-19.025611, 1.588819), tolerance = 1e-5)
})

###

test_that("BCADtocalBP returns expected calBP value", {
  result <- as.vector(BCADtocalBP(2000))
  expect_length(result, 1)
  expect_equal(result, c(-50), tolerance = 1e-5)
})

test_that("BCADtob2k returns expected b2k value", {
  result <- as.vector(BCADtob2k(2000))
  expect_length(result, 1)
  expect_equal(result, c(0), tolerance = 1e-5)
})

test_that("BCADtoC14 returns expected C14 age", {
  result <- as.vector(BCADtoC14(1000))
  expect_length(result, 2)
  expect_equal(result, c(1028, 10), tolerance = 1e-5)
})

test_that("BCADtoF14C returns expected F14C value", {
  result <- as.vector(BCADtoF14C(1000))
  expect_length(result, 2)
  expect_equal(result, c(0.8798779, 0.00109601), tolerance = 1e-5)
})

test_that("BCADtopMC returns expected pMC value", {
  result <- as.vector(BCADtopMC(1000))
  expect_length(result, 2)
  expect_equal(result, c(87.98779, 0.109601), tolerance = 1e-5)
})

test_that("BCADtoD14C returns expected D14C value", {
  result <- as.vector(BCADtoD14C(1000))
  expect_length(result, 2)
  expect_equal(result, c(-6.984748, 1.236938), tolerance = 1e-5)
})

###

test_that("b2ktocalBP returns expected calBP value", {
  result <- as.vector(b2ktocalBP(5000))
  expect_length(result, 1)
  expect_equal(result, c(4950), tolerance = 1e-5)
})

test_that("b2ktoBCAD returns expected BC/AD value", {
  result <- as.vector(b2ktoBCAD(5000))
  expect_length(result, 1)
  expect_equal(result, c(-3000), tolerance = 1e-5)
})

test_that("b2ktoC14 returns expected C14 age", {
  result <- as.vector(b2ktoC14(5000))
  expect_length(result, 2)
  expect_equal(result, c(4364, 16), tolerance = 1e-5)
})

test_that("b2ktoF14C returns expected F14C value", {
  result <- as.vector(b2ktoF14C(5000))
  expect_length(result, 2)
  expect_equal(result, c(0.58085213, 0.00115808), tolerance = 1e-5)
})

test_that("b2ktopMC returns expected pMC value", {
  result <- as.vector(b2ktopMC(5000))
  expect_length(result, 2)
  expect_equal(result, c(58.085213, 0.115808), tolerance = 1e-5)
})

test_that("b2ktoD14C returns expected D14C value", {
  result <- as.vector(b2ktoD14C(5000))
  expect_length(result, 2)
  expect_equal(result, c(63.489261, 2.120343), tolerance = 1e-5)
})

###

test_that("C14tocalBP returns expected cal BP values", {
  result <- as.vector(C14tocalBP(800))
  expect_length(result, 3)
  expect_equal(result, c(690.0040, 710.0000, 717.9961), tolerance = 1e-5)
})

test_that("C14toF14C returns expected F14C value", {
  result <- as.vector(C14toF14C(8000))
  expect_length(result, 1)
  expect_equal(result, c(0.3693938), tolerance = 1e-5)
})

test_that("C14toF14C returns expected F14C values", {
  result <- as.vector(C14toF14C(8000, 40))
  expect_length(result, 2)
  expect_equal(result, c(0.3693938, 0.00184397), tolerance = 1e-5)
})

test_that("C14toF14C returns expected F14C values", {
  result <- C14toF14C(c(8000, 8500), c(40,35))
  expect_length(as.vector(result), 4)
  expect_equal(result, cbind(c(0.3693938, 0.3471025), c(0.00184397, 0.0015163)), tolerance = 1e-5)
})

test_that("F14CtoC14 returns expected C14 age", {
  result <- as.vector(F14CtoC14(.5, .005))
  expect_length(result, 2)
  expect_equal(result, c(5568.051301, 79.931008), tolerance = 1e-5)
})

test_that("F14CtoD14C returns expected D14C values", {
  result <- as.vector(F14CtoD14C(.5, .005, t=5500))
  expect_length(result, 2)
  expect_equal(result, c(-27.467151, 9.725328), tolerance = 1e-5)
})

test_that("D14CtoC14 returns expected C14 ages", {
  result <- as.vector(D14CtoC14(20, .05, t=5500))
  expect_length(result, 2)
  expect_equal(result, c(5185.2462885, 0.3937649), tolerance = 1e-5)
})

test_that("calBPtoF14C returns expected F14C values", {
  # Example test input
  calBP <- 2000

  # Run the function with known input
  result <- calBPtoF14C(calBP)

  # Check that the result is a numeric vector of length 2 (F14C, uncertainty)
  expect_length(result, 2)

  # Check that the values are within expected bounds
  expect_gt(result[1], 0)     # F14C > 0

  # check with expected result
  expected <- matrix(c(0.7728371, 0.00134808), nrow=1)
  expect_equal(round(result,4), round(expected,4), tolerance = 1e-5)
})