.devel/testthat/test-rel-reflexivity.R

require('testthat')
require('Matrix')


test_that("rel_is_reflexive", {

   expect_true(rel_is_reflexive(matrix(c(TRUE), nrow=1)))
   expect_true(rel_is_reflexive(matrix(c(TRUE, FALSE, FALSE, TRUE), nrow=2)))
   expect_false(rel_is_reflexive(matrix(c(FALSE), nrow=1)))
   expect_true(rel_is_reflexive(matrix(c(1, 0, 0, 1), nrow=2)))
   expect_true(is.na(rel_is_reflexive(matrix(c(1, 0, 0, NA), nrow=2))))
   expect_false(rel_is_reflexive(matrix(c(1, 0, 0, 0), nrow=2)))
   expect_false(rel_is_reflexive(Matrix(c(1, 0, 0, 0), nrow=2)))
   expect_true(rel_is_reflexive(Matrix(c(999, 0, 0, 100), nrow=2)))

})



test_that("rel_closure_reflexive", {

   expect_equivalent(rel_closure_reflexive(
         matrix(0,ncol=4, nrow=4, byrow=TRUE)
      ),
      matrix(diag(4),ncol=4, byrow=TRUE))


   A <- structure(as.logical(diag(4)), dimnames=list(LETTERS[1:4], LETTERS[1:4]), dim=c(4,4))
   expect_identical(rel_closure_reflexive(A), A)

   # random graphs:
   expect_true(rel_is_reflexive(rel_closure_reflexive(
     matrix(runif(625)>0.1, ncol=25, byrow=TRUE))))

   expect_true(rel_is_reflexive(rel_closure_reflexive(
     matrix(runif(625)>0.3, ncol=25, byrow=TRUE))))

   expect_true(rel_is_reflexive(rel_closure_reflexive(
     matrix(runif(625)>0.5, ncol=25, byrow=TRUE))))

   expect_true(rel_is_reflexive(rel_closure_reflexive(
     matrix(runif(625)>0.7, ncol=25, byrow=TRUE))))

   expect_true(rel_is_reflexive(rel_closure_reflexive(
     matrix(runif(625)>0.9, ncol=25, byrow=TRUE))))
})


test_that("rel_reduction_reflexive", {

   expect_equivalent(rel_reduction_reflexive(
         matrix(diag(4),ncol=4, nrow=4, byrow=TRUE)
      ),
      matrix(FALSE, ncol=4, nrow=4, byrow=TRUE))

   A <- structure(as.logical(diag(4)), dimnames=list(LETTERS[1:4], LETTERS[1:4]), dim=c(4,4))
   expect_identical(rel_closure_reflexive(rel_reduction_reflexive(A)), A)

})
Rexamine/agop documentation built on Dec. 11, 2023, 10:02 p.m.