Simulation Experiment Recipe

Objectives

In Wang et al. (2023), high-throughput single-cell gene-silencing experiments were performed to assess the presence of epistasis (i.e,. non-additive interaction effects) in the genetic regulation of cardiac structure. More specifically, these single-cell experiments were designed to assess whether the genes pairs CCDC141 and TTN as well as CCDC141 and IGF1R interact non-additively to impact the cell diameter of human induced pluriopotent stem cell-derived cardiomyocytes. Mathematically, this boils down to studying the following regression problem: for a gene pair, say gene 1 and gene 2, consider

\[y_i = \beta_0 + \beta_1 x_{i1} + \beta_2 x_{i2} + \beta_{12} x_{i1} x_{i2} + \epsilon_i,\]

where \(y_i\) is the diameter (\(\mu m\)) of cell \(i\), \(x_{ij}\) is an indicator of whether cell \(i\) belonged to an experimental batch where gene \(j\) was silenced, and \(\epsilon_i\) is the random error or noise term for cell \(i\). Under this model, the quantity of interest is \(\beta_{12}\), or the coefficient corresponding to the non-additive interaction effect between genes 1 and 2, and we would like to assess the null hypothesis that \(\beta_{12} = 0\) (i.e., no epistatic effect) versus the alternative hypothesis that \(\beta_{12} \neq 0\) (i.e., a non-zero epistatic effect).

One possible complication in the statistical analysis of gene-silencing experiments, however, is that the gene-silencing efficiency may not be 100% and may vary across different experimental batches of cells. In this simulation study, we aim to assess whether variation in silencing efficiencies lead to spurious epistasis signals.

Overall, we found unsurprisingly that the signal-to-noise ratio plays a role in understanding the effect of varying silencing efficiencies. But notably, in the very low signal-to-noise regimes that we expect to be realistic in practice, we found that the differences in our observed gene-silencing efficiencies do not typically lead to high false positive rates. Furthermore, when the strength of the interaction effect is small but not 0, the resulting quantile regression p-values are frequently not significant. We realize however that this simulation models the gene-silencing efficiencies in a particular way, and model misspecification of not only the silencing efficiencies but the model form itself is rightfully a concern. We thus simply view this simulation study as a helpful starting point and acknowledge that additional simulations can be done to explore these avenues.

In what follows, we provide additional details regarding the simulation setup (see tabs above) as well as the results of the simulation study (see section tabs on the left side panel).

Data Generation

CCDC141-IGF1R, Diseased

To simulate data from batches with different gene-silencing efficiencies, we generated data via the following model:

\[y_i = \beta_0 + (\beta_1 x_{i1} + \beta_2 x_{i2} + \beta_{12} x_{i1} x_{i2}) \cdot z_i + \epsilon_i,\]

where \(y_i\) is the diameter (\(\mu m\)) of cell \(i\), \(x_{ij}\) is an indicator of whether cell \(i\) belonged to an experimental batch where gene \(j\) was silenced, \(z_i\) is an indicator of whether the gene(s) were successfully silenced in cell \(i\), generated via \(z_i \mid x_{i1} = j, x_{i2} = k \stackrel{ind.}{\sim} Bernoulli(\pi_{jk})\), and \(\epsilon_i \stackrel{iid}{\sim} N(0, \sigma^2)\) is the random error or noise term for cell \(i\).

We simulated \(n = 500\) cells for each of the gene-silencing conditions, that is,

  • 500 scrambled control cells (i.e., \(x_{i1} = x_{i2} = 0\)),
  • 500 cells wherein gene 1 was silenced (i.e., \(x_{i1} = 1\) and \(x_{i2} = 0\)),
  • 500 cells wherein gene 2 was silenced (i.e., \(x_{i1} = 0\) and \(x_{i2} = 1\)), and
  • 500 cells wherein gene 1 and gene 2 were silenced (i.e., \(x_{i1} = x_{i2} = 1\)).

Unless otherwise specified, simulations were run across 100 replicates, and the following simulation parameter values were used:

  • \(\beta_0 = 0\)
  • \(\beta_1 = -1\)
  • \(\beta_2 = -1\)
  • \(\beta_{12} = 0\)
  • \(\sigma^2 = 1\)

We then ran simulations that varied across the interaction effect strength \(\beta_{12} \in \{-1, -0.5, -0.2, -0.1, 0\}\) and the noise level \(\sigma \in \{0.1, 0.2, 0.5, 1, 1.5, 2\}\).

So that we may best connect to the observed experimental data in Wang et al. (2023), we took the observed gene-silencing efficiencies from the experiments and used these values for the simulated gene-silencing efficiencies \(\{\pi_{01}, \pi_{10}, \pi_{11}\}\). Since there were four scenarios of interest (i.e., the CCDC141-TTN interaction in the healthy and diseased (MYH7-R403Q) cell lines and the CCDC141-IGF1R interaction in the healthy and diseased (MYH7-R403Q) cell lines), we repeated the simulation study under four different choices of \(\{\pi_{01}, \pi_{10}, \pi_{11}\}\).

For this particular data-generating process, we set the simulated gene-silencing efficiencies to be those from the CCDC141-IGF1R experiment in the healthy cell line so that

  • \(\pi_{01} = 0.92\)
  • \(\pi_{10} = 0.83\)
  • \(\pi_{11} = 0.695\)

Remark: We note that the above model is a simplistic interpretation of how gene-silencing efficiencies may effect the response of interest and that there are other ways to incorporate gene-silencing efficiencies into the model that may be more realistic. We view this simulation study as a helpful starting point for understanding the effect of varying gene-silencing efficiencies on epistasis modeling.

Function

## function(n, beta0, beta1, beta2, beta12,
##                                eff1 = 1, eff2 = 1, int_eff = 1, 
##                                err_fun = rnorm, ...) {
##   eff1_vec <- rbinom(n, 1, eff1)
##   eff2_vec <- rbinom(n, 1, eff2)
##   int_eff_vec <- rbinom(n, 1, int_eff)
##   
##   y0 <- beta0 + err_fun(n, ...)
##   y1 <- beta0 + beta1 * eff1_vec + err_fun(n, ...)
##   y2 <- beta0 + beta2 * eff2_vec + err_fun(n, ...)
##   y12 <- beta0 + (beta1 + beta2 + beta12) * int_eff_vec + err_fun(n, ...)
##   
##   x <- matrix(0, nrow = 4 * n, ncol = 2)
##   x[(n + 1):(2 * n), 1] <- 1
##   x[(2 * n + 1):(3 * n), 2] <- 1
##   x[(3 * n + 1):(4 * n), 1:2] <- 1
##   y <- c(y0, y1, y2, y12)
##   
##   out <- list(
##     x = x,
##     y = y,
##     eff1_vec = eff1_vec,
##     eff2_vec = eff2_vec,
##     int_eff_vec = int_eff_vec
##   )
## }

Input Parameters

## $n
## [1] 500
## 
## $beta0
## [1] 0
## 
## $beta1
## [1] -1
## 
## $beta2
## [1] -1
## 
## $beta12
## [1] 0
## 
## $eff1
## [1] 0.92
## 
## $eff2
## [1] 0.83
## 
## $int_eff
## [1] 0.695
## 
## $err_fun
## function (n, mean = 0, sd = 1) 
## .Call(C_rnorm, n, mean, sd)
## <bytecode: 0x7feab78c3da0>
## <environment: namespace:stats>
## 
## $sd
## [1] 1

CCDC141-IGF1R, Healthy

In this data-generating process, we simulated data as described above but set the simulated gene-silencing efficiencies to be the observed efficiencies from the CCDC141-IGF1R experiment in the diseased cell line:

  • \(\pi_{01} = 0.66\)
  • \(\pi_{10} = 0.9\)
  • \(\pi_{11} = 0.725\)

Function

## function(n, beta0, beta1, beta2, beta12,
##                                eff1 = 1, eff2 = 1, int_eff = 1, 
##                                err_fun = rnorm, ...) {
##   eff1_vec <- rbinom(n, 1, eff1)
##   eff2_vec <- rbinom(n, 1, eff2)
##   int_eff_vec <- rbinom(n, 1, int_eff)
##   
##   y0 <- beta0 + err_fun(n, ...)
##   y1 <- beta0 + beta1 * eff1_vec + err_fun(n, ...)
##   y2 <- beta0 + beta2 * eff2_vec + err_fun(n, ...)
##   y12 <- beta0 + (beta1 + beta2 + beta12) * int_eff_vec + err_fun(n, ...)
##   
##   x <- matrix(0, nrow = 4 * n, ncol = 2)
##   x[(n + 1):(2 * n), 1] <- 1
##   x[(2 * n + 1):(3 * n), 2] <- 1
##   x[(3 * n + 1):(4 * n), 1:2] <- 1
##   y <- c(y0, y1, y2, y12)
##   
##   out <- list(
##     x = x,
##     y = y,
##     eff1_vec = eff1_vec,
##     eff2_vec = eff2_vec,
##     int_eff_vec = int_eff_vec
##   )
## }

Input Parameters

## $n
## [1] 500
## 
## $beta0
## [1] 0
## 
## $beta1
## [1] -1
## 
## $beta2
## [1] -1
## 
## $beta12
## [1] 0
## 
## $eff1
## [1] 0.66
## 
## $eff2
## [1] 0.9
## 
## $int_eff
## [1] 0.725
## 
## $err_fun
## function (n, mean = 0, sd = 1) 
## .Call(C_rnorm, n, mean, sd)
## <bytecode: 0x7feab72fd800>
## <environment: namespace:stats>
## 
## $sd
## [1] 1

CCDC141-TTN, Diseased

In this data-generating process, we simulated data as described above but set the simulated gene-silencing efficiencies to be the observed efficiencies from the CCDC141-TTN experiment in the diseased cell line:

  • \(\pi_{01} = 0.92\)
  • \(\pi_{10} = 0.99\)
  • \(\pi_{11} = 0.92\)

Function

## function(n, beta0, beta1, beta2, beta12,
##                                eff1 = 1, eff2 = 1, int_eff = 1, 
##                                err_fun = rnorm, ...) {
##   eff1_vec <- rbinom(n, 1, eff1)
##   eff2_vec <- rbinom(n, 1, eff2)
##   int_eff_vec <- rbinom(n, 1, int_eff)
##   
##   y0 <- beta0 + err_fun(n, ...)
##   y1 <- beta0 + beta1 * eff1_vec + err_fun(n, ...)
##   y2 <- beta0 + beta2 * eff2_vec + err_fun(n, ...)
##   y12 <- beta0 + (beta1 + beta2 + beta12) * int_eff_vec + err_fun(n, ...)
##   
##   x <- matrix(0, nrow = 4 * n, ncol = 2)
##   x[(n + 1):(2 * n), 1] <- 1
##   x[(2 * n + 1):(3 * n), 2] <- 1
##   x[(3 * n + 1):(4 * n), 1:2] <- 1
##   y <- c(y0, y1, y2, y12)
##   
##   out <- list(
##     x = x,
##     y = y,
##     eff1_vec = eff1_vec,
##     eff2_vec = eff2_vec,
##     int_eff_vec = int_eff_vec
##   )
## }

Input Parameters

## $n
## [1] 500
## 
## $beta0
## [1] 0
## 
## $beta1
## [1] -1
## 
## $beta2
## [1] -1
## 
## $beta12
## [1] 0
## 
## $eff1
## [1] 0.92
## 
## $eff2
## [1] 0.99
## 
## $int_eff
## [1] 0.92
## 
## $err_fun
## function (n, mean = 0, sd = 1) 
## .Call(C_rnorm, n, mean, sd)
## <bytecode: 0x7feab75fcfd8>
## <environment: namespace:stats>
## 
## $sd
## [1] 1

CCDC141-TTN, Healthy

In this data-generating process, we simulated data as described above but set the simulated gene-silencing efficiencies to be the observed efficiencies from the CCDC141-TTN experiment in the healthy cell line:

  • \(\pi_{01} = 0.66\)
  • \(\pi_{10} = 0.8\)
  • \(\pi_{11} = 0.93\)

Function

## function(n, beta0, beta1, beta2, beta12,
##                                eff1 = 1, eff2 = 1, int_eff = 1, 
##                                err_fun = rnorm, ...) {
##   eff1_vec <- rbinom(n, 1, eff1)
##   eff2_vec <- rbinom(n, 1, eff2)
##   int_eff_vec <- rbinom(n, 1, int_eff)
##   
##   y0 <- beta0 + err_fun(n, ...)
##   y1 <- beta0 + beta1 * eff1_vec + err_fun(n, ...)
##   y2 <- beta0 + beta2 * eff2_vec + err_fun(n, ...)
##   y12 <- beta0 + (beta1 + beta2 + beta12) * int_eff_vec + err_fun(n, ...)
##   
##   x <- matrix(0, nrow = 4 * n, ncol = 2)
##   x[(n + 1):(2 * n), 1] <- 1
##   x[(2 * n + 1):(3 * n), 2] <- 1
##   x[(3 * n + 1):(4 * n), 1:2] <- 1
##   y <- c(y0, y1, y2, y12)
##   
##   out <- list(
##     x = x,
##     y = y,
##     eff1_vec = eff1_vec,
##     eff2_vec = eff2_vec,
##     int_eff_vec = int_eff_vec
##   )
## }

Input Parameters

## $n
## [1] 500
## 
## $beta0
## [1] 0
## 
## $beta1
## [1] -1
## 
## $beta2
## [1] -1
## 
## $beta12
## [1] 0
## 
## $eff1
## [1] 0.66
## 
## $eff2
## [1] 0.8
## 
## $int_eff
## [1] 0.93
## 
## $err_fun
## function (n, mean = 0, sd = 1) 
## .Call(C_rnorm, n, mean, sd)
## <bytecode: 0x7feab7ce5bc0>
## <environment: namespace:stats>
## 
## $sd
## [1] 1

Efficiency = 1

In this data-generating process, we simulated data as described above but set the simulated gene-silencing efficiencies to be 1. That is, \(\pi_{01} = \pi_{10} = \pi_{11} = 1\). This data-generating process is akin to the setting where the gene-silencing efficiences are perfect, and there is no variation in gene-silencing efficiences across different batches/gene-silencing conditions.

Function

## function(n, beta0, beta1, beta2, beta12,
##                                eff1 = 1, eff2 = 1, int_eff = 1, 
##                                err_fun = rnorm, ...) {
##   eff1_vec <- rbinom(n, 1, eff1)
##   eff2_vec <- rbinom(n, 1, eff2)
##   int_eff_vec <- rbinom(n, 1, int_eff)
##   
##   y0 <- beta0 + err_fun(n, ...)
##   y1 <- beta0 + beta1 * eff1_vec + err_fun(n, ...)
##   y2 <- beta0 + beta2 * eff2_vec + err_fun(n, ...)
##   y12 <- beta0 + (beta1 + beta2 + beta12) * int_eff_vec + err_fun(n, ...)
##   
##   x <- matrix(0, nrow = 4 * n, ncol = 2)
##   x[(n + 1):(2 * n), 1] <- 1
##   x[(2 * n + 1):(3 * n), 2] <- 1
##   x[(3 * n + 1):(4 * n), 1:2] <- 1
##   y <- c(y0, y1, y2, y12)
##   
##   out <- list(
##     x = x,
##     y = y,
##     eff1_vec = eff1_vec,
##     eff2_vec = eff2_vec,
##     int_eff_vec = int_eff_vec
##   )
## }

Input Parameters

## $n
## [1] 500
## 
## $beta0
## [1] 0
## 
## $beta1
## [1] -1
## 
## $beta2
## [1] -1
## 
## $beta12
## [1] 0
## 
## $eff1
## [1] 1
## 
## $eff2
## [1] 1
## 
## $int_eff
## [1] 1
## 
## $err_fun
## function (n, mean = 0, sd = 1) 
## .Call(C_rnorm, n, mean, sd)
## <bytecode: 0x7feab664cf30>
## <environment: namespace:stats>
## 
## $sd
## [1] 1

Methods and Evaluation

Methods

Quantile Regression (tau = 0.5)

In Wang et al. (2023), epistasis was assessed in the experimental data by fitting a quantile regression model and testing whether the interaction effect term was significantly different from 0. In this simulation study, we thus follow suit. Given the simulated data, we fit a quantile regression model:

\[y_i = \beta_0 + \beta_1 x_{i1} + \beta_2 x_{i2} + \beta_{12} x_{i1} x_{i2} + \epsilon_i,\]

where \(y_i\) is the diameter (\(\mu m\)) of cell \(i\), \(x_{ij}\) is an indicator of whether cell \(i\) belonged to an experimental batch where gene \(j\) was silenced, and \(\epsilon_i\) is the random error or noise term for cell \(i\). Using this regression model, we performed a t-test to assess the null hypothesis that \(\beta_{12} = 0\) (i.e., no epistatic effect) versus the alternative hypothesis that \(\beta_{12} \neq 0\) (i.e., a non-zero epistatic effect) at the 0.5 quantile level (i.e., median).

We note that in Wang et al. (2023), percentile bootstrap hypothesis tests were conducted in addition to the t-tests, but as they gave similar conclusions in Wang et al. (2023), we omit the bootstrap hypothesis tests here to reduce the computational burden.

Function

## function(x, y, tau = 0.5, min_thr = 1e-16, ...) {
##   fit_df <- data.frame(.y = y, x)
##   fit <- quantreg::rq(formula = .y ~ X1 * X2, tau = tau, data = fit_df)
##   fit_summary <- summary(fit, se = "boot")$coefficients %>%
##     as.data.frame() %>%
##     tibble::rownames_to_column("Term") %>%
##     dplyr::rename(Coefficient = Value) %>%
##     dplyr::filter(Term == "X1:X2") %>%
##     dplyr::mutate(`Pr(>|t|)` = max(min_thr, `Pr(>|t|)`))
##   return(fit_summary)
##   # return(broom::tidy(fit_summary))
## }

Input Parameters

## list()

Evaluation

P-values Summary

N/A

Function

## function (fit_results, vary_params = NULL, nested_cols = NULL, 
##     feature_col, imp_col, group_cols = NULL, na_rm = FALSE, summary_funs = c("mean", 
##         "median", "min", "max", "sd", "raw"), custom_summary_funs = NULL, 
##     eval_id = "feature_importance") 
## {
##     group_vars <- c(".dgp_name", ".method_name", vary_params, 
##         group_cols, feature_col)
##     eval_tib <- eval_feature_importance(fit_results = fit_results, 
##         vary_params = vary_params, nested_cols = nested_cols, 
##         feature_col = feature_col, imp_col = imp_col, group_cols = group_cols) %>% 
##         dplyr::group_by(dplyr::across(tidyselect::all_of(group_vars)))
##     eval_summary <- eval_summarizer(eval_data = eval_tib, eval_id = eval_id, 
##         value_col = imp_col, summary_funs = summary_funs, custom_summary_funs = custom_summary_funs, 
##         na_rm = na_rm)
##     return(eval_summary)
## }
## <bytecode: 0x7feacc42f890>
## <environment: namespace:simChef>

Input Parameters

## $eval_id
## [1] "pval"
## 
## $nested_cols
## NULL
## 
## $feature_col
## [1] "Term"
## 
## $imp_col
## [1] "Pr(>|t|)"

Visualizations

P-values Plot

In this visualization, we plot the distribution of p-values corresponding to the interaction effect term (i.e., \(\beta_{12}\)) across 100 simulation replicates. The red line indicates the \(\alpha = 0.05\) significance level. The y-axis is on the log10 scale. When \(\beta_{12} = 0\), we would hope that the p-values are greater than the significance threshold.

CCDC141-IGF1R, Diseased Simulation - CCDC141-IGF1R, Diseased - Varying beta12

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab36244e0>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-IGF1R, Diseased Simulation - CCDC141-IGF1R, Diseased - Varying sd

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feace7e3600>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-IGF1R, Healthy Simulation - CCDC141-IGF1R, Healthy - Varying beta12

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab254bca0>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-IGF1R, Healthy Simulation - CCDC141-IGF1R, Healthy - Varying sd

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab4433fb8>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-TTN, Diseased Simulation - CCDC141-TTN, Diseased - Varying beta12

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab517f108>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-TTN, Diseased Simulation - CCDC141-TTN, Diseased - Varying sd

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab45e2240>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-TTN, Healthy Simulation - CCDC141-TTN, Healthy - Varying beta12

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab5405090>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
CCDC141-TTN, Healthy Simulation - CCDC141-TTN, Healthy - Varying sd

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab278b738>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
Efficiency = 1 Simulation - Efficiency = 1 - Varying beta12

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab5a03530>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>
Efficiency = 1 Simulation - Efficiency = 1 - Varying sd

Function

## function(fit_results = NULL, eval_results, vary_params = NULL,
##                          eval_name, eval_id, add_ggplot_layers = NULL) {
##   vary_params <- unique(vary_params)
##   plt_df <- eval_results[[eval_name]] %>%
##     dplyr::filter(!is.na(!!rlang::sym(vary_params)))
##   plt <- plt_df %>%
##     tidyr::unnest(cols = tidyselect::all_of(sprintf("raw_%s", eval_id))) %>%
##     vdocs::plot_boxplot(
##       x_str = vary_params, y_str = sprintf("raw_%s", eval_id)
##     ) +
##     ggplot2::labs(y = "P-value")
##   if (vary_params == "sd") {
##     plt <- plt + ggplot2::labs(x = bquote("SD of Noise" ~ (sigma)))
##   } else if (vary_params == "beta12") {
##     plt <- plt + ggplot2::labs(x = bquote("Interaction Effect" ~ (beta[12])))
##   }
##   n_dgps <- length(unique(plt_df$.dgp_name))
##   if (n_dgps > 1) {
##     plt <- plt + ggplot2::facet_wrap(~ .dgp_name)
##   }
##   if (!is.null(add_ggplot_layers)) {
##     for (ggplot_layer in add_ggplot_layers) {
##       plt <- plt + ggplot_layer
##     }
##   }
##   return(plt)
## }
## <bytecode: 0x7feab5a41ce0>

Input Parameters

## $eval_name
## [1] "P-values Summary"
## 
## $eval_id
## [1] "pval"
## 
## $add_ggplot_layers
## $add_ggplot_layers[[1]]
## <ScaleContinuousPosition>
##  Range:  
##  Limits:    0 --    1
## 
## $add_ggplot_layers[[2]]
## mapping: yintercept = ~yintercept 
## geom_hline: na.rm = FALSE
## stat_identity: na.rm = FALSE
## position_identity 
## 
## $add_ggplot_layers[[3]]
## <ggproto object: Class CoordCartesian, Coord, gg>
##     aspect: function
##     backtransform_range: function
##     clip: on
##     default: FALSE
##     distance: function
##     expand: TRUE
##     is_free: function
##     is_linear: function
##     labels: function
##     limits: list
##     modify_scales: function
##     range: function
##     render_axis_h: function
##     render_axis_v: function
##     render_bg: function
##     render_fg: function
##     setup_data: function
##     setup_layout: function
##     setup_panel_guides: function
##     setup_panel_params: function
##     setup_params: function
##     train_panel_guides: function
##     transform: function
##     super:  <ggproto object: Class CoordCartesian, Coord, gg>

CCDC141-IGF1R, Diseased Simulation

CCDC141-IGF1R, Diseased

Varying beta12

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-IGF1R, Diseased`
## $dgp$`CCDC141-IGF1R, Diseased`$beta12
## [1]  0.0 -0.1 -0.2 -0.5 -1.0
## 
## 
## 
## $method
## list()

Varying sd

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-IGF1R, Diseased`
## $dgp$`CCDC141-IGF1R, Diseased`$sd
## [1] 0.1 0.2 0.5 1.0 1.5 2.0
## 
## 
## 
## $method
## list()

CCDC141-IGF1R, Healthy Simulation

CCDC141-IGF1R, Healthy

Varying beta12

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-IGF1R, Healthy`
## $dgp$`CCDC141-IGF1R, Healthy`$beta12
## [1]  0.0 -0.1 -0.2 -0.5 -1.0
## 
## 
## 
## $method
## list()

Varying sd

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-IGF1R, Healthy`
## $dgp$`CCDC141-IGF1R, Healthy`$sd
## [1] 0.1 0.2 0.5 1.0 1.5 2.0
## 
## 
## 
## $method
## list()

CCDC141-TTN, Diseased Simulation

CCDC141-TTN, Diseased

Varying beta12

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-TTN, Diseased`
## $dgp$`CCDC141-TTN, Diseased`$beta12
## [1]  0.0 -0.1 -0.2 -0.5 -1.0
## 
## 
## 
## $method
## list()

Varying sd

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-TTN, Diseased`
## $dgp$`CCDC141-TTN, Diseased`$sd
## [1] 0.1 0.2 0.5 1.0 1.5 2.0
## 
## 
## 
## $method
## list()

CCDC141-TTN, Healthy Simulation

CCDC141-TTN, Healthy

Varying beta12

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-TTN, Healthy`
## $dgp$`CCDC141-TTN, Healthy`$beta12
## [1]  0.0 -0.1 -0.2 -0.5 -1.0
## 
## 
## 
## $method
## list()

Varying sd

P-values Plot

Parameter Values

## $dgp
## $dgp$`CCDC141-TTN, Healthy`
## $dgp$`CCDC141-TTN, Healthy`$sd
## [1] 0.1 0.2 0.5 1.0 1.5 2.0
## 
## 
## 
## $method
## list()

Efficiency = 1 Simulation

Efficiency = 1

Varying beta12

P-values Plot

Parameter Values

## $dgp
## $dgp$`Efficiency = 1`
## $dgp$`Efficiency = 1`$beta12
## [1]  0.0 -0.1 -0.2 -0.5 -1.0
## 
## 
## 
## $method
## list()

Varying sd

P-values Plot

Parameter Values

## $dgp
## $dgp$`Efficiency = 1`
## $dgp$`Efficiency = 1`$sd
## [1] 0.1 0.2 0.5 1.0 1.5 2.0
## 
## 
## 
## $method
## list()
LS0tCnRpdGxlOiAiYHIgcGFyYW1zJHNpbV9uYW1lYCIKYXV0aG9yOiAiYHIgcGFyYW1zJGF1dGhvcmAiCmRhdGU6ICJgciBmb3JtYXQoU3lzLnRpbWUoKSwgJyVCICVkLCAlWScpYCIKb3V0cHV0OiBybWFya2Rvd246Omh0bWxfZG9jdW1lbnQKcGFyYW1zOgogIGF1dGhvcjogCiAgICBsYWJlbDogIkF1dGhvcjoiCiAgICB2YWx1ZTogIiIKICBzaW1fbmFtZToKICAgIGxhYmVsOiAiU2ltdWxhdGlvbiBFeHBlcmltZW50IE5hbWU6IgogICAgdmFsdWU6ICIiCiAgc2ltX3BhdGg6CiAgICBsYWJlbDogIlBhdGggdG8gU2ltdWxhdGlvbiBFeHBlcmltZW50IEZvbGRlcjoiCiAgICB2YWx1ZTogIiIKICB3cml0ZV9maWxlbmFtZToKICAgIGxhYmVsOiAiT3V0cHV0IEZpbGU6IgogICAgdmFsdWU6ICIiCiAgc2hvd19jb2RlOgogICAgbGFiZWw6ICJTaG93IENvZGU6IgogICAgdmFsdWU6IFRSVUUKICBzaG93X2V2YWw6CiAgICBsYWJlbDogIlNob3cgRXZhbHVhdG9yczoiCiAgICB2YWx1ZTogVFJVRQogIHNob3dfdml6OgogICAgbGFiZWw6ICJTaG93IFZpc3VhbGl6ZXJzOiIKICAgIHZhbHVlOiBUUlVFCiAgZXZhbF9vcmRlcjoKICAgIGxhYmVsOiAiT3JkZXIgb2YgRXZhbHVhdG9yczoiCiAgICB2YWx1ZTogTlVMTAogIHZpel9vcmRlcjoKICAgIGxhYmVsOiAiT3JkZXIgb2YgVmlzdWFsaXplcnM6IgogICAgdmFsdWU6IE5VTEwKICB1c2VfaWNvbnM6CiAgICBsYWJlbDogIlVzZSBJY29uczoiCiAgICB2YWx1ZTogVFJVRQogIHVzZV92bW9kZXJuOgogICAgbGFiZWw6ICJVc2UgdnRoZW1lczo6dm1vZGVybjoiCiAgICB2YWx1ZTogVFJVRQogIHdyaXRlOgogICAgbGFiZWw6ICJXcml0ZSBGaWxlOiIKICAgIHZhbHVlOiBGQUxTRQogIHZlcmJvc2U6CiAgICBsYWJlbDogIlZlcmJvc2UgTGV2ZWw6IgogICAgdmFsdWU6IDIKLS0tCgpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KAogIGVjaG8gPSBGQUxTRSwKICB3YXJuaW5nID0gRkFMU0UsCiAgbWVzc2FnZSA9IEZBTFNFLAogIGNhY2hlID0gRkFMU0UsCiAgZmlnLmFsaWduID0gImNlbnRlciIsCiAgZmlnLnBvcyA9ICJIIiwKICBmaWcuaGVpZ2h0ID0gMTIsCiAgZmlnLndpZHRoID0gMTAKKQoKb3B0aW9ucygKICB3aWR0aCA9IDEwMDAwLAogIGtuaXRyLmthYmxlLk5BID0gIk5BIgopCgojIHNjcm9sbGFibGUgdGV4dCBvdXRwdXQKbG9jYWwoewogIGhvb2tfb3V0cHV0IDwtIGtuaXRyOjprbml0X2hvb2tzJGdldCgib3V0cHV0IikKICBrbml0cjo6a25pdF9ob29rcyRzZXQob3V0cHV0ID0gZnVuY3Rpb24oeCwgb3B0aW9ucykgewogICAgaWYgKCFpcy5udWxsKG9wdGlvbnMkbWF4LmhlaWdodCkpIHsKICAgICAgb3B0aW9ucyRhdHRyLm91dHB1dCA8LSBjKAogICAgICAgIG9wdGlvbnMkYXR0ci5vdXRwdXQsCiAgICAgICAgc3ByaW50Zignc3R5bGU9Im1heC1oZWlnaHQ6ICVzOyInLCBvcHRpb25zJG1heC5oZWlnaHQpCiAgICAgICkKICAgIH0KICAgIGhvb2tfb3V0cHV0KHgsIG9wdGlvbnMpCiAgfSkKfSkKCmNodW5rX2lkeCA8LSAxCmRvY19kaXIgPC0gZmlsZS5wYXRoKHBhcmFtcyRzaW1fcGF0aCwgImRvY3MiKQp3cml0ZV9maWxlbmFtZSA8LSBwYXJhbXMkd3JpdGVfZmlsZW5hbWUKYGBgCgpgYGB7ciBoZWxwZXItZnVuc30KCiMnIFdyYXAgdGV4dC9jb2RlIGluIGtuaXRyIGNvZGUgY2h1bmsgc3RyaW5nCiMnCiMnIEBwYXJhbSBjb2RlIFN0cmluZyBvZiBjb2RlIHRvIHdyYXAgaW4ga25pdHIgY29kZSBjaHVuawojJyBAcGFyYW0gY2h1bmtfYXJncyBTdHJpbmcgb2YgYXJndW1lbnRzIHRvIHBsYWNlIGluIHRoZSBrbml0ciBjb2RlIGNodW5rIGhlYWRlcgojJyBAcmV0dXJuIFN0cmluZyBvZiBjb2RlLCB3cmFwcGVkIGluc2lkZSB0aGUga25pdHIgY29kZSBjaHVuayBgYGAgbWFya2Vycwp3cml0ZV9jb2RlX2NodW5rIDwtIGZ1bmN0aW9uKGNvZGUgPSAiIiwgY2h1bmtfYXJncyA9ICIiKSB7CiAgc3ByaW50ZigiXG5gYGB7ciwgJXN9XG4lc1xuYGBgXG4iLCBjaHVua19hcmdzLCBjb2RlKQp9CgojJyBXcml0ZSB0ZXh0IHRvIHZlY3RvciAod3JpdGVfZmxhZyA9IFRSVUUpIG9yIHRvIGNvbnNvbGUgKHdyaXRlX2ZsYWcgPSBGQUxTRSkKIycKIycgQHBhcmFtIC4uLiBUZXh0IHRvIHdyaXRlIHRvIHZlY3RvciBvciB0byBjb25zb2xlCiMnIEBwYXJhbSBvbGRfdGV4dCBQcmV2aW91cyB0ZXh0IHRvIGFwcGVuZCB0byB3aGVuIHdyaXRpbmcgdG8gYSB2ZWN0b3IKIycgQHBhcmFtIHdyaXRlX2ZsYWcgQm9vbGVhbiBpbmRpY2F0aW5nIHdoZXRoZXIgdG8gd3JpdGUgdGV4dCB0byBhIHZlY3RvcgojJyAgICh3cml0ZV9mbGFnID0gVFJVRSkgb3IgdG8gY29uc29sZSAod3JpdGVfZmxhZyA9IEZBTFNFKQojJyBAcmV0dXJuIElmIHdyaXRlX2ZsYWcgPSBUUlVFLCByZXR1cm5zIHZlY3RvciBvZiB0ZXh0LiBPdGhlcndpc2UsIHRleHQgaXMKIycgICB3cml0dGVuIHRvIGNvbnNvbGUgdmlhIGBjYXQoKWAuCndyaXRlIDwtIGZ1bmN0aW9uKC4uLiwgb2xkX3RleHQgPSBOVUxMLCB3cml0ZV9mbGFnKSB7CiAgaWYgKHdyaXRlX2ZsYWcpIHsKICAgIHJldHVybihjKG9sZF90ZXh0LCAuLi4pKQogIH0gZWxzZSB7CiAgICBkb3RzX2xpc3QgPC0gbGlzdCguLi4pICU+JQogICAgICBwdXJycjo6bWFwKAogICAgICAgIGZ1bmN0aW9uKHgpIHsKICAgICAgICAgIGlmIChzdHJpbmdyOjpzdHJfZGV0ZWN0KHgsICJgciAuKmAiKSkgewogICAgICAgICAgICAjIHJ1biByIGNvZGUgYmVmb3JlIHByaW50aW5nIHJlc3VsdHMgaW4gY2F0KCkKICAgICAgICAgICAgb3V0IDwtIHN0cmluZ3I6OnN0cl9yZXBsYWNlKAogICAgICAgICAgICAgIHgsICJgciAuKmAiLAogICAgICAgICAgICAgIGV2YWwocGFyc2UodGV4dCA9IHN0cmluZ3I6OnN0cl9leHRyYWN0KHgsICIoPzw9YHIgKSguKj8pKD89YCkiKSkpCiAgICAgICAgICAgICkKICAgICAgICAgIH0gZWxzZSB7CiAgICAgICAgICAgIG91dCA8LSB4CiAgICAgICAgICB9CiAgICAgICAgICByZXR1cm4ob3V0KQogICAgICAgIH0KICAgICAgKQogICAgZG8uY2FsbChjYXQsIGFyZ3MgPSBjKGRvdHNfbGlzdCwgbGlzdChzZXAgPSAiIikpKQogIH0KfQoKIycgV3JpdGUgdGV4dCB0byBmaWxlCiMnCiMnIEBwYXJhbSBwYXRoIFBhdGggdG8gb3V0cHV0IGZpbGUKIycgQHBhcmFtIC4uLiBUZXh0IHRvIG91dHB1dCB0byBmaWxlCndyaXRlX3RvX2ZpbGUgPC0gZnVuY3Rpb24ocGF0aCwgLi4uKSB7CiAgc3RvcmVsaW5lcyA8LSByZWFkTGluZXMocGF0aCkKICBzdG9yZWxpbmVzIDwtIGMoc3RvcmVsaW5lcywgLi4uKQogIHdyaXRlTGluZXMoc3RvcmVsaW5lcywgcGF0aCkKfQoKIycgR2V0IG9yZGVyIG9mIG9iamVjdHMgdG8gZGlzcGxheQojJwojJyBAcGFyYW0gb2JqX25hbWVzIFZlY3RvciBvZiBhbGwgb2JqZWN0IG5hbWVzIHRoYXQgbmVlZCB0byBiZSBkaXNwbGF5ZWQuCiMnIEBwYXJhbSBvYmpfb3JkZXIgVmVjdG9yIG9mIG9iamVjdCBuYW1lcyBpbiB0aGUgZGVzaXJlZCBhcHBlYXJhbmNlIG9yZGVyLgojJyBAcmV0dXJuIFZlY3RvciBvZiBvYmplY3QgbmFtZXMgaW4gdGhlIG9yZGVyIGluIHdoaWNoIHRoZXkgd2lsbCBiZSBkaXNwbGF5ZWQuCmdldF9vYmplY3Rfb3JkZXIgPC0gZnVuY3Rpb24ob2JqX25hbWVzLCBvYmpfb3JkZXIgPSBOVUxMKSB7CiAgaWYgKGlzLm51bGwob2JqX29yZGVyKSkgewogICAgcmV0dXJuKG9ial9uYW1lcykKICB9IGVsc2UgewogICAgcmV0dXJuKGludGVyc2VjdChvYmpfb3JkZXIsIG9ial9uYW1lcykpCiAgfQp9CgojJyBHZXQgYWxsIGV4cGVyaW1lbnRzIHVuZGVyIGEgZ2l2ZW4gZGlyZWN0b3J5IG5hbWUKIycKIycgQHBhcmFtIGRpcl9uYW1lIG5hbWUgb2YgZGlyZWN0b3J5CiMnIEByZXR1cm4gbGlzdCBvZiBuYW1lZCBleHBlcmltZW50cwpnZXRfZGVzY2VuZGFudHMgPC0gZnVuY3Rpb24oZGlyX25hbWUpIHsKICBleHBlcmltZW50cyA8LSBsaXN0KCkKICBmb3IgKGQgaW4gbGlzdC5kaXJzKGRpcl9uYW1lKSkgewogICAgaWYgKGZpbGUuZXhpc3RzKGZpbGUucGF0aChkLCAiZXhwZXJpbWVudC5yZHMiKSkpIHsKICAgICAgaWYgKGlkZW50aWNhbChkLCBwYXJhbXMkc2ltX3BhdGgpKSB7CiAgICAgICAgZXhwX25hbWUgPC0gIkJhc2UiCiAgICAgIH0gZWxzZSB7CiAgICAgICAgZXhwX25hbWUgPC0gc3RyaW5ncjo6c3RyX3JlcGxhY2VfYWxsKAogICAgICAgICAgc3RyaW5ncjo6c3RyX3JlbW92ZShkLCBwYXN0ZTAocGFyYW1zJHNpbV9wYXRoLCAiLyIpKSwKICAgICAgICAgICIvIiwgIiAtICIKICAgICAgICApCiAgICAgIH0KICAgICAgZXhwZXJpbWVudHNbW2V4cF9uYW1lXV0gPC0gcmVhZFJEUyhmaWxlLnBhdGgoZCwgImV4cGVyaW1lbnQucmRzIikpCiAgICB9CiAgfQogIHJldHVybihleHBlcmltZW50cykKfQoKIycgQ2hlY2sgaWYgZXhwZXJpbWVudCBleGlzdHMKIycKIycgQHBhcmFtIGRpcl9uYW1lIG5hbWUgb2YgZGlyZWN0b3J5IG9yIHZlY3RvciB0aGVyZW9mCiMnIEBwYXJhbSByZWN1cnNpdmUgbG9naWNhbDsgaWYgVFJVRSwgY2hlY2tzIGlmIGV4cGVyaW1lbnQgZXhpc3RzIHVuZGVyIHRoZQojJyAgIGdpdmVuIGRpcmVjdG9yeShzKTsgaWYgRkFMU0UsIGNoZWNrcyBpZiBhbnkgZXhwZXJpbWVudCBleGlzdHMgdW5kZXIgdGhlCiMnICAgZGlyZWN0b3J5KHMpIGFuZCBpdHMgZGVzY2VuZGFudHMKIycgQHJldHVybiBUUlVFIGlmIGV4cGVyaW1lbnQgZXhpc3RzIGFuZCBGQUxTRSBvdGhlcndpc2UKZXhwZXJpbWVudF9leGlzdHMgPC0gZnVuY3Rpb24oZGlyX25hbWUsIHJlY3Vyc2l2ZSA9IEZBTFNFKSB7CiAgcmVzIDwtIHB1cnJyOjptYXBfbGdsKAogICAgZGlyX25hbWUsCiAgICBmdW5jdGlvbihkKSB7CiAgICAgIGlmICghcmVjdXJzaXZlKSB7CiAgICAgICAgZXhwX2ZuYW1lIDwtIGZpbGUucGF0aChkLCAiZXhwZXJpbWVudC5yZHMiKQogICAgICAgIHJldHVybihmaWxlLmV4aXN0cyhleHBfZm5hbWUpKQogICAgICB9IGVsc2UgewogICAgICAgIGRlc2NlbmRhbnRzIDwtIGdldF9kZXNjZW5kYW50cyhkKQogICAgICAgIHJldHVybihsZW5ndGgoZGVzY2VuZGFudHMpID4gMCkKICAgICAgfQogICAgfQogICkKICByZXR1cm4oYW55KHJlcykpCn0KCiMnIERpc3BsYXlzIGNvbnRlbnQgZm9yIHNwZWNpZmllZCBwYXJ0IG9mIHJlY2lwZQojJwojJyBAcGFyYW0gZmllbGRfbmFtZSBwYXJ0IG9mIHJlY2lwZSB0byBzaG93OyBtdXN0IGJlIG9uZSBvZiAiZGdwIiwgIm1ldGhvZCIsCiMnICAgImV2YWx1YXRvciIsIG9yICJ2aXN1YWxpemVyIgojJyBAcGFyYW0gd3JpdGVfZmxhZyBCb29sZWFuIGluZGljYXRpbmcgd2hldGhlciB0byB3cml0ZSB0ZXh0IHRvIGEgdmVjdG9yCiMnICAgKHdyaXRlX2ZsYWcgPSBUUlVFKSBvciB0byBjb25zb2xlICh3cml0ZV9mbGFnID0gRkFMU0UpCiMnIEByZXR1cm4gY29udGVudCBmb3IgcmVjaXBlCnNob3dfcmVjaXBlIDwtIGZ1bmN0aW9uKGZpZWxkX25hbWUgPSBjKAogICAgICAgICAgICAgICAgICAgICAgICAgICJkZ3AiLCAibWV0aG9kIiwgImV2YWx1YXRvciIsICJ2aXN1YWxpemVyIgogICAgICAgICAgICAgICAgICAgICAgICApLAogICAgICAgICAgICAgICAgICAgICAgICB3cml0ZV9mbGFnID0gRkFMU0UpIHsKICBmaWVsZF9uYW1lIDwtIG1hdGNoLmFyZyhmaWVsZF9uYW1lKQogIGZ1bmNfbmFtZSA8LSBkcGx5cjo6Y2FzZV93aGVuKAogICAgZmllbGRfbmFtZSA9PSAiZXZhbHVhdG9yIiB+ICJldmFsIiwKICAgIGZpZWxkX25hbWUgPT0gInZpc3VhbGl6ZXIiIH4gInZpeiIsCiAgICBUUlVFIH4gZmllbGRfbmFtZQogICkKICBkZXNjZW5kYW50cyA8LSBnZXRfZGVzY2VuZGFudHMoZGlyX25hbWUgPSBwYXJhbXMkc2ltX3BhdGgpCiAgb2JqcyA8LSBwdXJycjo6bWFwKGRlc2NlbmRhbnRzLCB+IC54W1twYXN0ZTAoImdldF8iLCBmaWVsZF9uYW1lLCAicyIpXV0oKSkKICBvYmpfbmFtZXMgPC0gdW5pcXVlKHB1cnJyOjpyZWR1Y2Uoc2FwcGx5KG9ianMsIG5hbWVzKSwgYykpCgogIGlmIChmaWVsZF9uYW1lICVpbiUgYygibWV0aG9kIiwgImV2YWx1YXRvciIpKSB7CiAgICBvYmpfaGVhZGVyIDwtICJcblxuIyMjIyAlcyB7LnRhYnNldCAudGFic2V0LXBpbGxzIC50YWJzZXQtY2lyY2xlIC50YWJzZXQtcmVjaXBlfVxuXG4iCiAgICBzaG93dHlwZV9oZWFkZXIgPC0gIlxuXG4jIyMjIyAlcyB7LnRhYnNldCAudGFic2V0LXBpbGxzfVxuXG4iCiAgICBleHBfaGVhZGVyIDwtICJcblxuIyMjIyMjICVzIFxuXG4iCiAgfSBlbHNlIHsKICAgIG9ial9oZWFkZXIgPC0gIlxuXG4jIyMgJXMgey50YWJzZXQgLnRhYnNldC1waWxscyAudGFic2V0LWNpcmNsZSAudGFic2V0LXJlY2lwZX1cblxuIgogICAgc2hvd3R5cGVfaGVhZGVyIDwtICJcblxuIyMjIyAlcyB7LnRhYnNldCAudGFic2V0LXBpbGxzfVxuXG4iCiAgICBleHBfaGVhZGVyIDwtICJcblxuIyMjIyMgJXMgXG5cbiIKICB9CgogIGlmIChwYXJhbXMkdXNlX2ljb25zKSB7CiAgICBpZiAocGFyYW1zJHVzZV92bW9kZXJuKSB7CiAgICAgIGRlc2NyaXB0aW9uX2xhYmVsIDwtICJgciBmb250YXdlc29tZTo6ZmEoJ3JlYWRtZScsIGZpbGwgPSAnd2hpdGUnKWAiCiAgICAgIGNvZGVfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgnY29kZScsIGZpbGwgPSAnd2hpdGUnKWAiCiAgICB9IGVsc2UgewogICAgICBkZXNjcmlwdGlvbl9sYWJlbCA8LSAiYHIgZm9udGF3ZXNvbWU6OmZhKCdyZWFkbWUnKWAiCiAgICAgIGNvZGVfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgnY29kZScpYCIKICAgIH0KICB9IGVsc2UgewogICAgZGVzY3JpcHRpb25fbGFiZWwgPC0gIkRlc2NyaXB0aW9uIgogICAgY29kZV9sYWJlbCA8LSAiQ29kZSIKICB9CgogIGlmIChhbGwoc2FwcGx5KG9ianMsIGxlbmd0aCkgPT0gMCkpIHsKICAgIGlmICh3cml0ZV9mbGFnKSB7CiAgICAgIHJldHVybigiTi9BIikKICAgIH0gZWxzZSB7CiAgICAgIHJldHVybihjYXQoIk4vQSIpKQogICAgfQogIH0KCiAgcmVjaXBlIDwtIGMoKQogIGZvciAoaWR4IGluIDE6bGVuZ3RoKG9ial9uYW1lcykpIHsKICAgIG9ial9uYW1lIDwtIG9ial9uYW1lc1tpZHhdCiAgICBkZXNjcmlwdGlvbl9mcGF0aCA8LSBmaWxlLnBhdGgoCiAgICAgIGRvY19kaXIsIHBhc3RlMChmaWVsZF9uYW1lLCAicyIpLCBwYXN0ZTAob2JqX25hbWUsICIubWQiKQogICAgKQogICAgCiAgICBpZiAocGFyYW1zJHVzZV92bW9kZXJuKSB7CiAgICAgIHJlY2lwZSA8LSB3cml0ZSgKICAgICAgICAiXG5cbjxkaXYgY2xhc3M9J3BhbmVsIHBhbmVsLWRlZmF1bHQgcGFkZGVkLXBhbmVsJz5cblxuIiwKICAgICAgICBvbGRfdGV4dCA9IHJlY2lwZSwgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICAgKQogICAgfQoKICAgIHJlY2lwZSA8LSB3cml0ZSgKICAgICAgc3ByaW50ZihvYmpfaGVhZGVyLCBvYmpfbmFtZSksCiAgICAgIHNwcmludGYoc2hvd3R5cGVfaGVhZGVyLCBkZXNjcmlwdGlvbl9sYWJlbCksCiAgICAgIHBhc3RlTWQoZGVzY3JpcHRpb25fZnBhdGgpLAogICAgICBvbGRfdGV4dCA9IHJlY2lwZSwgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICkKCiAgICBpZiAocGFyYW1zJHNob3dfY29kZSkgewogICAgICByZWNpcGUgPC0gd3JpdGUoCiAgICAgICAgc3ByaW50ZihzaG93dHlwZV9oZWFkZXIsIGNvZGVfbGFiZWwpLAogICAgICAgIG9sZF90ZXh0ID0gcmVjaXBlLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgICApCgogICAgICBrZWVwX29ianMgPC0gcHVycnI6OmNvbXBhY3QocHVycnI6Om1hcChvYmpzLCBvYmpfbmFtZSkpCiAgICAgIGlzX2lkZW50aWNhbCA8LSBhbGwoCiAgICAgICAgcHVycnI6Om1hcF9sZ2woa2VlcF9vYmpzLCB+IGlzVFJVRShjaGVja19lcXVhbCgueCwga2VlcF9vYmpzW1sxXV0pKSkKICAgICAgKQogICAgICBmb3IgKGV4cCBpbiBuYW1lcyhrZWVwX29ianMpKSB7CiAgICAgICAgb2JqIDwtIGtlZXBfb2Jqc1tbZXhwXV0KICAgICAgICBpZiAoIWlzX2lkZW50aWNhbCkgewogICAgICAgICAgcmVjaXBlIDwtIHdyaXRlKAogICAgICAgICAgICBzcHJpbnRmKGV4cF9oZWFkZXIsIGV4cCksCiAgICAgICAgICAgIG9sZF90ZXh0ID0gcmVjaXBlLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgICAgICAgKQogICAgICAgIH0KCiAgICAgICAgcmVjaXBlIDwtIHdyaXRlKAogICAgICAgICAgIlxuXG4qKkZ1bmN0aW9uKipcblxuIiwKICAgICAgICAgIG9sZF90ZXh0ID0gcmVjaXBlLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgICAgICkKICAgICAgICBpZiAod3JpdGVfZmxhZykgewogICAgICAgICAgcmVjaXBlIDwtIHNwcmludGYoCiAgICAgICAgICAgICJzaG93X3JlY2lwZSglc19vYmpzLCAnJXMnLCAnJXMnLCB3aGF0ID0gJ2Z1bmN0aW9uJykiLAogICAgICAgICAgICBmdW5jX25hbWUsIG9ial9uYW1lLCBleHAKICAgICAgICAgICkgJT4lCiAgICAgICAgICAgIHdyaXRlX2NvZGVfY2h1bmsoY2h1bmtfYXJncyA9ICJtYXguaGVpZ2h0PScyMDBweCciKSAlPiUKICAgICAgICAgICAgd3JpdGUob2xkX3RleHQgPSByZWNpcGUsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnKQogICAgICAgIH0gZWxzZSB7CiAgICAgICAgICB2dGhlbWVzOjpzdWJjaHVua2lmeSgKICAgICAgICAgICAgb2JqW1twYXN0ZTAoZnVuY19uYW1lLCAiX2Z1biIpXV0sIGNodW5rX2lkeCwKICAgICAgICAgICAgb3RoZXJfYXJncyA9ICJtYXguaGVpZ2h0PScyMDBweCciCiAgICAgICAgICApCiAgICAgICAgICBjaHVua19pZHggPDwtIGNodW5rX2lkeCArIDEKICAgICAgICB9CgogICAgICAgIHJlY2lwZSA8LSB3cml0ZSgKICAgICAgICAgICJcblxuKipJbnB1dCBQYXJhbWV0ZXJzKipcblxuIiwKICAgICAgICAgIG9sZF90ZXh0ID0gcmVjaXBlLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgICAgICkKICAgICAgICBpZiAod3JpdGVfZmxhZykgewogICAgICAgICAgcmVjaXBlIDwtIHNwcmludGYoCiAgICAgICAgICAgICJzaG93X3JlY2lwZSglc19vYmpzLCAnJXMnLCAnJXMnLCB3aGF0ID0gJ3BhcmFtZXRlcnMnKSIsCiAgICAgICAgICAgIGZ1bmNfbmFtZSwgb2JqX25hbWUsIGV4cAogICAgICAgICAgKSAlPiUKICAgICAgICAgICAgd3JpdGVfY29kZV9jaHVuayhjaHVua19hcmdzID0gIm1heC5oZWlnaHQ9JzIwMHB4JyIpICU+JQogICAgICAgICAgICB3cml0ZShvbGRfdGV4dCA9IHJlY2lwZSwgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcpCiAgICAgICAgfSBlbHNlIHsKICAgICAgICAgIHZ0aGVtZXM6OnN1YmNodW5raWZ5KAogICAgICAgICAgICBvYmpbW3Bhc3RlMChmdW5jX25hbWUsICJfcGFyYW1zIildXSwgY2h1bmtfaWR4LAogICAgICAgICAgICBvdGhlcl9hcmdzID0gIm1heC5oZWlnaHQ9JzIwMHB4JyIKICAgICAgICAgICkKICAgICAgICAgIGNodW5rX2lkeCA8PC0gY2h1bmtfaWR4ICsgMQogICAgICAgIH0KCiAgICAgICAgaWYgKGlzX2lkZW50aWNhbCkgewogICAgICAgICAgYnJlYWsKICAgICAgICB9CiAgICAgIH0KICAgIH0KCiAgICBpZiAocGFyYW1zJHVzZV92bW9kZXJuKSB7CiAgICAgIHJlY2lwZSA8LSB3cml0ZSgKICAgICAgICAiXG5cbjwvZGl2PlxuXG4iLCBvbGRfdGV4dCA9IHJlY2lwZSwgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICAgKQogICAgfQogIH0KCiAgcmV0dXJuKHJlY2lwZSkKfQoKIycgUmVhZHMgaW4gZmlsZSBpZiBpdCBleGlzdHMgYW5kIHJldHVybnMgTlVMTCBpZiB0aGUgZmlsZSBkb2VzIG5vdCBleGlzdAojJwojJyBAcGFyYW0gZmlsZW5hbWUgbmFtZSBvZiAucmRzIGZpbGUgdG8gdHJ5IHJlYWRpbmcgaW4KIycgQHJldHVybiBvdXRwdXQgb2YgZmlsZW5hbWUucmRzIGlmIHRoZSBmaWxlIGV4aXN0cyBhbmQgTlVMTCBvdGhlcndpc2UKZ2V0X3Jlc3VsdHMgPC0gZnVuY3Rpb24oZmlsZW5hbWUpIHsKICBpZiAoZmlsZS5leGlzdHMoZmlsZW5hbWUpKSB7CiAgICByZXN1bHRzIDwtIHJlYWRSRFMoZmlsZW5hbWUpCiAgfSBlbHNlIHsKICAgIHJlc3VsdHMgPC0gTlVMTAogIH0KICByZXR1cm4ocmVzdWx0cykKfQoKIycgRGlzcGxheXMgb3V0cHV0IChib3RoIGZyb20gZXZhbHVhdGUoKSBhbmQgdmlzdWFsaXplKCkpIGZyb20gc2F2ZWQgcmVzdWx0cyB1bmRlcgojJyBhIHNwZWNpZmllZCBkaXJlY3RvcnkKIycKIycgQHBhcmFtIGRpcl9uYW1lIG5hbWUgb2YgZGlyZWN0b3J5CiMnIEBwYXJhbSBkZXB0aCBpbnRlZ2VyOyBkZXB0aCBvZiBkaXJlY3RvcnkgZnJvbSBwYXJlbnQvYmFzZSBleHBlcmltZW50J3MgZm9sZGVyCiMnIEBwYXJhbSBiYXNlIGxvZ2ljYWw7IHdoZXRoZXIgb3Igbm90IHRoaXMgaXMgYSBiYXNlIGV4cGVyaW1lbnQKIycgQHBhcmFtIHNob3dfaGVhZGVyIGxvZ2ljYWw7IHdoZXRoZXIgb3Igbm90IHRvIHNob3cgc2VjdGlvbiBoZWFkZXIKIycgQHBhcmFtIHZlcmJvc2UgaW50ZWdlcjsgMCA9IG5vIG1lc3NhZ2VzOyAxID0gcHJpbnQgb3V0IGRpcmVjdG9yeSBuYW1lIG9ubHk7CiMnICAgMiA9IHByaW50IG91dCBkaXJlY3RvcnkgbmFtZSBhbmQgbmFtZSBvZiBldmFsdWF0b3JzL3Zpc3VhbGl6ZXJzCiMnIEBwYXJhbSB3cml0ZV9mbGFnIEJvb2xlYW4gaW5kaWNhdGluZyB3aGV0aGVyIHRvIHdyaXRlIHRleHQgdG8gYSB2ZWN0b3IKIycgICAod3JpdGVfZmxhZyA9IFRSVUUpIG9yIHRvIGNvbnNvbGUgKHdyaXRlX2ZsYWcgPSBGQUxTRSkKIycgQHJldHVybiBjb250ZW50IHJlc3VsdHMgZnJvbSBldmFsdWF0ZSgpIGFuZCB2aXN1YWxpemUoKSBmcm9tIHRoZSBleHBlcmltZW50CnNob3dfcmVzdWx0cyA8LSBmdW5jdGlvbihkaXJfbmFtZSwgZGVwdGgsIGJhc2UgPSBGQUxTRSwgc2hvd19oZWFkZXIgPSBUUlVFLAogICAgICAgICAgICAgICAgICAgICAgICAgdmVyYm9zZSA9IDEsIHdyaXRlX2ZsYWcgPSBGQUxTRSkgewogIGlmICh2ZXJib3NlID49IDEpIHsKICAgIGluZm9ybShwYXN0ZTAocGFzdGUocmVwKCIqIiwgZGVwdGgpLCBjb2xsYXBzZSA9ICIiKSwgYmFzZW5hbWUoZGlyX25hbWUpKSkKICB9CgogIGlmIChkZXB0aCA9PSAxKSB7CiAgICBoZWFkZXJfdGVtcGxhdGUgPC0gIlxuXG4lcyAlcyB7LnRhYnNldCAudGFic2V0LXBpbGxzIC50YWJzZXQtdm1vZGVybn1cblxuIgogIH0gZWxzZSB7CiAgICBpZiAoYmFzZSB8fCAhZXhwZXJpbWVudF9leGlzdHMoZGlyX25hbWUpKSB7CiAgICAgIGhlYWRlcl90ZW1wbGF0ZSA8LSAiXG5cbiVzICVzIHsudGFic2V0IC50YWJzZXQtcGlsbHN9XG5cbiIKICAgIH0gZWxzZSB7CiAgICAgIGhlYWRlcl90ZW1wbGF0ZSA8LSAiXG5cbiVzICVzIHsudGFic2V0IC50YWJzZXQtcGlsbHMgLnRhYnNldC1jaXJjbGV9XG5cbiIKICAgIH0KICB9CgogIHJlc3VsdHMgPC0gYygpCiAgaWYgKHNob3dfaGVhZGVyKSB7CiAgICByZXN1bHRzIDwtIHNwcmludGYoCiAgICAgIGhlYWRlcl90ZW1wbGF0ZSwKICAgICAgcGFzdGUocmVwKCIjIiwgZGVwdGgpLCBjb2xsYXBzZSA9ICIiKSwKICAgICAgYmFzZW5hbWUoZGlyX25hbWUpCiAgICApICU+JQogICAgICB3cml0ZShvbGRfdGV4dCA9IHJlc3VsdHMsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnKQogIH0KCiAgaWYgKGJhc2UpIHsKICAgIHJlc3VsdHMgPC0gc3ByaW50ZigKICAgICAgIlxuXG4lcyBCYXNlIC0gJXMgey50YWJzZXQgLnRhYnNldC1waWxscyAudGFic2V0LWNpcmNsZX1cblxuIiwKICAgICAgcGFzdGUocmVwKCIjIiwgZGVwdGggKyAxKSwgY29sbGFwc2UgPSAiIiksCiAgICAgIGJhc2VuYW1lKGRpcl9uYW1lKQogICAgKSAlPiUKICAgICAgd3JpdGUob2xkX3RleHQgPSByZXN1bHRzLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZykKICAgIGRlcHRoIDwtIGRlcHRoICsgMQogIH0KCiAgc2hvd3R5cGVfdGVtcGxhdGUgPC0gcGFzdGUwKAogICAgIlxuXG4iLCBwYXN0ZShyZXAoIiMiLCBkZXB0aCArIDEpLCBjb2xsYXBzZSA9ICIiKSwgIiAlc1xuXG4iCiAgKQogIGZpZ25hbWVfdGVtcGxhdGUgPC0gcGFzdGUwKAogICAgIlxuXG4iLCBwYXN0ZShyZXAoIiMiLCBkZXB0aCArIDIpLCBjb2xsYXBzZSA9ICIiKSwgIiAlc1xuXG4iCiAgKQogIGludmlzaWJsZV9oZWFkZXIgPC0gcGFzdGUwKAogICAgIlxuXG4iLCBwYXN0ZShyZXAoIiMiLCBkZXB0aCArIDMpLCBjb2xsYXBzZSA9ICIiKSwKICAgICIgey50YWJzZXQgLnRhYnNldC1waWxsc31cblxuIgogICkKICBwbHRfdGVtcGxhdGUgPC0gcGFzdGUwKAogICAgIlxuXG4iLCBwYXN0ZShyZXAoIiMiLCBkZXB0aCArIDQpLCBjb2xsYXBzZSA9ICIiKSwgIiAlc1xuXG4iCiAgKQoKICBpZiAocGFyYW1zJHVzZV9pY29ucykgewogICAgaWYgKHBhcmFtcyR1c2Vfdm1vZGVybikgewogICAgICBldmFsdWF0b3JfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgndGFibGUnLCBmaWxsID0gJ3doaXRlJylgIgogICAgICB2aXN1YWxpemVyX2xhYmVsIDwtICJgciBmb250YXdlc29tZTo6ZmEoJ2NoYXJ0LWJhcicsIGZpbGwgPSAnd2hpdGUnKWAiCiAgICAgIGNvZGVfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgnY29kZScsIGZpbGwgPSAnd2hpdGUnKWAiCiAgICB9IGVsc2UgewogICAgICBldmFsdWF0b3JfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgndGFibGUnKWAiCiAgICAgIHZpc3VhbGl6ZXJfbGFiZWwgPC0gImByIGZvbnRhd2Vzb21lOjpmYSgnY2hhcnQtYmFyJylgIgogICAgICBjb2RlX2xhYmVsIDwtICJgciBmb250YXdlc29tZTo6ZmEoJ2NvZGUnKWAiCiAgICB9CiAgfSBlbHNlIHsKICAgIGV2YWx1YXRvcl9sYWJlbCA8LSAiRXZhbHVhdG9ycyIKICAgIHZpc3VhbGl6ZXJfbGFiZWwgPC0gIlZpc3VhbGl6ZXJzIgogICAgY29kZV9sYWJlbCA8LSAiVmFyeWluZyBQYXJhbWV0ZXJzIgogIH0KCiAgZXhwX2ZuYW1lIDwtIGZpbGUucGF0aChkaXJfbmFtZSwgImV4cGVyaW1lbnQucmRzIikKICBldmFsX2ZuYW1lIDwtIGZpbGUucGF0aChkaXJfbmFtZSwgImV2YWxfcmVzdWx0cy5yZHMiKQogIHZpel9mbmFtZSA8LSBmaWxlLnBhdGgoZGlyX25hbWUsICJ2aXpfcmVzdWx0cy5yZHMiKQoKICBleHAgPC0gZ2V0X3Jlc3VsdHMoZXhwX2ZuYW1lKQogIGV2YWxfcmVzdWx0cyA8LSBnZXRfcmVzdWx0cyhldmFsX2ZuYW1lKQogIHZpel9yZXN1bHRzIDwtIGdldF9yZXN1bHRzKHZpel9mbmFtZSkKCiAgaWYgKCFpcy5udWxsKGV2YWxfcmVzdWx0cykgJiYgcGFyYW1zJHNob3dfZXZhbCkgewogICAgcmVzdWx0cyA8LSB3cml0ZSgKICAgICAgc3ByaW50ZihzaG93dHlwZV90ZW1wbGF0ZSwgZXZhbHVhdG9yX2xhYmVsKSwKICAgICAgb2xkX3RleHQgPSByZXN1bHRzLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgKQoKICAgIGV2YWxfbmFtZXMgPC0gZ2V0X29iamVjdF9vcmRlcihuYW1lcyhldmFsX3Jlc3VsdHMpLCBwYXJhbXMkZXZhbF9vcmRlcikKICAgIGZvciAoZXZhbF9uYW1lIGluIGV2YWxfbmFtZXMpIHsKICAgICAgZXZhbHVhdG9yIDwtIGV4cCRnZXRfZXZhbHVhdG9ycygpW1tldmFsX25hbWVdXQogICAgICBpZiAoZXZhbHVhdG9yJGRvY19zaG93KSB7CiAgICAgICAgaWYgKHZlcmJvc2UgPj0gMSkgewogICAgICAgICAgaW5mb3JtKHBhc3RlMChwYXN0ZShyZXAoIiAiLCBkZXB0aCArIDEpLCBjb2xsYXBzZSA9ICIiKSwgZXZhbF9uYW1lKSkKICAgICAgICB9CiAgICAgICAgcmVzdWx0cyA8LSB3cml0ZSgKICAgICAgICAgIHNwcmludGYoZmlnbmFtZV90ZW1wbGF0ZSwgZXZhbF9uYW1lKSwKICAgICAgICAgIG9sZF90ZXh0ID0gcmVzdWx0cywgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICAgICApCiAgICAgICAgaWYgKGlzLm51bGwoZXZhbHVhdG9yJGRvY19ucm93cykpIHsKICAgICAgICAgIGV2YWxfcmVzdWx0c19zaG93IDwtIGV2YWxfcmVzdWx0c1tbZXZhbF9uYW1lXV0KICAgICAgICB9IGVsc2UgewogICAgICAgICAga2VlcF9yb3dzIDwtIDE6bWluKGV2YWx1YXRvciRkb2NfbnJvd3MsIG5yb3coZXZhbF9yZXN1bHRzW1tldmFsX25hbWVdXSkpCiAgICAgICAgICBldmFsX3Jlc3VsdHNfc2hvdyA8LSBldmFsX3Jlc3VsdHNbW2V2YWxfbmFtZV1dW2tlZXBfcm93cywgXQogICAgICAgICAgaWYgKG5yb3coZXZhbF9yZXN1bHRzW1tldmFsX25hbWVdXSkgPiBldmFsdWF0b3IkZG9jX25yb3dzKSB7CiAgICAgICAgICAgIG9taXR0ZWRfbnJvd3MgPC0gbnJvdyhldmFsX3Jlc3VsdHNbW2V2YWxfbmFtZV1dKSAtIGV2YWx1YXRvciRkb2NfbnJvd3MKICAgICAgICAgICAgcmVzdWx0cyA8LSB3cml0ZSgKICAgICAgICAgICAgICBzcHJpbnRmKAogICAgICAgICAgICAgICAgIlNob3dpbmcgcHJldmlldyBvZiAlcyByZXN1bHRzLiAlcyByb3dzIGhhdmUgYmVlbiBvbWl0dGVkLlxuXG4iLAogICAgICAgICAgICAgICAgZXZhbF9uYW1lLCBvbWl0dGVkX25yb3dzCiAgICAgICAgICAgICAgKSwKICAgICAgICAgICAgICBvbGRfdGV4dCA9IHJlc3VsdHMsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnCiAgICAgICAgICAgICkKICAgICAgICAgIH0KICAgICAgICB9CiAgICAgICAgaWYgKHdyaXRlX2ZsYWcpIHsKICAgICAgICAgIHJlc3VsdHMgPC0gc3ByaW50ZigKICAgICAgICAgICAgInNob3dfcmVzdWx0cygnJXMnLCAnJXMnLCAnZXZhbHVhdG9yJykiLCBkaXJfbmFtZSwgZXZhbF9uYW1lCiAgICAgICAgICApICU+JQogICAgICAgICAgICB3cml0ZV9jb2RlX2NodW5rKGNodW5rX2FyZ3MgPSAicmVzdWx0cyA9ICdhc2lzJyIpICU+JQogICAgICAgICAgICB3cml0ZShvbGRfdGV4dCA9IHJlc3VsdHMsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnKQogICAgICAgIH0gZWxzZSB7CiAgICAgICAgICBkby5jYWxsKAogICAgICAgICAgICB2dGhlbWVzOjpwcmV0dHlfRFQsCiAgICAgICAgICAgIGMobGlzdChldmFsX3Jlc3VsdHNfc2hvdyksIGV2YWx1YXRvciRkb2Nfb3B0aW9ucykKICAgICAgICAgICkgJT4lCiAgICAgICAgICAgIHZ0aGVtZXM6OnN1YmNodW5raWZ5KGkgPSBjaHVua19pZHgpCiAgICAgICAgICBjaHVua19pZHggPDwtIGNodW5rX2lkeCArIDEKICAgICAgICB9CiAgICAgIH0KICAgIH0KICB9CgogIGlmICghaXMubnVsbCh2aXpfcmVzdWx0cykgJiYgcGFyYW1zJHNob3dfdml6KSB7CiAgICByZXN1bHRzIDwtIHdyaXRlKAogICAgICBzcHJpbnRmKHNob3d0eXBlX3RlbXBsYXRlLCB2aXN1YWxpemVyX2xhYmVsKSwKICAgICAgb2xkX3RleHQgPSByZXN1bHRzLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgKQoKICAgIHZpel9uYW1lcyA8LSBnZXRfb2JqZWN0X29yZGVyKG5hbWVzKHZpel9yZXN1bHRzKSwgcGFyYW1zJHZpel9vcmRlcikKICAgIGZvciAodml6X25hbWUgaW4gdml6X25hbWVzKSB7CiAgICAgIHZpc3VhbGl6ZXIgPC0gZXhwJGdldF92aXN1YWxpemVycygpW1t2aXpfbmFtZV1dCiAgICAgIGlmICh2aXN1YWxpemVyJGRvY19zaG93KSB7CiAgICAgICAgaWYgKHZlcmJvc2UgPj0gMSkgewogICAgICAgICAgaW5mb3JtKHBhc3RlMChwYXN0ZShyZXAoIiAiLCBkZXB0aCArIDEpLCBjb2xsYXBzZSA9ICIiKSwgdml6X25hbWUpKQogICAgICAgIH0KICAgICAgICByZXN1bHRzIDwtIHdyaXRlKAogICAgICAgICAgc3ByaW50ZihmaWduYW1lX3RlbXBsYXRlLCB2aXpfbmFtZSksCiAgICAgICAgICBpbnZpc2libGVfaGVhZGVyLAogICAgICAgICAgb2xkX3RleHQgPSByZXN1bHRzLCB3cml0ZV9mbGFnID0gd3JpdGVfZmxhZwogICAgICAgICkKICAgICAgICBwbHRzIDwtIHZpel9yZXN1bHRzW1t2aXpfbmFtZV1dCiAgICAgICAgaWYgKCFpbmhlcml0cyhwbHRzLCAibGlzdCIpKSB7CiAgICAgICAgICBwbHRzIDwtIGxpc3QocGx0ID0gcGx0cykKICAgICAgICB9CiAgICAgICAgaWYgKGlzLm51bGwobmFtZXMocGx0cykpKSB7CiAgICAgICAgICBuYW1lcyhwbHRzKSA8LSAxOmxlbmd0aChwbHRzKQogICAgICAgIH0KICAgICAgICBmb3IgKHBsdF9uYW1lIGluIG5hbWVzKHBsdHMpKSB7CiAgICAgICAgICBpZiAobGVuZ3RoKHBsdHMpICE9IDEpIHsKICAgICAgICAgICAgcmVzdWx0cyA8LSB3cml0ZSgKICAgICAgICAgICAgICBzcHJpbnRmKHBsdF90ZW1wbGF0ZSwgcGx0X25hbWUpLAogICAgICAgICAgICAgIG9sZF90ZXh0ID0gcmVzdWx0cywgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICAgICAgICAgKQogICAgICAgICAgfQogICAgICAgICAgcGx0IDwtIHBsdHNbW3BsdF9uYW1lXV0KICAgICAgICAgIGlzX3Bsb3QgPC0gaW5oZXJpdHMocGx0LCAicGxvdGx5IikgfHwgCiAgICAgICAgICAgIGluaGVyaXRzKHBsdCwgImdnIikgfHwgCiAgICAgICAgICAgIGluaGVyaXRzKHBsdCwgImdncGxvdCIpCiAgICAgICAgICAKICAgICAgICAgIGlmIChwYXJhbXMkdXNlX3Ztb2Rlcm4gJiYgaXNfcGxvdCkgewogICAgICAgICAgICBjaHVua19hcmdzIDwtICJmaWcuaGVpZ2h0ID0gJXMsIGZpZy53aWR0aCA9ICVzLCBvdXQud2lkdGggPSAnMTAwJSUnLCBhZGQucGFuZWwgPSBUUlVFIgogICAgICAgICAgICBhZGRfY2xhc3MgPC0gInBhbmVsIHBhbmVsLWRlZmF1bHQgcGFkZGVkLXBhbmVsIgogICAgICAgICAgfSBlbHNlIHsKICAgICAgICAgICAgY2h1bmtfYXJncyA8LSAiZmlnLmhlaWdodCA9ICVzLCBmaWcud2lkdGggPSAlcywgb3V0LndpZHRoID0gJzEwMCUlJyIKICAgICAgICAgICAgYWRkX2NsYXNzIDwtIE5VTEwKICAgICAgICAgIH0KICAgICAgICAgIAogICAgICAgICAgaWYgKHdyaXRlX2ZsYWcpIHsKICAgICAgICAgICAgcmVzdWx0cyA8LSBzcHJpbnRmKAogICAgICAgICAgICAgICJzaG93X3Jlc3VsdHMoJyVzJywgJyVzJywgJ3Zpc3VhbGl6ZXInKSIsIGRpcl9uYW1lLCB2aXpfbmFtZQogICAgICAgICAgICApICU+JQogICAgICAgICAgICAgIHdyaXRlX2NvZGVfY2h1bmsoCiAgICAgICAgICAgICAgICBjaHVua19hcmdzID0gc3ByaW50ZigKICAgICAgICAgICAgICAgICAgY2h1bmtfYXJncywKICAgICAgICAgICAgICAgICAgdmlzdWFsaXplciRkb2Nfb3B0aW9ucyRoZWlnaHQsIHZpc3VhbGl6ZXIkZG9jX29wdGlvbnMkd2lkdGgKICAgICAgICAgICAgICAgICkKICAgICAgICAgICAgICApICU+JQogICAgICAgICAgICAgIHdyaXRlKG9sZF90ZXh0ID0gcmVzdWx0cywgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcpCiAgICAgICAgICB9IGVsc2UgewogICAgICAgICAgICB2dGhlbWVzOjpzdWJjaHVua2lmeShwbHQsCiAgICAgICAgICAgICAgaSA9IGNodW5rX2lkeCwKICAgICAgICAgICAgICBmaWdfaGVpZ2h0ID0gdmlzdWFsaXplciRkb2Nfb3B0aW9ucyRoZWlnaHQsCiAgICAgICAgICAgICAgZmlnX3dpZHRoID0gdmlzdWFsaXplciRkb2Nfb3B0aW9ucyR3aWR0aCwKICAgICAgICAgICAgICBvdGhlcl9hcmdzID0gIm91dC53aWR0aCA9ICcxMDAlJyIsIAogICAgICAgICAgICAgIGFkZF9jbGFzcyA9IGFkZF9jbGFzcwogICAgICAgICAgICApCiAgICAgICAgICAgIGNodW5rX2lkeCA8PC0gY2h1bmtfaWR4ICsgMQogICAgICAgICAgfQogICAgICAgIH0KICAgICAgfQogICAgfQogIH0KCiAgaWYgKCFpcy5udWxsKGV4cCkgJiYgcGFyYW1zJHNob3dfY29kZSkgewogICAgaWYgKChsZW5ndGgoZXhwJGdldF92YXJ5X2Fjcm9zcygpJGRncCkgIT0gMCkgfHwKICAgICAgICAobGVuZ3RoKGV4cCRnZXRfdmFyeV9hY3Jvc3MoKSRtZXRob2QpICE9IDApKSB7CiAgICAgIHJlc3VsdHMgPC0gd3JpdGUoCiAgICAgICAgc3ByaW50ZihzaG93dHlwZV90ZW1wbGF0ZSwgY29kZV9sYWJlbCksCiAgICAgICAgIlxuXG4qKlBhcmFtZXRlciBWYWx1ZXMqKlxuXG4iLAogICAgICAgIG9sZF90ZXh0ID0gcmVzdWx0cywgd3JpdGVfZmxhZyA9IHdyaXRlX2ZsYWcKICAgICAgKQogICAgICBpZiAod3JpdGVfZmxhZykgewogICAgICAgIHJlc3VsdHMgPC0gc3ByaW50ZigKICAgICAgICAgICJzaG93X3Jlc3VsdHMoJyVzJywgTlVMTCwgJ3ZhcnlfcGFyYW1zJykiLCBkaXJfbmFtZQogICAgICAgICkgJT4lCiAgICAgICAgICB3cml0ZV9jb2RlX2NodW5rKGNodW5rX2FyZ3MgPSAibWF4LmhlaWdodD0nMjAwcHgnIikgJT4lCiAgICAgICAgICB3cml0ZShvbGRfdGV4dCA9IHJlc3VsdHMsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnKQogICAgICB9IGVsc2UgewogICAgICAgIHZ0aGVtZXM6OnN1YmNodW5raWZ5KGV4cCRnZXRfdmFyeV9hY3Jvc3MoKSwKICAgICAgICAgIGNodW5rX2lkeCwKICAgICAgICAgIG90aGVyX2FyZ3MgPSAibWF4LmhlaWdodD0nMjAwcHgnIgogICAgICAgICkKICAgICAgICBjaHVua19pZHggPDwtIGNodW5rX2lkeCArIDEKICAgICAgfQogICAgfQogIH0KCiAgcmV0dXJuKHJlc3VsdHMpCn0KCiMnIERpc3BsYXlzIG91dHB1dCBvZiBleHBlcmltZW50IGZvciBhbGwgb2YgaXRzIChzYXZlZCkgZGVzY2VuZGFudHMKIycKIycgQHBhcmFtIGRpcl9uYW1lIG5hbWUgb2YgcGFyZW50IGV4cGVyaW1lbnQgZGlyZWN0b3J5CiMnIEBwYXJhbSBkZXB0aCBwbGFjZWhvbGRlciBmb3IgcmVjdXJzaW9uOyBzaG91bGQgbm90IGJlIG1lc3NlZCB3aXRoCiMnIEBwYXJhbSB3cml0ZV9mbGFnIEJvb2xlYW4gaW5kaWNhdGluZyB3aGV0aGVyIHRvIHdyaXRlIHRleHQgdG8gYSBmaWxlCiMnICAgKHdyaXRlX2ZsYWcgPSBUUlVFKSBvciB0byBjb25zb2xlICh3cml0ZV9mbGFnID0gRkFMU0UpCiMnIEBwYXJhbSB3cml0ZV9maWxlbmFtZSBOYW1lIG9mIGZpbGUgdG8gd3JpdGUgdG8gaWYgd3JpdGVfZmxhZyA9IFRSVUUKIycgQHBhcmFtIC4uLiBvdGhlciBhcmd1bWVudHMgdG8gcGFzcyBpbnRvIHNob3dfcmVzdWx0cygpCnNob3dfZGVzY2VuZGFudF9yZXN1bHRzIDwtIGZ1bmN0aW9uKGRpcl9uYW1lLCBkZXB0aCA9IDEsIHdyaXRlX2ZsYWcgPSBGQUxTRSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgd3JpdGVfZmlsZW5hbWUgPSBOVUxMLCAuLi4pIHsKICBjaGlsZHJlbiA8LSBsaXN0LmRpcnMoZGlyX25hbWUsIHJlY3Vyc2l2ZSA9IEZBTFNFKQogIGlmIChsZW5ndGgoY2hpbGRyZW4pID09IDApIHsKICAgIHJldHVybigpCiAgfQogIGZvciAoY2hpbGRfaWR4IGluIDE6bGVuZ3RoKGNoaWxkcmVuKSkgewogICAgY2hpbGQgPC0gY2hpbGRyZW5bY2hpbGRfaWR4XQogICAgaWYgKCFleHBlcmltZW50X2V4aXN0cyhjaGlsZCwgcmVjdXJzaXZlID0gVFJVRSkpIHsKICAgICAgbmV4dAogICAgfQogICAgaWYgKGV4cGVyaW1lbnRfZXhpc3RzKGNoaWxkLCByZWN1cnNpdmUgPSBGQUxTRSkgJiYKICAgICAgKGV4cGVyaW1lbnRfZXhpc3RzKGxpc3QuZGlycyhjaGlsZCwgcmVjdXJzaXZlID0gVFJVRSlbLTFdKSB8fAogICAgICAgIChkZXB0aCA9PSAxKSkpIHsKICAgICAgYmFzZSA8LSBUUlVFCiAgICB9IGVsc2UgewogICAgICBiYXNlIDwtIEZBTFNFCiAgICB9CiAgICByZXN1bHRzIDwtIHNob3dfcmVzdWx0cyhjaGlsZCwgZGVwdGgsIGJhc2UsIHdyaXRlX2ZsYWcgPSB3cml0ZV9mbGFnLCAuLi4pCiAgICBpZiAod3JpdGVfZmxhZykgewogICAgICB3cml0ZV90b19maWxlKHBhdGggPSB3cml0ZV9maWxlbmFtZSwgcmVzdWx0cykKICAgIH0KICAgIHNob3dfZGVzY2VuZGFudF9yZXN1bHRzKGNoaWxkLCBkZXB0aCArIDEsIHdyaXRlX2ZsYWcsIHdyaXRlX2ZpbGVuYW1lLCAuLi4pCiAgfQp9CgojJyBDbGVhbiBvdXRwdXQgZmlsZSAoZS5nLiwgcmVtb3ZlIGV4Y2Vzc2l2ZSBibGFuayBsaW5lcykKIycKIycgQHBhcmFtIHBhdGggUGF0aCB0byBvdXRwdXQgZmlsZQpjbGVhbl9maWxlIDwtIGZ1bmN0aW9uKHBhdGgpIHsKICBzdG9yZWxpbmVzIDwtIHJlYWRMaW5lcyhwYXRoKQogIHJsZV9vdXQgPC0gcmxlKHN0b3JlbGluZXMgPT0gIiIpCiAgbGluZV9pZHMgPC0gd2hpY2goKHJsZV9vdXQkbGVuZ3RocyA+IDIpICYgcmxlX291dCR2YWx1ZXMpCiAga2VlcF9saW5lcyA8LSByZXAoVFJVRSwgbGVuZ3RoKHN0b3JlbGluZXMpKQogIGZvciAobGluZV9pZCBpbiBsaW5lX2lkcykgewogICAgbnVtX2JsYW5rIDwtIHJsZV9vdXQkbGVuZ3Roc1tsaW5lX2lkXQogICAgbGluZV9wdHIgPC0gc3VtKHJsZV9vdXQkbGVuZ3Roc1sxOmxpbmVfaWRdKQogICAgIyBvbmx5IGFsbG93IGZvciBtYXggb2YgdHdvIGNvbnNlY3V0aXZlIGJsYW5rIGxpbmVzCiAgICBrZWVwX2xpbmVzWyhsaW5lX3B0ciAtIG51bV9ibGFuayArIDMpOmxpbmVfcHRyXSA8LSBGQUxTRQogIH0KICB3cml0ZUxpbmVzKHN0b3JlbGluZXNba2VlcF9saW5lc10sIHBhdGgpCn0KCiMnIEluc2VydCBsaW5lcyB0byBhZGQgZXh0cmEgcmVzb3VyY2VzIChjc3MvanMpIGZvciBzaW1DaGVmIFIgTWFya2Rvd24gdGhlbWUKIycgCiMnIEBwYXJhbSBwYXRoIFBhdGggdG8gb3V0cHV0IGZpbGUKaW5zZXJ0X3NpbUNoZWZfcmVzb3VyY2VzIDwtIGZ1bmN0aW9uKHBhdGgpIHsKICBzdG9yZWxpbmVzIDwtIHJlYWRMaW5lcyhwYXRoKQogIHBhdHRlcm4gPC0gIjxJbnNlcnQgZXh0cmEgc2ltQ2hlZiByZXNvdXJjZXMgaGVyZT4iCiAgcmVwbGFjZSA8LSBzcHJpbnRmKAogICAgJzxzY3JpcHQgc3JjPSIlcyI+PC9zY3JpcHQ+XG5cbjxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iJXMiPicsIAogICAgc3lzdGVtLmZpbGUoInJtZCIsICJqcyIsICJzaW1jaGVmTmF2Q2xhc3MuanMiLAogICAgICAgICAgICAgICAgcGFja2FnZSA9IHV0aWxzOjpwYWNrYWdlTmFtZSgpKSwKICAgIHN5c3RlbS5maWxlKCJybWQiLCAiY3NzIiwgInNpbWNoZWYuY3NzIiwKICAgICAgICAgICAgICAgIHBhY2thZ2UgPSB1dGlsczo6cGFja2FnZU5hbWUoKSkKICApIAogIHN0b3JlbGluZXNbc3RvcmVsaW5lcyA9PSBwYXR0ZXJuXSA8LSByZXBsYWNlCiAgd3JpdGVMaW5lcyhzdG9yZWxpbmVzLCBwYXRoKQp9CgojJyBSZW1vdmUgbGluZXMgd2l0aCBzaW1DaGVmIFIgTWFya2Rvd24gdGhlbWUtc3BlY2lmaWMgY29kZQojJyAKIycgQHBhcmFtIHBhdGggUGF0aCB0byBvdXRwdXQgZmlsZQpyZW1vdmVfc2ltQ2hlZl9yZXNvdXJjZXMgPC0gZnVuY3Rpb24ocGF0aCkgewogIHN0b3JlbGluZXMgPC0gcmVhZExpbmVzKHBhdGgpCiAgCiAgcGF0dGVybiA8LSAiYWRkLnBhbmVsID0gZnVuY3Rpb24iCiAgbGluZV9pZCA8LSB3aGljaChzdHJpbmdyOjpzdHJfZGV0ZWN0KHN0b3JlbGluZXMsIHBhdHRlcm4pKQogIHJlbW92ZV9saW5lcyA8LSAobGluZV9pZCAtIDIpOihsaW5lX2lkICsgNSkKICBzdG9yZWxpbmVzIDwtIHN0b3JlbGluZXNbLXJlbW92ZV9saW5lc10KICAKICBwYXR0ZXJuIDwtICI8SW5zZXJ0IGV4dHJhIHNpbUNoZWYgcmVzb3VyY2VzIGhlcmU+IgogIHJlbW92ZV9saW5lcyA8LSB3aGljaChzdHJpbmdyOjpzdHJfZGV0ZWN0KHN0b3JlbGluZXMsIHBhdHRlcm4pKQogIHN0b3JlbGluZXMgPC0gc3RvcmVsaW5lc1stcmVtb3ZlX2xpbmVzXQogIAogIHdyaXRlTGluZXMoc3RvcmVsaW5lcywgcGF0aCkKfQoKYGBgCgpgYGB7cn0KaWYgKHBhcmFtcyR3cml0ZSkgewogIGlmIChwYXJhbXMkdXNlX3Ztb2Rlcm4pIHsKICAgIGluc2VydF9zaW1DaGVmX3Jlc291cmNlcyh3cml0ZV9maWxlbmFtZSkKICB9IGVsc2UgewogICAgcmVtb3ZlX3NpbUNoZWZfcmVzb3VyY2VzKHdyaXRlX2ZpbGVuYW1lKQogIH0KfSBlbHNlIHsKICBpZiAocGFyYW1zJHVzZV92bW9kZXJuKSB7CiAgICBodG1sdG9vbHM6OkhUTUwoJzxzY3JpcHQgc3JjPSJqcy9zaW1jaGVmTmF2Q2xhc3MuanMiPjwvc2NyaXB0PlxuXG48bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9ImNzcy9zaW1jaGVmLmNzcyI+JykKICB9Cn0KYGBgCgoKIyBTaW11bGF0aW9uIEV4cGVyaW1lbnQgUmVjaXBlIHsudGFic2V0IC50YWJzZXQtdm1vZGVybn0KCiMjIE9iamVjdGl2ZXMKCmBgYHtyIG9iamVjdGl2ZXMsIHJlc3VsdHMgPSAiYXNpcyJ9CmlmIChwYXJhbXMkdXNlX3Ztb2Rlcm4pIHsKICBvYmplY3RpdmVzIDwtIHdyaXRlKAogICAgIlxuXG48ZGl2IGNsYXNzPSdwYW5lbCBwYW5lbC1kZWZhdWx0IHBhZGRlZC1wYW5lbCc+XG5cbiIsCiAgICBwYXN0ZU1kKGZpbGUucGF0aChkb2NfZGlyLCAib2JqZWN0aXZlcy5tZCIpKSwKICAgICJcblxuPC9kaXY+XG5cbiIsCiAgICB3cml0ZV9mbGFnID0gcGFyYW1zJHdyaXRlCiAgKQp9IGVsc2UgewogIG9iamVjdGl2ZXMgPC0gd3JpdGUoCiAgICBwYXN0ZU1kKGZpbGUucGF0aChkb2NfZGlyLCAib2JqZWN0aXZlcy5tZCIpKSwKICAgIHdyaXRlX2ZsYWcgPSBwYXJhbXMkd3JpdGUKICApCn0KaWYgKHBhcmFtcyR3cml0ZSkgewogIHdyaXRlX3RvX2ZpbGUocGF0aCA9IHdyaXRlX2ZpbGVuYW1lLCAiXG5cbiMjIE9iamVjdGl2ZXNcblxuIiwgb2JqZWN0aXZlcykKfQpgYGAKCiMjIERhdGEgR2VuZXJhdGlvbgoKYGBge3IgZGdwcywgcmVzdWx0cyA9ICJhc2lzIn0KZGdwX3JlY2lwZSA8LSBzaG93X3JlY2lwZShmaWVsZF9uYW1lID0gImRncCIsIHdyaXRlX2ZsYWcgPSBwYXJhbXMkd3JpdGUpCmlmIChwYXJhbXMkd3JpdGUpIHsKICB3cml0ZV90b19maWxlKHBhdGggPSB3cml0ZV9maWxlbmFtZSwgIlxuXG4jIyBEYXRhIEdlbmVyYXRpb25cblxuIiwgZGdwX3JlY2lwZSkKfQpgYGAKCiMjIE1ldGhvZHMgYW5kIEV2YWx1YXRpb24KCiMjIyBNZXRob2RzCgpgYGB7ciBtZXRob2RzLCByZXN1bHRzID0gImFzaXMifQptZXRob2RfcmVjaXBlIDwtIHNob3dfcmVjaXBlKGZpZWxkX25hbWUgPSAibWV0aG9kIiwgd3JpdGVfZmxhZyA9IHBhcmFtcyR3cml0ZSkKaWYgKHBhcmFtcyR3cml0ZSkgewogIHdyaXRlX3RvX2ZpbGUoCiAgICBwYXRoID0gd3JpdGVfZmlsZW5hbWUsCiAgICAiXG5cbiMjIE1ldGhvZHMgYW5kIEV2YWx1YXRpb25cblxuIiwgIlxuXG4jIyMgTWV0aG9kc1xuXG4iLAogICAgbWV0aG9kX3JlY2lwZQogICkKfQpgYGAKCiMjIyBFdmFsdWF0aW9uCgpgYGB7ciBldmFsdWF0b3JzLCByZXN1bHRzID0gImFzaXMifQpldmFsX3JlY2lwZSA8LSBzaG93X3JlY2lwZShmaWVsZF9uYW1lID0gImV2YWx1YXRvciIsIHdyaXRlX2ZsYWcgPSBwYXJhbXMkd3JpdGUpCmlmIChwYXJhbXMkd3JpdGUpIHsKICB3cml0ZV90b19maWxlKHBhdGggPSB3cml0ZV9maWxlbmFtZSwgIlxuXG4jIyMgRXZhbHVhdGlvblxuXG4iLCBldmFsX3JlY2lwZSkKfQpgYGAKCiMjIFZpc3VhbGl6YXRpb25zCgpgYGB7ciB2aXN1YWxpemVycywgcmVzdWx0cyA9ICJhc2lzIn0Kdml6X3JlY2lwZSA8LSBzaG93X3JlY2lwZShmaWVsZF9uYW1lID0gInZpc3VhbGl6ZXIiLCB3cml0ZV9mbGFnID0gcGFyYW1zJHdyaXRlKQppZiAocGFyYW1zJHdyaXRlKSB7CiAgd3JpdGVfdG9fZmlsZShwYXRoID0gd3JpdGVfZmlsZW5hbWUsICJcblxuIyMgVmlzdWFsaXphdGlvbnNcblxuIiwgdml6X3JlY2lwZSkKfQpgYGAKCgoKYGBge3IgcmVzLCByZXN1bHRzID0gImFzaXMifQppZiAocGFyYW1zJHZlcmJvc2UgPiAwKSB7CiAgaW5mb3JtKHNwcmludGYoIkNyZWF0aW5nIFIgTWFya2Rvd24gcmVwb3J0IGZvciAlcy4uLiIsIHBhcmFtcyRzaW1fbmFtZSkpCn0KCiMgc2hvdyByZXN1bHRzCmlmIChleHBlcmltZW50X2V4aXN0cyhwYXJhbXMkc2ltX3BhdGgpKSB7CiAgYmFzZV9oZWFkZXIgPC0gd3JpdGUoCiAgICBzcHJpbnRmKCJcblxuIyBCYXNlICVzIFxuXG4iLCBwYXJhbXMkc2ltX25hbWUpLAogICAgIlxuXG4jIyB7LnRhYnNldCAudGFic2V0LXBpbGxzIC50YWJzZXQtY2lyY2xlfVxuXG4iLAogICAgd3JpdGVfZmxhZyA9IHBhcmFtcyR3cml0ZQogICkKICBiYXNlX3Jlc3VsdHMgPC0gc2hvd19yZXN1bHRzKAogICAgcGFyYW1zJHNpbV9wYXRoLAogICAgZGVwdGggPSAyLCBiYXNlID0gRkFMU0UsIHNob3dfaGVhZGVyID0gRkFMU0UsCiAgICB2ZXJib3NlID0gcGFyYW1zJHZlcmJvc2UsIHdyaXRlX2ZsYWcgPSBwYXJhbXMkd3JpdGUKICApCgogIGlmIChwYXJhbXMkd3JpdGUpIHsKICAgIHdyaXRlX3RvX2ZpbGUocGF0aCA9IHdyaXRlX2ZpbGVuYW1lLCBiYXNlX2hlYWRlciwgYmFzZV9yZXN1bHRzKQogIH0KfQoKc2hvd19kZXNjZW5kYW50X3Jlc3VsdHMoCiAgcGFyYW1zJHNpbV9wYXRoLAogIHZlcmJvc2UgPSBwYXJhbXMkdmVyYm9zZSwKICB3cml0ZV9mbGFnID0gcGFyYW1zJHdyaXRlLCAKICB3cml0ZV9maWxlbmFtZSA9IHdyaXRlX2ZpbGVuYW1lCikKYGBgCgoKCmBgYHtyfQppZiAocGFyYW1zJHdyaXRlKSB7CiAgY2xlYW5fZmlsZShwYXRoID0gd3JpdGVfZmlsZW5hbWUpCn0KYGBgCg==