# Create a unit test
test_that("spatmcacvall_rcpp returns the correct result", {
  # Set up test data (replace this with your actual test data)
  sxr <- matrix(1:20, nrow = 10, ncol = 2)
  syr <- matrix(21:40, nrow = 10, ncol = 2)
  Xr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  Yr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  M <- 5
  K <- 2
  tau1ur <- c(0.1, 0.5, 1.0)
  tau2ur <- c(0.2, 0.8)
  tau1vr <- c(0.1, 0.5, 1.0)
  tau2vr <- c(0.2, 0.8)
  nkr <- c(5, 10, 15)
  maxit <- 100
  tol <- 1e-6
  l2ur <- c(0.1, 0.5, 1.0)
  l2vr <- c(0.2, 0.8)
  
  # Call the Rcpp function
  result <- spatmcacvall_rcpp(sxr, syr, Xr, Yr, M, K, tau1ur, tau2ur, tau1vr, tau2vr, nkr, maxit, tol, l2ur, l2vr)
  
  # Check if the result is a list with expected components
  expect_type(result, "list")
  expect_true("cvall" %in% names(result))
  expect_true("Uest" %in% names(result))
  expect_true("Vest" %in% names(result))
  expect_true("Dest" %in% names(result))
  expect_true("cvtau1u" %in% names(result))
  expect_true("cvtau2u" %in% names(result))
  expect_true("cvtau1v" %in% names(result))
  expect_true("cvtau2v" %in% names(result))
})

test_that("tpm2 function produces correct output", {
  # Set up test data (replace this with your actual test data)
  z <- matrix(rnorm(100), nrow = 10, ncol = 10)
  P <- matrix(rnorm(100), nrow = 10, ncol = 10)
  Phi <- matrix(rnorm(100), nrow = 10, ncol = 10)
  
  # Call the tpm2 function
  result <- tpm2(z, P, Phi)
  
  expect_equal(dim(result), c(nrow(z), ncol(Phi)))
})


# Create a unit test
test_that("spatmcacv_rcpp function produces correct output", {
  # Set up test data (replace this with your actual test data)
  sxr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  syr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  Xr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  Yr <- matrix(rnorm(100), nrow = 10, ncol = 10)
  M <- 5
  K <- 3
  tau1ur <- seq(0.1, 1, length.out = 5)
  tau2ur <- seq(0.1, 1, length.out = 5)
  tau1vr <- seq(0.1, 1, length.out = 5)
  tau2vr <- seq(0.1, 1, length.out = 5)
  nkr <- seq(1, 5)
  maxit <- 100
  tol <- 1e-5
  l2ur <- seq(0.1, 1, length.out = 5)
  l2vr <- seq(0.1, 1, length.out = 5)
  
  # Call the spatmcacv_rcpp function
  result <- spatmcacv_rcpp(sxr, syr, Xr, Yr, M, K, tau1ur, tau2ur, tau1vr, tau2vr, nkr, maxit, tol, l2ur, l2vr)
  
  
  # Check if the result is a list
  expect_type(result, "list")
  
  # Check if the list contains expected components
  expect_true("cv1" %in% names(result))
  expect_true("cv2" %in% names(result))
  expect_true("Uest" %in% names(result))
  expect_true("Vest" %in% names(result))
  expect_true("Dest" %in% names(result))
  expect_true("cvtau1u" %in% names(result))
  expect_true("cvtau2u" %in% names(result))
  expect_true("cvtau1v" %in% names(result))
  expect_true("cvtau2v" %in% names(result))
})