library(BCodifSIS)

n <- 30
p <- 6
grid <- seq_len(n)
X <- vapply(seq_len(p), function(j) {
  sin(grid / (j + 1)) + cos(grid * j / 9)
}, numeric(n))
y <- X[, 1] + 0.25 * sin(grid / 4)

fast <- bcodif_score(X[, 1], y)
slow <- bcodif_score_slow(X[, 1], y)
stopifnot(isTRUE(all.equal(fast, slow, tolerance = 1e-10)))

fit <- bcodif_sis(X, y, d = 3)
stopifnot(inherits(fit, "bcodif_sis"))
stopifnot(length(fit$ix) == 3)
stopifnot(all(is.finite(fit$scores)))

qc <- quick_check()
stopifnot(is.list(qc))
stopifnot(qc$score_difference < 1e-10)