`fit` not equal to `fit2`.
Component "copula": Component "pair_copulas": Component 1: Component 2: Component 3: Mean relative difference: 0.4962467
Component "copula": Component "pair_copulas": Component 1: Component 4: Component 3: Mean relative difference: 0.4962467
Component "copula": Component "pair_copulas": Component 1: Component 6: Component 3: Mean relative difference: 0.4962467
Component "copula": Component "pair_copulas": Component 2: Component 1: Component 3: Mean relative difference: 0.1354241
Component "copula": Component "pair_copulas": Component 2: Component 2: Component 3: Mean relative difference: 0.01119977
Component "copula": Component "pair_copulas": Component 2: Component 3: Component 3: Mean relative difference: 0.1354241
Component "copula": Component "pair_copulas": Component 2: Component 4: Component 3: Mean relative difference: 0.01119977
Component "copula": Component "pair_copulas": Component 2: Component 5: Component 3: Mean relative difference: 0.1354241
Component "copula": Component "pair_copulas": Component 2: Component 6: Component 3: Mean relative difference: 0.01119977
...

test_that("creating custom models (1d)", {
lapply(1:3, function(j) replicate(4 - j, bc, simplify = FALSE)),
test_that("creating custom models (multivariate)", {
lapply(1:11, function(j) replicate(12 - j, bc, simplify = FALSE)),
expect_equal(unname(dim(svine_sim(10, 1, fit))), c(10, 1))
expect_equal(unname(dim(svine_sim(10, 1, fit, x))), c(10, 1))
expect_equal(unname(dim(svine_sim(10, 3, fit, x))), c(10, 1, 3))
fit <- svine(x, margin_families = c("exp", "norm"), p = 0)
expect_equal(unname(dim(svine_sim(10, 1, fit))), c(10, 2))
expect_equal(unname(dim(svine_sim(10, 1, fit, x))), c(10, 2))
expect_equal(unname(dim(svine_sim(10, 3, fit, x))), c(10, 2, 3))
fit <- svine(x, margin_families = c("exp", "norm"), p = 3)
expect_gt(min(svine_sim(10, 1, fit, x)), 0)
expect_gt(min(svine_sim(10, 5, fit, x)), 0)
expect_equal(fit, fit2, tolerance = 0.01)
expect_silent(svine(x, p = 1, type = "D", margin_families = c("exp", "norm")))
expect_silent(svine(x, p = 2, type = "M", margin_families = c("exp", "norm")))
expect_error(svine(x, p = 1, type = "R", margin_families = c("exp", "norm")))
fit <- svine(x, margin_families = c("std", "norm", "exp"), p = 0)
fit <- svine(x, margin_families = "empirical", p = 0)