diff --git a/code/analysis/investigating/invest_clust_perf.R b/code/analysis/investigating/invest_clust_perf.R new file mode 100644 index 0000000..649291c --- /dev/null +++ b/code/analysis/investigating/invest_clust_perf.R @@ -0,0 +1,47 @@ +library(colSBM) +library(microbenchmark) +library(here) + +set.seed(1234, "L'Ecuyer-CMRG") +data(dorebipartite) +str(dorebipartite) + +nb_cores = parallelly::availableCores() +concurrent_models <- min(length(seq(1,4)), nb_cores%/%2) +per_model <- (nb_cores - concurrent_models) %/% length(seq(1,4)) + +mb1 <- microbenchmark( + "fancy computation" = { + parallel::mclapply(seq(1, 4), function(idx) { + message("Start ", idx) + estimate_colBiSBM( + netlist = dorebipartite[1:4], + colsbm_model = "iid", + global_opts = list(nb_cores = per_model, verbosity = 0L) + ) + message("End ", idx) + }, + mc.cores = concurrent_models + ) + }, + times = 5L +) + +mb2 <- microbenchmark( + "Attributing all" = { + parallel::mclapply(seq(1, 4), function(idx) { + message("Start ", idx) + estimate_colBiSBM( + netlist = dorebipartite[1:4], + colsbm_model = "iid", + global_opts = list(nb_cores = parallelly::availableCores(omit = 1), verbosity = 1L) + ) + message("End ", idx) + }, + mc.cores = parallelly::availableCores(omit = 1) + ) + }, + times = 5L +) + +# Bilan : plus intéressant de bourriner \ No newline at end of file diff --git a/code/analysis/investigating/invest_inference_bernoulli_double_ARI.Rmd b/code/analysis/investigating/invest_inference_bernoulli_double_ARI.Rmd index 0a7e782..f8ed5a1 100644 --- a/code/analysis/investigating/invest_inference_bernoulli_double_ARI.Rmd +++ b/code/analysis/investigating/invest_inference_bernoulli_double_ARI.Rmd @@ -1,9 +1,179 @@ -```{r, echo = FALSE} -data <- readRDS("code/results//simulations/inference//bernoulli/bernoulli_inference_18-04-2024_09-41-45_1-972.Rds") - -prob_data <- data[which((data$pirho_double_row_ARI < 1 | data$pirho_double_col_ARI < 1)&(data$pirho_mean_row_ARI == 1 & data$pirho_mean_col_ARI == 1)),] +```{r libraries, echo = FALSE} +library(here) +library(dplyr) +library(colSBM) +library(aricode) +library(patchwork) +library(ggplot2) ``` -```{r , echo = FALSE} -knitr::kable(prob_data) +```{r plot-alpha, echo = FALSE} +plot_alpha <- function(alpha) { + p_alpha <- alpha[drop = FALSE] %>% + t() %>% + reshape2::melt() %>% + ggplot2::ggplot(ggplot2::aes(x = Var1, y = Var2, fill = value)) + + ggplot2::geom_tile() + + ggplot2::geom_text(ggplot2::aes(label = round(value, 2)), color = "black") + + ggplot2::scale_fill_gradient2("alpha", + low = "white", + high = "red", + limits = c( + 0, + 1 + ) + ) + + ggplot2::geom_hline(yintercept = seq(nrow(alpha)) + .5) + + ggplot2::geom_vline(xintercept = seq(ncol(alpha)) + .5) + + ggplot2::scale_y_reverse() + + ggplot2::theme_bw(base_size = 15, base_rect_size = 1, base_line_size = 1) + + ggplot2::xlab("") + + ggplot2::ylab("") + + ggplot2::coord_fixed(expand = FALSE) + return(p_alpha) +} +``` + +```{r, echo = FALSE} +data_folder <- file.path( + here(), "code", "results", + "investigating", "selection_error", "incorrect24-04-2024_17-18-42" +) + +file_list <- list.files(data_folder) +file <- file_list[[1]] +data <- readRDS(file.path(data_folder, file)) + +#  Extracting data to netlist and Z +netlist_memb <- data[["netlist"]] +netlist <- lapply(netlist_memb, function(full) full[["incidence_matrix"]]) +nr <- nrow(netlist[[1]]) +nc <- ncol(netlist[[1]]) + +row_blockmemberships <- lapply(netlist_memb, function(full) full[["row_blockmemberships"]]) +col_blockmemberships <- lapply(netlist_memb, function(full) full[["col_blockmemberships"]]) + +joined_row_memberships <- unlist(row_blockmemberships) +joined_col_memberships <- unlist(col_blockmemberships) + +row_taus <- lapply(row_blockmemberships, function(Z1_m) colSBM:::.one_hot(Z1_m, 4)) +col_taus <- lapply(col_blockmemberships, function(Z2_m) colSBM:::.one_hot(Z2_m, 4)) + + +``` + +```{r parameters, echo = FALSE} +set.seed(1234, "L'Ecuyer-CMRG") +fitted_bisbmpop_pirho <- estimate_colBiSBM( + netlist = netlist, + colsbm_model = "pirho", + nb_run = 3L, + distribution = "bernoulli", + global_opts = list( + verbosity = 3L, + plot_details = 0, + nb_cores = parallelly::availableCores(omit = 1) + ) +) +``` + +```{r, echo = FALSE} +bad_model <- fitted_bisbmpop_pirho[["best_fit"]]$clone() +good_model <- fitted_bisbmpop_pirho[["model_list"]][[4, 4]]$clone() +``` + +```{r ARI, echo = FALSE} +good_Z <- good_model[["Z"]] +good_row <- lapply(good_Z, function(model) model[["row"]]) +good_col <- lapply(good_Z, function(model) model[["col"]]) + +good_joined_row <- unlist(good_row) +good_joined_col <- unlist(good_col) + +sapply(seq_len(length(good_row)), function(m) { + ARI( + row_blockmemberships[[m]], + good_row[[m]] + ) +}) + +sapply(seq_len(length(good_row)), function(m) { + ARI( + col_blockmemberships[[m]], + good_col[[m]] + ) +}) + +ARI(joined_row_memberships, good_joined_row) +ARI(joined_col_memberships, good_joined_col) +``` + +```{r, echo = FALSE} +wrap_plots( + plot_alpha(data[["alpha"]])+ labs(caption = "Vrai alpha"), + plot(good_model, type = "meso", mixture = TRUE, values = TRUE) + labs(caption = "Q = (4,4)"), + plot(bad_model, type = "meso", mixture = TRUE, values = TRUE) + labs(caption = "Q = (4,5)") +) +``` + +```{r poisson?, echo = FALSE} +# set.seed(1234, "L'Ecuyer-CMRG") +# alpha_poisson <- alpha * 10 + +# netlist_poisson <- c(generate_bipartite_collection(nr = nr, nc = nc, pi = pi1, rho = rho1, alpha = alpha_poisson, M = 1, distribution = "poisson"), +# generate_bipartite_collection(nr = nr, nc = nc, pi = pi2, rho = rho2, alpha = alpha_poisson, M = 1, distribution = "poisson")) + +# fitted_bisbmpop_pirho_poisson <- estimate_colBiSBM( +# netlist = netlist_poisson, +# colsbm_model = "pirho", +# nb_run = 1, +# distribution = "poisson", +# global_opts = list( +# verbosity = 3L, +# plot_details = 0, +# nb_cores = parallelly::availableCores(omit = 1) +# ) +# ) +# wrap_plots( +# plot_alpha(alpha_poisson) + labs(caption = "Vrai alpha"), +# plot(fitted_bisbmpop_pirho_poisson$best_fit, type = "meso", values = TRUE, mixture = TRUE)) +``` + +Ici on clone le modèle et on lui donne les bons paramètres pour voir s'il fait +mieux que celui trouvé avant. + +```{r testing_true_params, echo = FALSE} +# Un clone +good_clone <- good_model$clone() + +# Les vrais paramètres +alpha <- data[["alpha"]] +taus <- lapply(seq_along(row_taus), function(m) { + list(row = row_taus[[m]], + col = col_taus[[m]]) +}) +# Pi params +pi1 <- as.vector(data[["pi1"]]) +pi1[1] <- pi1[1] + 1e-9 +rho1 <- rep(0.25, 4) +pi2 <- rep(0.25, 4) +rho2 <- as.vector(data[["rho2"]]) +rho2[3] <- rho2[3] + 1e-9 +pi <- list(list(pi1, rho1), list(pi2, rho2)) + +good_clone[["tau"]] <- taus +good_clone[["pi"]] <- pi +good_clone[["alpha"]] <- alpha + +result_BICL <- c(bad_BICL = bad_model$compute_BICL(), +good_model_BICL = good_model$compute_BICL(), +good_clone_BICL = good_clone$compute_BICL()) +print(result_BICL) +``` + +On vient mettre à la place du 4,4 le modèle avec les vrais paramètres. +Et on espère voir la bonne information se diffuser. +```{r adjusting, echo = FALSE} +fitted_bisbmpop_pirho$model_list[[4,4]] <- good_clone +fit_1_pass <- adjust_colBiSBM(fitted_bisbmpop_pirho, Q = c(4,4), depth = 1L, nb_pass = 1L) ``` \ No newline at end of file diff --git a/code/analysis/investigating/invest_inference_bernoulli_double_ARI.html b/code/analysis/investigating/invest_inference_bernoulli_double_ARI.html new file mode 100644 index 0000000..34b79f4 --- /dev/null +++ b/code/analysis/investigating/invest_inference_bernoulli_double_ARI.html @@ -0,0 +1,5115 @@ + + + + + + + + + + + + + +invest_inference_bernoulli_double_ARI.knit + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
+ + + + + + + +
## here() starts at /home/polarolouis/Documents/Temporary Projects/mia-stage-2024
+ +++++++++++++++++++++++++++++++++++++++++++++++ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
epsilon_alphapi1.1pi1.2pi1.3pi1.4rho2.1rho2.2rho2.3rho2.4repetitionsep_BICLsep_mean_row_ARIsep_mean_col_ARIsep_double_row_ARIsep_double_col_ARIiid_BICLiid_mean_row_ARIiid_mean_col_ARIiid_double_row_ARIiid_double_col_ARIiid_Q1iid_Q2pi_BICLpi_mean_row_ARIpi_mean_col_ARIpi_double_row_ARIpi_double_col_ARIpi_Q1pi_Q2rho_BICLrho_mean_row_ARIrho_mean_col_ARIrho_double_row_ARIrho_double_col_ARIrho_Q1rho_Q2pirho_BICLpirho_mean_row_ARIpirho_mean_col_ARIpirho_double_row_ARIpirho_double_col_ARIpirho_Q1pirho_Q2elapsed_secs
5550.150.00.40.20.40.33333330.33333330.00000000.33333333-17326.431.00000010.63540380.3475734-17372.740.89939830.95772380.46750790.408208445-17347.481.00000000.93345670.53107440.314500764-17326.570.91178110.98829420.51034700.437615345-17328.68110.41189670.4308905551027.1878 secs
5570.150.00.40.20.40.33333330.33333330.33333330.00000002-17661.351.00000010.56110580.3683275-17676.781.00000000.98888161.00000000.982912644-17655.461.00000000.98888161.00000000.982912644-17685.010.99685131.00000000.52023820.483686645-17651.53110.68896360.4588394451677.9255 secs
5680.150.20.00.40.40.00000000.33333330.33333330.33333331-17351.631.00000010.36142460.4932339-17421.390.98745340.95848190.38740320.319087154-17371.861.00000000.96473750.39951570.322912854-17366.990.90158781.00000000.29956100.525362246-17341.95110.45948110.6009543541095.9794 secs
5740.150.20.00.40.40.33333330.33333330.33333330.00000001-17289.301.00000011.00000000.5951142-17334.061.00000000.98561070.36287450.367454744-17292.151.00000001.00000001.00000001.000000044-17320.870.90783311.00000000.46629250.508691246-17291.33110.46304470.390044554877.5671 secs
5860.150.20.40.40.00.00000000.33333330.33333330.33333331-17761.271.00000010.59434570.3751656-17839.000.94460540.93824310.45167270.318922854-17759.481.00000000.99258471.00000000.989990344-17767.480.90699651.00000000.37931480.486338545-17753.26110.48363940.4851668541256.3595 secs
5980.150.40.00.20.40.33333330.33333330.00000000.33333331-17365.911.00000010.34230570.8046233-17381.451.00000001.00000001.00000001.000000044-17392.891.00000000.94928640.51873120.346666064-17392.650.91768140.97715770.40101620.517094146-17368.81110.47630360.4549627551002.0712 secs
6170.150.40.20.00.40.33333330.00000000.33333330.33333332-17275.181.00000010.35401240.6027146-17291.831.00000001.00000001.00000001.000000044-17276.661.00000001.00000001.00000001.000000044-17273.961.00000001.00000001.00000001.000000044-17279.25110.42884280.462798455666.5099 secs
6510.180.00.20.40.40.33333330.00000000.33333330.33333333-16658.301.00000010.36791780.3727995-16758.311.00000000.99300620.44316130.345637254-16685.241.00000000.97260710.52698640.354325064-16690.080.94696021.00000000.32880390.502030346-16663.02110.41433150.429937255918.0173 secs
6530.180.00.20.40.40.33333330.33333330.00000000.33333332-16525.571.00000010.75609130.3756357-16631.780.97240330.99559230.38516510.440941045-16542.311.00000000.88021410.55318960.297570764-16530.550.99376261.00000000.36564310.488640745-16530.44110.46153420.445801655753.9150 secs
6560.180.00.20.40.40.33333330.33333330.33333330.00000002-17111.041.00000010.39415030.3792336-17138.161.00000000.98758401.00000000.980982644-17111.601.00000000.98758401.00000000.980982644-17117.901.00000001.00000001.00000001.000000044-17101.76110.51916680.509608945990.8326 secs
6600.180.00.40.20.40.33333330.00000000.33333330.33333333-16736.341.00000010.78713400.3531989-16745.181.00000001.00000001.00000001.000000044-16724.641.00000001.00000001.00000001.000000044-16735.431.00000001.00000001.00000001.000000044-16728.42110.51840850.397606845881.5467 secs
6610.180.00.40.20.40.33333330.33333330.00000000.33333331-16750.091.00000010.69824250.6363895-16817.570.96746920.99351060.65307040.428250645-16781.250.99494490.98728540.51855190.340262264-16744.850.99151151.00000000.98989781.000000044-16741.69110.69543220.432268745857.1393 secs
6620.180.00.40.20.40.33333330.33333330.00000000.33333332-16361.051.00000010.63714020.5349127-16443.320.91475951.00000000.39927980.352666654-16382.491.00000000.96780370.52684240.329714064-16362.911.00000001.00000001.00000001.000000044-16354.91110.63014580.414378045902.6525 secs
6630.180.00.40.20.40.33333330.33333330.00000000.33333333-16779.091.00000010.60665310.5805976-16842.450.97637460.98720880.56444330.379449145-16808.591.00000000.96044330.47830270.317400064-16793.780.88048200.98586440.27760200.521723846-16773.39110.59993940.3904557451184.5061 secs
6760.180.20.00.40.40.00000000.33333330.33333330.33333331-16703.851.00000010.59263630.6435837-16774.470.93163280.99231110.41638250.511886554-16713.830.98738130.99231110.42950210.511886554-16691.691.00000001.00000001.00000001.000000044-16708.29110.42116880.415525555954.6656 secs
6770.180.20.00.40.40.00000000.33333330.33333330.33333332-16683.471.00000010.41075040.5398208-16773.161.00000000.99493990.41758130.344598054-16708.731.00000000.96780370.53148580.325523064-16707.290.99307341.00000000.36710740.477533846-16688.20110.41758130.405047655863.7348 secs
6780.180.20.00.40.40.00000000.33333330.33333330.33333333-16750.451.00000010.63803910.7694721-16825.110.94730761.00000000.47035310.562238354-16765.221.00000000.99377690.48265740.552841554-16770.780.95150871.00000000.34460630.537737946-16740.48110.48265740.562238354787.8910 secs
6820.180.20.00.40.40.33333330.33333330.33333330.00000001-17067.051.00000010.60839611.0000000-17084.631.00000001.00000001.00000001.000000044-17076.531.00000000.93848700.46748530.317003254-17069.521.00000001.00000001.00000001.000000044-17069.69110.46748530.4394105551251.0973 secs
6840.180.20.00.40.40.33333330.33333330.33333330.00000003-16359.261.00000010.37689920.6320516-16442.710.97910170.97170610.39020560.337594854-16375.461.00000000.93540930.50492170.303822764-16390.250.98768221.00000000.37210380.518180946-16361.35110.45352900.438474155953.8394 secs
6860.180.20.40.00.40.00000000.33333330.33333330.33333332-16594.751.00000010.49990680.3664218-16663.680.99316061.00000000.37724270.343818244-16601.540.97461410.93077200.40534750.306860554-16631.360.98903231.00000000.37026980.503467046-16586.38110.48037520.5050099541113.0697 secs
6880.180.20.40.00.40.33333330.00000000.33333330.33333331-16419.031.00000010.88995180.5834349-16526.611.00000000.93936630.37765290.379106945-16411.151.00000001.00000001.00000001.000000044-16445.390.98584761.00000000.36848490.507819846-16421.33110.46980170.447487555681.8163 secs
6890.180.20.40.00.40.33333330.00000000.33333330.33333332-16508.711.00000010.36466530.7076365-16611.830.92714821.00000000.41847840.372747554-16532.330.97916940.97067070.55462160.374394064-16542.910.86904521.00000000.28318760.644868747-16511.42110.49199700.434109155981.3183 secs
6940.180.20.40.40.00.00000000.33333330.33333330.33333331-16735.751.00000010.38011400.7539287-16739.760.99225121.00000000.98736051.000000044-16730.021.00000001.00000001.00000001.000000044-16754.940.96599431.00000000.31703610.535338446-16739.93110.39141420.4653847551037.5083 secs
6970.180.20.40.40.00.33333330.00000000.33333330.33333331-16546.251.00000010.34857680.6070932-16566.371.00000001.00000001.00000001.000000044-16551.771.00000001.00000001.00000001.000000044-16541.581.00000001.00000001.00000001.000000044-16549.37110.41128740.481441655694.4658 secs
7040.180.40.00.20.40.00000000.33333330.33333330.33333332-16741.691.00000010.34775851.0000000-16874.980.93184270.96375260.65384530.728867364-16772.680.97723100.96375260.70308560.728867364-16737.511.00000001.00000001.00000001.000000044-16750.15110.71473710.797345564740.7993 secs
7070.180.40.00.20.40.33333330.33333330.00000000.33333332-16464.821.00000010.37204300.3962905-16489.801.00000001.00000001.00000001.000000044-16495.871.00000000.91580780.66944890.386937074-16485.830.93809001.00000000.37807690.574273846-16468.01110.40765360.4862882551073.6763 secs
7090.180.40.00.20.40.33333330.33333330.33333330.00000001-16670.481.00000010.86884460.5464027-16689.421.00000001.00000001.00000001.000000044-16693.781.00000000.95197990.52862910.319872464-16676.461.00000001.00000001.00000001.000000044-16672.80110.50090210.457559355749.2522 secs
7160.180.40.00.40.20.33333330.33333330.00000000.33333332-16750.671.00000010.40196061.0000000-16840.691.00000000.96467190.37335690.410286745-16747.671.00000001.00000001.00000001.000000044-16776.240.99325301.00000000.37304850.530714446-16753.19110.48894050.451923155741.2067 secs
7170.180.40.00.40.20.33333330.33333330.00000000.33333333-16553.011.00000010.38078910.5751410-16568.121.00000001.00000001.00000001.000000044-16583.811.00000000.96362180.58553530.337290864-16552.831.00000001.00000001.00000001.000000044-16558.95110.47066000.450348855894.9349 secs
7190.180.40.00.40.20.33333330.33333330.33333330.00000002-16746.681.00000010.67557801.0000000-16760.251.00000001.00000001.00000001.000000044-16767.011.00000000.94372300.48150410.320858764-16779.300.96370971.00000000.39713960.528420446-16750.09110.42226290.4326288551067.0659 secs
7300.180.40.20.40.00.00000000.33333330.33333330.33333331-16706.871.00000010.34814290.3730593-16784.790.91189201.00000000.45906000.638201354-16724.211.00000000.98431570.46429420.599696854-16700.250.99253991.00000000.99146681.000000044-16697.16110.46429420.6382013541014.5758 secs
7390.180.40.40.00.20.00000000.33333330.33333330.33333331-16493.831.00000010.70954860.6138714-16598.891.00000000.94965010.49033070.344336554-16492.831.00000001.00000001.00000001.000000044-16519.250.98324621.00000000.36638030.536658146-16485.69110.46628200.4947522541246.3186 secs
7490.180.40.40.20.00.00000000.33333330.33333330.33333332-16717.471.00000010.60264020.5794619-16778.110.98062461.00000000.41298230.564861454-16731.011.00000000.98313090.41264210.537301354-16747.020.98796321.00000000.47323580.757413946-16709.14110.41264210.5648614541029.9068 secs
7530.180.40.40.20.00.33333330.00000000.33333330.33333333-16456.781.00000010.39204730.3930401-16477.151.00000001.00000001.00000001.000000044-16477.150.98960190.97941970.53101850.368538464-16474.420.91693431.00000000.32594640.551064046-16458.96110.48026300.4725277551174.7415 secs
7570.210.00.20.40.40.33333330.00000000.33333330.33333331-15626.071.00000010.37957960.3743909-15737.600.99730711.00000000.37739910.418716145-15655.531.00000000.94621010.52859250.370062964-15643.170.89858821.00000000.76450330.773580746-15619.20110.64431050.4694788451136.3157 secs
7580.210.00.20.40.40.33333330.00000000.33333330.33333332-15463.681.00000010.37864360.3829776-15603.100.95694770.97317630.66351310.706612546-15464.321.00000001.00000001.00000001.000000044-15482.450.92921451.00000000.33087580.543339646-15466.96110.42402390.4368330551323.6269 secs
7610.210.00.20.40.40.33333330.33333330.00000000.33333332-15747.531.00000010.48626120.6254415-15876.850.92086090.99318970.37203700.440315655-15778.521.00000000.97434140.49310170.359189564-15766.560.91464821.00000000.30722350.551753246-15751.86110.43910250.463566455875.4299 secs
7620.210.00.20.40.40.33333330.33333330.00000000.33333333-15593.911.00000010.47102420.3537180-15665.150.99190351.00000000.44663810.449063045-15622.111.00000000.97096410.48783750.349545864-15601.950.99190351.00000000.44663810.449063045-15584.32110.46107050.449063045836.4729 secs
7650.210.00.20.40.40.33333330.33333330.33333330.00000003-15862.041.00000010.70492920.3695618-15930.651.00000000.97510630.62402370.400163345-15856.631.00000001.00000001.00000001.000000044-15860.261.00000001.00000001.00000001.000000044-15854.26110.62402370.400163345983.0682 secs
7660.210.00.40.20.40.33333330.00000000.33333330.33333331-15804.861.00000010.63872060.8020132-15821.921.00000001.00000001.00000001.000000044-15806.321.00000001.00000001.00000001.000000044-15829.650.95058441.00000000.36980450.536873446-15807.62110.47276300.424336055934.6113 secs
7700.210.00.40.20.40.33333330.33333330.00000000.33333332-15530.411.00000010.57922780.3404878-15548.001.00000001.00000001.00000001.000000044-15546.901.00000000.93262710.55174630.302490864-15559.270.97562011.00000000.37870620.478451646-15534.05110.41775690.447387455792.1594 secs
7780.210.00.40.40.20.33333330.33333330.00000000.33333331-15768.661.00000011.00000000.3455340-15777.031.00000001.00000001.00000001.000000044-15758.761.00000001.00000001.00000001.000000044-15766.781.00000001.00000001.00000001.000000044-15771.88110.51993840.467736055671.9395 secs
7850.210.20.00.40.40.00000000.33333330.33333330.33333332-16078.521.00000010.36042910.3492673-16088.931.00000001.00000001.00000001.000000044-16072.541.00000001.00000001.00000001.000000044-16074.371.00000001.00000001.00000001.000000044-16084.76110.41112350.466658755542.3221 secs
7860.210.20.00.40.40.00000000.33333330.33333330.33333333-15832.091.00000010.36197540.3578411-15835.051.00000001.00000001.00000001.000000044-15821.971.00000001.00000001.00000001.000000044-15825.241.00000001.00000001.00000001.000000044-15835.53110.42424630.444655855914.9031 secs
7870.210.20.00.40.40.33333330.33333330.00000000.33333331-15212.441.00000010.38838650.3716715-15402.150.93577860.99434470.37348010.426077855-15210.071.00000000.99511211.00000000.992510744-15214.261.00000001.00000001.00000001.000000044-15215.99110.45318930.4651698551296.5766 secs
7890.210.20.00.40.40.33333330.33333330.00000000.33333333-15498.371.00000010.34845800.3643496-15627.350.95775650.99507230.36427260.491459846-15498.721.00000001.00000001.00000001.000000044-15528.250.95775651.00000000.36427260.496554446-15502.13110.43758270.443057455737.4774 secs
7960.210.20.40.00.40.33333330.00000000.33333330.33333331-15614.991.00000010.45296440.3751147-15736.831.00000000.87879090.39076290.406253046-15608.601.00000001.00000001.00000001.000000044-15618.931.00000001.00000001.00000001.000000044-15617.88110.53031320.442350855784.1966 secs
7980.210.20.40.00.40.33333330.00000000.33333330.33333333-15490.841.00000010.69149870.3588368-15613.380.86728681.00000000.43420480.364829264-15526.281.00000000.97853790.60209260.365670464-15494.271.00000001.00000001.00000001.000000044-15493.76110.51168150.460683855948.9699 secs
8070.210.20.40.40.00.33333330.00000000.33333330.33333333-15433.841.00000010.36630090.3608174-15534.910.88616320.95670200.30044370.448540746-15427.001.00000000.99393711.00000000.991209544-15433.301.00000001.00000001.00000001.000000044-15437.33110.48452280.428469955926.2239 secs
8110.210.40.00.20.40.00000000.33333330.33333330.33333331-15539.781.00000010.37178390.3481736-15549.631.00000001.00000001.00000001.000000044-15534.411.00000001.00000001.00000001.000000044-15566.220.94340591.00000000.33464110.497366246-15544.94110.47272900.4272450551009.2387 secs
8160.210.40.00.20.40.33333330.33333330.00000000.33333333-15286.791.00000010.36373970.6122522-15307.871.00000001.00000001.00000001.000000044-15287.671.00000001.00000001.00000001.000000044-15321.370.97627491.00000000.37893090.548292846-15290.40110.40897170.442852955788.0950 secs
8200.210.40.00.40.20.00000000.33333330.33333330.33333331-15494.311.00000010.36827140.6406356-15628.661.00000000.97903520.50389210.625057864-15492.201.00000001.00000001.00000001.000000044-15494.151.00000001.00000001.00000001.000000044-15485.49110.44124230.575656654762.1558 secs
8210.210.40.00.40.20.00000000.33333330.33333330.33333332-15383.171.00000010.40888730.5291564-15400.531.00000001.00000001.00000001.000000044-15415.141.00000000.94279100.61090270.332072364-15389.971.00000001.00000001.00000001.000000044-15387.64110.49752490.435491155526.9670 secs
8310.210.40.20.00.40.00000000.33333330.33333330.33333333-15412.601.00000010.39624630.5461138-15429.111.00000001.00000001.00000001.000000044-15406.121.00000001.00000001.00000001.000000044-15414.681.00000001.00000001.00000001.000000044-15415.31110.46865320.451543155656.4451 secs
8330.210.40.20.00.40.33333330.00000000.33333330.33333332-15589.521.00000010.40202180.5388090-15608.501.00000001.00000001.00000001.000000044-15587.051.00000001.00000001.00000001.000000044-15589.401.00000001.00000001.00000001.000000044-15592.05110.41251600.430098155324.3734 secs
8400.210.40.20.40.00.00000000.33333330.33333330.33333333-15632.501.00000010.40069620.3416710-15645.880.99287381.00000000.98933721.000000044-15626.501.00000001.00000001.00000001.000000044-15632.160.99287381.00000000.98933721.000000044-15625.87110.48048590.579365854824.6625 secs
8410.210.40.20.40.00.33333330.00000000.33333330.33333331-15123.661.00000010.41631710.5948807-15221.730.98290301.00000000.36122590.435338645-15115.081.00000001.00000001.00000001.000000044-15127.181.00000001.00000001.00000001.000000044-15128.29110.40451310.435338655407.0289 secs
8440.210.40.20.40.00.33333330.33333330.00000000.33333331-15388.291.00000010.60941820.7763398-15526.170.92565400.99331240.33411330.572087446-15390.541.00000001.00000001.00000001.000000044-15389.501.00000001.00000001.00000001.000000044-15392.24110.42445660.4937706551353.1661 secs
8500.210.40.40.00.20.33333330.00000000.33333330.33333331-15111.460.98942710.35509860.3559548-15121.431.00000001.00000001.00000001.000000044-15137.500.98942700.95856520.54096980.363623664-15107.551.00000001.00000001.00000001.000000044-15114.79110.44020790.432435155539.5275 secs
8510.210.40.40.00.20.33333330.00000000.33333330.33333332-15179.791.00000010.88957770.5178101-15190.921.00000001.00000001.00000001.000000044-15214.081.00000000.98825670.49821280.354706364-15212.390.96680351.00000000.34042150.510251146-15182.46110.40144130.429723655971.0546 secs
8570.210.40.40.20.00.00000000.33333330.33333330.33333332-15741.241.00000010.36341210.5862905-15818.560.98872291.00000000.41000920.575549254-15749.311.00000001.00000001.00000001.000000044-15740.061.00000001.00000001.00000001.000000044-15744.26110.46327190.486346355319.9372 secs
8580.210.40.40.20.00.00000000.33333330.33333330.33333333-15668.761.00000010.55343960.5643967-15756.930.98824291.00000000.50151900.594245254-15692.431.00000000.95567260.56772590.344554864-15691.990.97267901.00000000.35972560.534729846-15673.72110.49509170.447454755913.9304 secs
8600.210.40.40.20.00.33333330.00000000.33333330.33333332-15484.421.00000010.60033370.8155280-15497.491.00000001.00000001.00000001.000000044-15482.121.00000001.00000001.00000001.000000044-15482.101.00000001.00000001.00000001.000000044-15487.04110.45339690.435287355619.5330 secs
8610.210.40.40.20.00.33333330.00000000.33333330.33333333-15225.711.00000010.56063130.7916553-15351.520.93205311.00000000.48375450.341423864-15259.221.00000001.00000000.76929520.601626564-15257.140.98011001.00000000.35601590.478964146-15229.63110.45381110.4231780551124.2132 secs
8650.240.00.20.40.40.33333330.00000000.33333330.33333331-14412.391.00000010.69463270.3622749-14430.251.00000001.00000001.00000001.000000044-14448.871.00000000.95155440.63991430.321269374-14432.640.90770061.00000000.31517410.537168446-14406.22110.69463270.439828445994.6232 secs
8660.240.00.20.40.40.33333330.00000000.33333330.33333332-14513.011.00000010.34909300.5707612-14646.660.99616050.95008640.35226430.475256746-14503.781.00000001.00000001.00000001.000000044-14513.671.00000001.00000001.00000001.000000044-14519.69110.64247700.761830846616.2170 secs
8670.240.00.20.40.40.33333330.00000000.33333330.33333333-14336.091.00000010.36535480.6404642-14403.351.00000000.99215720.66591710.474065645-14352.031.00000000.93823560.48473890.336167864-14351.270.94052641.00000000.30387730.515723946-14338.69110.41967500.438852655366.3348 secs
8690.240.00.20.40.40.33333330.33333330.00000000.33333332-14245.511.00000010.55034540.5245335-14379.580.92508921.00000000.50271760.348036664-14261.211.00000000.92968690.58704900.293281264-14246.681.00000001.00000001.00000001.000000044-14249.24110.51144820.459229855778.0212 secs
8700.240.00.20.40.40.33333330.33333330.00000000.33333333-13950.581.00000010.49313730.5875910-14040.830.99324321.00000000.54020040.435724745-13941.351.00000001.00000001.00000001.000000044-13973.130.92063241.00000000.35327610.510059746-13957.76110.83095330.841971764885.8248 secs
8730.240.00.20.40.40.33333330.33333330.33333330.00000003-14170.601.00000010.55916070.3651567-14280.600.86458180.97037340.28437980.447444546-14165.121.00000001.00000001.00000001.000000044-14174.901.00000001.00000001.00000001.000000044-14173.76110.50979960.422623255690.2595 secs
8740.240.00.40.20.40.33333330.00000000.33333330.33333331-13797.501.00000010.38507330.6310183-13811.121.00000001.00000001.00000001.000000044-13794.621.00000001.00000001.00000001.000000044-13794.181.00000001.00000001.00000001.000000044-13812.35110.49579530.454296265475.9064 secs
8760.240.00.40.20.40.33333330.00000000.33333330.33333333-14090.191.00000010.37829180.6135664-14108.901.00000001.00000001.00000001.000000044-14087.841.00000001.00000001.00000001.000000044-14121.130.97554031.00000000.38009760.531460046-14095.42110.46816350.455231255822.1496 secs
8770.240.00.40.20.40.33333330.33333330.00000000.33333331-13727.871.00000010.63309890.6157026-13808.170.98577390.99489160.62842790.434012545-13758.561.00000000.98425710.53673580.346551364-13747.300.91234811.00000000.41752650.495767646-13719.44110.64920210.438169445853.9934 secs
8780.240.00.40.20.40.33333330.33333330.00000000.33333332-14097.321.00000011.00000000.5257670-14109.031.00000001.00000001.00000001.000000044-14091.301.00000001.00000001.00000001.000000044-14097.821.00000001.00000001.00000001.000000044-14101.34110.43273960.434294355702.5250 secs
8820.240.00.40.20.40.33333330.33333330.33333330.00000003-14594.461.00000010.38781130.3960076-14674.851.00000000.98948690.62457050.529518145-14598.851.00000001.00000001.00000001.000000044-14606.930.99216351.00000000.61495940.517833145-14583.47110.62457050.517833145678.1629 secs
8910.240.00.40.40.20.33333330.33333330.33333330.00000003-13873.801.00000010.62125290.3823466-13950.041.00000000.98462170.64079460.416066345-13888.691.00000000.93203710.53767890.303553464-13890.540.91783760.99030390.37978620.546471546-13877.64110.43425650.444985155415.2216 secs
8920.240.20.00.40.40.00000000.33333330.33333330.33333331-14347.321.00000010.35451090.5572280-14363.231.00000001.00000001.00000001.000000044-14372.671.00000000.98118240.50870400.361194064-14344.981.00000001.00000001.00000001.000000044-14362.05110.45805170.553829956574.2524 secs
8930.240.20.00.40.40.00000000.33333330.33333330.33333332-14362.291.00000010.42122920.6082431-14457.250.99283081.00000000.50416130.552090354-14384.171.00000001.00000000.50651160.552090354-14381.830.92098611.00000000.33323430.543872546-14364.87110.41148870.429634155746.0737 secs
8940.240.20.00.40.40.00000000.33333330.33333330.33333333-14894.771.00000010.47155830.5765483-14907.831.00000001.00000001.00000001.000000044-14923.571.00000000.94194620.52294590.338782264-14890.641.00000001.00000001.00000001.000000044-14897.55110.42390490.486941355415.8795 secs
8980.240.20.00.40.40.33333330.33333330.33333330.00000001-14067.021.00000010.65686190.3481523-14188.051.00000000.92168280.68459550.450638655-14093.631.00000000.97590470.81654870.611514164-14068.601.00000001.00000001.00000001.000000044-14070.15110.50111610.432090655804.7793 secs
8990.240.20.00.40.40.33333330.33333330.33333330.00000002-14124.491.00000010.35676720.3606921-14255.781.00000000.96327960.36060000.481774746-14162.571.00000000.93449150.63215610.313910574-14145.210.94416481.00000000.32594020.522866946-14126.72110.71416170.442779755645.6266 secs
9000.240.20.00.40.40.33333330.33333330.33333330.00000003-13501.011.00000010.37353790.3522930-13517.371.00000001.00000001.00000001.000000044-13494.981.00000001.00000001.00000001.000000044-13503.811.00000001.00000001.00000001.000000044-13506.90110.65607110.418807055688.9274 secs
9010.240.20.40.00.40.00000000.33333330.33333330.33333331-13955.521.00000010.86918970.3447676-14034.481.00000001.00000000.46600140.561945154-13972.171.00000001.00000000.46600140.561945154-13987.380.81545991.00000000.26958900.652524847-13945.65110.46600140.561945154695.1306 secs
9020.240.20.40.00.40.00000000.33333330.33333330.33333332-14380.771.00000010.37598700.3699391-14502.950.98367080.98851600.48282830.369585764-14396.991.00000000.94483180.54057720.316500064-14378.351.00000001.00000001.00000001.000000044-14384.48110.41986340.444472055261.8858 secs
9040.240.20.40.00.40.33333330.00000000.33333330.33333331-13685.151.00000010.59738280.3773899-13694.051.00000001.00000001.00000001.000000044-13676.991.00000001.00000001.00000001.000000044-13681.731.00000001.00000001.00000001.000000044-13692.67110.51781480.377068264255.6611 secs
9060.240.20.40.00.40.33333330.00000000.33333330.33333333-13811.661.00000010.50440900.3566067-13943.600.96666350.99529090.39935380.408585855-13804.051.00000001.00000001.00000001.000000044-13804.181.00000001.00000001.00000001.000000044-13814.99110.48900340.438873155530.4463 secs
9110.240.20.40.40.00.00000000.33333330.33333330.33333332-14535.641.00000010.35892340.3436176-14655.041.00000000.98282210.37248540.464780446-14542.531.00000000.98795250.49713290.549766054-14565.720.98735561.00000000.35123570.503367146-14530.05110.49713290.567989054196.2896 secs
9140.240.20.40.40.00.33333330.00000000.33333330.33333332-13778.371.00000010.77066310.5829787-13900.260.99303071.00000000.44637540.368228054-13810.491.00000001.00000000.53841410.368228064-13772.001.00000001.00000001.00000001.000000044-13782.02110.45360480.412960555705.7804 secs
9150.240.20.40.40.00.33333330.00000000.33333330.33333333-14332.291.00000010.40288630.5683112-14466.180.98575740.97707810.38839900.476344746-14349.871.00000000.94000120.54337140.312666364-14338.671.00000001.00000001.00000001.000000044-14336.31110.49729600.421776155794.6852 secs
9180.240.20.40.40.00.33333330.33333330.00000000.33333333-13927.611.00000010.54051130.5038732-14065.071.00000000.99550820.68880970.460370155-13928.621.00000001.00000001.00000001.000000044-13940.330.90659921.00000000.31799050.545598546-13932.66110.68880970.462705255622.3278 secs
9300.240.40.00.40.20.00000000.33333330.33333330.33333333-14314.401.00000010.38725220.5778284-14335.101.00000001.00000001.00000001.000000044-14332.191.00000001.00000000.41616440.565373454-14348.250.89711331.00000000.28584970.634025847-14328.90110.48926140.432672465294.2685 secs
9330.240.40.00.40.20.33333330.33333330.00000000.33333333-13642.031.00000010.36317990.7816660-13659.571.00000001.00000001.00000001.000000044-13637.381.00000001.00000001.00000001.000000044-13642.141.00000001.00000001.00000001.000000044-13644.44110.46866200.422879855377.8753 secs
9380.240.40.20.00.40.00000000.33333330.33333330.33333332-14040.521.00000010.48617690.3409453-14055.521.00000001.00000001.00000001.000000044-14030.541.00000001.00000001.00000001.000000044-14046.261.00000001.00000001.00000001.000000044-14043.66110.43410500.426810355632.0982 secs
9480.240.40.20.40.00.00000000.33333330.33333330.33333333-13950.391.00000010.74212380.3552372-13958.391.00000001.00000001.00000001.000000044-13964.601.00000000.99378790.46519830.606072254-13945.931.00000001.00000001.00000001.000000044-13953.01110.47556490.451303155468.0121 secs
9490.240.40.20.40.00.33333330.00000000.33333330.33333331-13325.041.00000010.60792650.7804317-13433.070.99491201.00000000.38766230.444557145-13362.201.00000000.95441250.66472980.437447974-13327.411.00000001.00000001.00000001.000000044-13327.10110.47916900.425704855371.5319 secs
9560.240.40.40.00.20.00000000.33333330.33333330.33333332-13711.721.00000010.38067390.3821680-13795.031.00000000.99430250.45909230.528906054-13717.601.00000001.00000001.00000001.000000044-13747.770.90606061.00000000.39951070.646393547-13716.69110.48233370.460605655417.2319 secs
9580.240.40.40.00.20.33333330.00000000.33333330.33333331-13455.861.00000010.37744410.5552071-13476.601.00000001.00000001.00000001.000000044-13493.401.00000000.91594510.64504600.417166174-13493.440.99281171.00000000.37565010.522823846-13460.29110.47214850.420902855240.1433 secs
9680.240.40.40.20.00.33333330.00000000.33333330.33333332-13987.921.00000010.37272970.7606305-14004.781.00000001.00000001.00000001.000000044-13982.961.00000001.00000001.00000001.000000044-14020.030.99585711.00000000.36302750.553501946-13993.70110.60900770.443843255654.0868 secs
+ + + + +
+ + + + + + + + + + + + + + + diff --git a/code/analysis/investigating/profiling_parallelization_speed.R b/code/analysis/investigating/profiling_parallelization_speed.R new file mode 100644 index 0000000..cee0216 --- /dev/null +++ b/code/analysis/investigating/profiling_parallelization_speed.R @@ -0,0 +1,42 @@ +library(colSBM) +library(here) + +set.seed(1234, "L'Ecuyer-CMRG") +prof_folder <- file.path(here(), "code", "results", "investigating", "profiling") +if (!dir.exists(prof_folder)) { + dir.create(prof_folder) +} +prof_file <- file.path(prof_folder, "no-parallel-asis.out") + +nr <- 70 +nc <- 50 + +pi <- c(0.3, 0.2, 0.5) +rho <- c(0.25, 0.15, 0.6) + +alpha <- matrix(c( + 0.75, 0.4, 0.1, + 0.2, 0.05, 0.05, + 0.1, 0.05, 0.3 +), nrow = 3) + +netlist <- generate_bipartite_collection( + nr = nr, nc = nc, + pi = pi, rho = rho, + alpha = alpha, + M = 4 +) + +Rprof(filename = prof_file, memory.profiling = TRUE) +fit <- estimate_colBiSBM( + netlist = netlist, + colsbm_model = "pirho", + global_opts = list(verbosity = 1L, backend = "no_mc") +) +Rprof(NULL) + +for (filename in list.files(prof_folder)) { + print(filename) + prof_sum <- summaryRprof(file.path(prof_folder, filename), memory = "both") + print(prof_sum[["by.total"]][1, ]) +} diff --git a/code/results/simulations/inference/bernoulli/bernoulli_inference_24-04-2024_17-14-46_1-1.Rds b/code/results/simulations/inference/bernoulli/bernoulli_inference_24-04-2024_17-14-46_1-1.Rds new file mode 100644 index 0000000..a76d99a Binary files /dev/null and b/code/results/simulations/inference/bernoulli/bernoulli_inference_24-04-2024_17-14-46_1-1.Rds differ