tests/testthat/test-02-conversions.R

In KopfLab/microbialkitchen: Culture Media Tools

context("Conversions")

test_that("Testing that standard vector operations are working", {
  
  # subsetting
  expect_equal(qty(1:5, "g")[2:3], qty(2:3, "g"))
  expect_equal(qty(1:5, "g")[-5], qty(1:4, "g"))
  expect_equal(qty(1:5, "g")[-c(1:5)], qty(numeric(0), "g"))
  
  # combining
  expect_equal(c_qty(qty(5, "g"), qty(c(10, 20), "mg")), qty(c(5000, 10, 20), "mg"))
  expect_equal(c(qty(5, "g"), qty(c(10, 20), "mg")), qty(c(5000, 10, 20), "mg"))
  expect_equal(c(qty(NA, "g"), qty(c(10, 20), "mg")), qty(c(NA, 10, 20), "mg"))
  expect_equal(c(qty(numeric(0), "g"), qty(c(10, 20), "mg")), qty(c(10, 20), "mg"))
  expect_error(c(qty(5, "g"), 5), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(c(qty(5, "g"), 5L), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(c(qty(5, "g"), qty(5, "L")), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_equal(c(5, qty(5, "L")), c(5, 5)) # problem with the vctrs dispatch
})

test_that("Testing that arithemtic calculations of units are working", {
  
  # comparisons
  expect_error(qty(2, "C") > qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(qty(2, "C") >= qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(qty(2, "C") < qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(qty(2, "C") <= qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(qty(2, "C") == qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_error(qty(2, "C") != qty(1, "g"), "Can't combine", class = "vctrs_error_incompatible_type")
  expect_true(qty(100, "mg") < qty(1, "g"))
  expect_true(qty(100, "mg") <= qty(0.1, "g"))
  expect_true(qty(100, "ng") > qty(10000, "pg"))
  expect_true(qty(10, "ng") >= qty(10000, "pg"))
  expect_true(qty(10, "ng") == qty(10000, "pg"))
  expect_false(qty(10, "ng") != qty(10000, "pg"))
  
  # addition / subtraction
  expect_error(qty(2, "C") + 1, "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(1 + qty(2, "C"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(qty(2, "C") + qty(1, "g"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(qty(2, "C") - 1, "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(1 - qty(2, "C"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(qty(2, "C") - qty(1, "g"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_equal({q <- qty(1, "kg") + qty(5, "g") - qty(0.0003, "Mg"); get_qty_units(q)}, "g")
  expect_equal(get_qty_value(q), 705) 
  expect_equal(qty(0.3, "mg") - qty(200, "µg") - qty(100, "µg"), qty(0, "g"))
  
  # division
  expect_error(qty(2, "C") / qty(1, "g"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_error(1000 / qty(1, "g"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_equal(qty(1, "g") / 1000, qty(1, "mg")) # auto-converts units
  expect_equal(qty(1, "g") / qty(1, "mg"), 1000) # should give plain number
  expect_equal(qty(5, "mM") / qty(1, "M"), 0.005) # should give plain number
  
  # multiplication
  expect_error(qty(1, "g") * qty(1, "g"), "not permitted", class = "vctrs_error_incompatible_op")
  expect_equal(qty(5, "g") * 100, qty(500, "g"))
  expect_equal(100 * qty(5, "g"), qty(500, "g"))
  
  # amount / volume = concentration
  expect_equal(qty(5, "pmol") / qty(1, "ml"), qty(5, "nM"))
  expect_equal(qty(5, "pmol") / qty(10, "mL"), qty(500, "pM"))
  # amount / molarity = volume
  expect_equal(qty(5, "nmol") / qty(10, "mM"), qty(500, "nL"))
  # volume * molarity = amount
  expect_equal(qty(1, "nl") * qty(10, "mM"), qty(10, "pmol"))
  expect_equal(qty(10, "mM") * qty(1, "nl"), qty(10, "pmol"))
  # chain: volume * molarity / volume = amount
  expect_equal(qty(1, "ml") * qty(10, "mM") / qty(1, "kL"), qty(10, "nM"))
  expect_equal(qty(1, "ml") * qty(0.5, "mM") / qty(1000, "mM"), qty(500, "nL"))
  
  # mass / MW = amount
  expect_equal(qty(5, "ng") / qty(100, "g/mol"), qty(50, "pmol"))
  # mass / amount = MW
  expect_equal(qty(2, "mg") / qty(1, "mmol"), qty(2, "g/mol"))
  # amount * MW = mass
  expect_equal(qty(10, "nmol") * qty(50, "g/mol"), qty(500, "ng"))
  expect_equal(qty(50, "g/mol") * qty(10, "nmol"), qty(500, "ng"))
  
  # molarity / pressure = solubility
  expect_equal(qty(10, "mM") / qty(200, "mbar"), qty(50, "mM/bar"))
  # molarity / solubility = pressure
  expect_equal(qty(10, "mM") / qty(50, "mM/bar"), qty(200, "mbar"))
  # solubility * pressure = molarity
  expect_equal(qty(50, "mM/bar") * qty (200, "mbar"), qty(10, "mM"))
  
})