set.seed(20250522) fcl <- function(need, covlist) distfreereg:::fill_covariance_list(need = need, covariance_list = covlist, matsqrt_tol = distfreereg:::default_distfreereg_tol()[["matsqrt_tol"]], solve_tol = distfreereg:::default_distfreereg_tol()[["solve_tol"]]) n <- 1e2 Sig <- rWishart(1, df = n, Sigma = diag(n))[,,1] P <- fcl(need = "P", list(Sigma = Sig))[["P"]] SqrtSigma <- fcl(need = "SqrtSigma", list(P = P))[["SqrtSigma"]] Q <- fcl(need = "Q", list(SqrtSigma = SqrtSigma))[["Q"]] Sig_2 <- fcl(need = "Sigma", list(SqrtSigma = SqrtSigma))[["Sigma"]] Sig_3 <- fcl(need = "Sigma", list(Q = Q))[["Sigma"]] P_2 <- fcl(need = "P", list(Q = Q))[["P"]] P_3 <- fcl(need = "P", list(SqrtSigma = SqrtSigma))[["P"]] tryCatch(distfreereg:::model2method("hello"), error = function(e) warning(e)) tryCatch(distfreereg::rejection("hello"), error = function(e) warning(e)) mats <- rWishart(n = 2, df = 5, Sigma = diag(5)) mat_list <- list(mats[,,1], mats[,,2]) mat_list_sq <- distfreereg:::covsquare(mat_list) stopifnot(isTRUE(all.equal(mat_list_sq, lapply(mat_list, crossprod)))) vec <- rexp(5) vec_sq <- distfreereg:::matsquare(vec) stopifnot(isTRUE(all.equal(vec_sq, vec * vec))) err_dist_args <- distfreereg:::create_err_dist_args_default( err_dist_fun_formals = c("Sigma", "SqrtSigma", "P", "Q"), cov_args = c("Sigma", "SqrtSigma", "P", "Q"), reps = 10, n = 1e2, true_covariance = list(Sigma = Sig, P = P, Q = Q, SqrtSigma = SqrtSigma)) names(err_dist_args) distfreereg:::matsqrt(matrix(4)) tryCatch(distfreereg:::covprod(cov = array(4), mat = matrix(1)), error = function(e) warning(e)) tryCatch(distfreereg:::scale_to_01(matrix(rep(1, 5)), one_value = "fail"), error = function(e) warning(e))