tests/testthat/test_oper_cps.R

In aphalo/photobiology: Photobiological Calculations

library("photobiology")
context("cps_spct")

test_that("constructor", {

  empty.spct <- cps_spct()
  expect_true(is.cps_spct(empty.spct))
  expect_true(is.any_spct(empty.spct))
  expect_named(empty.spct, c("w.length", "cps"))
  expect_equal(nrow(empty.spct), 0L)

  my.spct <- cps_spct(w.length = 400:409, cps = 1)

  expect_equal(my.spct[["cps"]], rep(1, length.out = 10))
  expect_equal(my.spct[["w.length"]], 400:409)
  expect_named(my.spct, c("w.length", "cps"))
  expect_true(is.cps_spct(my.spct))
  expect_true(is.any_spct(my.spct))
  expect_false(is.raw_spct(my.spct))
  expect_false(is.source_spct(my.spct))
  expect_false(is.filter_spct(my.spct))
  expect_false(is.reflector_spct(my.spct))
  expect_false(is.object_spct(my.spct))
  expect_false(is.response_spct(my.spct))
  expect_false(is.chroma_spct(my.spct))

  my.df <- data.frame(w.length = 400:409, cps = 1)
  my.spct <- as.cps_spct(my.df)

  expect_equal(my.spct[["cps"]], rep(1, length.out = 10))
  expect_equal(my.spct[["w.length"]], 400:409)
  expect_named(my.spct, c("w.length", "cps"))
  expect_true(is.cps_spct(my.spct))
  expect_true(is.any_spct(my.spct))
})


test_that("oper", {

  my.spct <- cps_spct(w.length = 400:409, cps = 1)
  my.2.spct <- cps_spct(w.length = 400:409, cps= 2)

  expect_equal(class(my.spct)[1:2], c("cps_spct", "generic_spct") )
  expect_equal(attr(my.spct, "spct.version", exact = TRUE), 2)

#  order of attribbutes is changed
#  expect_equal(my.spct + my.spct,  my.2.spct)
#  expect_equal(my.spct + my.spct,  my.spct * 2)
  expect_equal(my.spct * 2, my.2.spct)
  expect_equal(my.spct * 2L, my.2.spct)
  expect_equal(my.2.spct / 2, my.spct)
  expect_equal(my.2.spct / 2L, my.spct)
  expect_equal(-my.spct * -2, my.2.spct)
  expect_equal(-my.spct * -2L, my.2.spct)
  expect_equal(-my.2.spct / -2, my.spct)
  expect_equal(-my.2.spct / -2L, my.spct)
  expect_equal( 2 * my.spct, my.2.spct)
  expect_equal( 1 / (2 / my.2.spct), my.spct)
  expect_equal( 1 / my.spct, my.spct^-1)
  expect_equal(my.2.spct %/% 2L, my.spct)
  expect_equal(my.2.spct %% 2L, my.spct %% 1L)
})


test_that("math", {

  my.spct <- cps_spct(w.length = 400:409, cps = 1)

  expect_equal(log10(my.spct)[["cps"]],  rep(log10(1), length.out = 10))
  expect_equal(log(my.spct)[["cps"]],  rep(log(1), length.out = 10))
  expect_equal(log(my.spct, 2)[["cps"]],  rep(log(1, 2), length.out = 10))
  expect_equal(exp(my.spct)[["cps"]],  rep(exp(1), length.out = 10))
  expect_equal(sqrt(my.spct)[["cps"]],  rep(sqrt(1), length.out = 10))

})