tests/testthat/test-cost_function.R

In PoolDilutionR: Calculate Gross Biogeochemical Flux Rates from Isotope Pool Dilution Data

test_that("pdr_cost works", {
  times <- 0:5
  m <- c(10, 8, 6, 5, 4, 3)
  n <- c(1, 0.7, 0.6, 0.4, 0.3, 0.2)

  x <- pdr_cost(
    params = list(P = 0.5, k = 0.3),
    time = times,
    m = m, n = n, m_prec = 0.001, ap_prec = 0.01
  )
  expect_type(x, "double")
  expect_identical(length(x), 1L)

  # TODO: not sure what else to test

  # Log function called if supplied
  logfn_called <- FALSE
  logfn <- function(df) {
    # Cost functions should pass a data frame to the log function
    expect_s3_class(df, "data.frame")
    logfn_called <<- TRUE
  }
  pdr_cost(
    params = list(P = 0.5, k = 0.3),
    time = times,
    m = m, n = n, m_prec = 0.001, ap_prec = 0.01,
    log_progress = logfn
  )
  expect_true(logfn_called)
})