# Default-output baseline (frozen oracle, cairn/ORACLES.md).
#
# The fixture holds the loadings, weights, variance, factor_cor, and every
# edge matrix of five default fits, generated on master before the
# factor-correlation carry and the partialled edge columns were added
# (data-raw/baseline-m89.R names the commit). Each fit is refit here fresh
# (never through cached(): a memoized second call asserts nothing) and every
# stored value must reproduce within a platform tolerance. The fixture was
# generated on macOS (Accelerate BLAS). Linux and Windows link a different
# BLAS, which moves the pca/efa values by about 1e-12 (the polychoric fits
# most) and the seeded ESEM values by about 1e-6 at one level, because
# lavaan's optimizer takes a slightly different path from the same random
# starts (CI runs of 2026-09-16 and 2026-09-23). The tolerances are 1e-8 for
# pca/efa and 1e-5 for esem: loose enough to absorb that drift, tight enough
# that any method change (which moves values at 1e-3 or more) still fails.
# The ESEM fit is seeded, because lavaan's rotation draws random starts (see
# the generator's header).

baseline <- readRDS(test_path("fixtures", "baseline-m89.rds"))

.refit_baseline <- function(spec) {
  dat <- get(spec$data, envir = asNamespace("ackwards"))
  suppressWarnings(suppressMessages(
    ackwards(dat, k_max = spec$k_max, engine = spec$engine, cor = spec$cor, seed = spec$seed)
  ))
}

for (name in names(baseline)) {
  test_that(paste0("default output is unchanged against the baseline: ", name), {
    entry <- baseline[[name]]
    if (entry$spec$engine == "esem") skip_if_not_installed("lavaan")
    x <- .refit_baseline(entry$spec)
    tol <- if (entry$spec$engine == "esem") 1e-5 else 1e-8

    expect_identical(names(x$levels), names(entry$levels))
    for (ki in names(entry$levels)) {
      lev <- x$levels[[ki]]
      ref <- entry$levels[[ki]]
      expect_equal(lev$loadings, ref$loadings, tolerance = tol, label = paste(name, ki, "loadings"))
      expect_equal(lev$scoring$weights, ref$weights, tolerance = tol, label = paste(name, ki, "weights"))
      expect_equal(lev$variance, ref$variance, tolerance = tol, label = paste(name, ki, "variance"))
      # Values only: PCA and EFA factor_cor gained the level's labels after the
      # fixture was generated (ESEM's already carried them).
      expect_equal(unname(lev$factor_cor), unname(ref$factor_cor),
        tolerance = tol, label = paste(name, ki, "factor_cor")
      )
    }

    expect_identical(names(x$edges$matrices), names(entry$edges))
    for (key in names(entry$edges)) {
      expect_equal(x$edges$matrices[[key]], entry$edges[[key]], tolerance = tol, label = paste(name, "edge", key))
    }
  })
}