if TRUE, returns the Kullback-Leibler divergence with the independent case, if FALSE, the L2 distance. There is no re-normalization, contrary to what happens for simany.
nbsafe
the ratio between the number of sub-samples and the cardinality of the discretized copula.
nthreads
number of number of threads, assumed to be strictly positive. For "full throttle" computations, consider using parallel::detectCores()