test_that("mcca_fit returns correct structure", {

  set.seed(1)

  n <- 10
  T_len <- 6

  X <- list(
    lapply(seq_len(n), function(i) matrix(rnorm(T_len * 3), T_len, 3)),
    lapply(seq_len(n), function(i) matrix(rnorm(T_len * 2), T_len, 2))
  )

  fit <- mcca_fit(method = "kernel", X = X, ncomp = 2)

  expect_s3_class(fit, "mcca_fit")

  expect_true(is.list(fit$scores))
  expect_equal(length(fit$scores), 2)

  expect_equal(ncol(fit$scores[[1]]), 2)
  expect_equal(ncol(fit$scores[[2]]), 2)

})