tests/testthat/test-gsp_util.R

In gmonette/spida2: Collection of tools developed for the Summer Programme in Data Analysis 2000-2012

test_that("Polyshift works", {
  expect_equal(
    {
      coef <- 1:4
      x <- 4
      coef*(x^(0:3))
      sum((PolyShift(3,4) %*% coef)*(x-3)^(0:3) ) -
      sum(coef*(x^(0:3)))
    }, 0)
})