tests/testthat/test-constrLaplacianRank.R

In spectralGraphTopology: Learning Graphs from Data via Spectral Constraints

context("Jordan's constrained Laplacian rank")
library(testthat)
library(spectralGraphTopology)

set.seed(42)

test_that("pairwise row norm", {
  n <- sample(c(3:10), size = 1)
  m <- sample(c(3:10), size = 1)
  M <- matrix(runif(n * m), n, m)
  V <- pairwise_matrix_rownorm2(M)
  expect_equal(c(n, n), dim(V))
  expect_true(isSymmetric(V))
  expect_equal(diag(V), rep(0, n))

  M <- rbind(c(1, 2, 3),
             c(3, 2, 0),
             c(1, 1, 1))
  V <- pairwise_matrix_rownorm2(M)
  ANS <- rbind(c(0, 13, 5),
              c(13,  0, 6),
               c(5,  6, 0))
  expect_equal(V, ANS)
})

test_that("Jordan's initial smooth regular graph", {
  m <- sample(c(10:20), size = 1)
  p <- 5
  Y <- matrix(runif(p * m), p, m)
  LA <- learn_smooth_approx_graph(Y, m = 2)
  expect_equal(c(p, p), dim(LA))
  expect_equal(rowSums(LA), rep(0, p))
})