powered by
find.smooth.tau(dl, psi = exp(dl$hyper$mu_psi), omega = exp(dl$hyper$mu_omega), num.cores = default.num.cores(), num.tau.to.try = num.cores, num.tau.to.keep = num.cores, method = "metropolis", ...)