context("Comparison between the results of adjClust with sparse and dense matrices")
#Sys.setenv("OMP_THREAD_LIMIT" = 2)

mat <- matrix(c(1.0, 0.0, 0.0, 0.0, 0.0, 
                0.1, 1.0, 0.0, 0.0, 0.0, 
                0.5, 0.2, 1.0, 0.0, 0.0, 
                0.8, 0.6, 0.3, 1.0, 0.0, 
                0.0, 0.9, 0.7, 0.4, 1.0),
              nrow = 5)
mat <- mat + t(mat)

if (packageVersion("Matrix") < '1.5.0') {
  smat1 <- as(mat, "dgeMatrix")
  smat2 <- as(mat, "dgCMatrix")
  smat3 <- as(mat, "dgTMatrix")
  
  mat <- Matrix::forceSymmetric(mat)
  smat4 <- as(mat, "dsCMatrix")
  smat5 <- as(mat, "dsTMatrix")
  smat6 <- as(mat, "dsyMatrix")
} else {
  smat1 <- as(as(as(mat, "dMatrix"), "generalMatrix"), "unpackedMatrix")
  smat2 <- as(as(as(mat, "dMatrix"), "generalMatrix"), "CsparseMatrix")
  smat3 <- as(as(as(mat, "dMatrix"), "generalMatrix"), "TsparseMatrix")
  
  mat <- Matrix::forceSymmetric(mat)
  smat4 <- as(as(as(mat, "dMatrix"), "symmetricMatrix"), "CsparseMatrix")
  smat5 <- as(as(as(mat, "dMatrix"), "symmetricMatrix"), "TsparseMatrix")
  smat6 <- as(as(as(mat, "dMatrix"), "symmetricMatrix"), "unpackedMatrix")
}

mat <- as(mat, "matrix")
p <- nrow(mat)

test_that("test that adjClust gives identical results for sparse and dense matrices when h < p-1", {
  #Sys.setenv("OMP_THREAD_LIMIT" = 2)
  fit1 <- adjClust(mat, h = 2)
  fit2 <- adjClust(smat1, h = 2)
  fit3 <- adjClust(smat2, h = 2)
  fit4 <- adjClust(smat3, h = 2)
  fit5 <- adjClust(smat4, h = 2)
  fit6 <- adjClust(smat5, h = 2)
  fit7 <- adjClust(smat6, h = 2)
  
  expect_equal(fit1$merge, fit2$merge)
  expect_equal(fit1$height, fit2$height)
  
  expect_equal(fit1$merge, fit3$merge)
  expect_equal(fit1$height, fit3$height)
  
  expect_equal(fit1$merge, fit4$merge)
  expect_equal(fit1$height, fit4$height)
  
  expect_equal(fit1$merge, fit5$merge)
  expect_equal(fit1$height, fit5$height)
  
  expect_equal(fit1$merge, fit6$merge)
  expect_equal(fit1$height, fit6$height)
  
  expect_equal(fit1$merge, fit7$merge)
  expect_equal(fit1$height, fit7$height)
  
  # test that adjClust methods returns expected 'calls' for sparse matrices
  expect_identical(as.list(fit1$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit2$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit3$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit4$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit5$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit6$call)[[1]], as.symbol("adjClust"))
  expect_identical(as.list(fit7$call)[[1]], as.symbol("adjClust"))
  
})

test_that("test that adjClust gives identical results for sparse and dense matrices when h is p-1", {
  #Sys.setenv("OMP_THREAD_LIMIT" = 2)
  fit1 <- adjClust(mat)
  fit2 <- adjClust(smat1)
  fit3 <- adjClust(smat2)
  fit4 <- adjClust(smat3)
  fit5 <- adjClust(smat4)
  fit6 <- adjClust(smat5)
  fit7 <- adjClust(smat6)
  
  expect_equal(fit1$merge, fit2$merge)
  expect_equal(fit1$height, fit2$height)
  
  expect_equal(fit1$merge, fit3$merge)
  expect_equal(fit1$height, fit3$height)
  
  expect_equal(fit1$merge, fit4$merge)
  expect_equal(fit1$height, fit4$height)
  
  expect_equal(fit1$merge, fit5$merge)
  expect_equal(fit1$height, fit5$height)
  
  expect_equal(fit1$merge, fit6$merge)
  expect_equal(fit1$height, fit6$height)
  
  expect_equal(fit1$merge, fit7$merge)
  expect_equal(fit1$height, fit7$height)
})