testthat::skip_on_cran()

tolerance <- 1e-6

emb <- c(3,4)
tau <- c(6,10)
threshold <- c(1.2, 1.8)

test_that("AccRQA results match CRQA -- RR, euclidean",{
  input_data <- readRDS("testdata/input.rds")
  expected_RR <- readRDS("testdata/expected_CRQA_RR_euclidean.rds")
  
  results_RR <- accrqa_RR(input_data, 
                       tau_values = tau,
                       emb_values = emb,
                       threshold_values = threshold,
                       distance_type = "euclidean",
                       comp_platform = "cpu")
  

  expect_equal(results_RR, expected_RR, tolerance = tolerance)
})

test_that("AccRQA results match CRQA -- RR, maximal",{
  input_data <- readRDS("testdata/input.rds")
  expected_RR <- readRDS("testdata/expected_CRQA_RR_maximal.rds")

    results_RR <- accrqa_RR(input_data, 
                          tau_values = tau,
                          emb_values = emb,
                          threshold_values = threshold,
                          distance_type = "maximal",
                          comp_platform = "cpu")
  
  
  expect_equal(results_RR, expected_RR, tolerance = tolerance)
})

test_that("AccRQA results match CRQA -- DET, euclidean",{
  input_data <- readRDS("testdata/input.rds")
  expected_DET <- readRDS("testdata/expected_CRQA_DET_euclidean.rds")
  
  results_DET <- accrqa_DET(input_data, 
                            tau_values = tau,
                            emb_values = emb,
                            lmin_values = c(2),
                            threshold_values = threshold,
                            distance_type = "euclidean",
                            calculate_ENTR = TRUE,
                            comp_platform = "cpu")
  
  results_DET <- subset(results_DET, select = -c(Lmax))
  expect_equal(results_DET, expected_DET, tolerance = tolerance)
  
})

test_that("AccRQA results match CRQA -- DET, maximal",{
  input_data <- readRDS("testdata/input.rds")
  expected_DET <- readRDS("testdata/expected_CRQA_DET_maximal.rds")
  
  results_DET <- accrqa_DET(input_data, 
                            tau_values = tau,
                            emb_values = emb,
                            lmin_values = c(2),
                            threshold_values = threshold,
                            distance_type = "maximal",
                            calculate_ENTR = TRUE,
                            comp_platform = "cpu")
  
  results_DET <- subset(results_DET, select = -c(Lmax))
  expect_equal(results_DET, expected_DET, tolerance = tolerance)
  
})

test_that("AccRQA results match CRQA -- LAM, euclidean",{
  input_data <- readRDS("testdata/input.rds")
  expected_LAM <- readRDS("testdata/expected_CRQA_LAM_euclidean.rds")


  results_LAM <- accrqa_LAM(input_data, 
                            tau_values = tau,
                            emb_values = emb,
                            vmin_values = 2,
                            threshold_values = threshold,
                            distance_type = "euclidean",
                            calculate_ENTR = TRUE,
                            comp_platform = "cpu")
  
  results_LAM <- subset(results_LAM, select = -c(ENTR, TTmax))
  expect_equal(results_LAM, expected_LAM, tolerance = tolerance)
})


test_that("AccRQA results match CRQA -- LAM, maximal",{
  input_data <- readRDS("testdata/input.rds")
  expected_LAM <- readRDS("testdata/expected_CRQA_LAM_maximal.rds")
  
  results_LAM <- accrqa_LAM(input_data, 
                            tau_values = tau,
                            emb_values = emb,
                            vmin_values = 2,
                            threshold_values = threshold,
                            distance_type = "maximal",
                            calculate_ENTR = TRUE,
                            comp_platform = "cpu")
  
  results_LAM <- subset(results_LAM, select = -c(ENTR, TTmax))
  expect_equal(results_LAM, expected_LAM, tolerance = tolerance)
})