tests/testthat/test_ops_addition.R

In polyMatrix: Infrastructure for Manipulation Polynomial Matrices

context("ops addition")

test_that("with number", {
  pm <- parse.polyMatrix("  0,       1, 1 + x, x",
                         "x^2, 1 + x^2,     2, 0")
  expect_equal(pm + 11, parse.polyMatrix("      11,       12, 12 + x, 11 + x",
                                         "11 + x^2, 12 + x^2,     13, 11"))
  expect_equal(-2 + pm, parse.polyMatrix("      -2,       -1, -1 + x, -2 + x",
                                         "-2 + x^2, -1 + x^2,      0, -2"))
})

test_that("with polynomial", {
  pm <- parse.polyMatrix("  0, x",
                         "x^2, 3",
                         "1 + x, 1 + x^2")
  expect_equal(pm + parse.polynomial("1 + x"),
               parse.polyMatrix("      1 + x, 1 + 2x",
                                "1 + x + x^2, 4 + x",
                                "     2 + 2x, 2 + x+ x^2"))
  expect_equal(parse.polynomial("1 - x^2") + pm,
               parse.polyMatrix("    1 - x^2, 1 + x - x^2",
                                "          1, 4 - x^2",
                                "2 + x - x^2, 2"))
})

test_that("with matrix", {
  pm <- parse.polyMatrix("  0, x",
                         "x^2, 3",
                         "1 + x, 1 + x^2")

  expect_equal(pm + matrix(1:6, 3, 2, byrow = TRUE),
               parse.polyMatrix("      1, 2 + x",
                                "3 + x^2, 7",
                                "  6 + x, 7 + x^2"))
  expect_equal(matrix(1:6, 3, 2) + pm,
               parse.polyMatrix("1, 4 + x",
                                "2 + x^2, 8",
                                "4 + x, 7+x^2"))
  expect_error(pm + matrix(0, 1, 2), "non-conformable")
})

test_that("with polyMatrix", {
  first <- parse.polyMatrix("1 + x, 1 + x^2",
                            "    2, 2 + 3x")
  second <- parse.polyMatrix("1, 2 + x",
                             "-3, -4 - 3x")
  expect_equal(first + second, parse.polyMatrix("2 + x,   3 + x + x^2",
                                                "    -1, -2"))
  expect_equal(second + first, parse.polyMatrix("2 + x,   3 + x + x^2",
                                                "    -1, -2"))

  expect_error(first + polyMatrix(0, 1, 2), "non-conformable")
})