R Under development (unstable) (2025-06-30 r88369 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> 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))
Warning message:
In distfreereg:::model2method("hello") :
  invalid model type found: character
> 
> tryCatch(distfreereg::rejection("hello"), error = function(e) warning(e))
Warning message:
In validate_args_rejection(object = object, alpha = alpha, stat = stat) :
  'object' must have class 'compare'
> 
> 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)
[1] "Sigma"     "SqrtSigma" "P"         "Q"        
> 
> 
> distfreereg:::matsqrt(matrix(4))
     [,1]
[1,]    2
> 
> tryCatch(distfreereg:::covprod(cov = array(4), mat = matrix(1)),
+          error = function(e) warning(e))
Warning message:
In distfreereg:::covprod(cov = array(4), mat = matrix(1)) :
  invalid covariance type
> 
> tryCatch(distfreereg:::scale_to_01(matrix(rep(1, 5)), one_value = "fail"),
+          error = function(e) warning(e))
Warning message:
In distfreereg:::scale_to_01(matrix(rep(1, 5)), one_value = "fail") :
  At least one column of sorting matrix has only one value
> 
> proc.time()
   user  system elapsed 
   1.73    0.18    1.90