tests/testthat/test.fem.R

In rSPDE: Rational Approximations of Fractional Stochastic Partial Differential Equations

context("FEM calculations")

test_that("FEM matrix construction", {
  testthat::skip_on_cran()
  local_rspde_safe_inla()

  x <- c(0,1,1.75,2)
  fem <- rSPDE.fem1d(x)
  mesh <- fmesher::fm_mesh_1d(x)
  fem_inla <- fmesher::fm_fem(mesh)

  # Pr multiplication
  expect_equal(fem$C, fem_inla$c1, tolerance = 1e-10)
  expect_equal(fem$G, fem_inla$g1, tolerance = 1e-10)

})

test_that("A matrix construction", {
  testthat::skip_on_cran()
  local_rspde_safe_inla()

  x <- c(0,1,1.75,2)
  fem <- rSPDE.fem1d(x)
  loc <- c(0.6,1.1,2)
  A <- rSPDE.A1d(x,loc)
  mesh <- fmesher::fm_mesh_1d(x)
  A_inla <- fmesher::fm_basis(mesh,loc)

  # Pr multiplication
  # From fmesher 0.7.0, fm_basis removes redundant explicit zeros, so
  # cannot use expect_equal(A, B), as that would trigger an error due to
  # mismatching internal matrix storage.
  expect_equal(sum(abs(A - A_inla)), 0, tolerance = 1e-10)

})